R：RasterLayerから緯度と経度を取得する方法は？

私は地理データの絶対的な初心者なので、質問が適切でない場合はご容赦ください。

NCDC NARRからデータをダウンロードし、rasterパッケージを使用してRにロードすることに成功しました。緯度、経度、値のリストを取得したいのですが。私はそれ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」を使用して、実際にラスタを地理的（10進度）投影に再投影する必要があります。目的の投影文字列を指定するCRS sp定義も参照してください。「projectRaster」のヘルプの例は、これを行う方法が非常に明確です。

ラスターデータをSpatialPointsDataFrameオブジェクトに強制する場合、「spTransform」を使用して@coordinatesスロットから座標を取得し、@ dataスロットのdata.frameに追加します。これがどのように見えるかの例を次に示します。

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)

ラスタをベクトルオブジェクトクラスに変換することはお勧めできません。メモリセーフな処理を提供するラスタパッケージの利点は無効になります。多くの場合、問題について本当に考えて、問題に正しく近づいているかどうかを評価するのが賢明です。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での地理座標系の変換