PHPで潮位を計算する

プログラムを実行する

1. サンプル・プログラムの実行例
2. サンプル・プログラムのダウンロード
3. 準備：外部クラスなど
4. 全体の流れ
5. 潮汐と分潮
6. 天文引数の計算
7. 補正係数 f_i、u_i の計算
8. 表-C.1、表-C.2
9. Z_i、H_i、u_i の入手
10. ある日の干潮、満潮時刻と潮位を求める
11. 解説：各種定数
12. 解説：地点セレクタ
13. 解説：潮位表の計算と作成
14. 解説：潮位グラフの計算と作成
15. 解説：メインプログラム
16. 参考サイト

ダウンロード（PHP8対応）

  43: //jqPlotのあるフォルダ
  44: define('JQPLOT', '../../../../common/jqplot/');
  45: 
  46: //月の満ち欠け画像ファイルの場所
  47: // 画像作成方法→https://www.pahoo.org/e-soul/webtech/phpgd/phpgd-17-01.shtm
  48: define('MOONAGE', './moon/');
  49: 
  50: //潮位表：計算期間（日）の初期値
  51: define('INTERVAL_DEF', 7);
  52: 
  53: //潮位グラフ：計算期間（日）の初期値
  54: define('INTERVAL_GRAPH_DEF', 5);
  55: 
  56: //潮位グラフのプロット間隔（日）
  57: define('INTERVAL_GRAPH', 0.02);
  58: 
  59: //潮位グラフの名前
  60: define('GRAPH_TIDE', 'jqPlot_tide');
  61: 
  62: //地図描画サービスの選択
  63: //    0：Google
  64: //    2：地理院地図・OSM
  65: define('MAPSERVICE', 0);
  66: 
  67: //住所検索サービスの選択
  68: //    0：Google
  69: //    1：Yahoo!JAPAN
  70: //   11：HeartRails Geo API
  71: define('GEOSERVICE', 1);
  72: 
  73: //逆ジオコーディングサービスの選択
  74: //    0：Google
  75: //    1：Yahoo!JAPAN
  76: //   11：HeartRails Geo API
  77: //   21：簡易ジオコーディングサービス
  78: define('REVGEOSERVICE', 21);
  79: 
  80: //指定できる西暦年の範囲
  81: define('MIN_YEAR', 1901);
  82: define('MAX_YEAR', 2099);
  83: 
  84: //表示言語（jp:日本語, en:英語, en3:英語略記）
  85: define('LANGUAGE', 'jp');
  86: 
  87: //世界時からの時差（日本標準時）
  88: define('UTCDIFF', +9.0);
  89: 
  90: //観測点のデフォルト値（東京駅）
  91: define('DEF_LOCATION', 'TK');   //東京
  92: define('DEF_LATITUDE',  35.681111);
  93: define('DEF_LONGITUDE', 139.766667);
  94: define('DEF_HEIGHT',    4.0);
  95: 
  96: //月齢計算時刻
  97: define('MOONAGE_HOUR', 21);
  98: 
  99: //マップID
 100: define('MAPID', 'map_id');
 101: //初期値
 102: define('MAP_WIDTH',     800);           //表示幅
 103: define('MAP_HEIGHT',    300);           //表示高
 104: define('DEF_TYPE',      'roadmap');     //デフォルトのマップタイプ
 105: define('DEF_ZOOM',      8);             //デフォルトの拡大率
 106: define('DEF_CAT',       'address');     //カテゴリ
 107: 
 108: //require_once()で呼ぶファイルはinclude_pathが通っているフォルダに配置すること．
 109: //暦計算クラス
 110: require_once('pahooCalendar.php');
 111: 
 112: //住所・緯度・経度に関わるクラス
 113: require_once('pahooGeoCode.php');
 114: 
 115: //データ入力に関わる関数群
 116: require_once('pahooInputData.php');

  37: class pahooGeoCode {
  38:     var $items;     //検索結果格納用
  39:     var $error;     //エラー・フラグ
  40:     var $errmsg;    //エラー・メッセージ
  41:     var $hits;      //検索ヒット件数
  42:     var $webapi;    //直前に呼び出したWebAPI URL
  43: 
  44:     //Google Cloud Platform APIキー
  45:     //https://cloud.google.com/maps-platform/
  46:     //※Google Maps APIを利用しないのなら登録不要
  47:     var $GOOGLE_API_KEY_1 = '**************************';   //HTTPリファラ用
  48:     var $GOOGLE_API_KEY_2 = '**************************';   //IP制限用
  49: 
  50:     //Yahoo! JAPAN Webサービス アプリケーションID
  51:     //https://e.developer.yahoo.co.jp/register
  52:     //※Yahoo! JAPAN Webサービスを利用しないのなら登録不要
  53:     var $YAHOO_APPLICATION_ID = '*****************************';

  48: //地図描画サービスの選択
  49: //    0：Google
  50: //    2：地理院地図・OSM
  51: define('MAPSERVICE', 2);
  52: 
  53: //住所検索サービスの選択
  54: //    0：Google
  55: //    1：Yahoo!JAPAN
  56: //   11：HeartRails Geo API
  57: //   12：OSM Nominatim Search API
  58: //   13：国土地理院ジオコーディングAPI
  59: define('GEOSERVICE', 12);
  60: 
  61: //逆ジオコーディングサービスの選択
  62: //    0：Google
  63: //    1：Yahoo!JAPAN
  64: //   11：HeartRails Geo API
  65: //   21：簡易ジオコーディングサービス
  66: define('REVGEOSERVICE', 21);

2644: // 潮位計算 ===============================================================
2645: // @参考URL https://www.pahoo.org/e-soul/webtech/php02/php02-51-01.shtm
2646: class pahooTide {
2647:     const FILE_ZIPNAME = 'tideData.zip';    //計算用パラメータ格納ZIPファイル
2648:     const FILE_C1      = '_C1.txt';         //表-C.1
2649:     const FILE_C2      = '_C2.txt';         //表-C.2
2650:     const FILE_INDEX   = '_index.txt';      //地点一覧
2651:     const FILE_EXT     = '.txt';            //拡張子
2652:     const JST          = -9;                //日本時の時差
2653:     const COMMENT = '#';                    //コメント文字
2654:     var $error, $errmsg;                    //エラーフラグ，エラーメッセージ
2655:     var $_s, $_h, $_p, $_N;                 //天文引数
2656: 
2657:     var $c1, $c2;                           //表-C.1, C.2
2658:     var $index;                             //地点一覧
2659:     var $port;                              //ある地点の分潮一覧表
2660: 
2661: //分潮記号
2662: 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');

2747: /**
2748:  * 指定した年月日から天文引数を求める．
2749:  * @param   int $year, $month, $day  グレゴリオ暦による年月日
2750:  * @return  なし
2751: */
2752: function sun_moon($year, $month, $day) {
2753:     $y = (float)$year - 2000.0;
2754:     $L = (float)floor(($year + 3.0) / 4.0) - 500.0;
2755:     $d = (float)date('z', mktime(0, 0, 0, $month, $day, $year)) + $L;
2756: 
2757:     $this->_s = (float)211.728 + 129.38471 * $y + 13.176396 * $d;
2758:     $this->_h = (float)279.974 -   0.23871 * $y +  0.985647 * $d;
2759:     $this->_p = (float) 83.298 +  40.66229 * $y +  0.111404 * $d;
2760:     $this->_N = (float)125.071 -  19.32812 * $y -  0.052954 * $d;
2761: 
2762:     $this->_s = $this->_stdegree($this->_s);
2763:     $this->_h = $this->_stdegree($this->_h);
2764:     $this->_p = $this->_stdegree($this->_p);
2765:     $this->_N = $this->_stdegree($this->_N);
2766: }

2768: /**
2769:  * L2, M1分潮のfi, uiを求める．
2770:  * @param   string $i  分潮記号 'L2' または 'M1'
2771:  * @param   string $fn 'fi' または 'ui'
2772:  * @return  float 計算結果
2773: */
2774: function fiui_L2M1($i, $fn) {
2775:     //L2
2776:     if ($i == 'L2') {
2777:         $fcosu = (float)1 - 0.2505 * cos(deg2rad(2 * $this->_p))
2778:                 - 0.1102 * cos(deg2rad(2 * $this->_p - $this->_N))
2779:                 - 0.0156 * cos(deg2rad(2 * $this->_p - 2 * $this->_N))
2780:                 - 0.0370 * cos(deg2rad($this->_N));
2781:         $fsinu = (float)-0.2505 * sin(deg2rad(2 * $this->_p))
2782:                 - 0.1102 * sin(deg2rad(2 * $this->_p - $this->_N))
2783:                 - 0.0156 * sin(deg2rad(2 * $this->_p - 2 * $this->_N))
2784:                 - 0.0370 * sin(deg2rad($this->_N));
2785:     //M1
2786:     } else {
2787:         $fcosu = (float)2 * cos(deg2rad($this->_p))
2788:                 + 0.4 * cos(deg2rad($this->_p - $this->_N));
2789:         $fsinu = (float)sin(deg2rad($this->_p))
2790:                  + 0.2 * sin(deg2rad($this->_p - $this->_N));
2791:     }
2792: 
2793:     //fi
2794:     if ($fn == 'fi') {
2795:         $res = (float)sqrt($fcosu * $fcosu + $fsinu * $fsinu);
2796:     //ui
2797:     } else {
2798:         $res = (float)atan($fcosu / $fsinu);
2799:     }
2800: 
2801:     return $res;
2802: }

2804: /**
2805:  * fi, uiを求める．
2806:  * @param   string $i  分潮記号
2807:  * @param   string $fn 'fi' または 'ui'
2808:  * @return  float 計算結果
2809: */
2810: function fiui($i, $fn) {
2811:     //L2 or M1
2812:     if (($i == 'L2') || ($i == 'M1')) {
2813:         $res = (float)$this->fiui_L2M1($i, $fn);
2814: 
2815:     //fi
2816:     } else if ($fn == 'fi') {
2817:         //表-C.1による計算
2818:         if (isset($this->c1[$i])) {
2819:             $res = (float)$this->c1[$i][0]
2820:             + (float)$this->c1[$i][1] * cos(deg2rad($this->_N))
2821:             + (float)$this->c1[$i][2] * cos(2 * deg2rad($this->_N))
2822:             + (float)$this->c1[$i][3] * cos(3 * deg2rad($this->_N));
2823:         //表-C.2の参照
2824:         } else {
2825:             $res = $this->c2[$i][5];
2826:             if (! is_numeric($res)) $res = (float)$this->fiui($res, $fn);
2827:         }
2828:     //ui
2829:     } else {
2830:         //表-C.1による計算
2831:         if (isset($this->c1[$i])) {
2832:             $res = (float)$this->c1[$i][4] * sin(deg2rad($this->_N))
2833:             + (float)$this->c1[$i][5] * sin(2 * deg2rad($this->_N))
2834:             + (float)$this->c1[$i][6] * sin(3 * deg2rad($this->_N));
2835:         //表-C.2の参照
2836:         } else {
2837:             $res = (float)$this->c2[$i][6];
2838:             if (! is_numeric($res)) $res = (float)$this->fiui($res, $fn);
2839:         }
2840:     }
2841:     return $res;
2842: }

2713: /**
2714:  * 表から配列へパラメータを格納する．
2715:  * @param   string $str パラメータ表
2716:  * @param   array  $arr パラメータを格納する配列
2717:  * @return  なし
2718: */
2719: function str2array($str, &$arr) {
2720:     $tok = strtok($str, "\n");
2721:     while ($tok != FALSE) {
2722:         $ss = trim($tok);
2723:         if ($ss == "")      continue;
2724:         $cols = preg_split("/\t/iu", $ss);
2725:         $key = trim($cols[0]);
2726:         if (mb_substr($key, 0, 1) != self::COMMENT) {
2727:             for ($i = 1; $i < count($cols); $i++) {
2728:                 if (mb_substr($cols[$i], 0, 1) == self::COMMENT)    break;      //コメントから行末まで無視
2729:                 $arr[$key][$i - 1] = trim($cols[$i]);
2730:             }
2731:         }
2732:         $tok = strtok("\n");
2733:     }
2734: }

2961: /**
2962:  * 指定地点のある日の干潮，満潮時刻と潮位を求める．
2963:  * @param   string $code 地点記号
2964:  * @param   int $year, $month, $day グレゴリオ暦による年月日
2965:  * @param   array $items 結果を格納する配列
2966:  *              array('high'=>array('hh:mm',潮位), 'low'=>array('hh:mm',潮位))
2967:  * @return  TRUE/FALSE
2968: */
2969: function tide_day($code, $year, $month, $day, &$items) {
2970:     $this->setLocation($code);
2971:     if ($this->error)   return FALSE;
2972: 
2973:     $this->sun_moon($year, $month, $day);
2974: 
2975:     $interval = 15.0;           //第1段階刻み（分）
2976:     $td0 = $this->tide_level($year, $month, $day, self::JST, -$interval);
2977:     $flag = 0;
2978:     $cnt_high = 0;
2979:     $cnt_low = 0;
2980:     for ($i = 0; $i <= 24 * 60; $i += $interval) {
2981:         //第1段階刻み
2982:         $hh = floor($i / 60.0);
2983:         $mm = $i - $hh * 60.0;
2984:         $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
2985:         if ($flag == 0) {
2986:             if ($td1 > $td0)    $flag = +1;     //上昇
2987:             else                $flag = -1;     //下降
2988:         } else if ($flag < 0) {
2989:             if ($td1 > $td0) {                      //上昇へ転じた
2990:                 //第2段階刻み
2991:                 $td0 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm - $interval - 1);
2992:                 for ($j = -$interval; $j < $interval; $j++) {
2993:                     $hh = floor(($i + $j) / 60.0);
2994:                     $mm = ($i + $j) - $hh * 60.0;
2995:                     $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
2996:                     if ($td1 > $td0) {              //上昇へ転じた
2997:                         $items['low'][$cnt_low]['hhmm'] = sprintf("%02d:%02d", $hh, $mm);
2998:                         $items['low'][$cnt_low]['lev'] = (int)$td0;
2999:                         $cnt_low++;
3000:                         $flag = +1;
3001:                         break;
3002:                     }
3003:                     $td0 = $td1;
3004:                 }
3005:             }
3006:         } else {
3007:             if ($td1 < $td0) {                       //下降へ転じた
3008:                 //第2段階刻み
3009:                 $td0 = $this->tide_level($year, $month, $day, self::JST, $mm - $interval - 1);
3010:                 for ($j = -$interval; $j < $interval; $j++) {
3011:                     $hh = floor(($i + $j) / 60.0);
3012:                     $mm = ($i + $j) - $hh * 60.0;
3013:                     $td1 = $this->tide_level($year, $month, $day, $hh + self::JST, $mm);
3014:                     if ($td1 < $td0) {               //下降へ転じた
3015:                         $items['high'][$cnt_high]['hhmm'] = sprintf("%02d:%02d", $hh, $mm);
3016:                         $items['high'][$cnt_high]['lev'] = (int)$td0;
3017:                         $cnt_high++;
3018:                         $flag = -1;
3019:                         break;
3020:                     }
3021:                     $td0 = $td1;
3022:                 }
3023:             }
3024:         }
3025:         $td0 = $td1;
3026:     }
3027:     return TRUE;
3028: }

  43: //jqPlotのあるフォルダ
  44: define('JQPLOT', '../../../../common/jqplot/');
  45: 
  46: //月の満ち欠け画像ファイルの場所
  47: // 画像作成方法→https://www.pahoo.org/e-soul/webtech/phpgd/phpgd-17-01.shtm
  48: define('MOONAGE', './moon/');
  49: 
  50: //潮位表：計算期間（日）の初期値
  51: define('INTERVAL_DEF', 7);
  52: 
  53: //潮位グラフ：計算期間（日）の初期値
  54: define('INTERVAL_GRAPH_DEF', 5);
  55: 
  56: //潮位グラフのプロット間隔（日）
  57: define('INTERVAL_GRAPH', 0.02);
  58: 
  59: //潮位グラフの名前
  60: define('GRAPH_TIDE', 'jqPlot_tide');

 362: /**
 363:  * 地点セレクタを作成する．
 364:  * @param   pahooTide $pt  潮位計算クラス
 365:  * @param   string $code 地点記号
 366:  * @param   string $errmsg エラーメッセージ格納用；エラーなければ空文字
 367:  * @return  string HTML／FALSE：失敗
 368: */
 369: function makeJSlocation($pt, $code) {
 370:     static $prefs = array(
 371: '--選択--','北海道', '青森県', '岩手県', '宮城県', '秋田県', '山形県', '福島県', 
 372: '茨城県', '栃木県', '群馬県', '埼玉県', '千葉県', '東京都', '神奈川県', 
 373: '新潟県', '富山県', '石川県', '福井県', '山梨県', '長野県', '岐阜県', '静岡県', 
 374: '愛知県', '三重県', '滋賀県', '京都府', '大阪府', '兵庫県', '奈良県',
 375: '和歌山県', '鳥取県', '島根県', '岡山県', '広島県', '山口県', '徳島県',
 376: '香川県', '愛媛県', '高知県', '福岡県', '佐賀県', '長崎県', '熊本県', '大分県', 
 377: '宮崎県', '鹿児島県', '沖縄県');
 378: 
 379:     if ($code == '') {
 380:         $pp = '';
 381:     } else if (! isset($pt->index[$code])) {
 382:         $errmsg = "not exist location '{$code}'";
 383:         return FALSE;
 384:     } else {
 385:         $pp = $pt->index[$code][4];
 386:     }
 387: 
 388:     //都道府県
 389:     $js = "var pref = [\n";
 390:     foreach ($prefs as $key=>$val) {
 391:         $selected = ($val == $pp) ? 'selected' : '';
 392:         $js .= "{'type':{$key}, 'text':'{$val}', 'value':'{$val}', 'selected':'{$selected}'},\n";
 393:     }
 394:     //地点
 395:     $js .= "];\nvar loc = [\n";
 396:     foreach ($prefs as $key=>$val) {
 397:         $js .= "{'type':{$key}, 'text':'--選択--', 'value':''},\n";
 398:         foreach ($pt->index as $cd=>$arr) {
 399:             if ($arr[4] == $val) {
 400:                 $selected = ($cd == $code) ? 'selected' : '';
 401:                 $js .= "{'type':{$key}, 'text':'{$arr[0]}', 'value':'{$cd}', 'selected':'{$selected}'},\n";
 402:             }
 403:         }
 404:     }
 405: 
 406:     //スクリプト
 407:     $js .=<<< EOT
 408: ];
 409: 
 410: $(function() {
 411:     var html = '';
 412:     for (var i in pref) {
 413:         html += '<option value="'+ pref[i].value + '" ' + pref[i].selected + '>'+ pref[i].text +'</option>';
 414:     }
 415:     $('#pref').append(html);
 416: 
 417:     var html = '';
 418:     var index = $('#pref option:selected').index();
 419:     for (var i in loc) {
 420:         if (loc[i].type == index) {
 421:             html += '<option value="'+ loc[i].value + '" ' + loc[i].selected + '>' + loc[i].text +'</option>';
 422:         }
 423:     }
 424:     $('#code').append(html);
 425: 
 426:     //都道府県を選択した時の挙動
 427:     $('#pref').change(function() {
 428:         var html = '';
 429:         var index = $('#pref option:selected').index();
 430:         //地点の内容を削除
 431:         $('#code').empty();
 432:         //対応する地点の内容を表示
 433:         for (var i in loc) {
 434:             if (loc[i].type == index) {
 435:                 html += '<option value="'+ loc[i].value +'">'+ loc[i].text +'</option>';
 436:             }
 437:         }
 438:         $('#code').append(html);
 439:     });
 440: });
 441: 
 442: EOT;
 443:     return $js;
 444: }

 447: /**
 448:  * 潮位表：潮位・月齢を計算し，配列へ格納する．
 449:  * @param   pahooTide $pt  潮位計算クラス
 450:  * @param   pahooCalendar $pc  暦計算クラス
 451:  * @param   array  $inputs 計算用パラメータ
 452:  * @param   array  $items  計算結果格納用配列
 453:  * @param   array  $locs   地点情報格納用配列
 454:  * @param   string $errmsg エラーメッセージ格納用；エラーなければ空文字
 455:  * @return  bool TRUE：計算成功／FALSE：失敗
 456: */
 457: function calcTideTable($pt, $pc, $params, &$items, &$locs, &$errmsg) {
 458:     //地点を設定
 459:     $pt->setLocation($params['code']);
 460:     if ($pt->iserror()) return FALSE;
 461: 
 462:     //地点を取得
 463:     $pt->getLocation($params['code'], $locs);
 464:     if ($pt->iserror()) return FALSE;
 465: 
 466:     //月齢・潮位を計算
 467:     $res = TRUE;
 468:     $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], 0, 0, 0);
 469:     for ($i = 0; $i < $params['interval']; $i++) {
 470:         $arr = array();
 471:         list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
 472:         $yb = $pc->getWeekString($year, $month, $day);
 473:         $items[$i]['dt'] = sprintf("%04d-%02d-%02d(%s)", $year, $month, $day, $yb);
 474:         $items[$i]['moonage'] = $pc->moon_age($year, $month, $day, MOONAGE_HOUR, 0, 0);
 475:         $moonmeridian = $pc->moon_time(2, $params['longitude'], $params['latitude'], $params['height'], $year, $month, $day);
 476:         $moonmeridian = ($moonmeridian == FALSE) ? '---' : $pc->day2hhmm($moonmeridian);
 477:         $items[$i]['tide'] = (preg_match('/([0-9]+)\:([0-9]+)/', $moonmeridian, $arr) > 0) ? $pc->tide($year, $month, $day, $arr[1], $arr[2], 0) : '';
 478:         unset($arr);
 479: 
 480:         $arr = array();
 481:         $pt->tide_day($params['code'], $year, $month, $day, $arr);
 482:         //満潮
 483:         for ($j = 0; $j < 4; $j++) {
 484:             if (isset($arr['high'][$j])) {
 485:                 $items[$i]['high'][$j]['hhmm']  = $arr['high'][$j]['hhmm'];
 486:                 $items[$i]['high'][$j]['lev']   = sprintf("%d", $arr['high'][$j]['lev']);
 487:                 $items[$i]['high'][$j]['align'] = 'text-align:right;';
 488:             } else {
 489:                 $items[$i]['high'][$j]['hhmm']  = '*';
 490:                 $items[$i]['high'][$j]['lev']   = '*';
 491:                 $items[$i]['high'][$j]['align'] = 'text-align:center;';
 492:             }
 493:         }
 494:         //干潮
 495:         for ($j = 0; $j < 4; $j++) {
 496:             if (isset($arr['low'][$j])) {
 497:                 $items[$i]['low'][$j]['hhmm']  = $arr['low'][$j]['hhmm'];
 498:                 $items[$i]['low'][$j]['lev']   = sprintf("%d", $arr['low'][$j]['lev']);
 499:                 $items[$i]['low'][$j]['align'] = 'text-align:right;';
 500:             } else {
 501:                 $items[$i]['low'][$j]['hhmm']  = '*';
 502:                 $items[$i]['low'][$j]['lev']   = '*';
 503:                 $items[$i]['low'][$j]['align'] = 'text-align:center;';
 504:             }
 505:         }
 506:         $jd++;
 507:     }
 508: 
 509:     return $res;
 510: }

 512: /**
 513:  * 潮位表：表示用HTMLを作成する．
 514:  * @param   pahooTide     $pt      潮位計算クラス
 515:  * @param   pahooCalendar $pc      暦計算クラス
 516:  * @param   array         $params  計算用パラメータ
 517:  * @param   string        $js      スクリプト格納用
 518:  * @param   string        $html    HTML格納用
 519:  * @param   string        $errmsg  エラーメッセージ格納用
 520:  * @return  bool TRUE：成功／FALSE：失敗
 521:  * @return  表示用HTML
 522: */
 523: function makeTable($pt, $pc, $params, &$js, &$html, &$errmsg) {
 524:     $js = $html = '';
 525:     $items = array();
 526:     $locs  = array();
 527:     if (calcTideTable($pt, $pc, $params, $items, $locs, $errmsg) == FALSE)  return FALSE;
 528: 
 529:     $hour = MOONAGE_HOUR;
 530:     $html =<<< EOT
 531: <table class="tide">
 532: <caption>{$locs['title']}（{$locs['prefecture']}{$locs['address']}）</caption>
 533: <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>
 534: 
 535: EOT;
 536:     $html .= "<tr>";
 537:     for ($j = 0; $j < 8; $j++) {
 538:         $html .= "<th class=\"index\">時刻</th><th class=\"index\">潮位<br />(cm)</th>";
 539:     }
 540:     $html .= "</tr>\n";
 541:     foreach ($items as $item) {
 542:         $html .= sprintf("<tr><td>%s</td><td>%.1f</td><td>%s</td>", $item['dt'], $item['moonage'], $item['tide']);
 543:         for ($j = 0; $j < 4; $j++) {
 544:             $html .= sprintf("<td>%s</td><td style=\"%s\">%s</td>", $item['high'][$j]['hhmm'], $item['high'][$j]['align'], $item['high'][$j]['lev']);
 545:         }
 546:         for ($j = 0; $j < 4; $j++) {
 547:             $html .= sprintf("<td>%s</td><td style=\"%s\">%s</td>", $item['low'][$j]['hhmm'], $item['low'][$j]['align'], $item['low'][$j]['lev']);
 548:         }
 549:         $html .= "</tr>\n";
 550:     }
 551:     $html .= "</table>\n";
 552: 
 553:     return TRUE;
 554: }

 556: /**
 557:  * 潮位グラフ：潮位・月齢を計算し，配列に格納する．
 558:  * @param   pahooTide $pt  潮位計算クラス
 559:  * @param   pahooCalendar $pc  暦計算クラス
 560:  * @param   array  $params 計算用パラメータ
 561:  * @param   array  $items  潮位計算結果格納用配列
 562:  * @param   array  $moons  月齢計算結果格納用配列
 563:  * @param   array  $locs   地点情報格納用配列
 564:  * @param   string $errmsg エラーメッセージ格納用；エラーなければ空文字
 565:  * @return  bool TRUE：計算成功／FALSE：失敗
 566: */
 567: function calcTideGraph($pt, $pc, $params, &$items, &$moons, &$locs, &$errmsg) {
 568:     //地点を設定
 569:     $pt->setLocation($params['code']);
 570:     if ($pt->iserror()) return FALSE;
 571: 
 572:     //地点を取得
 573:     $pt->getLocation($params['code'], $locs);
 574:     if ($pt->iserror()) return FALSE;
 575: 
 576:     $res = TRUE;
 577: 
 578:     //潮位を計算
 579:     $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], 0, 0);
 580:     $n = $params['interval'] / INTERVAL_GRAPH;
 581:     for ($i = 0; $i < $n; $i++) {
 582:         list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
 583:         $items[$i]['dt'] = sprintf("%04d/%02d/%02d %02d:%02d", $year, $month, $day, $hour, $min);
 584:         list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
 585:         $pt->sun_moon($year, $month, $day);
 586:         $items[$i]['tlevel'] = $pt->tide_level($year, $month, $day, $hour, $min);
 587:         $jd += INTERVAL_GRAPH;
 588:     }
 589: 
 590:     //月齢を計算
 591:     $height = 0;        //標高0メートルを仮定
 592:     $n = $params['interval'];
 593:     $jd = $pc->Gregorian2JD($params['year'], $params['month'], $params['day'], 0, 0, 0);
 594:     $i = $d = 0;
 595:     while ($i < $n) {
 596:         //月の南中時刻
 597:         list($year, $month, $day, $hour, $min, $sec) = $pc->JD2Gregorian($jd);
 598:         $moonmeridian = $pc->moon_time(2, $params['longitude'], $params['latitude'], $height, $year, $month, $day);
 599:         if ($moonmeridian != FALSE) {
 600:             $x = (double)$i + $moonmeridian;
 601:             $moons[$i]['t']   = $d + $moonmeridian;
 602:             $moons[$i]['age'] = $pc->moon_age($year, $month, $day, $moonmeridian * 24, 0, 0);
 603:             $i++;
 604:         }
 605:         $jd++;
 606:         $d++;
 607:     }
 608: 
 609:     return $res;
 610: }

 612: /**
 613:  * 潮位グラフ：表示用スクリプトおよびHTMLを求める．
 614:  * @param   pahooTide     $pt      潮位計算クラス
 615:  * @param   pahooCalendar $pc      暦計算クラス
 616:  * @param   array         $params  計算用パラメータ
 617:  * @param   string        $js      スクリプト格納用
 618:  * @param   string        $html    HTML格納用
 619:  * @param   string        $errmsg  エラーメッセージ格納用
 620:  * @return  bool TRUE：成功／FALSE：失敗
 621: */
 622: function makeGraph($pt, $pc, $params, &$js, &$html, &$errmsg) {
 623:     $js = $html = '';
 624:     $items = array();
 625:     $moons = array();
 626:     $locs  = array();
 627:     if (calcTideGraph($pt, $pc, $params, $items, $moons, $locs, $errmsg) == FALSE)  return FALSE;
 628: 
 629:     $title = "{$locs['title']}（{$locs['prefecture']}{$locs['address']}）";
 630:     $name   = GRAPH_TIDE;
 631:     $width  = MAP_WIDTH;
 632:     $height = MAP_HEIGHT;
 633: 
 634:     //潮位プロット・データ作成
 635:     $data = '';
 636:     foreach ($items as $item) {
 637:         $data .= sprintf("['%s', %f], ", $item['dt'], $item['tlevel']);
 638:     }
 639: 
 640:     //月齢データ作成
 641:     $moon = '';
 642:     $path = MOONAGE;
 643:     foreach ($moons as $val) {
 644:         $x = round(30 + (MAP_WIDTH - 70) / $params['interval'] * $val['t']);
 645:         $fname = sprintf("{$path}moon_%02d.png", round($val['age']));
 646:         $moon .=<<< EOT
 647: <img style="position:absolute; top:35px; left:{$x}px; width:40px; height:40px; z-index:999;" src="{$fname}" />
 648: 
 649: EOT;
 650:     }
 651: 
 652:     $js =<<< EOT
 653: $(function() {
 654:     jQuery.jqplot('{$name}',
 655:     [
 656:         [ {$data} ]
 657:     ],
 658:     {
 659:         //タイトル
 660:         title: {
 661:             text: '{$title}',
 662:             show: true,
 663:             fontSize: '16px',
 664:             textAlign: 'center',
 665:             textColor: 'black'
 666:         },
 667:         //背景
 668:         grid: {
 669:             background: '#EEFFFF'
 670:         },
 671:         //グラフ
 672:         seriesDefaults: {
 673:             showLine: true,
 674:             rendererOptions: { smooth: false },
 675:             markerOptions: { size: 0 },
 676:             color: 'blue',
 677:         },
 678:         //軸ラベル
 679:         axes: {
 680:             xaxis: {
 681:                 renderer: $.jqplot.DateAxisRenderer,
 682:                 tickOptions: { formatString: '%m/%d' },
 683:                 tickInterval: '1 days'
 684:             },
 685:             yaxis: {
 686:                 label: '潮位(cm)'
 687:             }
 688:         }
 689:     }
 690:     );
 691: });
 692: 
 693: EOT;
 694:     $html =<<< EOT
 695: <div id="{$name}" style="width:{$width}px; height:{$height}px;">
 696: {$moon}
 697: </div>
 698: </body>
 699: 
 700: EOT;
 701: 
 702:     return TRUE;
 703: }

• 各種WebAPIの登録方法：ぱふぅ家のホームページ
• PHPで日出没・月出没・月齢・潮を計算：ぱふぅ家のホームページ
• PHPでNHK政治意識月例調査をグラフ表示：ぱふぅ家のホームページ
• PHPで月の満ち欠けを描画：ぱふぅ家のホームページ
• 潮位表掲載地点一覧表：気象庁