आर में रेखापुंज ग्रिड ग्रिड


10

मैं आर में दो रेखापुंज ग्रिड को संरेखित करने की कोशिश कर रहा हूं। एक बार संरेखित करने के बाद मैं उन्हें एक साथ जोड़ना चाहूंगा।

मैंने यह जांचने की कोशिश की है कि क्या stackकाम करना होगा:

grid_snap <- stack(habi_sdw, Pop_sdw)

और मुझे निम्नलिखित त्रुटि मिलती है:

तुलनात्मकता (x) में त्रुटि: अलग-अलग सीमा

रेखापुंज ग्रिड में निम्नलिखित गुण होते हैं:

show(habi_sdw)
# class       : RasterLayer 
# dimensions  : 9187, 9717, 89270079  (nrow, ncol, ncell)
# resolution  : 0.00892857, 0.00892857  (x, y)
# extent      : -28.83706, 57.92186, -36.02464, 46.00214  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs 
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_172902_12860_17067.grd 
# names       : layer 
# values      : 0, 333707.6  (min, max)

show(Pop_sdw)
# class       : RasterLayer 
# dimensions  : 10143, 8858, 89846694  (nrow, ncol, ncell)
# resolution  : 0.008333333, 0.008333333  (x, y)
# extent      : -17.53524, 56.28143, -46.97893, 37.54607  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_170421_12860_12760.grd 
# names       : pop2010ppp 
# values      : 0, 128925.9  (min, max)

alignExtent()रेखापुंज पैकेज में उपयोग करना सही दृष्टिकोण नहीं है।

क्या मुझे फिर से तैयार करने की आवश्यकता है क्योंकि संकल्प थोड़े अलग हैं?

(0.00892857 x 0.00892857) बनाम (0.008333333 बनाम 0.008333333)

जवाबों:


15

यह प्रश्न समान है: डेटा निष्कर्षण और रिज़ॉल्यूशन परिवर्तन के साथ रेखापुंज द्वारा क्लिप रेखापुंज , लेकिन एक अलग कोण से आ रहा है। हालाँकि, मुझे लगता है कि उत्तर की संभावना समान है। सबसे पहले, यह चुनें कि आप कौन सा रेखापुंज निश्चित करना चाहते हैं। मैं अपना पिछला उत्तर आसानी से यहाँ दूंगा:

आवश्यक लाइब्रेरी लोड करें:

library(raster)
library(rgdal)

पढ़ें आपदाओं:

r1 = raster("./dir/r1.tif")
r2 = raster("./dir/r2.tif")

उसी ग्रिड में फिर से भरना:

r.new = resample(r1, r2, "bilinear")

यदि आवश्यक हो (मास्किंग के लिए), मिलान करने के लिए extents सेट करें:

ex = extent(r1)
r2 = crop(r2, ex)

हटाए गए डेटा जो कि आपदाओं में से एक के बाहर आते हैं (यदि आपको आवश्यकता है):

r.new = mask(r.new, r2)

आपके चीर हरण अब मेल खाते हैं।


0

हाँ। उनके आकार के समान होने के लिए आपको अपने चूहों को फिर से आकार देने की आवश्यकता है। आर खुद के साथ सौदा नहीं करता है। यह देखते हुए कि आपके किसी भी रैस्टर में पूरी तरह से दूसरे नहीं हैं, आपको अपने पसंदीदा रिज़ॉल्यूशन के साथ न्यूनतम-सीमा का रैस्टर बनाने पर विचार करना चाहिए, और फिर दूसरों से मिलान करने के लिए फिर से क्रॉप करना चाहिए।


2
कृपया अपना उत्तर विस्तृत करें, उदाहरण के लिए कुछ नमूना कोड प्रदान करके।
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.