การนับจำนวนจาก Spatial Join
จะใช้ชุดข้อมูลใบอนุญาตก่อสร้างที่ออกในเมือง Nashville ในปี 2017 โดย DataFrame ชื่อ permits ถูกโหลดไว้ให้แล้ว พร้อมกับ GeoDataFrame ชื่อ council_districts
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การแสดงผลข้อมูลเชิงภูมิศาสตร์ด้วย Python
คำแนะนำการฝึกหัด
- สร้างคอลัมน์
geometryในpermitsจากlatและlng - สร้าง
permits_geoเป็น GeoDataFrame โดยใช้permits, ค่า crs จากcouncil_districtsและgeometryจากpermits - ใช้ spatial join เพื่อค้นหาใบอนุญาตที่อยู่
withinแต่ละเขต (council district) แล้วพิมพ์ 2 แถวแรก - สร้าง
permit_countsเพื่อแสดงจำนวนใบอนุญาตในแต่ละเขต โดยใช้groupby()และ.size()จากนั้นพิมพ์permit_counts
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
import geopandas as gpd
# Create a geometry column in permits from lat and lng
permits[____] = gpd.points_from_xy(____.____ , ____.____)
# Build a GeoDataFrame: permits_geo
permits_geo = gpd.GeoDataFrame(____, crs = ____.crs, geometry = ____.geometry)
# Spatial join of permits_geo and council_districts
permits_by_district = gpd.____(permits_geo, council_districts, predicate = ____)
print(permits_by_district.head(2))
# Create permit_counts
permit_counts = permits_by_district.groupby([____]).____()
print(permit_counts)