พื้นที่เป็นตารางกิโลเมตร, Geometry เป็น Decimal Degrees

ในแบบฝึกหัดนี้ จะเริ่มต้นด้วย GeoDataFrame council_districts และ DataFrame permits อีกครั้ง โดยจะเปลี่ยน council_districts ให้ใช้ระบบอ้างอิงพิกัด EPSG 3857 ก่อนสร้างคอลัมน์ area จากนั้นเมื่อสร้างคอลัมน์ area เรียบร้อยแล้ว ให้เปลี่ยน CRS กลับเป็น EPSG 4326 เพื่อให้ geometry อยู่ในรูปแบบ decimal degrees

แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร

การแสดงผลข้อมูลเชิงภูมิศาสตร์ด้วย Python

ดูคอร์ส

คำแนะนำการฝึกหัด

เปลี่ยนระบบอ้างอิงพิกัดของ council_districts เป็น EPSG 3857 แล้วพิมพ์ crs และสองแถวแรกอีกครั้ง
สร้างคอลัมน์ชื่อ area โดยหารพื้นที่ของแต่ละ polygon ด้วย sqm_to_sqkm เพื่อให้ได้พื้นที่เป็นตารางกิโลเมตร
เปลี่ยนระบบอ้างอิงพิกัดของ council_districts กลับเป็น EPSG 4326 แล้วพิมพ์ crs และสองแถวแรก

แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ

ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์

# 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())

แก้ไขและรันโค้ด