ascii - In R create maximum raster from mulitple raster files -

-

i kind of new r maybe stupid question cant figure out myself.

this problem. have multiple asc files same gridsize , cover same area. want maximum value each grid asc files. tried multiple things: 

   (i in 1:144){     asc0<-rastertopoints(raster(asc0))         asc1<-rastertopoints(raster(asc[i]))         asc0[,3] <-pmax(asc0[,3], asc1[,3])     }

this 1 fails when loop threw files, becasue leaves out na, asc0 (my base file) different size next file asc1[2].

does know way this? have loops ready go through files, there 13x144 files. cant figure out way maximum value, store it, , compare next file.

thanks, appriciate help!!!

use stack , max:

r1 <- r2 <- r3 <- raster(nrows=10, ncols=10, xmn=0, xmx=10, ymn=0, ymx=10) r1[] <- 1   # raster 1 r2[] <- 2   # raster 2 r3[] <- 3   # raster 3 s1 <- stack(r1, r2, r3) s1 #class       : rasterstack  #dimensions  : 10, 10, 100, 3  (nrow, ncol, ncell, nlayers) #resolution  : 1, 1  (x, y) #extent      : 0, 10, 0, 10  (xmin, xmax, ymin, ymax) #coord. ref. : +proj=longlat +datum=wgs84 +ellps=wgs84 +towgs84=0,0,0  #names       : layer.1, layer.2, layer.3  #min values  :       1,       2,       3    <- values between 1 , 3 #max values  :       1,       2,       3    <- values between 1 , 3  max(s1) #class       : rasterlayer  #dimensions  : 10, 10, 100  (nrow, ncol, ncell) #resolution  : 1, 1  (x, y) #extent      : 0, 10, 0, 10  (xmin, xmax, ymin, ymax) #coord. ref. : +proj=longlat +datum=wgs84 +ellps=wgs84 +towgs84=0,0,0  #data source : in memory #names       : layer  #values      : 3, 3  (min, max)   <- 3

Comments

Post a Comment

Popular posts from this blog

python - No exponential form of the z-axis in matplotlib-3D-plots -

-
Image
i have similar problem described in how prevent numbers being changed exponential form in python matplotlib figure : i don't want (in special case) weird scientific formatting of axis. problem different have problem @ z -axis. 2-d plots can use ax.get_yaxis().get_major_formatter().set_useoffset(false) . , there no function ax.get_zaxis() what use format z -axis same way? edit: example: from mpl_toolkits.mplot3d import axes3d import numpy np import sys import matplotlib import matplotlib.pyplot pyplot def func(xi, ti): res = 10e3 + np.cos(ti) * np.sin(xi) return res if __name__ == '__main__': timespacing = 20 timestart = 0 timeend = 1 time = np.linspace(timestart, timeend, timespacing) widthspacing = 50 widthstart = 0 widthend = 3 width = np.linspace(widthstart, widthend, widthspacing) reslist = [] matplotlib.rcparams['legend.fontsize'] = 10 fig = pyplot.figure() ax = fig.gca(projection = '3d...
Read more

php - Best Light server (Linux + Web server + Database) for Raspberry Pi -

-
i install web server database on raspberry pi (little computer). computer has 1gb ram. i want know best combination: linux distribution , web server , dbms run local server multiple users minimal latency, use php on server. , best settings performance , not have bugs (memory usage, disable plugin, disable service, etc)? i thought light debian , lighttpd server , sqlite database. solution? i think lighttpd + sqlite great choice. linux distro, debian @ centos or tiny core linux , although i'm not sure of compatibility pi. obviously, can't go wrong raspbian if want use in production , more stable performance, few more pi's , set them in cluster .
Read more

c# - "Newtonsoft.Json.JsonSerializationException unable to find constructor to use for types" error when deserializing class -

-
i developing app in xamarin android. have splash screen serialize class , pass mainactivity using intent . when try deserialize in mainactivity error message : "serializationexception unable find constructor use types" serialization : void loaddata() { currency currency = new currency(this); intent = new intent(this, typeof(mainactivity)); intent.putextra("currency", newtonsoft.json.jsonconvert.serializeobject(currency)); startactivity(intent); } deserialization : cur = newtonsoft.json.jsonconvert.deserializeobject<currency> (intent.getstringextra("currency")); // error here class constructor : public currency(context context) { thiscontext = context; } i began experiencing problem after added parameter context context constructor. has caused this? how can solve it? possible serialize/deserialize class has parameters? you need have empty constructor in currency class in order deserialize it. ...
Read more