サンプル・プログラムの実行例
目次
サンプル・プログラム
| earthquake.php | サンプル・プログラム本体。 |
| pahooGeoCode.php | 住所・緯度・経度に関わるクラス pahooGeoCode。 使い方は「PHPで住所・ランドマークから最寄り駅を求める」などを参照。include_path が通ったディレクトリに配置すること。 |
| pahooCache.php | キャッシュ処理に関わるクラス pahooCache。 キャッシュ処理に関わるクラスの使い方は「PHPで天気予報を求める」を参照。include_path が通ったディレクトリに配置すること。 |
プログラムの方針
正規表現を使って、このページから
- 発生日時分
- 震源地
- 震源の位置
- 震源の深さ
- 規模
- 最大震度
気象庁防災情報XMLフォーマット
http://www.data.jma.go.jp/developer/xml/data/yyyymmddhhmmss_番号_電文コード_連番.xml
震源・震度に関する情報は、長期フィードの地震火山の中にある電文コード VXSE53 に入っている。そこで、フィードから配信日時 yyyymmddhhmmss が最も大きく、VXSE53 を含む URL を取り出せば、それが目指す情報 XML となる。
VXSE53の構造
ここから必要な情報を取り出す。
準備：pahooGeoCode クラス
0037: class pahooGeoCode {
0038: var $items; //検索結果格納用
0039: var $error; //エラーフラグ
0040: var $hits; //検索ヒット件数
0041: var $webapi; //直前に呼び出したWebAPI URL
0042:
0043: //Google Cloud Platform APIキー
0044: //https://cloud.google.com/maps-platform/
0045: //※Google Maps APIを利用しないのなら登録不要
0046: var $GOOGLE_API_KEY_1 = '**************************'; //HTTPリファラ用
0047: var $GOOGLE_API_KEY_2 = '**************************'; //IP制限用
0048:
0049: //Yahoo! JAPAN Webサービス アプリケーションID
0050: //https://e.developer.yahoo.co.jp/register
0051: //※Yahoo! JAPAN Webサービスを利用しないのなら登録不要
0052: var $YAHOO_APPLICATION_ID = '*****************************';
クラスについては「PHP でクラスを使ってテキストの読みやすさを調べる」を参照されたい。
地図として Google マップを利用するのであれば、Google Cloud Platform API キー が必要で、その入手方法は「Google Cloud Platform - WebAPI の登録方法」を参照されたい。
準備：地図サービスの選択
0039: //地図描画サービスの選択
0040: // 0：Google
0041: // 2：地理院地図・OSM
0042: define('MAPSERVICE', 2);
準備：キャッシュ・システム
0060: //キャッシュ保持時間（分） 0：キャッシュしない
0061: //気象庁へのアクセス負荷軽減のため，適切な時間を設定してください．
0062: define('LIFE_CACHE_FEED', 5); //高頻度 - 地震火山フィードに対して
0063: define('LIFE_CACHE_FEED_L', 120); //長期 - 地震火山フィードに対して
0064: define('LIFE_CACHE_DATA', 720); //地震情報に対して
0065:
0066: //キャッシュ・ディレクトリ
0067: //書き込み可能で，外部からアクセスされないディレクトリを指定してください．
0068: //天気予報系のプログラムとは別のディレクトリを指定してください．
0069: define('DIR_CACHE_FEED', './pcache1/');
0070: define('DIR_CACHE_FEED_L', './pcache2/');
0071: define('DIR_CACHE_DATA', './pcache3/');
キャッシュ保持時間、キャッシュ・ディレクトリともに、自サイトの環境に応じて変更してほしい。天気予報系プログラムと別のキャッシュ・ディレクトリにした方が、お互いのキャッシュ保持時間の干渉を受けなくなる。
配布ファイルは、新しい地震情報が入ってくることを考え、フィードの方を短く、地震情報の方は長くキャッシュ保持時間を設定してある。
準備：各種定数
0044: //マップの表示サイズ（単位：ピクセル）
0045: define('MAP_WIDTH', 600);
0046: define('MAP_HEIGHT', 400);
0047: //マップID
0048: define('MAPID', 'map_id');
0049: //初期値
0050: define('DEF_TYPE', 'HYBRID'); //マップタイプ
0051: define('DEF_ZOOM', 6); //ズーム
0052: define('INFO_WIDTH', (MAP_WIDTH * 0.75)); //情報ウィンドウの最大幅
0053:
0054: //Spinner - jQuery UI を使用するかどうか
0055: define('USESPINNER', TRUE);
0056:
0057: define('DEF_NUMBER', 1); //求めたい震源の数（初期値）
0058: define('MAX_NUMBER', 50); //求めたい震源の数（最大）
解説：気象庁防災情報XMLから最新の震源・震度に関する情報URLを取得
0249: /**
0250: * 気象庁防災情報XMLから震源・震度に関する情報URLを取得
0251: * @param int $mode 0：最新1件のみ取得，1：すべての情報URL取得
0252: * @param array $urls URL格納配列
0253: * @param string $errmsg エラーメッセージ格納用
0254: * @return bool TRUE：取得成功／FALSE：取得失敗
0255: */
0256: function jma_getLastEarthquakeURLs($mode, &$urls, &$errmsg) {
0257: //URLパターン
0258: $vxse53 = '/http\:\/\/www\.data\.jma\.go\.jp\/developer\/xml\/data\/([0-9\_]+)VXSE53\_[0-9]+\.xml/ui';
0259:
0260: //随時フィードの解析
0261: $cnt = 0;
0262: $pcc = new pahooCache(LIFE_CACHE_FEED, DIR_CACHE_FEED);
0263: $xml = $pcc->simplexml_load(FEED);
0264: //レスポンス・チェック
0265: if ($pcc->iserror() || !isset($xml->entry)) {
0266: $errmsg = '気象庁防災情報XMLにアクセスできません';
0267: return FALSE;
0268: }
0269: foreach ($xml->entry as $node) {
0270: //URLを取得
0271: if (preg_match($vxse53, $node->id, $arr) > 0) {
0272: $urls[$cnt] = $arr[0];
0273: $cnt++;
0274: }
0275: }
0276: $pcc = NULL;
0277:
0278: //長期フィードの解析
0279: $pcc = new pahooCache(LIFE_CACHE_FEED_L, DIR_CACHE_FEED_L);
0280: $xml = $pcc->simplexml_load(FEED_L);
0281: //レスポンス・チェック
0282: if ($pcc->iserror() || !isset($xml->entry)) {
0283: $errmsg = '気象庁防災情報XMLにアクセスできません';
0284: return FALSE;
0285: }
0286: foreach ($xml->entry as $node) {
0287: //URLを取得
0288: if (preg_match($vxse53, $node->id, $arr) > 0) {
0289: if (array_search($arr[0], $urls) === FALSE) {
0290: $urls[$cnt] = $arr[0];
0291: $cnt++;
0292: }
0293: }
0294: }
0295: $pcc = NULL;
0296:
0297: //エラー・チェック
0298: if ($cnt == 0) {
0299: $errmsg = '直近の地震情報はありません';
0300: return FALSE;
0301: }
0302:
0303: //URLを日時の新しい順にソート
0304: rsort($urls);
0305:
0306: return TRUE;
0307: }
VXSE53 を含む URL を配列 $urls に格納していき、最後に配列 $urls を大きい順にソートすることで、ひづけんの新しい順にソートしたことになる。
解説：地震情報の取り出し
0309: /**
0310: * 地震情報取得（気象庁防災情報XMLから）
0311: * @param object $pgc pahooGeoCodeオブジェクト
0312: * @param array $items 地震情報を格納する配列
0313: * @param string $urls 情報XMLのURLを格納する配列
0314: * @param string $errmsg エラーメッセージ格納用
0315: * @param int $count 取得件数（省略時=1）
0316: * @return bool TRUE：取得成功／FALSE：失敗
0317: */
0318: function get_earthquake($pgc, &$items, &$urls, &$errmsg, $count=1) {
0319: //名前空間
0320: define('JMX_EB', 'http://xml.kishou.go.jp/jmaxml1/elementBasis1/');
0321: //マッチングパターン
0322: $pat1 = '/([0-9]+)\-([0-9]+)\-([0-9]+)T([0-9]+)\:([0-9]+)/ui'; //年月日時分
0323: $pat2 = '/([\+\-][0-9\.]+)([\+\-][0-9\.]+)[\-\+]([0-9]+)/ui'; //緯度・経度・深さ
0324:
0325: //オブジェクト生成
0326: $pcc = new pahooCache(LIFE_CACHE_DATA, DIR_CACHE_DATA);
0327:
0328: //最新の震源・震度に関する情報URLを取得
0329: $urls = array();
0330: $mode = ($count == 1) ? 0 : 1;
0331: jma_getLastEarthquakeURLs($mode, $urls, $errmsg);
0332: if ($errmsg != '') return FALSE;
0333:
0334: foreach ($urls as $key=>$vxse53) {
0335: //取得件数が上限だったらループ脱出
0336: if ($key >= $count) break;
0337:
0338: //震源情報の取得
0339: $xml = $pcc->simplexml_load($vxse53);
0340: //レスポンス・チェック
0341: if ($pcc->iserror() || !isset($xml->Body->Earthquake)) {
0342: $errmsg = '気象庁防災情報XMLから最新の震源・震度情報を取得できません';
0343: return FALSE;
0344: }
0345:
0346: //地震発生日時の取得
0347: if ((preg_match($pat1, (string)$xml->Body->Earthquake->OriginTime, $arr) > 0) && isset($arr[5])) {
0348: $items[$key]['year'] = (int)$arr[1];
0349: $items[$key]['month'] = (int)$arr[2];
0350: $items[$key]['day'] = (int)$arr[3];
0351: $items[$key]['hour'] = (int)$arr[4];
0352: $items[$key]['minuite'] = (int)$arr[5];
0353: } else {
0354: $errmsg = '気象庁防災情報XMLから最新の震源・震度情報を取得できません';
0355: return FALSE;
0356: }
0357:
0358: //震源地の取得
0359: if (isset($xml->Body->Earthquake->Hypocenter->Area->Name)) {
0360: $items[$key]['location'] = $xml->Body->Earthquake->Hypocenter->Area->Name;
0361: } else {
0362: $items[$key]['location'] = '不明';
0363: }
0364: if ($items[$key]['location'] == '') {
0365: $items[$key]['location'] = '不明';
0366: }
0367: $node = $xml->Body->Earthquake->Hypocenter->Area->children(JMX_EB);
0368: if (preg_match($pat2, (string)$node->Coordinate, $arr)) {
0369: if (isset($arr[1]) && isset($arr[2])) {
0370: list($items[$key]['longitude'], $items[$key]['latitude']) = $pgc->tokyo_wgs84((float)$arr[2], (float)$arr[1]);
0371: } else {
0372: $items[$key]['latitude'] = $items[$key]['longitude'] = '不明';
0373: }
0374: if (isset($arr[3])) {
0375: $items[$key]['depth'] = (float)$arr[3];
0376: } else {
0377: $items[$key]['depth'] = '不明';
0378: }
0379: }
0380:
0381: //マグニチュードの取得
0382: $node = $xml->Body->Earthquake->children(JMX_EB);
0383: if (isset($node->Magnitude)) {
0384: $items[$key]['magnitude'] = (float)$node->Magnitude;
0385: } else {
0386: $items[$key]['magnitude'] = '不明';
0387: }
0388:
0389: //震度を取得
0390: if (isset($xml->Body->Intensity->Observation->MaxInt)) {
0391: $items[$key]['maxintensity'] = (int)$xml->Body->Intensity->Observation->MaxInt;
0392: } else {
0393: $items[$key]['maxintensity'] = '不明';
0394: }
0395: }
0396:
0397: //オブジェクト解法
0398: $pcc = NULL;
0399:
0400: return TRUE;
0401: }
0574: /**
0575: * 日本測地系を世界測地系に変換する
0576: * @param float $long 経度（日本測地系）
0577: * @param float $lat 緯度（日本測地系）
0578: * @return float array(経度,緯度)（世界測地系）
0579: */
0580: function tokyo_wgs84($long, $lat) {
0581: $glong = $long - $lat * 0.000046038 - $long * 0.000083043 + 0.010040;
0582: $glat = $lat - $lat * 0.00010695 + $long * 0.000017464 + 0.0046017;
0583: return array($glong, $glat);
0584: }
解説：地震情報をマッピング情報に変換
0414: /**
0415: * 地震情報をマッピング情報に変換
0416: * @param array $items 地震情報を格納した配列
0417: * @param array $points マッピング情報を格納する配列
0418: * @return int 変換したマッピング情報の件数
0419: */
0420: function info2points(&$items, &$points) {
0421: $i = 0;
0422: $j = 0;
0423: foreach ($items as $i=>$item) {
0424: //同じマッピング位置があるかどうか
0425: $flag = FALSE;
0426: for ($k = 0; $k < $j; $k++) {
0427: if (($items[$i]['latitude'] == $points[$k + 1]['latitude']) && ($items[$i]['longitude'] == $points[$k + 1]['longitude'])) {
0428: $flag = TRUE;
0429: $points[$k + 1]['description'] .= '<br /><hr />';
0430: break;
0431: }
0432: }
0433: //新規のマッピング位置
0434: if (! $flag) {
0435: $k = $j;
0436: $points[$k + 1]['latitude'] = $items[$i]['latitude'];
0437: $points[$k + 1]['longitude'] = $items[$i]['longitude'];
0438: $points[$k + 1]['title'] = '';
0439: $points[$k + 1]['description'] = '';
0440: $j++;
0441: }
0442: $items[$i]['id'] = num2alpha($k + 1);
0443: $dt = makeDateTime($items[$i]);
0444: $latitude = sprintf('%.1f', $items[$i]['latitude']);
0445: $longitude = sprintf('%.1f', $items[$i]['longitude']);
0446: if (is_numeric($items[$i]['depth'])) {
0447: if ($items[$i]['depth'] == 0) {
0448: $depth = 'ごく浅い';
0449: } else {
0450: $depth = sprintf('%dkm', $items[$i]['depth'] / 1000);
0451: }
0452: } else {
0453: $depth = $items[$i]['depth'];
0454: }
0455: $points[$k + 1]['description'] .=<<< EOT
0456: 発生日時：{$dt}<br />震源地：{$items[$i]['location']}<br />震源の位置：北緯 {$latitude}度，東経 {$longitude}度<br />震源の深さ：{$depth}<br />地震の規模：マグニチュード {$items[$i]['magnitude']}<br />最大震度：{$items[$i]['maxintensity']}
0457: EOT;
0458: //打ち切り条件
0459: if ($j >= 26) break;
0460: }
0461:
0462: return $j;
0463: }
そこで、ユーザー関数 info2points を使って、震源 1 つに対して 1 つの要素が対応する配列 [$items] に対し、同じ緯度・経度に対して 1 つの要素が対応する配列 [$points] に情報を変換してやる。
こうすることで、同一震源にマップ・アイコンを 1 つだけ立てて、情報ウィンドウに複数回の地震情報を羅列するようにできる。
解説：地図描画について
解説：表示とURLパラメータ
0691: //パラメータ
0692: $number = getParam('number', FALSE, DEF_NUMBER);
0693: $number = trim($number);
0694: if (preg_match('/^[0-9]+$/i', $number) > 0) {
0695: $errmsg = '';
0696: if (($number < 1) || ($number > MAX_NUMBER)) {
0697: $errmsg = sprintf('判定できる整数の範囲は，1から%dまでです', MAX_NUMBER);
0698: }
0699: } else {
0700: $errmsg = '数値は正の整数（自然数）を指定してください';
0701: }
earthquake.php?number=10
のようにすることで、number で指定した値が求めたい震源の数になる。
活用例
参考サイト
- WebAPI の登録方法：ぱふぅ家のホームページ
- PHP で住所・ランドマークから最寄り駅を求める：ぱふぅ家のホームページ
- PHP で緯度・経度から住所を求める：ぱふぅ家のホームページ
- 日本の直近の地震情報：みんなの知識 ちょっと便利帳
そこで今回は、気象庁防災情報 XMLから、直近の地震情報を取得する PHP スクリプトを作成してみることにする。
余震が多いことが分かりやすいように、求めたい震源の数を任意に指定できるように改良した。
（2021 年 4 月 17 日）27 件目以降の情報を表示しない不具合を修正．
（2021 年 4 月 8 日）キャッシュ・システム導入：pahooCache クラス．
（2021 年 3 月 27 日）求めたい震源の数を変更できるようにした。プログラム・タイトル変更。
（2021 年 3 月 3 日）震源の深さが「ごく浅い」ときにエラーが出る不具合を修正した。
（2021 年 2 月 27 日）気象庁サイト・リニューアルに対応。気象庁防災情報 XML から情報取得するよう変更した。