Samuele E. Locatelli
410 lines
14 KiB
PHP
410 lines
14 KiB
PHP
<?php
|
|
//JMA 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);
|
|
}
|
|
|
|
mysql_set_charset("utf8");
|
|
|
|
$ctx = stream_context_create(array(
|
|
'http' => array(
|
|
'timeout' => 15
|
|
)
|
|
)
|
|
);
|
|
|
|
$data = file_get_contents("http://www.jma.go.jp/jp/quake/quake_local_index.html",0,$ctx);
|
|
if ($data!=false)
|
|
{
|
|
//extract the maximum date
|
|
$sql=mysql_query("select date from last_notifications where provider='JMA'");
|
|
if (mysql_num_rows($sql)){
|
|
$max_date=mysql_result($sql,0,"date");
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_jma_temp");
|
|
|
|
$data = str_replace("\n", '', $data); // remove new lines
|
|
$data = str_replace("\r", '', $data); // remove carriage returns
|
|
preg_match_all('|<a href=\..*?html>|',$data,$result);
|
|
$result=$result[0];
|
|
|
|
date_default_timezone_set("Asia/Tokyo");
|
|
$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=0;$i<count($result);$i++){
|
|
$address=$result[$i];
|
|
$address=substr($address,9,strlen($address)-10);
|
|
|
|
$url="http://www.jma.go.jp/jp/quake".$address;
|
|
|
|
$data2 = file_get_contents($url,0,$ctx);
|
|
$data2 = str_replace("\n", '', $data2); // remove new lines
|
|
$data2 = str_replace("\r", '', $data2); // remove carriage returns
|
|
|
|
$idx1=strpos($data2,'平成')+6;
|
|
$idx2=strpos($data2,'年');
|
|
$year=substr($data2,$idx1,$idx2-$idx1)+1988;
|
|
|
|
$idx1=strpos($data2,'年')+3;
|
|
$idx2=$idx1+2;
|
|
$month=substr($data2,$idx1,$idx2-$idx1);
|
|
|
|
//rimuove la parte iniziale di $data2
|
|
$idx1=strpos($data2,'気象庁発表<br>');
|
|
$data2=substr($data2,$idx1,1000);
|
|
|
|
$idx1=strpos($data2,'気象庁発表<br>')+19;
|
|
$idx2=strpos($data2,'日');
|
|
$day=substr($data2,$idx1,$idx2-$idx1);
|
|
|
|
$idx1=strpos($data2,'気象庁発表<br>');
|
|
$hour=substr($data2,$idx1+24,2);
|
|
$minute=substr($data2,$idx1+29,2);
|
|
|
|
$timestamp = strtotime($year.'-'.$month.'-'.$day.' '.$hour.':'.$minute.':00');
|
|
|
|
$stamp2=strtotime($delta,$timestamp);
|
|
$stamp=date('Y-m-d H:i:s',$stamp2);
|
|
|
|
$idx1=strpos($data2,'北緯')+6;
|
|
$idx2=strpos($data2,'度、東');
|
|
$lat=substr($data2,$idx1,$idx2-$idx1);
|
|
|
|
$idx1=strpos($data2,'度、東経')+12;
|
|
$idx2=strpos($data2,'度',$idx1);
|
|
$lon=substr($data2,$idx1,$idx2-$idx1);
|
|
|
|
$idx1=strpos($data2,'約')+3;
|
|
$idx2=strpos($data2,'km',$idx1);
|
|
if ($idx2==FALSE){
|
|
$depth="0";
|
|
}
|
|
else{
|
|
$depth=substr($data2,$idx1,$idx2-$idx1);
|
|
}
|
|
|
|
$idx1=strpos($data2,'）は')+6;
|
|
$idx2=strpos($data2,'と推');
|
|
$magnitude=substr($data2,$idx1,$idx2-$idx1);
|
|
$magnitude_type='ML';
|
|
|
|
$lat=$lat+$magnitude/500;
|
|
|
|
$idx1=strpos($data2,'震源地は')+12;
|
|
$idx2=strpos($data2,'（北緯');
|
|
$place=substr($data2,$idx1,$idx2-$idx1);
|
|
$code=hash('ripemd160',$place);
|
|
|
|
//rimuove la parte iniziale di $data2
|
|
$idx1=strpos($data2,'<pre>');
|
|
$data2=substr($data2,$idx1,100);
|
|
|
|
$idx1=strpos($data2,'震度')+6;
|
|
$intensity=substr($data2,$idx1,3);
|
|
$intensity=(string)bin2hex(iconv('UTF-8', 'UCS-2',$intensity));
|
|
$intensity=intval(substr($intensity,1,1));
|
|
|
|
$tsunami=-1;
|
|
|
|
$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_jma_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,'JMA',$weather_code,$temperature,$humidity,$pressure,$wind_speed,$wind_dir,$clouds,'$icon',$picture_count)");
|
|
}
|
|
else{
|
|
$sql=mysql_query("insert into automatic_jma_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,'JMA',$picture_count)");
|
|
}
|
|
}
|
|
|
|
//empty the table
|
|
$sql=mysql_query("TRUNCATE automatic_jma");
|
|
$sql=mysql_query("INSERT INTO automatic_jma SELECT * FROM automatic_jma_temp");
|
|
|
|
$sql=mysql_query("DELETE from automatic_all WHERE provider like 'JMA'");
|
|
$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_jma_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_jma_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_jma (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"),
|
|
'data'=>mysql_result($sql,$j,"data"),
|
|
'difference'=>mysql_result($sql,$j,"difference"),
|
|
'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_jma=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='JMA'");
|
|
$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='JMA'");
|
|
}
|
|
}
|
|
$sql_z=mysql_query("select max(ID) as max_id FROM delay_jma");
|
|
$max_id=mysql_result($sql_z,0,"max_id");
|
|
$sql_z=mysql_query("delete from delay_jma where ID<$max_id-40");
|
|
}
|
|
}
|
|
else{
|
|
print("Timeout");
|
|
}
|
|
|
|
mysql_close();
|
|
?>
|