แปลง CRS ให้ตรงกันและบันทึกลงไฟล์
จากแบบฝึกหัดก่อนหน้า เราพบว่าชุดข้อมูลทั้งสองใช้ Coordinate Reference System (CRS) ที่แตกต่างกัน ซึ่งเห็นได้ชัดจากกราฟแรกในแบบฝึกหัดนี้ (โค้ดสำหรับกราฟนี้เตรียมไว้ให้แล้วในสคริปต์) เนื่องจากชุดข้อมูลทั้งสองครอบคลุมพื้นที่เดียวกัน พิกัดควรจะซ้อนทับกัน แต่กลับไม่เป็นเช่นนั้น
เพื่อให้การวิเคราะห์ในบทนี้ดำเนินต่อไปได้ เราจะแปลงชุดข้อมูลทั้งสองให้ใช้ CRS เดียวกัน แล้วบันทึกลงไฟล์ใหม่ เพื่อให้สามารถคำนวณระยะทางได้ จะแปลงไปยัง projected CRS ได้แก่ UTM zone 35 ในท้องถิ่น ซึ่งระบุด้วย EPSG:32735 (https://epsg.io/32735)
ชุดข้อมูลแหล่งขุดแร่ (mining_sites) และอุทยานแห่งชาติ (national_parks) โหลดไว้แล้ว รวมถึง GeoPandas และ matplotlib ก็ถูก import มาแล้วเช่นกัน
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การทำงานกับข้อมูลเชิงพื้นที่ใน Python
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Plot the natural parks and mining site data
ax = national_parks.plot()
mining_sites.plot(ax=ax, color='red')
plt.show()
# Convert both datasets to UTM projection
mining_sites_utm = ____
national_parks_utm = ____
# Plot the converted data again
ax = national_parks_utm.plot()
mining_sites_utm.plot(ax=ax, color='red')
plt.show()