อ่านไฟล์ raster เข้าสู่ R
วิธีที่ง่ายที่สุดในการอ่านไฟล์ raster เข้าสู่ R คือการใช้ฟังก์ชัน raster() จากแพ็กเกจ raster เพียงส่งชื่อไฟล์ (รวมนามสกุล) ของ raster เป็นอาร์กิวเมนต์แรก x
ฟังก์ชัน raster() ใช้ฟังก์ชันภายในของแพ็กเกจ raster สำหรับการอ่านไฟล์บางประเภท (โดยพิจารณาจากนามสกุลไฟล์) และส่งต่อการอ่านไฟล์ประเภทอื่นให้ readGDAL() จากแพ็กเกจ rgdal ข้อดีของการไม่ใช้ readGDAL() โดยตรงคือ raster() จะคืนค่าเป็น object ประเภท RasterLayer
ไฟล์ raster ที่พบบ่อยคือ GeoTIFF ซึ่งมีนามสกุล .tif หรือ .tiff เราได้ดาวน์โหลดข้อมูล raster รายได้เฉลี่ยจาก US census และวางไว้ในไดเรกทอรีทำงานของคุณแล้ว
มาดูและอ่านข้อมูลนี้กัน
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การสร้างภาพข้อมูลภูมิสารสนเทศใน R
คำแนะนำการฝึกหัด
- ใช้
dir()เพื่อดูไฟล์ในไดเรกทอรีทำงาน - ใช้
dir()อีกครั้งเพื่อดูภายในไดเรกทอรีnyc_grid_data - ใช้
raster()เพื่ออ่านข้อมูล raster รายได้เฉลี่ยเก็บไว้ในตัวแปรincome_gridโดยส่ง path แบบเต็มไปยังไฟล์.tif - ใช้
summary()เพื่อตรวจสอบว่า raster ถูกเก็บในรูปแบบRasterLayer - ใช้
plot()เพื่อตรวจสอบเนื้อหาของ raster
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
library(raster)
# Call dir()
# Call dir() on the directory
# Use raster() with file path: income_grid
# Call summary() on income_grid
# Call plot() on income_grid