স্ট্যাকড রাস্টার থেকে সর্বাধিক কক্ষের মান

Question

আমি কীভাবে সজ্জিত রাস্টার থেকে সর্বাধিক সেল মান খুঁজে পেতে পারি।

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

ভাল কাজ করে কিন্তু

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

এনএ দেয়।

অবশ্যই কোনও স্ট্যাকড রাস্টারে নেই।

আমার কোডটি এখানে:

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, সমস্ত স্তরের সর্বাধিক সেল মান, বৃদ্ধি) এর মতো একই স্কেল পেতে আমি সমস্ত স্তরগুলির সর্বাধিক সন্ধান করছি। ধন্যবাদ, নাহম

— নাহম

আমি এটি সেলস্ট্যাট দিয়ে পেয়েছি # জিওস্ট্যাট- কোর্স. org / সিস্টেমে / ফাইলস / লেউইস_টিউটর.এমপিডিএফ Rad1998.max <- সেলস্ট্যাটস (RAD1998.all, 'সর্বোচ্চ') Rad1998.all.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)) স্ট্যাকের সমস্ত রেস্টারদের সর্বোচ্চ মান সরবরাহ করবে।