Files
EQN/distquake_automatic_csi.php
T
2017-09-22 15:02:55 +02:00

359 lines
13 KiB
PHP

<?php
//CSI 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' => 35
)
)
);
$data = file_get_contents("http://www.csi.ac.cn/publish/main/1/100170/index.html",0,$ctx);
if ($data!=false)
{
//extract the maximum date
$sql=mysql_query("select date from last_notifications where provider='CSI'");
if (mysql_num_rows($sql)){
$max_date=mysql_result($sql,0,"date");
}
//empty the table
$sql=mysql_query("TRUNCATE automatic_csi_temp");
$data = str_replace("\n", '', $data); // remove new lines
$data = str_replace("\r", '', $data); // remove carriage returns
preg_match_all('|<tr onmouseover.*?</tr>|',$data,$result);
$result=$result[0];
date_default_timezone_set("Asia/Shanghai");
$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++){
preg_match_all('|<td.*?</td>|',$result[$i],$item);
$item=$item[0];
$idx=strpos($item[0],'/');
$url=substr($item[0],$idx,strlen($item[0])-$idx-32);
$url="http://www.csi.ac.cn".$url;
$date=substr($item[1],9,18);
$date=str_replace('/','-',$date);
$timestamp = strtotime($date);
$stamp2=strtotime($delta,$timestamp);
$stamp=date('Y-m-d H:i:s',$stamp2);
$lat=substr($item[2],9,strlen($item[2])-20);
$lon=substr($item[3],9,strlen($item[3])-20);
$idx1=strpos($item[4],'"');
$idx2=strrpos($item[4],'"');
$depth=substr($item[4],$idx1+1,$idx2-$idx1-1);
$magnitude=substr($item[5],6,strlen($item[5])-15);
$magnitude_type='ML';
$idx1=strpos($item[6],'self">');
$idx2=strrpos($item[6],'</a>');
$place=substr($item[6],$idx1+6,$idx2-$idx1-6);
$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_csi_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,'CSI',$weather_code,$temperature,$humidity,$pressure,$wind_speed,$wind_dir,$clouds,'$icon',$picture_count)");
}
else{
$sql=mysql_query("insert into automatic_csi_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,'CSI',$picture_count)");
}
}
$sql=mysql_query("select * from automatic_csi_temp");
if (mysql_num_rows($sql)){
//empty the table
$sql=mysql_query("TRUNCATE automatic_csi");
$sql=mysql_query("INSERT INTO automatic_csi SELECT * FROM automatic_csi_temp");
$sql=mysql_query("DELETE from automatic_all WHERE provider like 'CSI'");
$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_csi_temp");
//generate the file
$filename="seismic_free.txt";
$sql=mysql_query('SET CHARACTER SET utf8');
$sql=mysql_query("SELECT latitude,longitude,depth,magnitude,magitude_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_csi_temp where TIMESTAMPDIFF(MINUTE,'$max_date',data)>1 order by difference 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_csi (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"),
'pop100'=>mysql_result($sql,$j,"pop100"),
'difference'=>mysql_result($sql,$j,"difference"),
'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_csi=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='CSI'");
$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='CSI'");
}
}
$sql_z=mysql_query("select max(ID) as max_id FROM delay_csi");
$max_id=mysql_result($sql_z,0,"max_id");
$sql_z=mysql_query("delete from delay_csi where ID<$max_id-40");
}
}
}
else{
print("Timeout");
}
mysql_close();
?>