平方キロメートルの面積、小数度のジオメトリ
この演習では、council_districts のGeoDataFrameと permits のDataFrameから再開します。area 列を作成する前に、council_districts をEPSG 3857の座標参照系に変更します。area 列を作成したら、ジオメトリを小数度にするためにCRSをEPSG 4326に戻します。
この演習はコースの一部です
Pythonで可視化する地理空間データ
演習の手順
council_districtsの座標参照系をEPSG 3857に変更し、crsと先頭2行をもう一度表示します。areaという列を作成します。各ポリゴンの面積をsqm_to_sqkmで割って、面積を平方キロメートルで求めます。council_districtsの座標参照系をEPSG 4326に戻します。crsと先頭2行を表示します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Change council_districts crs to epsg 3857
council_districts = council_districts.____(____ = ____)
print(council_districts.crs)
print(council_districts.head())
# Create area in square km
sqm_to_sqkm = 10**6
council_districts['area'] = council_districts.____ / sqm_to_sqkm
# Change council_districts crs back to epsg 4326
council_districts = council_districts.____(____ = ____)
print(council_districts.crs)
print(council_districts.head())