Cách gán CRS cho RasterLayer trong R

Tôi có một lớp với coord. ref. : NA.

Làm thế nào tôi có thể gán coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0đến coord. ref. : NA?

class       : RasterLayer 
dimensions  : 127, 146, 18542  (nrow, ncol, ncell)
resolution  : 0.01, 0.01  (x, y)
extent      : -122.39, -120.93, 36.72, 37.99  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif 
names       : X2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc 
values      : 0, 0.0782172  (min, max)


class       : RasterLayer 
dimensions  : 180, 240, 43200  (nrow, ncol, ncell)
resolution  : 0.01, 0.01  (x, y)
extent      : -123.305, -120.905, 36.705, 38.505  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +to wgs84=0,0,0 
data source : C:\Users\nahmlee\Documents\1998Test.tif 
names       : X1998Test 
values      : 0, 0.07884484  (min, max)

— Không
nguồn

Sử dụng crschức năng. Nếu rlà raster của bạn:

crs(r) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"

Tất nhiên, bạn cần chắc chắn rằng đó là hình chiếu chính xác cho raster của bạn.

Kết quả là như nhau, nhưng phiên bản này dễ thích nghi hơn với các hệ thống tham chiếu khác: crs(r) <- CRS('+init=EPSG:4326')nếu bạn biết số EPSG.

— Julian

Điều này nên được CRS("+init=epsg:28992")cho rgdal.

— mace