![SRSの定義
GEOGCS["JGD2011"
,DATUM["Japanese Geodetic Datum 2011",
SPHEROID["GRS 1980",
6378137,
298.257222101,
AUTHORITY["EPSG","7019"]]
,AUTHORITY["EPSG","1128"]]
,PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]]
,UNIT["degree",0.017453292519943278,
AUTHORITY["EPSG","9122"]]
,AXIS["Lat",NORTH]
,AXIS["Lon",EAST]
,AUTHORITY["EPSG","6668"]]
? MySQLでは SRSの定義を INFORMATION_SCHEMA の以下のテーブル
で参照できる
ST_SPATIAL_REFERENCE_SYSTEMS
地理座標系の定義(JGD2011)：](https://image.slidesharecdn.com/osc2018tokyofallmysqlgis-181030105309/85/MySQL8-0-GIS-22-320.jpg)
![PROJCS["JGD2011 / Japan Plane Rectangular CS IX",
GEOGCS["JGD2011",
....]
,PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]]
,PARAMETER["Latitude of natural origin" , 36
,AUTHORITY["EPSG","8801"]]
,PARAMETER["Longitude of natural origin",139.833333333333
,AUTHORITY["EPSG","8802"]]
,PARAMETER["Scale factor at natural origin",0.9999
,AUTHORITY["EPSG","8805"]]
,PARAMETER["False easting" ,0 ,AUTHORITY["EPSG","8806"]]
,PARAMETER["False northing",0 ,AUTHORITY["EPSG","8807"]]
,UNIT["metre",1 ,AUTHORITY["EPSG","9001"]]
,AXIS["X",NORTH]
,AXIS["Y",EAST]
,AUTHORITY["EPSG","6677"]
]
投影座標系の定義(JGD2011 平面直角座標系 IX(東京付近))：](https://image.slidesharecdn.com/osc2018tokyofallmysqlgis-181030105309/85/MySQL8-0-GIS-23-320.jpg)
惭测厂蚕尝8.0の「地理情报」を理解しよう入门～いまからはじめる骋滨厂
- 1. MySQL 8.0 の新機能 「地理情報」 を理解しよう入門 ～いまからはじめるGIS～ オープンソースカンファレンス2018 Tokyo/Fall 2018/10/28 日本MySQLユーザ会 坂井 恵 （@sakaik） ハッシュタグ： #mysql_jp #osc18tk
- 2. 速報 MySQL 8.0.13 が出ました！(2018/10/24) そして、MySQL Cluster 8.0.13 dmr も！ 「なぜ MySQL 5.7 の次が 8.0？」の疑問が、いま明らかに！ 「６」は、昔、一旦出たバージョンだからスキップするよ → うんうん 「７」はMySQL Cluster で使っているからスキップするよ → ファ？ ※感想には個人差があります。諸説あります。 MySQLとCluster のバージョンを棲み分けしたいのではなく、すべての MySQLファミリーのソフトウェア群のバージョンをそろえたかったんだ！
- 3. フルバージョンは slideshare の資料を 見てください（ClubMySQL #4 の資料 が一番濃厚です） /sakaik
- 4. データベースは楽しい！ ? データベースは、現実世界の射影 ? 手に取るように世界が見える ? 数字、文字列、日付を中心とした情報 ? 例：売上、売れ筋商品、顧客、学生の成績、etc ? →そして、そこにいま「場所」が加わる GIS機能を楽しもう！
- 5. MySQLのGISの歴史 ? MySQL 4.1.0 (2003): MyISAM に "Spatial Extensions" 実装 http://ftp.nchu.edu.tw/MySQL/doc/refman/4.1/en/spatial-extensions.html ? MySQL 5.0.16 (2005): InnoDB にも実装 ? MySQL 5.7.1 (2013): spatialデータ専用の型を実装（これまでは BLOBに保管） Boostライブラリを利用して InnoDBに再実装 ? MySQL 8.0.0 (2016): 関数名をST_* 系に統一したり関数増やした り 地球が丸いことを知ったり、良い感じに 2018/04/19 8.0.11でGAに！ 2018/10/24 8.0.13 リリース ※MySQLでは、「GIS」ではなく「Spatial（空間情報）」というキーワード
- 6. MySQLのGIS機能習得への道 習得へのポイント 1. 平面での「座標」の扱い方の理解 2. 地球（球体）では違うぞ、と認識 3. どう违うのかの理解?表现方法の习得
- 7. 平面での座標の扱い方 座標の基本： （ｘ，ｙ） 表現できる「形」： ? 点： ひとつの点 ? 線： 線分の接続集合 ? 面： 線分の接続集合で囲まれたエリア （始点と終点が同じ位置となる）
- 8. 型に関する知識 空間情報(Spatial) をあらわす３つの形態 ?点 (POINT) ?線 (LINESTRING) ?面 (POLYGON) そのほか、これらの集合である MULTIPOINT MULTILINESTRING MULTIPOLYGON GEOMETRYCOLLECTION もあります。 ひっくるめて 幾何(形状） (GEOMETRY)
- 9. Spatialの ３つの表現 テキスト／バイナリ／MySQL内部バイナリ表現 ? テキスト：Well known text (WKT) ? バイナリ：Well known binary (WKB) ? 内部バイナリ： WKB にSRID情報を追加した もの ※ここからは基本的に、WKT を使って話を進めます。
- 10. 点?線?面のあらわし方 (まず覚えよう) ? 点： POINT (3 5) ? 線： LINESTRING (1 2, 5 3, 7 9, 3 12) ? 面： POLYGON((1 1, 5 1, 5 8, 1 8, 1 1)) WKT (Well known text)： コンマではないことに注意 ? 複数点： MULTIPOINT (3 5, 1 2, 8 3) MULTIPOINT ((3 5), (1 2), (8 3)) ? 複数線： MULTILINESTRING ((1 2, 5 3, 7 9, 3 12), (4 1, 3 3, -1 4)) ? 複数面： MULTIPOLYGON (略) 必ず元の点にもどる POLYGON((1 1, 9 1, 9 9, 1 9, 1 1), (3 3, 5 3, 5 5, 3 5, 3 3)) 中をくり抜く 参考：
- 11. WKT をMySQL内部表現に変換するには、 ST_GeomFromText() 関数を使用する mysql> SELECT ST_GeomFromText('POINT(3 5)') mygeo; +---------------------------+ | mygeo | +---------------------------+ | @ @ | +---------------------------+ mysql> SELECT HEX(ST_GeomFromText('POINT(3 5)')) mygeo; +----------------------------------------------------+ | mygeo | +----------------------------------------------------+ | 00000000010100000000000000000008400000000000001440 | +----------------------------------------------------+ mysql> SELECT ST_AsText(ST_GeomFromText('POINT(3 5)')) mygeo; +------------+ | mygeo | +------------+ | POINT(3 5) | +------------+
- 12. 点?線?面のあらわし方 WKT → MySQLの内部バイナリ： ? ST_GeomFromText() 関数を使用する。 MySQLの内部バイナリ → WKT： ? ST_AsText() 関数を使用する。 mysql> SELECT ST_AsText(ST_GeomFromText('LINESTRING (1 2, 5 3, 7 9, 3 12)')) geom; +------------------------------+ | geom | +------------------------------+ | LINESTRING(1 2,5 3,7 9,3 12) | +------------------------------+ Geom の代わりに Point, LineString, Polygon などの関数が 個別にそれぞれあります
- 13. 様々な関数 個別の関数名の紹介は今日はしませんが??? ? この２つの線(面)は交わってますか？ ? このエリア内に線(点、面)は含まれていますか？ ? この点、線、面を含む最小の矩形を返して ? この線の長さは？ ? この２つの点のあいだの距離は？ のような事ができる関数群が用意されています。
- 14. 平面での座標の扱い(おさらい） ? 点、線、面、それらの上位である Geometryがある ? 人間が見やすく書きやすいのは「ＷＫＴ」 ? MySQL内部で扱うのは「MySQL内部バイナリ」 ? WKTとMySQL内部バイナリの変換は、 ST_GeomFromText() と ST_AsText() を使用する
- 17. いろいろな地球のモデル 楕円体のままで位置を特定する「地理座標系」 ? 長半径と扁平率の定義が異なる様々な地球の形がある ? それぞれにID(SRS_ID;SRID)と SRS名が付けられている ? 日本で使われる主な 地理座標系の SRS SRS_NAME SRS_ID 備考 長半径 扁平率 の逆数 Tokyo 4301 今はほとんど使われない。 WGS84との差も大きい 6,377,397.155 m 299.1528128 JGD2000 4612 WGS84との差は小さい 6,378,137m 298.257222101 JGD2011 6668 東日本大震災の地殻変化 に対応 6,378,137m 298.257222101 WGS 84 4326 世界でよく使われる 6,378,137m 298.257223563 ※とりあえず JGD2011 (6668)、WGS84 (4326) の２つの SRS_ID（数字）は覚えておくとよさげ
- 18. いろいろな地球 平面の地図に落とし込む「投影座標系」 ? 回転楕円体を平面の地図に落とし込む(＝投影する)ルール ? 地図の中心となる緯度、経度や投影方法などの情報を持つ ? それぞれにID(SRS_ID)と SRS名が付けられている mysql> SELECT SRS_NAME, SRS_ID, ORGANIZATION -> FROM ST_SPATIAL_REFERENCE_SYSTEMS WHERE SRS_NAME LIKE 'JGD2011%'; +--------------------------------------------+--------+--------------+ | SRS_NAME | SRS_ID | ORGANIZATION | +--------------------------------------------+--------+--------------+ | JGD2011 | 6668 | EPSG | | JGD2011 / Japan Plane Rectangular CS I | 6669 | EPSG | | JGD2011 / Japan Plane Rectangular CS II | 6670 | EPSG | | JGD2011 / Japan Plane Rectangular CS III | 6671 | EPSG | | JGD2011 / Japan Plane Rectangular CS IV | 6672 | EPSG | | JGD2011 / Japan Plane Rectangular CS V | 6673 | EPSG | | JGD2011 / Japan Plane Rectangular CS VI | 6674 | EPSG | | JGD2011 / Japan Plane Rectangular CS VII | 6675 | EPSG | | JGD2011 / Japan Plane Rectangular CS VIII | 6676 | EPSG | | JGD2011 / Japan Plane Rectangular CS IX | 6677 | EPSG | | JGD2011 / Japan Plane Rectangular CS X | 6678 | EPSG | | JGD2011 / Japan Plane Rectangular CS XI | 6679 | EPSG | | JGD2011 / Japan Plane Rectangular CS XII | 6680 | EPSG | | JGD2011 / Japan Plane Rectangular CS XIII | 6681 | EPSG | | JGD2011 / Japan Plane Rectangular CS XIV | 6682 | EPSG | | JGD2011 / Japan Plane Rectangular CS XV | 6683 | EPSG | | JGD2011 / Japan Plane Rectangular CS XVI | 6684 | EPSG | | JGD2011 / Japan Plane Rectangular CS XVII | 6685 | EPSG | | JGD2011 / Japan Plane Rectangular CS XVIII | 6686 | EPSG | | JGD2011 / Japan Plane Rectangular CS XIX | 6687 | EPSG | | JGD2011 / UTM zone 51N | 6688 | EPSG | | JGD2011 / UTM zone 52N | 6689 | EPSG | | JGD2011 / UTM zone 53N | 6690 | EPSG | | JGD2011 / UTM zone 54N | 6691 | EPSG | | JGD2011 / UTM zone 55N | 6692 | EPSG | +--------------------------------------------+--------+--------------+ 25 rows in set (0.00 sec) JGD2011だけでも、こんなに！ 各中心点(原点)は 測量法付随の国交省告示で http://www.gsi.go.jp/LAW/heimencho.html
- 19. JGD2011 「日本測地系2011」(学術系) 「世界測地系」(法律等) 「測地成果2011」(国土地理院)を元に制定 ? １つの回転楕円体（地理座標系；緯度経度）と ? 19の平面直角座標系と ? 5つの鲍罢惭
- 20. Japan Plane Rectangular CS ? 「日本の平面直角座標系」 ? 国土交通省告示第9号で定義 ? 全国を19のエリアに分割 ? 位置は、各エリアの原点か らの北、東方向への距離で 表す ? 各エリアは自治体界を境界 線とするよう定義されてい る（機械的に長方形で区切 られているわけではない） 国土地理院サイトより http://www.gsi.go.jp/sokuchikijun/jpc.html
- 21. UTM zone ? 日本近郊は zone 51N～55Nでカバーされてい る ? Universal Transverse Mercator ? 経度6度ごとに区切った細長い帯に 1～60の番 号を付与 ? 西経180度から東回りに 1から順に附番 ? 本来は緯度8度ごとに更に分割するが、MySQL では(?)、北側と南側の２つにのみ分割して 扱っている模様（だから zone 51は、51N と 51S の2 種類のみが存在する）
- 22. SRSの定義 GEOGCS["JGD2011" ,DATUM["Japanese Geodetic Datum 2011", SPHEROID["GRS 1980", 6378137, 298.257222101, AUTHORITY["EPSG","7019"]] ,AUTHORITY["EPSG","1128"]] ,PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]] ,UNIT["degree",0.017453292519943278, AUTHORITY["EPSG","9122"]] ,AXIS["Lat",NORTH] ,AXIS["Lon",EAST] ,AUTHORITY["EPSG","6668"]] ? MySQLでは SRSの定義を INFORMATION_SCHEMA の以下のテーブル で参照できる ST_SPATIAL_REFERENCE_SYSTEMS 地理座標系の定義(JGD2011)：
- 23. PROJCS["JGD2011 / Japan Plane Rectangular CS IX", GEOGCS["JGD2011", ....] ,PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]] ,PARAMETER["Latitude of natural origin" , 36 ,AUTHORITY["EPSG","8801"]] ,PARAMETER["Longitude of natural origin",139.833333333333 ,AUTHORITY["EPSG","8802"]] ,PARAMETER["Scale factor at natural origin",0.9999 ,AUTHORITY["EPSG","8805"]] ,PARAMETER["False easting" ,0 ,AUTHORITY["EPSG","8806"]] ,PARAMETER["False northing",0 ,AUTHORITY["EPSG","8807"]] ,UNIT["metre",1 ,AUTHORITY["EPSG","9001"]] ,AXIS["X",NORTH] ,AXIS["Y",EAST] ,AUTHORITY["EPSG","6677"] ] 投影座標系の定義(JGD2011 平面直角座標系 IX(東京付近))：
- 24. MySQLでの地理座標系の座標の表し方 ■テーブルの作成。 扱う座標系が決まっている場合はSRIDを指定する。 CREATE TABLE mygeo3 (id INTEGER, pnt POINT SRID 4326); INSERT INTO mygeo3 VALUES (1, ST_GeomFromText('POINT(43.057265 141.389053)', 4326)); ■データの挿入。SRIDを指定する。 ■SRIDがマッチしていないとエラーになる mysql> INSERT INTO mygeo3 VALUES (1, ST_GeomFromText('POINT(43.057265 141.389053)')); ERROR 3643 (HY000): The SRID of the geometry does not match the SRID of the column 'pnt'. The SRID of the geometry is 0, but the SRID of the column is 4326. Consider changing the SRID of the geometry or the SRID property of the column.
- 25. その他のトピック
- 26. 緯度、経度の 表し方 ? 度分秒の表し方と、度(小 数)の表し方がある ? MySQLでは、度(小数)を使 用する ? この部屋(28号館)の位置は、 およそ、東経139度24分25.6 秒、北緯35度38分42.7秒 度 分 秒 35 38 42.7 ↓ x1 ↓÷60 ↓÷60÷60 合計 35 0.633333 0.011861111 35.64519 ※図は、iOS用アプリ「FieldAccess2」より 変換方法の例：
- 27. 轴の向き
- 31. 人丸前駅ホーム
- 33. 地理情報に興味を持ったら ? FOSS4G 2018 Okayama ? 2018/11/04-05 ? https://www.osgeo.jp/events/foss4g-2018/foss4g- 2018-okayama ? FOSS4G 2018 Tokyo ? 2018/11/09-10 ? https://www.osgeo.jp/events/foss4g-2018/foss4g- 2018-tokyo
- 34. ? 引き続き12時からは、Oracle社(MySQL Community Team)の 「MySQL InnoDB Cluster」のお話があります。 ? ブースにもお越しください（日本MySQLユーザ会：3階）
- 35. ～ おしまい ～