サンプル・プログラムの実行例
目次
サンプル・プログラムのダウンロード
| tide.php | サンプル・プログラム本体。 |
| pahooCalendar.php | 暦計算クラス pahooCalendar、潮位計算クラス pahooTide。 include_pathが通ったディレクトリに配置すること。 |
| tideData.zip | 計算パラメータ。pahooTideが解凍せずに参照する圧縮ファイル。 |
| pahooGeoCode.php | 住所・緯度・経度に関わるクラス pahooGeoCode。 使い方は「PHPで住所・ランドマークから最寄り駅を求める」「PHPで住所・ランドマークから緯度・経度を求める」などを参照。include_path が通ったディレクトリに配置すること。 |
| moon/moon_00.png～moon_30.png | 月齢画像。 |
| makeTideDataFiles.php | tideData.zip を自動生成する。 サンプル・プログラム実行時には不要。 |
| table_C1C2.xlsx | 表-C.1, C.2の元データ。 サンプル・プログラム実行時には不要。 |
| dispC1C2table.php | 表-C.1, C.2をTABLE表示する。 サンプル・プログラム実行時には不要。 |
準備：外部クラスなど
Google マップの表示や住所検索には、ユーザークラス "pahooGeoCode" を利用する。クラスファイル "pahooGeoCode.php" を include_path が通ったディレクトリに配置すること。
また、ユーザークラス "pahooTide" が参照するデータファイル "tideData.zip" はクラスファイルと同じディレクトリに配置すること。
準備：pahooGeoCode クラス
0035: class pahooGeoCode {
0036: var $items; //検索結果格納用
0037: var $error; //エラーフラグ
0038: var $hits; //検索ヒット件数
0039: var $webapi; //直前に呼び出したWebAPI URL
0040:
0041: //Google Cloud Platform APIキー
0042: //https://cloud.google.com/maps-platform/
0043: //※Google Maps APIを利用しないのなら登録不要
0044: var $GOOGLE_API_KEY_1 = '**************************'; //HTTPリファラ用
0045: var $GOOGLE_API_KEY_2 = '**************************'; //IP制限用
0046:
0047: //Yahoo! JAPAN Webサービス アプリケーションID
0048: //https://e.developer.yahoo.co.jp/register
0049: //※Yahoo! JAPAN Webサービスを利用しないのなら登録不要
0050: var $YAHOO_APPLICATION_ID = '*****************************';
クラスについては「PHP でクラスを使ってテキストの読みやすさを調べる」を参照されたい。
地図として Google マップを利用するのであれば、Google Cloud Platform API キー が必要で、その入手方法は「Google Cloud Platform - WebAPI の登録方法」を、YOLP マップを利用するのであれば、Yahoo! JAPAN Web サービス アプリケーション IDが必要で、その入手方法は「Google Cloud Platform - WebAPI の登録方法」を、それぞれ参照されたい。
準備：地図サービス（WebAPI）の選択
0046: //地図描画サービスの選択
0047: // 0：Google
0048: // 1：Yahoo!JAPAN
0049: // 2：地理院地図・OSM
0050: define('MAPSERVICE', 0);
0051:
0052: //住所検索サービスの選択
0053: // 0：Google
0054: // 1：Yahoo!JAPAN
0055: // 11：HeartRails Geo API
0056: define('GEOSERVICE', 0);
0057:
0058: //逆ジオコーディングサービスの選択
0059: // 0：Google
0060: // 1：Yahoo!JAPAN
0061: // 11：HeartRails Geo API
0062: // 21：簡易ジオコーディングサービス
0063: define('REVGEOSERVICE', 21);
あらかじめ、定数 GEOSERVICE に値を設定すること。
全体の流れ
潮汐は複数の分潮の合成であり、下記の式で求めることができる。
‥‥(C.1)
| η | 潮位 |
| Z0 | 観測値の平均水面 |
| fi | 分潮iの振幅に対する補正係数 |
| Hi | 観測地において計算された振幅 |
| Vi | 時間とともに変化する分潮引数 |
| ui | 時間とともに変化する位相の補正係数 |
| ti | 観測地において計算された位相遅れ |
、 
、 
は観測地に固有の値であり、後述するが、気象庁が公開している各観測地点における分潮一覧表から入手することができる。
1631: // 潮位計算 ===============================================================
1632: // @参考URL https://www.pahoo.org/e-soul/webtech/php02/php02-51-01.shtm
1633: class pahooTide {
1634: const FILE_ZIPNAME = 'tideData.zip'; //計算用パラメータ格納ZIPファイル
1635: const FILE_C1 = '_C1.txt'; //表-C.1
1636: const FILE_C2 = '_C2.txt'; //表-C.2
1637: const FILE_INDEX = '_index.txt'; //地点一覧
1638: const FILE_EXT = '.txt'; //拡張子
1639: const JST = -9; //日本時の時差
1640: const COMMENT = '#'; //コメント文字
1641: var $error, $errmsg; //エラーフラグ，エラーメッセージ
1642: var $_s, $_h, $_p, $_N; //天文引数
1643:
1644: var $c1, $c2; //表-C.1, C.2
1645: var $index; //地点一覧
1646: var $port; //ある地点の分潮一覧表
1647:
1648: //分潮記号
1649: var $keys = array('Sa', 'Ssa', 'Mm', 'MSf', 'Mf', '2Q1', 'σ1', 'Q1', 'ρ1', 'O1', 'MP1', 'M1', 'χ1', 'π1', 'P1', 'S1', 'K1', 'ψ1', 'φ1', 'θ1', 'J1', 'SO1', 'OO1', 'OQ2', 'MNS2', '2N2', 'μ2', 'N2', 'ν2', 'OP2', 'M2', 'MKS2', 'λ2', 'L2', 'T2', 'S2', 'R2', 'K2', 'MSN2', 'KJ2', '2SM2', 'MO3', 'M3', 'SO3', 'MK3', 'SK3', 'MN4', 'M4', 'SN4', 'MS4', 'MK4', 'S4', 'SK4', '2MN6', 'M6', 'MSN6', '2MS6', '2MK6', '2SM6', 'MSK6');
1650:
1651: /**
1652: * コンストラクタ
1653: * @param なし
1654: * @return なし
1655: */
1656: function __construct() {
1657: $this->error = FALSE;
1658: $this->errmsg = '';
1659:
1660: //表-C.1, C.2, 地点一覧を配列に格納
1661: $zipfname = __DIR__ . '/'. self::FILE_ZIPNAME;
1662: $zip = new ZipArchive;
1663: if ($zip->open($zipfname) == TRUE) {
1664: $this->str2array($zip->getFromName(self::FILE_C1), $this->c1);
1665: $this->str2array($zip->getFromName(self::FILE_C2), $this->c2);
1666: $this->str2array($zip->getFromName(self::FILE_INDEX), $this->index);
1667: $zip->close();
1668: } else {
1669: $this->error = TRUE;
1670: $this->errmsg = "cannot read \"{$zipfname}\"";
1671: }
1672: $zip = NULL;
1673: }
1674:
1675: /**
1676: * デストラクタ
1677: * @return なし
1678: */
1679: function __destruct() {
1680: unset($this->items);
1681: }
1682:
潮汐と分潮
1 つ 1 つの周期を分潮と呼び、60種類が定義されている。このうち影響の大きい下記の 4 つを、主要4 分潮を呼ぶ。
| M2 | 月の引力によるもので、周期は約12時間25分。 |
| S2 | 太陽の引力によるもので、周期は約12時間。 |
| O1 | 月の引力によるもので、周期は約25時間49分。 |
| K1 | 観月と太陽の合成引力によるもので、周期は約23時間56分。 |
天文引数の計算
‥‥(C.2)
‥‥(C.3)
‥‥(C.4)
‥‥(C.5)
‥‥(C.6)
‥‥(C.7)
| Y | 西暦 |
| D | その年の1月1日からの経過日数 |
| L | 西暦2000年からその年の年初までの閏日の数 |
{ }は中の数値を整数に丸めることを意味する。L は 2000 年以前では負数とする。
1734: /**
1735: * 天文引数を求める
1736: * @param int $year, $month, $day グレゴリオ暦による年月日
1737: * @return なし
1738: */
1739: function sun_moon($year, $month, $day) {
1740: $y = (double)$year - 2000;
1741: $L = (double)floor(($year + 3) / 4) - 500;
1742: $d = (double)date('z', mktime(0, 0, 0, $month, $day, $year)) + $L;
1743:
1744: $this->_s = (double)211.728 + 129.38471 * $y + 13.176396 * $d;
1745: $this->_h = (double)279.974 - 0.23871 * $y + 0.985647 * $d;
1746: $this->_p = (double) 83.298 + 40.66229 * $y + 0.111404 * $d;
1747: $this->_N = (double)125.071 - 19.32812 * $y - 0.052954 * $d;
1748:
1749: $this->_s = $this->_stdegree($this->_s);
1750: $this->_h = $this->_stdegree($this->_h);
1751: $this->_p = $this->_stdegree($this->_p);
1752: $this->_N = $this->_stdegree($this->_N);
1753: }
分潮引数 Vi の計算
は次のように表せる。
‥‥(C.9)
| V0i | 世界時0時における分潮引数[tex:V_i] |
| σi | 分潮iの角速度で、後述する表-C.2から取得する。 |
| n | 1日周期は1、半日周期は2、1/4周期は4‥‥となる定数で、後述する表-C.2から取得する。 |
| tl | ローカル時間（時）。 |
| ts | 世界時との時差（時）。日本標準時の場合は -9。 |
| Lat | 観測地点の経度。東経はプラス、西経はマイナス。 |
基準となる
は次のようにして計算できる。
‥‥(C.10)
、 
、 
、 
は、後述する表-C.2 から取得できる。
補正係数 fi、ui の計算
、 
の計算方法は 3種類に分かれる。
1）L2、M1分潮
‥‥(C.17)
‥‥(C.18)
‥‥(C.19)
‥‥(C.20)
‥‥(C.21)
‥‥(C.22)
1755: /**
1756: * L2, M1分潮のfi, uiの計算
1757: * @param string $i 分潮記号 'L2' または 'M1'
1758: * @param string $fn 'fi' または 'ui'
1759: * @return double 計算結果
1760: */
1761: function fiui_L2M1($i, $fn) {
1762: //L2
1763: if ($i == 'L2') {
1764: $fcosu = (double)1 - 0.2505 * cos(deg2rad(2 * $this->_p))
1765: - 0.1102 * cos(deg2rad(2 * $this->_p - $this->_N))
1766: - 0.0156 * cos(deg2rad(2 * $this->_p - 2 * $this->_N))
1767: - 0.0370 * cos(deg2rad($this->_N));
1768: $fsinu = (double)-0.2505 * sin(deg2rad(2 * $this->_p))
1769: - 0.1102 * sin(deg2rad(2 * $this->_p - $this->_N))
1770: - 0.0156 * sin(deg2rad(2 * $this->_p - 2 * $this->_N))
1771: - 0.0370 * sin(deg2rad($this->_N));
1772: //M1
1773: } else {
1774: $fcosu = (double)2 * cos(deg2rad($this->_p))
1775: + 0.4 * cos(deg2rad($this->_p - $this->_N));
1776: $fsinu = (double)sin(deg2rad($this->_p))
1777: + 0.2 * sin(deg2rad($this->_p - $this->_N));
1778: }
1779:
1780: //fi
1781: if ($fn == 'fi') {
1782: $res = (double)sqrt($fcosu * $fcosu + $fsinu * $fsinu);
1783: //ui
1784: } else {
1785: $res = (double)atan($fcosu / $fsinu);
1786: }
1787:
1788: return $res;
1789: }
‥‥(C.11)
‥‥(C.12)
3）1，2 以外の分潮
表-C.2 の列に記載した計算式によって求める。
1791: /**
1792: * fi, uiの計算
1793: * @param string $i 分潮記号
1794: * @param string $fn 'fi' または 'ui'
1795: * @return double 計算結果
1796: */
1797: function fiui($i, $fn) {
1798: //L2 or M1
1799: if (($i == 'L2') || ($i == 'M1')) {
1800: $res = (double)$this->fiui_L2M1($i, $fn);
1801:
1802: //fi
1803: } else if ($fn == 'fi') {
1804: //表-C.1による計算
1805: if (isset($this->c1[$i])) {
1806: $res = (double)$this->c1[$i][0]
1807: + (double)$this->c1[$i][1] * cos(deg2rad($this->_N))
1808: + (double)$this->c1[$i][2] * cos(2 * deg2rad($this->_N))
1809: + (double)$this->c1[$i][3] * cos(3 * deg2rad($this->_N));
1810: //表-C.2の参照
1811: } else {
1812: $res = $this->c2[$i][5];
1813: if (! is_numeric($res)) $res = (double)$this->fiui($res, $fn);
1814: }
1815: //ui
1816: } else {
1817: //表-C.1による計算
1818: if (isset($this->c1[$i])) {
1819: $res = (double)$this->c1[$i][4] * sin(deg2rad($this->_N))
1820: + (double)$this->c1[$i][5] * sin(2 * deg2rad($this->_N))
1821: + (double)$this->c1[$i][6] * sin(3 * deg2rad($this->_N));
1822: //表-C.2の参照
1823: } else {
1824: $res = (double)$this->c2[$i][6];
1825: if (! is_numeric($res)) $res = (double)$this->fiui($res, $fn);
1826: }
1827: }
1828: return $res;
1829: }
表-C.1、表-C.2
1700: /**
1701: * 表から配列へパラメータを格納する
1702: * @param string $str パラメータ表
1703: * @param array $arr パラメータを格納する配列
1704: * @return なし
1705: */
1706: function str2array($str, &$arr) {
1707: $tok = strtok($str, "\n");
1708: while ($tok != FALSE) {
1709: $ss = trim($tok);
1710: if ($ss == '') continue;
1711: $cols = preg_split("/\t/iu", $ss);
1712: $key = trim($cols[0]);
1713: if (mb_substr($key, 0, 1) != self::COMMENT) {
1714: for ($i = 1; $i < count($cols); $i++) {
1715: if (mb_substr($cols[$i], 0, 1) == self::COMMENT) break; //コメントから行末まで無視
1716: $arr[$key][$i - 1] = trim($cols[$i]);
1717: }
1718: }
1719: $tok = strtok("\n");
1720: }
1721: }
| 分潮記号 | fi | ui | |||||
|---|---|---|---|---|---|---|---|
| 1 | cos N | cos 2N | cos 3N | sin N | sin 2N | sin 3N | |
| Mm | 1.0000 | -0.1300 | 0.0013 | 0.0000 | 0.00 | 0.00 | 0.00 |
| Mf | 1.0429 | 0.4135 | -0.0040 | 0.0000 | -23.74 | 2.68 | -0.38 |
| O1 | 1.0089 | 0.1871 | -0.0147 | 0.0014 | 10.80 | -1.34 | 0.19 |
| K1 | 1.0060 | 0.1150 | -0.0088 | 0.0006 | -8.86 | 0.68 | -0.07 |
| J1 | 1.0129 | 0.1676 | -0.0170 | 0.0016 | -12.94 | 1.34 | -0.19 |
| OO1 | 1.1027 | 0.6504 | 0.0317 | -0.0014 | -36.68 | 4.02 | -0.57 |
| M2 | 1.0004 | -0.0373 | 0.0002 | 0.0000 | -2.14 | 0.00 | 0.00 |
| K2 | 1.0241 | 0.2863 | 0.0083 | -0.0015 | -17.74 | 0.68 | -0.04 |
| 分潮記号 | s | h | p | c | n | fi | ui | σ (deg/hour) | period (hour) |
|---|---|---|---|---|---|---|---|---|---|
| Sa | 0 | 1 | 0 | 0 | 0 | 1 | 0 | 0.0410686 | 8765.8211 |
| Ssa | 0 | 2 | 0 | 0 | 0 | 1 | 0 | 0.0821373 | 4382.9052 |
| Mm | 1 | 0 | -1 | 0 | 0 | Mm | Mm | 0.5443747 | 661.3092 |
| MSf | 2 | -2 | 0 | 0 | 0 | M2 | -M2 | 1.0158958 | 354.3671 |
| Mf | 2 | 0 | 0 | 0 | 0 | Mf | Mf | 1.0980331 | 327.8590 |
| 2Q1 | -4 | 1 | 2 | 270 | 1 | O1 | O1 | 12.8542862 | 28.0062 |
| σ1 | -4 | 3 | 0 | 270 | 1 | O1 | O1 | 12.9271398 | 27.8484 |
| Q1 | -3 | 1 | 1 | 270 | 1 | O1 | O1 | 12.3986609 | 26.8684 |
| ρ1 | -3 | 3 | -1 | 270 | 1 | O1 | O1 | 13.4715145 | 26.7231 |
| O1 | -2 | 1 | 0 | 270 | 1 | O1 | O1 | 13.9430356 | 25.8193 |
| MP1 | -2 | 3 | 0 | 90 | 1 | M2 | M2 | 14.0251729 | 25.6681 |
| M1 | -1 | 1 | 0 | 90 | 1 | M1 | M1 | 14.4920521 | 24.8412 |
| χ1 | -1 | 3 | -1 | 90 | 1 | J1 | J1 | 14.5695476 | 24.7091 |
| π1 | 0 | -2 | 0 | 193 | 1 | 1 | 0 | 14.9178647 | 24.1321 |
| P1 | 0 | -1 | 0 | 270 | 1 | 1 | 0 | 14.9589314 | 24.0659 |
| S1 | 0 | 0 | 0 | 180 | 1 | 1 | 0 | 15.0000000 | 24.0000 |
| K1 | 0 | 1 | 0 | 90 | 1 | K1 | K1 | 15.0410686 | 23.9345 |
| ψ1 | 0 | 2 | 0 | 167 | 1 | 1 | 0 | 15.0821353 | 23.8693 |
| φ1 | 0 | 3 | 0 | 90 | 1 | 1 | 0 | 15.1232059 | 23.8045 |
| θ1 | 1 | -1 | 1 | 90 | 1 | J1 | J1 | 15.5125897 | 23.2070 |
| J1 | 1 | 1 | -1 | 90 | 1 | J1 | J1 | 15.5854433 | 23.0985 |
| SO1 | 2 | -1 | 0 | 90 | 1 | J1 | J1 | 16.0569644 | 22.4202 |
| OO1 | 2 | 1 | 0 | 90 | 1 | OO1 | OO1 | 16.1391017 | 22.3061 |
| OQ2 | -5 | 2 | 1 | 180 | 2 | O1*Q1 | O1+Q1 | 27.3416964 | 13.1667 |
| MNS2 | -5 | 4 | 1 | 0 | 2 | M2^2 | 2*M2 | 27.4238337 | 13.1273 |
| 2N2 | -4 | 2 | 2 | 0 | 2 | M2 | M2 | 27.8953548 | 12.9054 |
| μ2 | -4 | 4 | 0 | 0 | 2 | M2 | M2 | 27.9682084 | 12.8718 |
| N2 | -3 | 2 | 1 | 0 | 2 | M2 | M2 | 28.4397295 | 12.6583 |
| ν2 | -3 | 4 | -1 | 0 | 2 | M2 | M2 | 28.5125831 | 12.6260 |
| OP2 | -2 | 0 | 0 | 180 | 2 | O1*P1 | O1+P1 | 28.9019669 | 12.4559 |
| M2 | -2 | 2 | 0 | 0 | 2 | M2 | M2 | 28.9841042 | 12.4206 |
| MKS2 | -2 | 4 | 0 | 0 | 2 | M2*K2 | M2+K2 | 29.0662415 | 12.3855 |
| λ2 | -1 | 0 | 1 | 180 | 2 | M2 | M2 | 29.4556253 | 12.2218 |
| L2 | -1 | 2 | -1 | 180 | 2 | L2 | L2 | 29.5284789 | 12.1916 |
| T2 | 0 | -1 | 0 | 283 | 2 | 1 | 0 | 29.9589333 | 12.0164 |
| S2 | 0 | 0 | 0 | 0 | 2 | 1 | 0 | 30.0000000 | 12.0000 |
| R2 | 0 | 1 | 0 | 257 | 2 | 1 | 0 | 30.0410667 | 11.9836 |
| K2 | 0 | 2 | 0 | 0 | 2 | K2 | K2 | 30.0821373 | 11.9672 |
| MSN2 | 1 | 0 | -1 | 0 | 2 | M2^2 | 2*M2 | 30.5443747 | 11.7861 |
| KJ2 | 1 | 2 | -1 | 180 | 2 | K1*J1 | K1+J1 | 30.6265120 | 11.7545 |
| 2SM2 | 2 | -2 | 0 | 0 | 2 | M2 | -M2 | 31.0158958 | 11.6070 |
| MO3 | -4 | 3 | 0 | 270 | 3 | M2*O1 | M2+O1 | 42.9271398 | 8.3863 |
| M3 | -3 | 3 | 0 | 180 | 3 | M2^1.5 | 1.5*M2 | 43.4761563 | 8.2804 |
| SO3 | -2 | 1 | 0 | 270 | 3 | O1 | O1 | 43.9430356 | 8.1924 |
| MK3 | -2 | 3 | 0 | 90 | 3 | M2*K1 | M2+K1 | 44.0251729 | 8.1771 |
| SK3 | 0 | 1 | 0 | 90 | 3 | K1 | K1 | 45.0410686 | 7.9927 |
| MN4 | -5 | 4 | 1 | 0 | 4 | M2^2 | 2*M2 | 57.4238337 | 6.2692 |
| M4 | -4 | 4 | 0 | 0 | 4 | M2^2 | 2*M2 | 57.9682084 | 6.2103 |
| SN4 | -3 | 2 | 1 | 0 | 4 | M2 | M2 | 58.4397295 | 6.1602 |
| MS4 | -2 | 2 | 0 | 0 | 4 | M2 | M2 | 58.9841042 | 6.1033 |
| MK4 | -2 | 4 | 0 | 0 | 4 | M2*K2 | M2+K2 | 59.0662415 | 6.0949 |
| S4 | 0 | 0 | 0 | 0 | 4 | 1 | 0 | 60.0000000 | 6.0000 |
| SK4 | 0 | 2 | 0 | 0 | 4 | K2 | K2 | 60.0821373 | 5.9918 |
| 2MN6 | -7 | 6 | 1 | 0 | 6 | M2^3 | 3*M2 | 86.4079380 | 4.1663 |
| M6 | -6 | 6 | 0 | 0 | 6 | M2^3 | 3*M2 | 86.9523127 | 4.1402 |
| MSN6 | -5 | 4 | 1 | 0 | 6 | M2^2 | 2*M2 | 87.4238337 | 4.1179 |
| 2MS6 | -4 | 4 | 0 | 0 | 6 | M2^2 | 2*M2 | 87.9682084 | 4.0924 |
| 2MK6 | -4 | 6 | 0 | 0 | 6 | M2^2*K2 | 2*M2+K2 | 88.0503457 | 4.0886 |
| 2SM6 | -2 | 2 | 0 | 0 | 6 | M2 | M2 | 88.9841042 | 4.0457 |
| MSK6 | -2 | 4 | 0 | 0 | 6 | M2*K2 | M2+K2 | 89.0662415 | 4.0419 |
Zi、Hi、ui の入手
、 
、 
は観測地点の固有の値である。これらの値は、気象庁の潮位表掲載地点一覧表から取得することができる。対象は、一覧表のうち、分潮一覧表の列のリンク先にコンテンツがあるもののみである。
これらのコンテンツを自動的に読み込んでゆき、各地点の分潮一覧表を ZIP 圧縮ファイルに格納していくプログラムが "makeTideDataFiles.php" である。
ある日の干潮、満潮時刻と潮位を求める
干潮、満潮は、1 日に各々最大 4 回起きる。1 分ごとに潮位を計算し、上昇から下降へ転じる点が満潮時刻、その逆が干潮時刻である。だが、1 分ごとに計算すると、1 日で 60×24＝1440 回も計算しなければならない。
計算回数を減らすために、まず 15 分ごとに計算することで変化点を検出し、そこから±15 分間を 1 分ごとに計算することで正確な時刻を求める。計算回数は、1 日で 4×24＋15×8
＝216 回に減らすことができる。
1948: /**
1949: * 指定地点のある日の干潮、満潮時刻と潮位を求める
1950: * @param string $code 地点記号
1951: * @param int $year, $month, $day グレゴリオ暦による年月日
1952: * @param array $items 結果を格納する配列
1953: * array('high'=>array('hh:mm',潮位), 'low'=>array('hh:mm',潮位))
1954: * @return TRUE/FALSE
1955: */
1956: function tide_day($code, $year, $month, $day, &$items) {
1957: $this->setLocation($code);
1958: if ($this->error) return FALSE;
1959:
1960: $this->sun_moon($year, $month, $day);
1961:
1962: $interval = 15; //第1段階刻み（分）
1963: $td0 = $this->tide_level($year, $month, $day, self::JST, -$interval);
1964: $flag = 0;
1965: $cnt_high = 0;
1966: $cnt_low = 0;
1967: for ($i = 0; $i <= 24 * 60; $i += $interval) {
1968: //第1段階刻み
1969: $hh = floor($i / 60);
1970: $mm = $i - $hh * 60;;
1971: $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
1972: if ($flag == 0) {
1973: if ($td1 > $td0) $flag = +1; //上昇
1974: else $flag = -1; //下降
1975: } else if ($flag < 0) {
1976: if ($td1 > $td0) { //上昇へ転じた
1977: //第2段階刻み
1978: $td0 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm - $interval - 1);
1979: for ($j = -$interval; $j < $interval; $j++) {
1980: $hh = floor(($i + $j) / 60);
1981: $mm = ($i + $j) - $hh * 60;;
1982: $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
1983: if ($td1 > $td0) { //上昇へ転じた
1984: $items['low'][$cnt_low]['hhmm'] = sprintf("%02d:%02d", $hh, $mm);
1985: $items['low'][$cnt_low]['lev'] = (int)$td0;
1986: $cnt_low++;
1987: $flag = +1;
1988: break;
1989: }
1990: $td0 = $td1;
1991: }
1992: }
1993: } else {
1994: if ($td1 < $td0) { //下降へ転じた
1995: //第2段階刻み
1996: $td0 = $this->tide_level($year, $month, $day, self::JST, $mm - $interval - 1);
1997: for ($j = -$interval; $j < $interval; $j++) {
1998: $hh = floor(($i + $j) / 60);
1999: $mm = ($i + $j) - $hh * 60;;
2000: $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
2001: if ($td1 < $td0) { //下降へ転じた
2002: $items['high'][$cnt_high]['hhmm'] = sprintf("%02d:%02d", $hh, $mm);
2003: $items['high'][$cnt_high]['lev'] = (int)$td0;
2004: $cnt_high++;
2005: $flag = -1;
2006: break;
2007: }
2008: $td0 = $td1;
2009: }
2010: }
2011: }
2012: $td0 = $td1;
2013: }
2014: return TRUE;
2015: }
解説：各種定数
0017: //jqPlotのあるフォルダ
0018: define('JQPLOT', '../../../../common/jqplot/');
0019:
0020: //月の満ち欠け画像ファイルの場所
0021: // 画像作成方法→http://www.pahoo.org/e-soul/webtech/phpgd/phpgd-17-01.shtm
0022: define('MOONAGE', './moon/');
0023:
0024: //プログラム・タイトル
0025: define('TITLE', '潮位・月齢を計算');
0026:
0027: //リリース・フラグ（公開時にはTRUEにすること）
0028: define('FLAG_RELEASE', TRUE);
0029:
0030: //リファラ・チェック（直リン防止用；空文字ならチェックしない）
0031: define('REFER_ON', 'www.pahoo.org');
0032: //define('REFER_ON', '');
0033:
0034: //潮位表：計算期間（日）の初期値
0035: define('INTERVAL_DEF', 7);
0036:
0037: //潮位グラフ：計算期間（日）の初期値
0038: define('INTERVAL_GRAPH_DEF', 5);
0039:
0040: //潮位グラフのプロット間隔（日）
0041: define('INTERVAL_GRAPH', 0.02);
0042:
0043: //潮位グラフの名前
0044: define('GRAPH_TIDE', 'jqPlot_tide');
ただし、冒頭に述べたとおり計算量が大きいので、INTERVAL_DEF や INTERVAL_GRAPH_DEF は、あまり大きな値にしない方がいいだろう。
潮位グラフを描くために、jQuery プラグイン「jqPlot」を利用する。使い方については、「PHP で NHK 政治意識月例調査をグラフ表示」を参照のこと。定数 JQPLOT の示すディレクトリに配置する。
また、月齢を示す月の画像は、定数 MOONAGE の示すディレクトリに配置する。圧縮ファイルに含まれている画像は、「PHP で月の満ち欠けを描画」で作成したモノだ。
解説：地点セレクタ
0358: /**
0359: * 地点セレクタを作成する
0360: * @param pahooTide $pt 潮位計算クラス
0361: * @param string $code 地点記号
0362: * @param string $errmsg エラーメッセージ格納用；エラーなければ空文字
0363: * @return string HTML／FALSE：失敗
0364: */
0365: function makeJSlocation($pt, $code) {
0366: static $prefs = array(
0367: '--選択--','北海道', '青森県', '岩手県', '宮城県', '秋田県', '山形県', '福島県',
0368: '茨城県', '栃木県', '群馬県', '埼玉県', '千葉県', '東京都', '神奈川県',
0369: '新潟県', '富山県', '石川県', '福井県', '山梨県', '長野県', '岐阜県', '静岡県',
0370: '愛知県', '三重県', '滋賀県', '京都府', '大阪府', '兵庫県', '奈良県',
0371: '和歌山県', '鳥取県', '島根県', '岡山県', '広島県', '山口県', '徳島県',
0372: '香川県', '愛媛県', '高知県', '福岡県', '佐賀県', '長崎県', '熊本県', '大分県',
0373: '宮崎県', '鹿児島県', '沖縄県');
0374:
0375: if ($code == '') {
0376: $pp = '';
0377: } else if (! isset($pt->index[$code])) {
0378: $errmsg = "not exist location '{$code}'";
0379: return FALSE;
0380: } else {
0381: $pp = $pt->index[$code][4];
0382: }
0383:
0384: //都道府県
0385: $js = "var pref = [\n";
0386: foreach ($prefs as $key=>$val) {
0387: $selected = ($val == $pp) ? 'selected' : '';
0388: $js .= "{'type':{$key}, 'text':'{$val}', 'value':'{$val}', 'selected':'{$selected}'},\n";
0389: }
0390: //地点
0391: $js .= "];\nvar loc = [\n";
0392: foreach ($prefs as $key=>$val) {
0393: $js .= "{'type':{$key}, 'text':'--選択--', 'value':''},\n";
0394: foreach ($pt->index as $cd=>$arr) {
0395: if ($arr[4] == $val) {
0396: $selected = ($cd == $code) ? 'selected' : '';
0397: $js .= "{'type':{$key}, 'text':'{$arr[0]}', 'value':'{$cd}', 'selected':'{$selected}'},\n";
0398: }
0399: }
0400: }
0401:
0402: //スクリプト
0403: $js .=<<< EOD
0404: ];
0405:
0406: $(function() {
0407: var html = '';
0408: for (var i in pref) {
0409: html += '<option value="'+ pref[i].value + '" ' + pref[i].selected + '>'+ pref[i].text +'</option>';
0410: }
0411: $('#pref').append(html);
0412:
0413: var html = '';
0414: var index = $('#pref option:selected').index();
0415: for (var i in loc) {
0416: if (loc[i].type == index) {
0417: html += '<option value="'+ loc[i].value + '" ' + loc[i].selected + '>' + loc[i].text +'</option>';
0418: }
0419: }
0420: $('#code').append(html);
0421:
0422: //都道府県を選択した時の挙動
0423: $('#pref').change(function() {
0424: var html = '';
0425: var index = $('#pref option:selected').index();
0426: //地点の内容を削除
0427: $('#code').empty();
0428: //対応する地点の内容を表示
0429: for (var i in loc) {
0430: if (loc[i].type == index) {
0431: html += '<option value="'+ loc[i].value +'">'+ loc[i].text +'</option>';
0432: }
0433: }
0434: $('#code').append(html);
0435: });
0436: });
0437:
0438: EOD;
0439: return $js;
0440: }
観測地点セレクタを都道府県セレクタに連動させるために、jQuery を利用した。
解説：潮位表の計算と作成
0443: /**
0444: * 潮位表：潮位・月齢を計算する
0445: * @param pahooTide $pt 潮位計算クラス
0446: * @param pahooCalendar $pc 暦計算クラス
0447: * @param array $inputs 計算用パラメータ
0448: * @param array $items 計算結果格納用配列
0449: * @param array $locs 地点情報格納用配列
0450: * @param string $errmsg エラーメッセージ格納用；エラーなければ空文字
0451: * @return bool TRUE：計算成功／FALSE：失敗
0452: */
0453: function calcTideTable($pt, $pc, $params, &$items, &$locs, &$errmsg) {
0454: //地点を設定
0455: $pt->setLocation($params['code']);
0456: if ($pt->iserror()) return FALSE;
0457:
0458: //地点を取得
0459: $pt->getLocation($params['code'], $locs);
0460: if ($pt->iserror()) return FALSE;
0461:
0462: //月齢・潮位を計算
0463: $res = TRUE;
0464: $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], 0, 0, 0);
0465: for ($i = 0; $i < $params['interval']; $i++) {
0466: $arr = array();
0467: list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
0468: $yb = $pc->getWeekString($year, $month, $day);
0469: $items[$i]['dt'] = sprintf("%04d-%02d-%02d(%s)", $year, $month, $day, $yb);
0470: $items[$i]['moonage'] = $pc->moon_age($year, $month, $day, MOONAGE_HOUR, 0, 0);
0471: $moonmeridian = $pc->moon_time(2, $params['longitude'], $params['latitude'], $params['height'], $year, $month, $day);
0472: $moonmeridian = ($moonmeridian == FALSE) ? '---' : $pc->day2hhmm($moonmeridian);
0473: $items[$i]['tide'] = (preg_match('/([0-9]+)\:([0-9]+)/', $moonmeridian, $arr) > 0) ? $pc->tide($year, $month, $day, $arr[1], $arr[2], 0) : '';
0474: unset($arr);
0475:
0476: $arr = array();
0477: $pt->tide_day($params['code'], $year, $month, $day, $arr);
0478: //満潮
0479: for ($j = 0; $j < 4; $j++) {
0480: if (isset($arr['high'][$j])) {
0481: $items[$i]['high'][$j]['hhmm'] = $arr['high'][$j]['hhmm'];
0482: $items[$i]['high'][$j]['lev'] = sprintf("%d", $arr['high'][$j]['lev']);
0483: $items[$i]['high'][$j]['align'] = 'text-align:right;';
0484: } else {
0485: $items[$i]['high'][$j]['hhmm'] = '*';
0486: $items[$i]['high'][$j]['lev'] = '*';
0487: $items[$i]['high'][$j]['align'] = 'text-align:center;';
0488: }
0489: }
0490: //干潮
0491: for ($j = 0; $j < 4; $j++) {
0492: if (isset($arr['low'][$j])) {
0493: $items[$i]['low'][$j]['hhmm'] = $arr['low'][$j]['hhmm'];
0494: $items[$i]['low'][$j]['lev'] = sprintf("%d", $arr['low'][$j]['lev']);
0495: $items[$i]['low'][$j]['align'] = 'text-align:right;';
0496: } else {
0497: $items[$i]['low'][$j]['hhmm'] = '*';
0498: $items[$i]['low'][$j]['lev'] = '*';
0499: $items[$i]['low'][$j]['align'] = 'text-align:center;';
0500: }
0501: }
0502: $jd++;
0503: }
0504:
0505: return $res;
0506: }
月齢・潮の計算は、「PHP で日出没・月出没・月齢・潮を計算」で紹介したとおりである。
0508: /**
0509: * 潮位表：表示用HTMLを作成
0510: * @param pahooTide $pt 潮位計算クラス
0511: * @param pahooCalendar $pc 暦計算クラス
0512: * @param array $params 計算用パラメータ
0513: * @param string $js スクリプト格納用
0514: * @param string $html HTML格納用
0515: * @param string $errmsg エラーメッセージ格納用
0516: * @return bool TRUE：成功／FALSE：失敗
0517: * @return 表示用HTML
0518: */
0519: function makeTable($pt, $pc, $params, &$js, &$html, &$errmsg) {
0520: $js = $html = '';
0521: $items = array();
0522: $locs = array();
0523: if (calcTideTable($pt, $pc, $params, $items, $locs, $errmsg) == FALSE) return FALSE;
0524:
0525: $hour = MOONAGE_HOUR;
0526: $html =<<< EOD
0527: <table class="tide">
0528: <caption>{$locs['title']}（{$locs['prefecture']}{$locs['address']}）</caption>
0529: <tr><th rowspan="2">年月日</th><th rowspan="2">月齢<br /><span style="font-size:small; font-weight:normal;">{$hour}時</span></th><th rowspan="2">潮</th><th colspan="8">満潮</th><th colspan="8">干潮</th></tr>
0530:
0531: EOD;
0532: $html .= "<tr>";
0533: for ($j = 0; $j < 8; $j++) {
0534: $html .= "<th class=\"index\">時刻</th><th class=\"index\">潮位<br />(cm)</th>";
0535: }
0536: $html .= "</tr>\n";
0537: foreach ($items as $item) {
0538: $html .= sprintf("<tr><td>%s</td><td>%.1f</td><td>%s</td>", $item['dt'], $item['moonage'], $item['tide']);
0539: for ($j = 0; $j < 4; $j++) {
0540: $html .= sprintf("<td>%s</td><td style=\"%s\">%s</td>", $item['high'][$j]['hhmm'], $item['high'][$j]['align'], $item['high'][$j]['lev']);
0541: }
0542: for ($j = 0; $j < 4; $j++) {
0543: $html .= sprintf("<td>%s</td><td style=\"%s\">%s</td>", $item['low'][$j]['hhmm'], $item['low'][$j]['align'], $item['low'][$j]['lev']);
0544: }
0545: $html .= "</tr>\n";
0546: }
0547: $html .= "</table>\n";
0548:
0549: return TRUE;
0550: }
解説：潮位グラフの計算と作成
0552: /**
0553: * 潮位グラフ：潮位・月齢を計算する
0554: * @param pahooTide $pt 潮位計算クラス
0555: * @param pahooCalendar $pc 暦計算クラス
0556: * @param array $params 計算用パラメータ
0557: * @param array $items 潮位計算結果格納用配列
0558: * @param array $moons 月齢計算結果格納用配列
0559: * @param array $locs 地点情報格納用配列
0560: * @param string $errmsg エラーメッセージ格納用；エラーなければ空文字
0561: * @return bool TRUE：計算成功／FALSE：失敗
0562: */
0563: function calcTideGraph($pt, $pc, $params, &$items, &$moons, &$locs, &$errmsg) {
0564: //地点を設定
0565: $pt->setLocation($params['code']);
0566: if ($pt->iserror()) return FALSE;
0567:
0568: //地点を取得
0569: $pt->getLocation($params['code'], $locs);
0570: if ($pt->iserror()) return FALSE;
0571:
0572: $res = TRUE;
0573:
0574: //潮位を計算
0575: $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], -JST, 0, 0);
0576: $n = $params['interval'] / INTERVAL_GRAPH;
0577: for ($i = 0; $i < $n; $i++) {
0578: list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd + JST / 24);
0579: $items[$i]['dt'] = sprintf("%04d/%02d/%02d %02d:%02d", $year, $month, $day, $hour, $min);
0580: list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
0581: $pt->sun_moon($year, $month, $day);
0582: $items[$i]['tlevel'] = $pt->tide_level($year, $month, $day, $hour, $min);
0583: $jd += INTERVAL_GRAPH;
0584: }
0585:
0586: //月齢を計算
0587: $height = 0; //標高0メートルを仮定
0588: $n = $params['interval'];
0589: $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], 0, 0, 0);
0590: $i = $d = 0;
0591: while ($i < $n) {
0592: //月の南中時刻
0593: list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
0594: $moonmeridian = $pc->moon_time(2, $params['longitude'], $params['latitude'], $height, $year, $month, $day);
0595: if ($moonmeridian != FALSE) {
0596: $x = (double)$i + $moonmeridian;
0597: $moons[$i]['t'] = $d + $moonmeridian;
0598: $moons[$i]['age'] = $pc->moon_age($year, $month, $day, $moonmeridian * 24, 0, 0);
0599: $i++;
0600: }
0601: $jd++;
0602: $d++;
0603: }
0604:
0605: return $res;
0606: }
月齢は、月の画像をグラフ上にオーバーライトすることを目指す。プロットする場所は、月の南中時刻とすべく、メソッド pahooCalendar::moon_time によって南中時刻を求める。
0608: /**
0609: * 潮位グラフ：表示用スクリプトおよびHTMLを作成
0610: * @param pahooTide $pt 潮位計算クラス
0611: * @param pahooCalendar $pc 暦計算クラス
0612: * @param array $params 計算用パラメータ
0613: * @param string $js スクリプト格納用
0614: * @param string $html HTML格納用
0615: * @param string $errmsg エラーメッセージ格納用
0616: * @return bool TRUE：成功／FALSE：失敗
0617: */
0618: function makeGraph($pt, $pc, $params, &$js, &$html, &$errmsg) {
0619: $js = $html = '';
0620: $items = array();
0621: $moons = array();
0622: $locs = array();
0623: if (calcTideGraph($pt, $pc, $params, $items, $moons, $locs, $errmsg) == FALSE) return FALSE;
0624:
0625: $title = "{$locs['title']}（{$locs['prefecture']}{$locs['address']}）";
0626: $name = GRAPH_TIDE;
0627: $width = MAP_WIDTH;
0628: $height = MAP_HEIGHT;
0629:
0630: //潮位プロット・データ作成
0631: $data = '';
0632: foreach ($items as $item) {
0633: $data .= sprintf("['%s', %f], ", $item['dt'], $item['tlevel']);
0634: }
0635:
0636: //月齢データ作成
0637: $moon = '';
0638: $path = MOONAGE;
0639: foreach ($moons as $val) {
0640: $x = round(30 + (MAP_WIDTH - 70) / $params['interval'] * $val['t']);
0641: $fname = sprintf("{$path}moon_%02d.png", round($val['age']));
0642: $moon .=<<< EOD
0643: <img style="position:absolute; top:35px; left:{$x}px; width:40px; height:40px; z-index:999;" src="{$fname}" />
0644:
0645: EOD;
0646: }
0647:
0648: $js =<<< EOD
0649: $(function() {
0650: jQuery.jqplot('{$name}',
0651: [
0652: [ {$data} ]
0653: ],
0654: {
0655: //タイトル
0656: title: {
0657: text: '{$title}',
0658: show: true,
0659: fontSize: '16px',
0660: textAlign: 'center',
0661: textColor: 'black'
0662: },
0663: //背景
0664: grid: {
0665: background: '#EEFFFF'
0666: },
0667: //グラフ
0668: seriesDefaults: {
0669: showLine: true,
0670: rendererOptions: { smooth: false },
0671: markerOptions: { size: 0 },
0672: color: 'blue',
0673: },
0674: //軸ラベル
0675: axes: {
0676: xaxis: {
0677: renderer: $.jqplot.DateAxisRenderer,
0678: tickOptions: { formatString: '%m/%d' },
0679: tickInterval: '1 days'
0680: },
0681: yaxis: {
0682: label: '潮位(cm)'
0683: }
0684: }
0685: }
0686: );
0687: });
0688:
0689: EOD;
0690: $html =<<< EOD
0691: <div id="{$name}" style="width:{$width}px; height:{$height}px;">
0692: {$moon}
0693: </div>
0694:
0695: EOD;
0696:
0697: return TRUE;
0698: }
jqPlot に渡すスクリプトと、表示する HTML は別々に生成する。jqPlot で画像を扱えないため、無理をして HTML 側でグラフ上にオーバーライトしている。
解説：メインプログラム
Google マップの表示や住所検索については「PHP で Google等を利用して住所から緯度・経度を求める」を、UI：Datepicker を利用した日付入力については「PHP で日付入力：カレンダーから選択」を参照してほしい。
参考サイト
- 各種WebAPI の登録方法：ぱふぅ家のホームページ
- PHP で日出没・月出没・月齢・潮を計算：ぱふぅ家のホームページ
- PHP で NHK 政治意識月例調査をグラフ表示：ぱふぅ家のホームページ
- PHP で月の満ち欠けを描画：ぱふぅ家のホームページ
- 潮位表掲載地点一覧表：気象庁
さて、「PHP で日出没・月出没・月齢・潮を計算」では潮の計算を行うプログラムを紹介したが、潮位を計算するには、はるかに複雑な計算を行う必要がある。幸い、コンピュータの計算能力が飛躍的に向上したおかげで、潮位の計算もストレスなく行えるようになった。
そこで、これまで作ってきた "pahooCalendar.php" に、あらた潮位計算クラス "pahooTide" を追加し、Google マップに表示されている地点の近くの海岸において、指定日から 1週間の毎日の満潮・干潮の時刻と潮位を一覧およびグラフ表示するサンプル・プログラムを作ってみることにする。
（2019 年 6 月 24 日）地理院地図、OpenStreetMap も利用できるようにした。緯度・経度から住所を検索する逆ジオコーディングサービスの選択肢を増やした。地図検索，リセット・ボタンを追加した。
（2019 年 4 月 20 日）Google マップと Yahoo!マップを切り替えて使えるようにした。