PHPで地図上に恵方を表示

「PHP で二十四節気・七十二候一覧を作成」では、毎年の節分の日を計算するプログラムを作成した。節分の日に恵方の方角を向いて恵方巻きを無言で食べると良いとされるが、この恵方は、年によって方角が変化する。そこで今回は、PHP を使い、指定した年の恵方の方角を地図上に矢印で描く PHP プログラムを作ってみることにする。

プログラムを実行する

1. サンプル・プログラムの実行例
2. サンプル・プログラムのダウンロード
3. 恵方とは
4. 解説：pahooCalendar クラス
5. 準備：pahooGeoCode クラス
6. 準備：地図サービス（WebAPI）の選択
7. 準備：住所検索サービス（WebAPI）の選択
8. 準備：各種パラメータ
9. 解説：恵方の方位角を求める
10. 解説：地図上に矢印を描く
11. 解説：方位角と距離を指定した地点の緯度・経度を求める
12. 解説：Google JavaScript マップ用スクリプト生成（直線描画）
13. 解説：Leaflet用スクリプト生成（直線描画）
14. 解説：Leaflet用スクリプト（直線描画）
15. 参考サイト

ダウンロード（PHP8対応）

圧縮ファイルの内容

eho.php	サンプル・プログラム本体
pahooGeoCode.php	住所・緯度・経度に関わるクラス pahooGeoCode。 使い方は「PHPで住所・ランドマークから最寄り駅を求める」などを参照。include_path が通ったディレクトリに配置すること。
pahooCalendar.php	暦計算クラス pahooCalendar。 暦計算クラスの使い方は「PHPで日出没・月出没・月齢・潮を計算」を参照。include_path が通ったディレクトリに配置すること。

恵方は、その年の十干 （じっかん） によって決まる。十干は、甲・乙・丙・丁・戊・己・庚・辛・壬・癸の 10 の要素が割り当てられている。つまり、年号を 10 で割った余りで、その年の十干を知ることができる。
恵方は 4種類しかなく、たとえば甲の年と己の年の恵方は同じである。
十干と恵方の対応を整理すると、下表のようになる。

十干	西暦の下1桁	恵方
甲	4	東北東微東	75度
乙	5	西南西微西	255度
丙	6	南南東微南	165度
丁	7	北北西微北	345度
戊	8	南南東微南	165度
己	9	東北東微東	75度
庚	0	西南西微西	255度
辛	1	南南東微南	165度
壬	2	北北西微北	345度
癸	3	南南東微南	165度

0076: //暦計算クラス：include_pathが通ったディレクトリに配置
0077: require_once('pahooCalendar.php');

0351: //pahooCalendarクラス
0352: $pcl = new pahooCalendar();
0353: $pcl->setLanguage('jp');

恵方の方位角の計算は、ユーザークラス "pahooCalendar" に分離している。
まず、クラスファイル "pahooCalendar.php" を  require_once  し、オブジェクトを生成する。

0073: //住所・緯度・経度に関わるクラス：include_pathが通ったディレクトリに配置
0074: require_once('pahooGeoCode.php');

0349: //pahooGeoCodeクラス
0350: $pgc = new pahooGeoCode();

地図サービスを利用するために、クラスファイル "pahooGeoCode.php" を使用する。組み込み関数  require_once  を使って読めるディレクトリに配置する。ディレクトリは、設定ファイル php.ini に記述されているオプション設定 include_path に設定しておく。
クラスについては「PHP でクラスを使ってテキストの読みやすさを調べる」を参照されたい。

Yahoo! JAPAN Web サービスを利用するには Yahoo! JAPAN Web サービス アプリケーション IDが必要で、その入手方法は「Yahoo!JAPAN デベロッパーネットワーク - WebAPI の登録方法」を参照されたい。
また、地図として Google マップを利用するのであれば、Google Cloud Platform API キー が必要で、その入手方法は「Google Cloud Platform - WebAPI の登録方法」を参照されたい。また、Google マップで経路探索を行う場合は、DirectionsService APIを有効にする必要がある。課金率が高いので留意されたい。

0034: //地図サービス（WebAPI）の選択
0035: //    0：Google
0036: //    2：地理院地図・OSM
0037: define('MAPSERVICE', 2);

表示する地図は、Google マップと地理院地図・ OSM から選べる。
あらかじめ、定数 MAPSERVICE に値を設定すること。

0039: //住所検索サービスの選択
0040: //    0：Google
0041: //    1：Yahoo!JAPAN
0042: //   11：HeartRails Geo API
0043: //   12：OSM Nominatim Search API
0044: define('GEOSERVICE', 0);

住所やワンドマークから緯度・軽度を検索する WebAPI を選ぶことができる。
あらかじめ、定数 GEOSERVICE に値を設定すること。

0053: //マップの表示サイズ（単位：ピクセル）
0054: define('MAP_WIDTH',  600);
0055: define('MAP_HEIGHT', 400);
0056: //マップID
0057: define('MAPID', 'map_id');
0058: 
0059: //初期値
0060: define('DEF_LONGITUDE', 139.766667); //地図中心（経度）
0061: define('DEF_LATITUDE',  35.681111);      //　　　　（緯度）
0062: define('DEF_QUERY',        '東京駅');     //検索クエリ
0063: define('DEF_TYPE',     'roadmap');     //マップタイプ
0064: define('DEF_ZOOM',     12);         //ズーム
0065: define('DEF_CAT',      'address');     //カテゴリ（Yahoo!）
0066: define('LINE_COLOR',   '#FF0000');     //恵方の色
0067: define('LINE_OPACITY', 0.6);         //恵方の透明度
0068: define('LINE_WEIGHT',  5);              //恵方の太さ
0069: define('LINE_LENGTH',  5000);           //恵方の全長（メートル）
0070: define('ARROW_ANGLE',  20);         //恵方：矢印の角度
0071: define('ARROW_LENGTH', 0.1);         //恵方：矢の長さ（全長に対する比率）

地図や矢印の表示レイアウトを決めるための各種パラメータである。変更不可としているもの以外は、任意の値を設定できる。

0740: /**
0741:  * 恵方の方位角を求める
0742:  * @param   int $year 西暦年
0743:  * @return  float 恵方の方位角（北を0度として時計回り）
0744: */
0745: function eho($year) {
0746:     static $table = array(
0747:         255,        //庚年（西南西微西）
0748:         165,        //辛年（南南東微南）
0749:         345,        //壬年（北北西微北）
0750:         165,        //癸年（南南東微南）
0751:          75,        //甲年（東北東微東）
0752:         255,        //乙年（西南西微西）
0753:         165,        //丙年（南南東微南）
0754:         345,        //丁年（北北西微北）
0755:         165         //戊年（南南東微南）
0756:     );
0757:     $i = $year % 10;
0758: 
0759:     return $table[$i];
0760: }

恵方の方位角を求めるには、前述の「恵方とは」で紹介した一覧を配列として用意しておき、西暦年を 10 で割った剰余と紐付けてやる。

0202: /**
0203:  * 地図上に矢印を描くスクリプトを生成
0204:  * @param   object $pgc       pahooGeoCodeクラス
0205:  * @param   float  $lat, $lng 開始点の緯度・経度
0206:  * @param   float  $angle     方位角（北を0度として時計回り）
0207:  * @param   int    $api   0：Google Maps JavaScript（省略時）
0208:  *                       11：地理院地図・OSM（Leaflet使用）
0209:  * @return  string JavaScript
0210: */
0211: function jsArrow($pgc, $latitude, $longitude, $angle, $api=0) {
0212:     //本線
0213:     list($lat1, $lng1) = $pgc->getPointAngle($longitude, $latitude, $angle, LINE_LENGTH);
0214:     $points[0]['latitude']  = $latitude;
0215:     $points[0]['longitude'] = $longitude;
0216:     $points[1]['latitude']  = $lat1;
0217:     $points[1]['longitude'] = $lng1;
0218:     $js = $pgc->jsLine($points, LINE_COLOR, LINE_OPACITY, LINE_WEIGHT, MAPSERVICE);
0219: 
0220:     //矢印
0221:     $a = $angle - 180 - ARROW_ANGLE;
0222:     if ($a < 0)     $a += 360;
0223:     list($lat2, $lng2) = $pgc->getPointAngle($lng1, $lat1, $a, LINE_LENGTH * ARROW_LENGTH);
0224:     $points[0]['latitude']  = $lat1;
0225:     $points[0]['longitude'] = $lng1;
0226:     $points[1]['latitude']  = $lat2;
0227:     $points[1]['longitude'] = $lng2;
0228:     $js .= $pgc->jsLine($points, LINE_COLOR, LINE_OPACITY, LINE_WEIGHT, MAPSERVICE);
0229: 
0230:     $a = $angle - 180 + ARROW_ANGLE;
0231:     if ($a < 0)     $a += 360;
0232:     list($lat3, $lng3) = $pgc->getPointAngle($lng1, $lat1, $a, LINE_LENGTH * ARROW_LENGTH);
0233:     $points[1]['latitude']  = $lat3;
0234:     $points[1]['longitude'] = $lng3;
0235:     $js .= $pgc->jsLine($points, LINE_COLOR, LINE_OPACITY, LINE_WEIGHT, MAPSERVICE);
0236: 
0237:     return $js;
0238: }

恵方の方角へ向かって、地図上に矢印を描く。
矢印を 3 本の直線に分解して描くことを考える。方位角を示す本線と、その終点にぶら下がる矢印（2 本の短い線分）である。

本線は、eho 関数で得られた方位角へ向かって直線を描けばいい。始点と終点は緯度・経度で指定する必要があるので、始点から、方位角と距離を指定して終点を算出するユーザー関数 getPointAngle を用意した（後述）。

次に矢印部分であるが、これは、終点から見て、角度 ＋ARROW_ANGLE 方向、－ARROW_ANGLE 方向へ短い線分を描けばいい。線分の長さは、本線との比率 ARROW_LENGTH で計算している。

0648: /**
0649:  * ある地点から方位角と距離を指定した地点の緯度・経度を求める
0650:  * @param   float $longitude 経度（世界測地系）
0651:  * @param   float $latitude  緯度（世界測地系）
0652:  * @param   float $angle     方位角（度）
0653:  * @param   float $distance  距離（メートル）
0654:  * @return  float array(緯度,経度)
0655: */
0656: function getPointAngle($longitude, $latitude, $angle, $distance) {
0657:     $rad = 6378137.0;         //地球の半径（メートル）
0658:     $e2  = 6.69437999019758E-03;
0659: 
0660:     $wt  = sqrt(1.0 - $e2 * pow(sin($latitude * pi() / 180.0), 2));
0661:     $mt  = $rad * (1.0 - $e2) / pow($wt, 3);
0662:     $dit = $distance * cos($angle * pi() / 180.0) / $mt;
0663:     $i   = $latitude * pi() / 180.0 + $dit / 2;
0664:     $w   = sqrt(1.0 - $e2 * pow(sin($i), 2));
0665: 
0666:     //緯度
0667:     $m   = $rad * (1 - $e2) / pow($w, 3);
0668:     $di  = $distance * cos($angle * pi() / 180) / $m;
0669:     $lat = $latitude + $di * 180 / pi();
0670: 
0671:     //経度
0672:     $n   = $rad / $w;
0673:     $dk  = $distance * sin($angle * pi() / 180) / ($n * cos($i));
0674:     $lng = $longitude + $dk * 180 / pi();
0675: 
0676:     return array($lat, $lng);
0677: }

ある地点から方位角と距離を指定した地点の緯度・経度を求めるユーザー関数が getPointAngle である。球面三角法の方程式を解くことで計算している。

1010: /**
1011:  * 直線描画スクリプト：Googleマップ用
1012:  * @param   array  $points  直線の座標配列
1013:  *                      [$n]['longitude'] 経度（世界測地系）
1014:  *                      [$n]['latitude']  緯度（世界測地系）
1015:  * @param   string $color   描画色（省略時=#FF0000）
1016:  * @param   float  $opacity 透明度（省略時=1）
1017:  * @param   int    $weight  線の太さ（省略時=1）
1018:  * @return  string JavaScript
1019: */
1020: function jsLine_Gmap($points, $color='#FF0000', $opacity=1.0, $weight=1) {
1021:     $ss = '';
1022:     $cnt = 0;
1023:     foreach ($points as $pt) {
1024:         if ($cnt > 0)   $ss .= ",\n";
1025: //      $ss .= "\t\tnew google.maps.LatLng({$pt['latitude']}, {$pt['longitude']})";
1026:         $ss .= "\t\t{ lat: {$pt['latitude']}, lng: {$pt['longitude']} }";
1027:         $cnt++;
1028:     }
1029: 
1030: $js =<<< EOT
1031:     var pt = [
1032: {$ss}
1033:     ];
1034:     var lines = new google.maps.Polyline({
1035:         map: map,
1036:         path: pt,
1037:         strokeColor: '{$color}',
1038:         strokeOpacity: {$opacity},
1039:         strokeWeight: {$weight}
1040:     });
1041:     lines.setMap(map);
1042: 
1043: EOT;
1044:     return $js;
1045: }

Google マップ上に直線（折れ線）を描くスクリプトを生成する。緯度・経度座標は配列で渡す。

1955: /**
1956:  * 直線描画スクリプト：Leaflet用
1957:  * @param   array  $points  直線の座標配列
1958:  *                      [$n]['longitude'] 経度（世界測地系）
1959:  *                      [$n]['latitude']  緯度（世界測地系）
1960:  * @param   string $color   描画色（省略時=#FF0000）
1961:  * @param   float  $opacity 透明度（省略時=1）
1962:  * @param   int    $weight  線の太さ（省略時=1）
1963:  * @return  string JavaScript
1964: */
1965: function jsLine_Leaflet($points, $color='#FF0000', $opacity='1.0', $weight=1) {
1966:     $ss = '';
1967:     $cnt = 0;
1968:     foreach ($points as $pt) {
1969:         if ($cnt > 0)   $ss .= ', ';
1970:         $ss .= "[{$pt['latitude']}, {$pt['longitude']}]";
1971:         $cnt++;
1972:     }
1973: 
1974: $js =<<< EOT
1975:     L.polyline(
1976:         [{$ss}], {
1977:         'color':  '{$color}',
1978:         'opacity': {$opacity},
1979:         'weight':  {$weight}
1980:     }).addTo(map);
1981: 
1982: EOT;
1983:     return $js;
1984: }

Google マップ上に直線（折れ線）を描くスクリプトを生成する。緯度・経度座標は配列で渡す。

2369: /**
2370:  * マップ上に直線を描くスクリプト生成
2371:  * @param   array  $points  直線の座標配列
2372:  *                      [$n]['longitude'] 経度（世界測地系）
2373:  *                      [$n]['latitude']  緯度（世界測地系）
2374:  * @param   string $color   描画色（省略時=#FF0000）
2375:  * @param   float  $opacity 透明度（省略時=1）
2376:  * @param   int    $weight  線の太さ（省略時=1）
2377:  * @param   int    $api   0：Google Maps JavaScript（省略時）
2378:  *                       11：地理院地図・OSM（Leaflet使用）
2379:  * @return  string JavaScript
2380: */
2381: function jsLine($points, $color='#FF0000', $opacity='1.0', $weight=1, $api=0) {
2382:     switch ($api) {
2383:     //Google Maps JavaScript
2384:     case 0;
2385:         $js = $this->jsLine_Gmap($points, $color, $opacity, $weight, $api);
2386:         break;
2387:     //地理院地図・OSM（Leaflet使用）
2388:     case 2:
2389:         $js = $this->jsLine_Leaflet($points, $color, $opacity, $weight, $api);
2390:         break;
2391:     }
2392: 
2393:     return $js;
2394: }

jsLine_Gmap または jsLine_Leaflet を呼び出し、マップ上に直線を描くスクリプトを生成する。緯度・経度座標は配列で渡す。
#jsLine_Leaflet:title=jsLine_Leaflet] で生成したスクリプトは、drawJSmap メソッドの引数 $call に渡してやる。

• 各種WebAPI の登録方法：ぱふぅ家のホームページ
• PHP で二十四節気・七十二候一覧を作成
• PHP で大圏航路を描く：ぱふぅ家のホームページ