R: วิธีรับละติจูดและลองจิจูดจาก RasterLayer

ฉันเป็นผู้เริ่มต้นที่สมบูรณ์ของข้อมูลทางภูมิศาสตร์ดังนั้นโปรดยกโทษให้ฉันหากคำถามไม่เหมาะสม

ฉันดาวน์โหลดข้อมูลจาก NCDC NARR และจัดการโหลดลงใน R โดยใช้rasterแพ็คเกจ ฉันต้องการรับรายการที่มีละติจูดลองจิจูดและค่า ฉันเข้าใจว่าrasterToPoints()ควรทำสิ่งที่ฉันต้องการอย่างแน่นอนอย่างไรก็ตามค่าละติจูดและลองจิจูดของฉันดูแปลก ๆ :

r <- raster(myfile)
data_matrix <- rasterToPoints(r)
head(data_matrix)
            x       y value
[1,] -5405401 4347242    70
[2,] -5372938 4347242    88
[3,] -5340475 4347242    76
[4,] -5308012 4347242    85
[5,] -5275549 4347242    87
[6,] -5243086 4347242    88

ฉันคิดว่าฉันควรทำอะไรกับการฉายภาพซึ่งปัจจุบันเป็น Lambert Conformal Conic (LCC) นี่คือข้อมูลเพิ่มเติมเกี่ยวกับแรสเตอร์

> r
class       : RasterLayer 
dimensions  : 277, 349, 96673  (nrow, ncol, ncell)
resolution  : 32463, 32463  (x, y)
extent      : -5648874, 5680713, -4628777, 4363474  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=lcc +lat_1=50 +lat_2=50 +lat_0=50 +lon_0=-107 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs 
data source : mypath-to-file
names       : value

ฉันต้องทำอย่างไรเพื่อให้ได้ค่าละติจูดและลองจิจูดในสหรัฐอเมริกาที่แท้จริง

คุณต้องปฏิเสธภาพแรสเตอร์ลงในโปรเจคชันทางภูมิศาสตร์ (องศาทศนิยม) โดยใช้ "projectRaster" หรือ "spTransform" ดูคำจำกัดความ sp CRS ที่ระบุสตริงการฉายที่คุณต้องการ ตัวอย่างในวิธีใช้สำหรับ "projectRaster" ค่อนข้างชัดเจนในวิธีการทำเช่นนี้

หากคุณบีบบังคับข้อมูลแรสเตอร์ของคุณให้เป็นวัตถุ SpatialPointsDataFrame คุณจะต้องใช้ "spTransform" และดึงพิกัดจากช่อง @coordinates และเพิ่มลงใน data.frame ในช่อง @data นี่คือตัวอย่างของสิ่งที่จะมีลักษณะ

library(raster)
library(rgdal) # for spTransform

# Create data
r <- raster(ncols=100, nrows=100)
  r[] <- runif(ncell(r))
  crs(r) <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84"
  projection(r)

# Convert raster to SpatialPointsDataFrame
r.pts <- rasterToPoints(r, spatial=TRUE)
  proj4string(r.pts)

# reproject sp object
geo.prj <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0" 
r.pts <- spTransform(r.pts, CRS(geo.prj)) 
  proj4string(r.pts)

# Assign coordinates to @data slot, display first 6 rows of data.frame
r.pts@data <- data.frame(r.pts@data, long=coordinates(r.pts)[,1],
                         lat=coordinates(r.pts)[,2])                         
head(r.pts@data)

ฉันควรทราบว่ามันไม่ใช่วิธีปฏิบัติที่ดีในการแปลง rasters เป็นคลาสเวกเตอร์วัตถุและคัดค้านข้อดีของแพ็คเกจ raster ที่ให้การประมวลผลที่ปลอดภัยของหน่วยความจำ บ่อยครั้งที่รอบคอบที่จะคิดถึงปัญหาของคุณและประเมินว่าคุณกำลังเข้าใกล้อย่างถูกต้องหรือไม่ หาก OP ได้ให้บริบทว่าทำไมพวกเขาต้องการพิกัด [x, y] สำหรับทุกเซลล์ชุมชนฟอรัมอาจให้ทางเลือกในการคำนวณที่จะทำให้เกิดปัญหาในสภาพแวดล้อมแบบแรสเตอร์

รับพิกัดของศูนย์เซลล์และสร้างวัตถุเชิงพื้นที่:

spts <- rasterToPoints(r, spatial = TRUE)

เปลี่ยนคะแนนเป็นเป้าหมายที่คุณต้องการ:

library(rgdal)
llprj <-  "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"
llpts <- spTransform(spts, CRS(llprj))

ค่าจะถูกคัดลอกเป็นคอลัมน์ใน SpatialPointsDataFrame นี้

print(llpts)

ตอนนี้ให้เสร็จเพื่อรับ data.frame:

x <- as.data.frame(llpts)

มีการใช้งานทั่วไปของสิ่งนี้ในแพ็คเกจ SGAT ดูฟังก์ชั่นlonlatFromCellที่นี่:

https://github.com/SWotherspoon/SGAT/blob/master/R/Raster.R

ดูเหมือนว่าคุณมีพิกัดที่คาดการณ์ไว้ที่นั่น (ไม่ใช่ละติจูด / ลองจิจูดหรือที่รู้จักกันในพิกัด GCS) อาจไม่ชัดเจนสำหรับคุณที่เป็นปัญหา ดูโพสต์นี้ การแปลงระบบพิกัดทางภูมิศาสตร์ใน R