Samuele E. Locatelli
383 lines
14 KiB
PHP
383 lines
14 KiB
PHP
<?php
|
|
//USGS earthquakes
|
|
$ini = parse_ini_file("conf.ini", true);
|
|
mysql_connect($ini['db']['host'],$ini['db']['user'],$ini['db']['pass']);
|
|
mysql_select_db($ini['db']['db']);
|
|
|
|
if ($ini['db']['log']=="1"){
|
|
$arrival_time=microtime(true);
|
|
$filename=basename(__FILE__);
|
|
$filename="log_".$filename.".txt";
|
|
$chunck=strval($arrival_time)." ".$_SERVER['REQUEST_URI']."\r\n";
|
|
file_put_contents($filename,$chunck,FILE_APPEND);
|
|
}
|
|
|
|
$url="http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/2.5_week.geojson";
|
|
$json = file_get_contents($url);
|
|
$data = json_decode($json, TRUE);
|
|
|
|
//extract the maximum date
|
|
$sql=mysql_query("select date from last_notifications where provider='USGS'");
|
|
if (mysql_num_rows($sql)){
|
|
$max_date=mysql_result($sql,0,"date");
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_usgs_temp");
|
|
|
|
foreach($data['features'] as $quake)
|
|
{
|
|
$place=mysql_real_escape_string($quake['properties']['place']);
|
|
$stamp2=$quake['properties']['time']/1000;
|
|
$stamp=date('Y-m-d H:i:s',$stamp2);
|
|
$felt=$quake['properties']['felt'];
|
|
$cdi=$quake['properties']['cdi'];
|
|
$mmi=$quake['properties']['mmi'];
|
|
$tsunami=$quake['properties']['tsunami'];
|
|
$url=$quake['properties']['url'];
|
|
if ($felt==null){
|
|
$felt=0;
|
|
}
|
|
if ($cdi==null){
|
|
$cdi=0;
|
|
}
|
|
if ($mmi==null){
|
|
$mmi=0;
|
|
}
|
|
|
|
$cdi=max($cdi,$mmi);
|
|
|
|
if ($tsunami==null){
|
|
$tsunami=-1;
|
|
}
|
|
|
|
$magnitude=(float)$quake['properties']['mag'];
|
|
$magnitude_type=$quake['properties']['magType'];
|
|
if (strcmp($magnitude_type,'mb_lg')==0){
|
|
$magnitude_type='mbLg';
|
|
}
|
|
if (strcmp($magnitude_type,'ml')==0){
|
|
$magnitude_type='ML';
|
|
}
|
|
|
|
$lat=(float)$quake['geometry']['coordinates'][1];
|
|
$lon=(float)$quake['geometry']['coordinates'][0];
|
|
$depth=(float)$quake['geometry']['coordinates'][2];
|
|
$code=$quake['properties']['code'];
|
|
|
|
$lat_idx=$lat;
|
|
$lat_idx=(round($lat_idx*10)/10+0.025)*1000;
|
|
if ($lat_idx>90000){
|
|
$lat_idx=90000-25;
|
|
}
|
|
$lon_idx=$lon;
|
|
$lon_idx=(round($lon_idx*10)/10+0.025)*1000;
|
|
if ($lon_idx>180000){
|
|
$lon_idx=$lon_idx-360000;
|
|
}
|
|
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop_recent where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
if (mysql_num_rows($sql)){
|
|
$pop25=mysql_result($sql,0,"pop25");
|
|
$pop100=mysql_result($sql,0,"pop100");
|
|
$pop200=mysql_result($sql,0,"pop200");
|
|
}
|
|
else{
|
|
if ($lon_idx<=-120025){
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop1 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
else{
|
|
if ($lon_idx<=-60025){
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop2 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
else{
|
|
if ($lon_idx<=-25){
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop3 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
else{
|
|
if ($lon_idx<=59975){
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop4 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
else{
|
|
if ($lon_idx<=119975){
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop5 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
else{
|
|
$sql=mysql_query("select pop25,pop100,pop200 from worldpop6 where latitude=$lat_idx and longitude=$lon_idx LIMIT 1");
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
if (mysql_num_rows($sql))
|
|
{
|
|
$pop25=mysql_result($sql,0,"pop25");
|
|
$pop100=mysql_result($sql,0,"pop100");
|
|
$pop200=mysql_result($sql,0,"pop200");
|
|
}
|
|
else
|
|
{
|
|
$pop25=0;
|
|
$pop100=0;
|
|
$pop200=0;
|
|
}
|
|
$sql=mysql_query("insert into worldpop_recent (latitude,longitude,pop25,pop100,pop200) values ($lat_idx,$lon_idx,$pop25,$pop100,$pop200)");
|
|
}
|
|
|
|
$ok=0;
|
|
$sql=mysql_query("select * from weather_recent where abs(latitude-$lat)<0.01 and abs(longitude-$lon)<0.01 and abs(TIMESTAMPDIFF(MINUTE,earthquake_date,'$stamp'))<5 LIMIT 1");
|
|
if (mysql_num_rows($sql)){
|
|
$weather_code=mysql_result($sql,0,"weather_code");
|
|
$temperature=mysql_result($sql,0,"temperature");
|
|
$humidity=mysql_result($sql,0,"humidity");
|
|
$pressure=mysql_result($sql,0,"pressure");
|
|
$wind_speed=mysql_result($sql,0,"wind_speed");
|
|
$wind_dir=mysql_result($sql,0,"wind_dir");
|
|
$clouds=mysql_result($sql,0,"clouds");
|
|
$icon=mysql_result($sql,0,"icon");
|
|
$ok=1;
|
|
}
|
|
else{
|
|
$current_time=time();
|
|
if ($current_time-$stamp2<3600){
|
|
$url_weather="http://api.openweathermap.org/data/2.5/weather?lat=".$lat."&lon=".$lon."&appid=2de02dac92c8a20f8571feb5ad5bdcaa";
|
|
$json = file_get_contents($url_weather);
|
|
$data_weather = json_decode($json, TRUE);
|
|
|
|
$weather_code=$data_weather["weather"][0]["id"];
|
|
$temperature=$data_weather["main"]["temp"];
|
|
$humidity=$data_weather["main"]["humidity"];
|
|
$pressure=$data_weather["main"]["pressure"];
|
|
if (array_key_exists('speed', $data_weather["wind"])){
|
|
$wind_speed=$data_weather["wind"]["speed"];
|
|
}
|
|
else{
|
|
$wind_speed=0;
|
|
}
|
|
if (array_key_exists('deg', $data_weather["wind"])){
|
|
$wind_dir=$data_weather["wind"]["deg"];
|
|
}
|
|
else{
|
|
$wind_dir=0;
|
|
}
|
|
$clouds=$data_weather["clouds"]["all"];
|
|
$icon=$data_weather["weather"][0]["icon"];
|
|
|
|
if (!is_null($weather_code)){
|
|
$sql=mysql_query("insert into weather_recent (latitude,longitude,earthquake_date,weather_code,temperature,humidity,pressure,wind_speed,wind_dir,clouds,icon) values ($lat,$lon,'$stamp',$weather_code,$temperature,$humidity,$pressure,$wind_speed,$wind_dir,$clouds,'$icon')");
|
|
$ok=1;
|
|
}
|
|
}
|
|
else{
|
|
$ok=0;
|
|
}
|
|
}
|
|
|
|
if ($magnitude>=5.5){
|
|
if ($magnitude>=5.5 and $magnitude<6.0)
|
|
{
|
|
$max_distance=40;
|
|
}
|
|
else{
|
|
if ($magnitude>=6.0 and $magnitude<6.5){
|
|
$max_distance=100;
|
|
}
|
|
else{
|
|
if ($magnitude>=6.5 and $magnitude<7.0){
|
|
$max_distance=200;
|
|
}
|
|
else{
|
|
if ($magnitude>=7.0 and $magnitude<7.5){
|
|
$max_distance=500;
|
|
}
|
|
else{
|
|
if ($magnitude>=7.5 and $magnitude<8.0){
|
|
$max_distance=800;
|
|
}
|
|
else{
|
|
if ($magnitude>=8.0 and $magnitude<8.5){
|
|
$max_distance=1000;
|
|
}
|
|
else{
|
|
if ($magnitude>=8.5 and $magnitude<9.0){
|
|
$max_distance=2000;
|
|
}
|
|
else{
|
|
if ($magnitude>=9.0){
|
|
$max_distance=3000;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
$max_distance_deg=$max_distance/111;
|
|
$sql=mysql_query("select count(*) as count from manual_picture where abs(latitude-$lat)<$max_distance_deg and abs(longitude-$lon)<$max_distance_deg and abs(TIMESTAMPDIFF(MINUTE,date,'$stamp'))<1440");
|
|
if (mysql_num_rows($sql)){
|
|
$picture_count=mysql_result($sql,0,"count");
|
|
}
|
|
else{
|
|
$picture_count=0;
|
|
}
|
|
}
|
|
else{
|
|
$picture_count=0;
|
|
}
|
|
|
|
if ($ok==1){
|
|
$sql=mysql_query("insert into automatic_usgs_temp (latitude,longitude,depth,magnitude,magnitude_type,pop25,pop100,pop200,place,data,code,tsunami,provider,weather_code,temperature,humidity,pressure,wind_speed,wind_dir,clouds,icon,picture_count) values ($lat,$lon,$depth,$magnitude,'$magnitude_type',$pop25,$pop100,$pop200,'$place','$stamp','$code',$tsunami,'USGS',$weather_code,$temperature,$humidity,$pressure,$wind_speed,$wind_dir,$clouds,'$icon',$picture_count)");
|
|
}
|
|
else{
|
|
$sql=mysql_query("insert into automatic_usgs_temp (latitude,longitude,depth,magnitude,magnitude_type,pop25,pop100,pop200,place,data,code,tsunami,provider,picture_count) values ($lat,$lon,$depth,$magnitude,'$magnitude_type',$pop25,$pop100,$pop200,'$place','$stamp','$code',$tsunami,'USGS',$picture_count)");
|
|
}
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_usgs");
|
|
$sql=mysql_query("INSERT INTO automatic_usgs SELECT * FROM automatic_usgs_temp");
|
|
$sql=mysql_query("DELETE from automatic_all WHERE provider like 'USGS'");
|
|
$sql=mysql_query("INSERT INTO automatic_all (latitude, longitude, depth, magnitude, magnitude_type, felt, cdi, pop25, pop100, pop200, data, code, place, tsunami, url, provider, weather_code, temperature, humidity, pressure, wind_speed, wind_dir, clouds, icon, picture_count) SELECT latitude, longitude, depth, magnitude, magnitude_type, felt, cdi, pop25, pop100, pop200, data, code, place, tsunami, url, provider, weather_code, temperature, humidity, pressure, wind_speed, wind_dir, clouds, icon, picture_count FROM automatic_usgs_temp");
|
|
|
|
//generate the file
|
|
$filename="seismic_free.txt";
|
|
|
|
$sql=mysql_query('SET CHARACTER SET utf8');
|
|
$sql=mysql_query("SELECT latitude,longitude,depth,magnitude,magnitude_type,pop25,pop100,pop200,cdi AS intensity,place,tsunami,url,provider,data,TIMESTAMPDIFF(MINUTE,data,NOW()) AS difference FROM (select * from automatic_all) as a WHERE a.data > date_sub(NOW(), interval 1440 minute) AND a.magnitude>=2 ORDER BY a.data DESC");
|
|
if (mysql_num_rows($sql)){
|
|
while($e=mysql_fetch_assoc($sql))
|
|
$output_file[]=$e;
|
|
$json_file = json_encode($output_file);
|
|
file_put_contents($filename,$json_file);
|
|
}
|
|
|
|
//extract the new earthquakes
|
|
$sql=mysql_query("select latitude,longitude,magnitude,magnitude_type,place,provider,depth,data,pop100,TIMESTAMPDIFF(MINUTE,data,NOW()) AS difference from automatic_usgs_temp where TIMESTAMPDIFF(MINUTE,'$max_date',data)>1 order by TIMESTAMPDIFF(MINUTE,data,NOW()) desc");
|
|
if (mysql_num_rows($sql))
|
|
{
|
|
$num_rows = mysql_num_rows($sql);
|
|
for ($j = 0; $j < $num_rows; $j++)
|
|
{
|
|
$delay=mysql_result($sql,$j,"difference");
|
|
$sql_z=mysql_query("INSERT INTO delay_usgs (delay) values ($delay)");
|
|
|
|
$dataArray = array('latitude'=>mysql_result($sql,$j,"latitude"),
|
|
'longitude'=>mysql_result($sql,$j,"longitude"),
|
|
'magnitude'=>mysql_result($sql,$j,"magnitude"),
|
|
'magnitude_type'=>mysql_result($sql,$j,"magnitude"),
|
|
'place'=>mysql_result($sql,$j,"place"),
|
|
'provider'=>mysql_result($sql,$j,"provider"),
|
|
'depth'=>mysql_result($sql,$j,"depth"),
|
|
'difference'=>mysql_result($sql,$j,"difference"),
|
|
'data'=>mysql_result($sql,$j,"data"),
|
|
'pop100'=>mysql_result($sql,$j,"pop100"),
|
|
'type'=>'official');
|
|
|
|
//extract users for notifications
|
|
|
|
$magnitude_notification=mysql_result($sql,$j,"magnitude");
|
|
$current_date=mysql_result($sql,$j,"data");
|
|
|
|
$sql2=mysql_query("select user_id from gcm order by user_id limit 1");
|
|
$user_id_ref = mysql_result($sql2, 0, "user_id");
|
|
$user_id_ref=$user_id_ref-1;
|
|
$row_limit=250000;
|
|
$keep_going=1;
|
|
|
|
while ($keep_going){
|
|
$sql2=mysql_query("SELECT user_id,reg_id FROM gcm where gcm.reg_id<>'' and gcm.notification_official_usgs=1 and ($magnitude_notification>=gcm.notification_official_minmag or gcm.notification_near=1) and gcm.notification_official=1 and user_id>$user_id_ref order by user_id limit $row_limit");
|
|
$num_rows2 = mysql_num_rows($sql2);
|
|
|
|
if ($num_rows2>0){
|
|
while($data = mysql_fetch_assoc($sql2)){
|
|
$regidArray[] = $data['reg_id'];
|
|
$IDArray[] = $data['user_id'];
|
|
}
|
|
$user_id_ref=end($IDArray);
|
|
|
|
$block_size=999;
|
|
$iter=ceil($num_rows2/$block_size);
|
|
|
|
$folder_eqn=uniqid();
|
|
mkdir($folder_eqn);
|
|
$filename_data=$folder_eqn."/dataarray.json";
|
|
file_put_contents($filename_data,json_encode($dataArray));
|
|
|
|
for ($i=0; $i<$iter; $i++){
|
|
$index_start=$i*$block_size;
|
|
$index_end=($i+1)*$block_size-1;
|
|
if ($index_end>$num_rows2-1){
|
|
$index_end=$num_rows2-1;
|
|
}
|
|
$regidArrayGCM = array_slice($regidArray,$index_start,$index_end-$index_start+1);
|
|
|
|
$filename=$folder_eqn."/gcmid".$i.".json";
|
|
file_put_contents($filename,json_encode($regidArrayGCM));
|
|
|
|
exec("php /var/www/earthquakenetwork.it/mysql/distquake_gcm_official_call.php $filename_data $filename"."> /dev/null 2>/dev/null &");
|
|
sleep(0.1);
|
|
}
|
|
sleep(1);
|
|
$files = glob($folder_eqn.'/*',GLOB_MARK);
|
|
foreach ($files as $file) {
|
|
unlink($file);
|
|
}
|
|
rmdir($folder_eqn.'/');
|
|
unset($regidArray);
|
|
unset($IDArray);
|
|
}
|
|
else{
|
|
$keep_going=0;
|
|
}
|
|
}
|
|
|
|
$sql_z=mysql_query("select date FROM last_notifications where provider='USGS'");
|
|
$last_date=mysql_result($sql_z,0,"date");
|
|
if ($current_date>$last_date){
|
|
$sql_z=mysql_query("update last_notifications set date='$current_date' where provider='USGS'");
|
|
}
|
|
}
|
|
$sql_z=mysql_query("select max(ID) as max_id FROM delay_usgs");
|
|
$max_id=mysql_result($sql_z,0,"max_id");
|
|
$sql_z=mysql_query("delete from delay_usgs where ID<$max_id-40");
|
|
}
|
|
|
|
//pressure
|
|
|
|
$sql=mysql_query("select count(*) as mag_45 from automatic_usgs_temp where magnitude>4.5 and data > date_sub(now(), interval 3 day)");
|
|
$mag_45=mysql_result($sql,0,"mag_45");
|
|
$sql=mysql_query("select count(*) as mag_55 from automatic_usgs_temp where magnitude>5.5 and data > date_sub(now(), interval 3 day)");
|
|
$mag_55=mysql_result($sql,0,"mag_55");
|
|
$sql=mysql_query("select count(*) as mag_65 from automatic_usgs_temp where magnitude>6.5 and data > date_sub(now(), interval 3 day)");
|
|
$mag_65=mysql_result($sql,0,"mag_65");
|
|
|
|
$sql=mysql_query("select count(*) as mer_40 from automatic_usgs_temp where cdi>4.0 and data > date_sub(now(), interval 3 day)");
|
|
$mer_40=mysql_result($sql,0,"mer_40");
|
|
$sql=mysql_query("select count(*) as mer_60 from automatic_usgs_temp where cdi>6.0 and data > date_sub(now(), interval 3 day)");
|
|
$mer_60=mysql_result($sql,0,"mer_60");
|
|
$sql=mysql_query("select count(*) as mer_80 from automatic_usgs_temp where cdi>8.0 and data > date_sub(now(), interval 3 day)");
|
|
$mer_80=mysql_result($sql,0,"mer_80");
|
|
|
|
$sql=mysql_query("select magnitude,cdi,TIMESTAMPDIFF(MINUTE,data, NOW()) AS difference from automatic_usgs_temp where data > date_sub(now(), interval 3 day)");
|
|
|
|
$num_rows = mysql_num_rows($sql);
|
|
$geo_index=0;
|
|
$pop_index=0;
|
|
for ($i=0;$i<$num_rows;$i++){
|
|
$magnitude=mysql_result($sql,$i,"magnitude");
|
|
$cdi=mysql_result($sql,$i,"cdi");
|
|
$difference=mysql_result($sql,$i,"difference");
|
|
$factor_time=pow(1-$difference/(3*1440),0.3);
|
|
|
|
$geo_index=$geo_index+exp($magnitude)*$factor_time;
|
|
if ($cdi>0){
|
|
$pop_index=$pop_index+exp($cdi)*$factor_time;
|
|
}
|
|
}
|
|
|
|
$sql=mysql_query("insert into pressure (mag_45,mag_55,mag_65,mer_40,mer_60,mer_80,geo_index,pop_index) values ($mag_45,$mag_55,$mag_65,$mer_40,$mer_60,$mer_80,$geo_index,$pop_index)");
|
|
|
|
mysql_close();
|
|
?>
|