Dear all, I'm working with data in the form of an ascii grid (British Ordnance survey grid). I've plotted the data using image.plot with success and would like to draw a boundary around the region so that values which are zero and the white plotting region can be easily distinguished. Has anyone any experience with this? Appreciate some guidance... Thanks Muhammad
AsciiGrid - draw boundary
2 messages · Muhammad Rahiz, Zia Ahmed
Hi all, After CSGS realization, I like to calculate proportion of simulated values exceeding a particular thresh values.(Suppose in a prediction grid, out of ten realizations, in five realizations- simulated value exceed a target threshold value, here proportion is 0.5). I want all calculate this proportion in all by prediction grid. If some help me write the code - it will help me a lot. Here meuse data example: Thanks library(gstat) data(meuse) data(meuse.grid) coordinates(meuse) = ~x + y coordinates(meuse.grid) = ~x + y # Variogram log Zn lzn.vgm = variogram(log(zinc) ~ 1, meuse) lzn.fit = fit.variogram(lzn.vgm, model = vgm(1, "Sph", 900, 1)) #Conditional simulation nsim = 10 lzn.sim = krige(log(zinc) ~ 1, meuse, meuse.grid, model = lzn.fit, nmax = 30, nsim = nsim) # target values. ########### target.zn <- log(100) -------------- next part -------------- A non-text attachment was scrubbed... Name: zua3.vcf Type: text/x-vcard Size: 281 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100319/f201c892/attachment.vcf>