สร้างแรสเตอร์ด้วยข้อมูลอ้างอิงทางภูมิศาสตร์ใน R

อ้างถึงคำถามก่อนหน้าของฉันฉันจะ "geocode" / แปลง - ใช้ R - ภาพแรสเตอร์ที่สร้างขึ้นในรูปแบบที่สามารถแสดงเป็นเลเยอร์บนแผนที่ได้อย่างไร ในระยะสั้นฉันตั้งใจจะสร้างผลลัพธ์ของซอฟต์แวร์ SpatialKey ( ภาพหน้าจอ ) อีกครั้ง

ฉันใช้อินสแตนซ์ GeoServer เพื่อโฮสต์ข้อมูลสำหรับแผนที่ แต่ฉันจะรวมการฉายภาพและการประสานงานข้อมูลที่ต้องการลงในภาพได้อย่างไรเพื่อที่จะจัดตำแหน่งให้ถูกต้องกับตำแหน่งทางภูมิศาสตร์ที่แท้จริง Heatmap ความหนาแน่นถูกสร้างขึ้นโดยใช้ข้อมูลเชิงพื้นที่

แพ็คเกจแรสเตอร์ช่วยให้คุณสร้างแรสเตอร์แบบสุ่ม (ขนาดการฉายภาพ) หรือใช้วัตถุที่มีอยู่เช่นเมทริกซ์ จากนั้นคุณสามารถมอบหมายการฉาย & ขอบเขต image()ฟังก์ชั่นของคุณใช้อาร์กิวเมนต์เมทริกซ์ซึ่งสามารถใช้โดยตรง:

library(raster)    
## Create a matrix with random data & use image()
xy <- matrix(rnorm(400),20,20)
image(xy)

# Turn the matrix into a raster
rast <- raster(xy)
# Give it lat/lon coords for 36-37°E, 3-2°S
extent(rast) <- c(36,37,-3,-2)
# ... and assign a projection
projection(rast) <- CRS("+proj=longlat +datum=WGS84")
plot(rast)

จากนั้นคุณสามารถใช้writeRaster()กับรูปแบบจำนวนเท่าใดก็ได้เช่น

writeRaster(rast, "~/myraster.asc", format = "ascii")

มีหลายวิธีที่แตกต่างกัน แต่การแทงที่คุณต้องการใช้เครื่องมือโอเพ่นซอร์สบน Windows ฉันสามารถแนะนำ 3 ตัวเลือก:

1. ใช้เครื่องมือ gdal_translate ในGDAL (ไลบรารีการจัดการทางภูมิศาสตร์แบบโอเพนซอร์ส) คุณอาจต้องการ-gcp pixel line easting northing elevationพารามิเตอร์ที่พิกเซลและเส้นเป็น x / y ใน raster ( http://www.gdal.org/gdal_datamodel.html ) และ easting / northing พิกัดทางภูมิศาสตร์ คุณจะต้องระบุการฉายด้วย-a_srs srs_def
2. ใช้การโยง rเพื่อ GDAL และทำสิ่งเดียวกัน
3. วิธี "GUI" ที่ง่ายขึ้นอาจจะเป็นไปตามบทช่วยสอนนี้และใช้ QGIS กับ GDAL ข้อเสียของสิ่งนี้คือคุณต้องตั้งค่าทั้ง QGIS และ GDAL หากคุณยังไม่ได้ตั้งค่า

ฉันจะดูที่การพึ่งพาของแต่ละบิตของซอฟต์แวร์ที่คุณต้องการอย่างรอบคอบก่อนที่จะลงเส้นทางใดก็ตามที่คุณเลือก