การหาจุดตัดระหว่าง GeoDataFrame กับ Polygon
เมื่อนำชุดข้อมูลการใช้ที่ดินและชุดข้อมูลเขตมารวมกัน เราสามารถวิเคราะห์ได้ว่าเขตหนึ่ง ๆ มีการใช้ที่ดินในลักษณะใดบ้าง
ขั้นแรก เราต้องหาจุดตัดระหว่างชุดข้อมูลการใช้ที่ดินกับเขตที่กำหนด โดยจะใช้เขต Muette เป็นตัวอย่าง
ชุดข้อมูลการใช้ที่ดินและเขตถูกนำเข้ามาแล้วในชื่อ land_use และ districts และได้แยกเขต Muette ออกมาเป็น shapely polygon ชื่อ muette พร้อมกับนำเข้า GeoPandas และ matplotlib แล้ว
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การทำงานกับข้อมูลเชิงพื้นที่ใน Python
คำแนะนำการฝึกหัด
- คำนวณจุดตัดระหว่าง polygon ใน
land_useกับ polygonmuetteและเก็บผลลัพธ์ไว้ในชื่อland_use_muette - สร้างกราฟแสดงผลจุดตัดนี้ โดยส่ง
edgecolor='black'เพื่อให้เห็นขอบเขตของแต่ละ polygon ได้ชัดเจนขึ้น - แสดง 5 แถวแรกของ
land_use_muette
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)