การดึงข้อมูลโดยใช้แพ็กเกจ
การอ่านข้อมูลเชิงพื้นที่จากไฟล์เป็นวิธีหนึ่งในการนำข้อมูลเข้าสู่ R แต่ยังมีแพ็กเกจบางส่วนที่ให้ข้อมูลเชิงพื้นที่ที่ใช้กันทั่วไปด้วย ตัวอย่างเช่น แพ็กเกจ rnaturalearth ให้ข้อมูลจาก Natural Earth ซึ่งเป็นแหล่งข้อมูลแผนที่โลกความละเอียดสูง ครอบคลุมแนวชายฝั่ง รัฐ และสถานที่ที่มีประชากรอาศัยอยู่ แหล่งข้อมูลนี้คือแหล่งเดียวกับที่ใช้ในบทที่ 2
เราจะวิเคราะห์รายได้มัธยฐานในระดับ census tract ของ New York County (หรือที่รู้จักในชื่อ Borough of Manhattan) แต่ก่อนอื่นต้องรู้ขอบเขตของ census tract เสียก่อน แพ็กเกจ tigris ใน R ช่วยให้ดาวน์โหลดและนำเข้า shapefile ตามภูมิศาสตร์ของ US Census ได้อย่างสะดวก โดยจะใช้ฟังก์ชัน tracts() เพื่อดาวน์โหลดขอบเขต tract นอกจากนี้ tigris ยังมีฟังก์ชัน states(), counties(), places() และอื่น ๆ อีกมากที่สอดคล้องกับระดับต่าง ๆ ของหน่วยทางภูมิศาสตร์ที่ Census กำหนดไว้
มาดึงข้อมูลเชิงพื้นที่ของ tract กันเลย
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การสร้างภาพข้อมูลภูมิสารสนเทศใน R
คำแนะนำการฝึกหัด
- เรียกใช้
tracts()โดยกำหนดstate = "NY",county = "New York", และcb = TRUEแล้วเก็บผลลัพธ์ไว้ในตัวแปรnyc_tracts - ใช้
summary()กับnyc_tractsเพื่อตรวจสอบว่าออบเจกต์ที่ได้รับคืนเป็นSpatialPolygonsDataFrame - Plot
nyc_tractsเพื่อตรวจสอบเนื้อหาด้วยplot()
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
library(sp)
library(tigris)
# Call tracts(): nyc_tracts
# Call summary() on nyc_tracts
# Plot nyc_tracts