国土数値情報の点データをmysqlにぶち込むプログラム

書いた。

$xml = simplexml_load_file(".xmlのパス");
$array = array();

if( count($xml->GM_Point) == count($xml->FG01) ){

  for( $i=0;$i<count($xml->GM_Point);$i++ ){

    //XMLから緯度経度情報を取得する
    $pointArr = (array)$xml->GM_Point[$i];
    $pointArr = (array)$pointArr['GM_Point.position'];
    $pointArr = (array)$pointArr['DirectPosition'];
    $point = $pointArr['DirectPosition.coordinate'];

    //mysqlのジオメトリ型に対応するため緯度と経度を入れ替える
    ereg('(.+) (.+)',$point, $regs);
    $point = $regs[2].' '.$regs[1];

    //住所とカテゴリを取得する
    $adsArr = (array)$xml->FG01[$i];
    $address = (string)$adsArr['ADS'];
    $category = (string)$adsArr['FSC'];

    //配列にブッ込む
    $array[$i] = array();
    $array[$i][0] = $point;
    $array[$i][1] = $address;
    $array[$i][2] = $category;

  }//for end

  //DB接続
  $con = mysql_connect("mysql223.db.sakura.ne.jp","あかうんと","ぱすわーど")
  or die("接続に失敗しました");

  $db_selected = mysql_select_db('DB名', $con);
  print(mysql_error());

  //データ挿入

  //utf8でブッ込む
  $result = mysql_query("set names utf8");

  //こんなやり方もあるらしいが今回は使わない
  //$result = mysql_query("SET @g = GeomFromText('POINT(139.48053 35.299015)')");
  //$result = mysql_query("insert into sh_spot (point,name,category) values (@g,'enoshima',1)");

  for($i=0;$i<count($array);$i++){

    $result = mysql_query("insert into sh_spot (point,address,category) values ( GeomFromText('POINT(".$array[$i][0].")'),'".$array[$i][1]."','".$array[$i][2]."')");

  }//for end
  print(mysql_error());
}//if end

ちょっと必要になったので、国土数値情報の燃料給油所データ（一都五県で約１万件）からデータベースを構築した。DBはmysql。今回、mysqlの空間データ型”Geometry”の存在におそばせながら気づき、驚愕した。
距離計算とか簡単な空間分析ならmysqlでも既に出来るみたい。パネェ！
googlemapとの相性もイイらしい。

あとはXMLデータの扱いに若干苦労した。なんやら”SimpleXMLElement”とかゆうめんどくさいのがでてきたのを、(array)とかつけて解決。ちなみにXMLの中身はあらかじめ余分なデータを削除して整形しておいた。

配列作らず直でinsertしていけばループ減って早くなるだろうが、全然おkだった。


参考URL
・国土数値情報ダウンロードサービス
・MySQL 4.1 リファレンスマニュアル - 10.2.2. Geometry クラス
・MySQL における空間情報の機能
・位置情報サービスのはじめ方 - 第5回　位置情報を保存しよう（前編）
・システム開発の備忘録 - SimpleXML関数の扱い方
・PHP Manual - 例1 正規表現の例