ตรวจสอบผลลัพธ์ของ overlay
หลังจากสร้าง overlay ระหว่างชุดข้อมูลการใช้ที่ดินและเขตการปกครองแล้ว เราสามารถวิเคราะห์การใช้ที่ดินในแต่ละเขตได้ง่ายขึ้น มาดูตัวอย่างของเขต Muette และตรวจสอบการใช้ที่ดินในเขตนั้นกัน
นำเข้า GeoPandas และ Matplotlib ไว้แล้ว ผลลัพธ์ของฟังก์ชัน overlay() จากแบบฝึกหัดก่อนหน้าพร้อมใช้งานในชื่อ combined
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การทำงานกับข้อมูลเชิงพื้นที่ใน Python
คำแนะนำการฝึกหัด
- เพิ่มคอลัมน์ใหม่ชื่อ
'area'ที่เก็บค่าพื้นที่ของแต่ละ polygon ลงใน GeoDataFramecombined - สร้างชุดข้อมูลย่อยชื่อ
land_use_muetteโดยกรองเฉพาะแถวที่'district_name'เท่ากับ "Muette" - สร้างกราฟของ
land_use_muetteโดยใช้คอลัมน์'class'เพื่อระบายสี polygon - คำนวณพื้นที่รวมของแต่ละ
'class'ในland_use_muetteโดยใช้เมธอดgroupby()แล้วแสดงผลลัพธ์
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Print the first rows of the overlay result
print(combined.head())
# Add the area as a column
____ = ____
# Take a subset for the Muette district
land_use_muette = combined[____]
# Visualize the land use of the Muette district
land_use_muette.____(____)
plt.show()
# Calculate the total area for each land use class
print(land_use_muette.____(____)['area'].____() / 1000**2)