1. 学ぶ
    2. /
  2. コース
    3. /
  3. Pythonで扱う地理空間データ
Connected

演習

ポイントの投影変換

前のチャプターでは、エッフェル塔の位置を扱いました。距離を計算できるように、今回も投影座標系で座標を用意しました。ここではもう一度この象徴的なランドマークに戻り、地理座標で表してみましょう。48°51′29.6″N, 2°17′40.2″E、10 進表記では緯度 48.8584、経度 2.2945 です。

Shapely のジオメトリオブジェクト自体には CRS の概念がないため、直接ほかの CRS に変換することはできません。そこで GeoPandas を使って、エッフェル塔のポイント位置を別の CRS に変換します。単一点を GeoSeries に入れ、to_crs() メソッドで変換し、再びポイントを取り出します。

GeoPandas はすでにインポートされています。

指示

100 XP
  • エッフェル塔の座標で shapely の Point オブジェクトを作成し、eiffel_tower という変数に代入します。
  • エッフェル塔のみを要素に持つ GeoSeries（s_eiffel_tower）を作成し、CRS に EPSG:4326 を指定します。
  • s_eiffel_tower を EPSG:2154 に変換し、結果を s_eiffel_tower_projected として保存します。