[R-sig-Geo] Help: PROBABILITY INTERVAL maps in gstat

I like to create a probability interval map of soil salinity (for example from 2 to 4 dS/cm) using IK or from ordinary kriging, How do I wrote the code in gstat? Help will be appreciated. Thanks Zia <>___ R-sig-Geo mailing list R-sig-Geo@s

[R-sig-Geo] Help: control font size of spplot of names.attr

Is it possible to change the font size of attribute name in sppot? also the background color? Thanks Zia spplot(zn, c("a", "b", "c", "d"), names.attr = c("ordinary kriging", "universal kriging with dist to river", "stratified kriging with flood freq", "inverse distance"), as.table = TRUE

[R-sig-Geo] Help: CDF plot

I want to write a R- function to make a CDF plot of all realizations (sim1, sim2, sim3) with CDF of observed data (suppose zinc of muse data). Help will be appreciated Thanks Zia library(gstat) data(meuse) data(meuse.grid) coordinates(meuse) = ~x + y coordinates(meuse.grid) = ~x

[R-sig-Geo] Help: Merging several geotiff files

Hi, I have several DEM sheets as geotiff format. I like to merge/mosaic them to create a grid file using a R-function. Help will be appreciated. Zia <>___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig

[R-sig-Geo] How do I calculate the proportion of simulated values exceeding a particular thresh values?

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 cal

Re: [R-sig-Geo] Help: R-function to plot all CSGS realization variograms in single a plot

, lzn.fit) } plot(do.call(rbind, out), m, layout=c(10,10), skip = FALSE, scales = list(y = list(relation = "same"))) Zia Ahmed wrote: Hi, I am trying to plot 100 variograms from 100 CSGS realizations in a single plot with the observed variogram model. I used following codes to p

[R-sig-Geo] Help: R-function to plot all CSGS realization variograms in single a plot

Hi, I am trying to plot 100 variograms from 100 CSGS realizations in a single plot with the observed variogram model. I used following codes to plot them together. But problem is this code is: I have to create variogram for all realizations, ie 100 times I have to repeat this code! Thi

Re: [R-sig-Geo] Box-cox transformation and Conditional Sequential Gaussian Simulation!

wrote: Yes, you can. It is less flexible than the normal score transform -- try for instance transforming a bimodal distribution -- but has the advantage that the back-transform is straightforward. How do you for example use the inverse normal score transformation simulated values outside the ra

[R-sig-Geo] Box-cox transformation and Conditional Sequential Gaussian Simulation!

I have a basic question related to the data transformation CSG simulation. Can I use Box-cox transformation of the data instead of normal score transformation before performing CSG simulation? Suggestion will be appreciated. Thanks ZIa ___ R-sig-G

[R-sig-Geo] transparent polygon outline Re: Shows error when I try to write shape file-Re: Combing attribute table with polygon??

length(m)) colnames(m) = c('REGION', 'v') # merge x = merge(d, m, by='REGION', all.x=TRUE) # sort x = x[order(x$sort),] # put the data back wrld_si...@data = x # now it works col = terrain.colors(nrow(m)) plot(wrld_simpl, col=col[wrld_simpl$v]) On Wed, Jan 27,

[R-sig-Geo] Shows error when I try to write shape file-Re: Combing attribute table with polygon??

quot; [11] "SAR""CaCO3" "SOM""BD" "AWC" [16] "KSAT" "Kf" "Kw" "LL" "LEP" >polyShp <- "D:/test/mapunit_poly.shp" >

Re: [R-sig-Geo] Combing attribute table with polygon??

Thank you so much Roberts! it works. Zia Robert J. Hijmans wrote: Zia, something like this might work as long as MUKEY has no NA values: poly.data = merge(poly.data, data.table, by="MUKEY" ,all.x = TRUE) Robert On Wed, Jan 27, 2010 at 1:58 PM, Zia Ahmed wrote: I have create

[R-sig-Geo] Combing attribute table with polygon??

I have created one data table and one polygon in R. Both have a similar type of column. I want to combine or join them with this column or field ("MUKEY"). Then I will convert this polygon as a ESRI shape file. But when I used spCbind, I got error. I think i miss something here. Is

Re: [R-sig-Geo] Helping with write a function. Merging shape file with all fields?

,files[i])) n <- length(slot(temp.data, "polygons")) temp.data <- spChFIDs(temp.data, as.character(uid:(uid+n-1))) uid <- uid + n poly.data <- spRbind(poly.data,temp.data) } writeOGR(poly.data, dsn=combinedShp, layer="combined", driver="ESRI

Re: [R-sig-Geo] Helping with write a function. Merging shape file with all fields?

ns")) temp.data <- spChFIDs(temp.data, as.character(uid:(uid+n-1))) uid <- uid + n poly.data <- spRbind(poly.data,temp.data) } writeOGR(poly.data, dsn=combinedShp, layer="combined", driver="ESRI Shapefile") Matt Beard Purdue University On Fri, Jan

[R-sig-Geo] Helping with write a function. Merging shape file with all fields?

Hi Matt: With help of Jim Holtman (my R Guru!), I am able to modified your code for selecting particular type of shape files from different folders and able to combine them together. Thank you so much! However, final polygon does not contain any filed (suppose county name, FIPS code etc

Re: [R-sig-Geo] Helping with write a function. Merging shape files?

the R 'plot' function order the polygons before displaying them? Is it possible to view the code for a function like plot(SpatialPolygons)? Matt On Wed, Jan 20, 2010 at 4:14 PM, Zia Ahmed <mailto:z...@cornell.edu>> wrote: Thanks! I am trying to write a

Re: [R-sig-Geo] Helping with write a function. Merging shape files?

s,slot(data.temp, "polygons")) the SoutheasternStates polygons completely covered up the Florida polygons in QGis, but the Florida polygons showed up when plotting in R: plot(spatialPolygons) Does the R 'plot' function order the polygons before displaying them? Is it possible to vi

[R-sig-Geo] Helping with write a function. Merging shape files?

one helps me to write correct code for solving this it will be great. Zia Hannes Reuter wrote: Dear Zia, 1)check GDAL (www.gdal.org) using the ogr2ogr function -append as well as the RGDAL library out. 2) 3) -->gdal_rasterize --> check as well the RGDAL library out. make sure you

[R-sig-Geo] Merging shape files?

Hi: I like use R for data mining, organization, manipulation and mapping of SSURO data of USA. I have already done first two; I need help for mapping these data.in R It will be appreciated if someone helps me to solve following three issues: *1. How do I merge several shape files?* I

[R-sig-Geo] Resent: colorkey: Placement of labels and ticks in center positions of "color key" of a levelplot ?

Is it possible to place labels and ticks in center positions of each cut of "color key" of levelplot ? Here is my code that I used to produce a map (see attachment). The two labels and ticks were placed in left sides each segment of the color key. I want to show them center of the cut

[R-sig-Geo] Risk mapping: Parametric vs Non-Parametric Conditional Simulation approaches?

Hi, I like to compare areas not exceeding of particular threshold values between parametric and non-parametric approaches of 1000 realizations. For non-parametric approach, I did IK for 1000 realizations (conditional simulation) and calculated the area not exceeding (0.5 probability) of the

Re: [R-sig-Geo] plot shape file using lattice xyplot

o use xyplot, you can use the panel function that sp defines: sp.polygon You can read in the shapefile with the RGDAL package: shp <- readOGR(...) 2009/12/17 Zia Ahmed : Hi, I am trying to plot a shape file in R. I can do it easily using plot or spplot function. But I want plot this

[R-sig-Geo] plot shape file using lattice xyplot

Hi, I am trying to plot a shape file in R. I can do it easily using plot or spplot function. But I want plot this map with lattice xyplot function. Any one has any idea? I do not want use "maps" package. Help will be appreciated... Thanks Zia library(maptools) library(lattice) bound<-r

Re: [R-sig-Geo] Moving window correlation- Further

-0.3807374 -0.3576051 -0.2260482 0.05467783 -0.2882067 -0.2124191 -0.08551635 3 -0.3943537 -0.3750882 -0.2435564 0.03511977 -0.3081105 -0.2300869 -0.10426607 4 0.324 0.349 0.442 0.4815 0.047 0.104 0.2347 Jaime R. Garcia M. wrote: On Tue, 2009-12-08 at 21:54 +01

[R-sig-Geo] Moving Windows Correlation----???

.99972521 2 -0.3807374 -0.3576051 -0.2260482 0.05467783 -0.2882067 -0.2124191 -0.08551635 3 -0.3943537 -0.3750882 -0.2435564 0.03511977 -0.3081105 -0.2300869 -0.10426607 4 0.324 0.349 0.442 0.4815 0.047 0.104 0.2347 Jaime R. Garcia M. wrote: On Tue, 2009-12-0

[R-sig-Geo] Moving window correlation

I have two grid maps (100 m cell size). I want to a create another raster gird of correlation coefficient ( r value) of two maps using moving windows approach. Has someone any idea how to do it in R? Thanks Zia <>___ R-sig-Geo mailing list

[R-sig-Geo] Memory issue in R gstat: Windows vs Linux (Ubuntu)

If someone facing problem to analyze large data set in gstat in Windows due to memory issue, they can try in Ubuntu using "R-gedit". It is very easy to install Ubuntu in a Windows machine (as a dual booting system- not as a guest OS!). After installing r with all necessary packages user

Re: [R-sig-Geo] Risks as hazardous (false positive) and safe (false negative)

willing to accept the indicator kriging values as estimates of these probabilities, you're essentially done. Of course, the kriged values can be outside [0, 1], so you need to deal with that. Best regards, -- Edzer Zia Ahmed wrote: Dear Edzer, Thanks for your replay. I have done IK o

[R-sig-Geo] How do change projection from geographic coordinate to projected coordinate system in R?

Hi all, I want project several several raster(.tif) and vector (.shp) files in R using rgdal driver. One set of files are in geographic coordinate system (WGS84) other are in LCC-Bangladesh coordinate system. .I want to re-project them to EPSG:3104 (Gulshan 303/ BTM). If some help me t

Re: [R-sig-Geo] Overlay of SpatialGrid/Raster and polygons

I am trying to install "raster" package, but it is not available in this site. http://R-Forge.R-project.org . Zia Robert J. Hijmans wrote: Hi Ned, Here is an approach to get values from a RasterStack to all cells in each polygon: library(raster) # a polygon data(meuse.riv) pol <- Spatial

Re: [R-sig-Geo] Risks as hazardous (false positive) and safe (false negative)

these probabilities, you're essentially done. Of course, the kriged values can be outside [0, 1], so you need to deal with that. Best regards, -- Edzer Zia Ahmed wrote: Dear Edzer, Thanks for your replay. I have done IK of water arsenic using quantiles of the data. I want to do miss classifi

[R-sig-Geo] Risks as hazardous (false positive) and safe (false negative)

Hi all, Is it possible in gstat to get misclassification of risks as hazardous (false positive) and safe (false negative) as describe by Goovaerts (1997) from parametric or non-parametric probability mapping?. help will appreciated. Thanks Zia <>

Re: [R-sig-Geo] bbox 4-dimensional objects - suitable method?

Hi all, Is it possible in gstat to get misclassification of risks as hazardous (false positive) and safe (false negative) as describe by Goovaerts (1997) from parametric or non-parametric probability mapping?. help will appreciated. Thanks Zia <>_

[R-sig-Geo] Parametric/Non parametric Probability Kriging: get misclassification of risks as hazardous (false positive) and safe (false negative) in gstat?

variance of each simulations? I have limited knowledge in geostatistic as well as R. Thanks Zia Ahmed ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo