Samuele E. Locatelli
402 lines
16 KiB
PHP
402 lines
16 KiB
PHP
<?php
|
|
//INGV 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);
|
|
}
|
|
|
|
$ctx = stream_context_create(array(
|
|
'http' => array(
|
|
'timeout' => 15
|
|
)
|
|
)
|
|
);
|
|
|
|
$current_date=date('Y-m-d', time());
|
|
$first_date=strtotime("-1 day", strtotime($current_date));
|
|
$first_date=date("Y-m-d", $first_date);
|
|
|
|
$web="http://webservices.rm.ingv.it/fdsnws/event/1/query?starttime=".$first_date."+00%3A00%3A00&endtime=".$current_date."+23%3A59%3A59&last_nd=7&minmag=2&maxmag=10&minlat=35&maxlat=48&minlon=6&maxlon=19&limit=50&page=1&orderby=time&radius_lat=&radius_lon=&radius_km=-1&format=text&limit=250";
|
|
|
|
$data = file_get_contents($web,0,$ctx);
|
|
|
|
if ($data!=false)
|
|
{
|
|
//extract the maximum date
|
|
$sql=mysql_query("select date from last_notifications where provider='INGV'");
|
|
if (mysql_num_rows($sql)){
|
|
$max_date=mysql_result($sql,0,"date");
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_ingv_temp");
|
|
$result = explode(PHP_EOL, $data);
|
|
|
|
date_default_timezone_set("UTC");
|
|
$date1=date("Y-m-d H:i:s", time());
|
|
date_default_timezone_set("Europe/Rome");
|
|
$date2=date("Y-m-d H:i:s", time());
|
|
$hourdiff = round((strtotime($date2) - strtotime($date1))/3600, 1);
|
|
$delta='+'.$hourdiff.' hour';
|
|
|
|
for ($i=1;$i<count($result);$i++){
|
|
$result_single = explode("|", $result[$i]);
|
|
if (count($result_single)>1){
|
|
$code=$result_single[0];
|
|
|
|
$sql_place=mysql_query("select place from automatic_ingv where code=$code limit 1");
|
|
|
|
if (mysql_num_rows($sql_place))
|
|
{
|
|
$place=mysql_result($sql_place,0,"place");
|
|
$place=mysql_real_escape_string($place);
|
|
}
|
|
else
|
|
{
|
|
$location_data = file_get_contents("http://cnt.rm.ingv.it/event/".$code,0,$ctx);
|
|
$location_data = str_replace("\n", '', $location_data); // remove new lines
|
|
$location_data = str_replace("\r", '', $location_data); // remove carriage returns
|
|
preg_match_all('|<tr class.*?</tr>|',$location_data,$result_location);
|
|
if (!empty($result_location[0])){
|
|
preg_match_all('|<td>.*?</td>|',$result_location[0][0],$result_location2);
|
|
$result_location2=$result_location2[0];
|
|
$town=$result_location2[0];
|
|
$town=substr($town,4,strlen($town)-9);
|
|
$distance_km=$result_location2[2];
|
|
$distance_km=substr($distance_km,4,strlen($distance_km)-9);
|
|
$place=$result_single[12];
|
|
$place=$distance_km."km da ".$town.", Prov. ".$place;
|
|
}
|
|
else{
|
|
$place=$result_single[12];
|
|
}
|
|
$place=str_replace(''',' ',$place);
|
|
$place=str_replace('ò','o',$place);
|
|
$place=mysql_real_escape_string($place);
|
|
}
|
|
|
|
$lat=$result_single[2];
|
|
$lon=$result_single[3];
|
|
$depth=$result_single[4];
|
|
$magnitude_type=$result_single[9];
|
|
$magnitude=$result_single[10];
|
|
$date=$result_single[1];
|
|
$date=str_replace('T',' ',$date);
|
|
$timestamp = strtotime($date);
|
|
$stamp2=strtotime($delta,$timestamp);
|
|
$stamp=date('Y-m-d H:i:s',$stamp2);
|
|
$tsunami=-1;
|
|
|
|
if ($magnitude>=2){
|
|
$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_ingv_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,'INGV',$weather_code,$temperature,$humidity,$pressure,$wind_speed,$wind_dir,$clouds,'$icon',$picture_count)");
|
|
}
|
|
else{
|
|
$sql=mysql_query("insert into automatic_ingv_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,'INGV',$picture_count)");
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_ingv");
|
|
$sql=mysql_query("INSERT INTO automatic_ingv SELECT * FROM automatic_ingv_temp");
|
|
$sql=mysql_query("DELETE from automatic_all WHERE provider like 'INGV'");
|
|
$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_ingv_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_ingv_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_ingv (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_type"),
|
|
'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_ingv=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='INGV'");
|
|
$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='INGV'");
|
|
}
|
|
}
|
|
$sql_z=mysql_query("select max(ID) as max_id FROM delay_ingv");
|
|
$max_id=mysql_result($sql_z,0,"max_id");
|
|
$sql_z=mysql_query("delete from delay_ingv where ID<$max_id-40");
|
|
}
|
|
|
|
$sql=mysql_query("select min(abs(timestampdiff(minute,date,now()))) as delta_min from italian_swarm");
|
|
$delta_min=mysql_result($sql,0,"delta_min");
|
|
if (is_null($delta_min)){
|
|
$delta_min=9999;
|
|
}
|
|
if ($delta_min>14){
|
|
$sql=mysql_query("SET @prevValue:=(SELECT data FROM automatic_ingv where latitude>42 and latitude<43 and longitude>12 and longitude<14 and timestampdiff(MINUTE,data,now())<480 order by data limit 1)");
|
|
$sql=mysql_query("SELECT avg(delta) as mean_time from (SELECT abs(timestampdiff(MINUTE,data,@prevValue)) AS delta,@prevValue:=data FROM automatic_ingv where latitude>42 and latitude<43.5 and longitude>12 and longitude<14 and timestampdiff(MINUTE,data,now())<480 order by data) as res");
|
|
$mean_time=mysql_result($sql,0,"mean_time");
|
|
if ($mean_time>0){
|
|
$sql=mysql_query("SELECT avg(magnitude) as mean_magnitude from automatic_ingv where latitude>42 and latitude<43.5 and longitude>12 and longitude<14 and timestampdiff(MINUTE,data,now())<480");
|
|
$mean_magnitude=mysql_result($sql,0,"mean_magnitude");
|
|
$sql=mysql_query("SELECT max(magnitude) as max_magnitude from automatic_ingv where latitude>42 and latitude<43.5 and longitude>12 and longitude<14 and timestampdiff(MINUTE,data,now())<480");
|
|
$max_magnitude=mysql_result($sql,0,"max_magnitude");
|
|
|
|
$sql=mysql_query("INSERT INTO italian_swarm (mean_time,mean_magnitude,max_magnitude) values ($mean_time,$mean_magnitude,$max_magnitude)");
|
|
}
|
|
}
|
|
}
|
|
else{
|
|
print("Timeout");
|
|
}
|
|
|
|
mysql_close();
|
|
?>
|