座標系の変換
緯度経度は度になっているが
面積、距離などを求める時には変換が必要
https://sorabatake.jp/20670/
も参考にする
Name: WGS 84 – Lat[north]: Geodetic latitude (degree) – Lon[east]: Geodetic longitude (degree)
がよく使われる
これはGPSで使われる
WGS 84
これを
EPSGにすると4326
https://epsg.io/4326
単位は度
これは世界
日本近辺に特化したものが
EPSG:6668
https://epsg.io/6668
単位は度
日本付近に絞るならこちらを使う方が良さそう
EPSG:6691
https://epsg.io/6691
これは
単位がメートルになっている
こちらも日本周辺がメイン
Area of use: Japan – between 138°E and 144°E, onshore and offshore.
となっているので
それ以外のエリアの精度が正しくないことがあるので注意
とりあえず品川の情報を表示
shinagawa N03_001 N03_002 N03_003 N03_004 N03_007 geometry 87 東京都 None 品川区 None 13109 POLYGON ((139.75501 35.61771, 139.75507 35.617... 88 東京都 None 品川区 None 13109 POLYGON ((139.77231 35.62188, 139.77302 35.621... 89 東京都 None 品川区 None 13109 POLYGON ((139.75953 35.625, 139.75966 35.6216,... 90 東京都 None 品川区 None 13109 POLYGON ((139.71943 35.64167, 139.71953 35.641...
座標情報を見るなら crs で表示できる
shinagawa.crs
結果は
<Geographic 2D CRS: EPSG:6668> Name: JGD2011 Axis Info [ellipsoidal]: - Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree) Area of Use: - name: Japan - onshore and offshore. - bounds: (122.38, 17.09, 157.65, 46.05) Datum: Japanese Geodetic Datum 2011 - Ellipsoid: GRS 1980 - Prime Meridian: Greenwich
この中で注目するのは
EPSG:6668
と
- Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree)
から
単位は度であることがわかる
同様に路線図の crs も確認する
gdf.crs
結果は
<Geographic 2D CRS: EPSG:4326> Name: WGS 84 Axis Info [ellipsoidal]: - Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree) Area of Use: - name: World. - bounds: (-180.0, -90.0, 180.0, 90.0) Datum: World Geodetic System 1984 ensemble - Ellipsoid: WGS 84 - Prime Meridian: Greenwich
こちらも
EPSG:4326
と
- Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree)
から単位は度であることがわかる
この２つはEPSGコードが異なるため
このままだと一緒に使えないので変換をして合わせる
今回は路線図のEPSGコードを6668に変換するので
to_crs を使う
gdf.to_crs(epsg=6668,inplace=True) gdf.crs
これで
<Geographic 2D CRS: EPSG:6668> Name: JGD2011 Axis Info [ellipsoidal]: - Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree) Area of Use: - name: Japan - onshore and offshore. - bounds: (122.38, 17.09, 157.65, 46.05) Datum: Japanese Geodetic Datum 2011 - Ellipsoid: GRS 1980 - Prime Meridian: Greenwich
となる