R में RasterLayer को CRS कैसे असाइन करें

Question

मेरे पास एक परत है coord. ref. : NA।

मैं कैसे प्रदान कर सकते हैं coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0करने के लिए 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)

Answer 1

crsफ़ंक्शन का उपयोग करें । यदि rआपका रेखापुंज है:

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

बेशक, आपको यह सुनिश्चित करने की आवश्यकता है कि यह आपके रास्टर के लिए सही प्रक्षेपण है।

Comments

परिणाम समान है, लेकिन यह संस्करण अन्य संदर्भ प्रणालियों के लिए आसान अनुकूल है: crs(r) <- CRS('+init=EPSG:4326')यदि आप ईपीएसजी नंबर जानते हैं।

— जूलियन

यह होना चाहिए CRS("+init=epsg:28992")के लिए rgdal।

— गदा