स्टैक्ड रैस्टर से मैक्स सेल वैल्यू

Question

स्टैक्ड रैस्टर से मैं अधिकतम सेल वैल्यू कैसे पा सकता हूं।

Rmax <- maxValue(RAD1998.all[[1]]) 

ठीक काम करता है लेकिन

Rmax <- maxValue(RAD1998.all[[2]]) 

NA देता है।

बेशक एक स्टैक्ड रैस्टर में नहीं।

यहाँ मेरा कोड है:

RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = ""))
list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) 
RAD1998.all <- stack(list.ras)

Comments

क्या आप सभी परतों की अधिकतम या प्रत्येक परत की अधिकतम के लिए देख रहे हैं? वैसे भी, आप maxValueसही तरीके का उपयोग नहीं कर रहे हैं । सहायता पृष्ठ के अनुसार, आप बेहतर तर्क का उपयोग करते हैं... Additional argument: layer number (for RasterStack or RasterBrick objects)

मैं सभी परतों की अधिकतम के लिए देख रहा हूँ जैसे my.at <- seq (0, सभी परतों का अधिकतम सेल मान, वेतन वृद्धि)। धन्यवाद, नहम

— नह्म

मैंने इसे cellStats # geostat-course.org/system/files/lewis_tutorAM.pdf Rad1998.max <- cellStats (RAD1998.all, 'max' - Rad1998.all.max <- max (Rad1998.max) Rad1998.all के साथ मिला। .मैक्स

— नह्म

Answer 1

निम्नलिखित उदाहरण एक स्टैक में अधिकतम रेखापुंज मूल्य पर प्राप्त करने के दो तरीके दिखाता है। पहला उपयोग max()जो आपको अन्य उपयोगी जानकारी का एक मेजबान भी देता है। दूसरी विधि का उपयोग करता है maxValue(), जो स्टैक में दोनों चूहों का अधिकतम मूल्य देता है

library(raster)  

# Generate some georeferenced raster data
x = matrix(rnorm(400),20,20)
rast = raster(x)
extent(rast) = c(36,37,-3,-2)
projection(rast) = CRS("+proj=longlat +datum=WGS84")

y = matrix(rnorm(400),20,20)
rast2 = raster(y)
extent(rast2) = c(36,37,-3,-2)
projection(rast2) = CRS("+proj=longlat +datum=WGS84")

raster = stack(rast, rast2)

# Now run the statistics
max(raster) # Provides min, max and additional details  # Example 1

maxValue(raster)  # Gives both values                   # Example 2...
maxValue(raster)[[1]] # Gives first in stack max value
maxValue(raster)[[2]] # Gives second in stack max value

---

> maxValue(raster)  # Gives both values
[1] 2.688376 2.971443
> maxValue(raster)[[1]] # Gives first in stack max value
[1] 2.688376
> maxValue(raster)[[2]] # Gives second in stack max value
[1] 2.971443
> 
> max(raster) # Provides min, max and additional details
class       : RasterLayer 
dimensions  : 20, 20, 400  (nrow, ncol, ncell)
resolution  : 0.05, 0.05  (x, y)
extent      : 36, 37, -3, -2  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : -1.457908, 2.971443  (min, max)

Answer 2

max(MaxValue(raster)) स्टैक में सभी रैस्टर्स का अधिकतम मूल्य प्रदान करेगा।