[R-sig-Geo] Use estimateAnisotropy by intamap package

Dear All, I tried to inspect the geometric anisotropy by using the function of estimateAnisotropy with parameter "formulaString" in intamap, but never worked. Any help will be highly appreciated. Cheers, Yong -Original Message- From: r-sig-geo-boun...@stat.math.ethz.ch 代表 r-sig-geo-r

[R-sig-Geo] Carry out non-isotropic kriging

Dear ALL, I wondering that any expert here could show me some examples to compute directional sample variograms and fit/plot non-isotropic models and to do kriging using non-isotropic models by using gstat package. I can do the calculation of sample variograms at different directions, but not y

[R-sig-Geo] How to generate the probability map of a variable under a threshold using gstat?

Hi All, I want to produce a probability map of a soil variable under a threshold. I know I can make it in ArcGIS, but not sure using gstat. Any tips are highly appreciated. Regards, Yong Li ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https

[R-sig-Geo] FW: Interpolation option: IDW or OK?

al already. I appreciate you could give me your response as soon as possible. Regards Yong Li -Original Message- From: r-sig-geo-boun...@stat.math.ethz.ch [mailto:r-sig-geo-boun...@stat.math.ethz.ch] On Behalf Of r-sig-geo-requ...@stat.math.ethz.ch Sent: Tuesday, 10 February 2009 10:0

[R-sig-Geo] Interpolcation option: IDW or OK?

Many thanks for the message, Edzer. -Original Message- From: Edzer Pebesma [mailto:edzer.pebe...@uni-muenster.de] Sent: Monday, 9 February 2009 7:08 PM To: Yong Li Cc: r-sig-geo@stat.math.ethz.ch Subject: Re: FW: [R-sig-Geo] Interpolcation option: IDW or OK? Yong Li wrote: > Hi Ed

[R-sig-Geo] FW: Interpolcation option: IDW or OK?

: Edzer Pebesma [mailto:edzer.pebe...@uni-muenster.de] Sent: Monday, 9 February 2009 9:41 AM To: Yong Li Cc: r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] Interpolcation option: IDW or OK? Yong Li wrote: > Hi ALL, > > I have been with the attached dataset and R code to use OK to >

[R-sig-Geo] How to evaluate cross-validaiton of interpolation?

d and how and why? Cheers Yong Li [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] How to interpolate a raster-based point data?

appreciated. Regards Yong Li [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] FW: How to interpolate a points raster file (gridded)?

. Regards Yong Li ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] How to pick up cell values of grids using points?

Hi ALL, I want to pick up cells values of a grid (say a RS image) using points' positions. For example, I have an image and a point shape file with the same coordinate system, then I try to overlay points on the image, and want to get the cell values of the on-point cell and 8 surrounding cells

Re: [R-sig-Geo] how to load sp, gstate or other packages into R using R Dcom for standardalone applications?

Hi ALL, I had a problem in using r Dcom to develop applications for applying geostatistics using gstat. I transferred all code, which can run in R successfully, into VB, but I could not make them running. However, data management and simple graphing can be done without any problem. Any help and ti

[R-sig-Geo] How to estimate variables of a given vgm model?

Hi ALL, How can I calculate the corresponding semivariance or distance if a vgm model (model, nug, psill, range are known) and a value of either semivariance or distance are given? I am using gstat R package. Regards Yong ___ R-sig-Geo mailing list R-

Re: [R-sig-Geo] RE rgdal: Memory error!

ly 2008 12:30 AM To: Yong Li Cc: r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] RE rgdal: Memory error! Hi Yong, What operating system are you using (probably windows) and is it 32 bit of 64 bit and after what command did your computer crash. Without this kind of information it is hard

[R-sig-Geo] RE rgdal: Memory error!

Hi All, I ran the following code and got the below error: "Error: cannot allocate vector of size 30.3 Mb". Actually my laptop has 4096 MB ram, and I increased the limit already. Pls advise any tips to let me pass through. I guess Dr Hengl experienced this before. Regards Yong #--

Re: [R-sig-Geo] How to subset a spatial grid?

: Yong Li Cc: r-sig-geo@stat.math.ethz.ch Subject: Re: [R-sig-Geo] How to subset a spatial grid? Have you tried Map1 <- Map[1:100,1:100] ? Note that this works differently for SpatialGridDataFrames and SpatialPixelsDataFrames; the first knows about [rows,cols,bands], the second about [pix

[R-sig-Geo] How to subset a spatial grid?

<- Map$band1 Map$band1 <- NULL Map1 <- a subset of Map$va based on selected cols and rows. ????.. Cheers Yong Li ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Problem with universal kriging using gstat

Hi ALL, Can any expert to see my mistake in the following R script when I try the gstat universal kriging? The error occurs at the last step. It sounds I missed the setup for the names of coordinates of the SpatialGrid "data.grid". Regards Yong > memory.size() [1] 47930616 > memory.size(T

[R-sig-Geo] How to produce a R2 for semivariance model fitting

. Yong Li [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] How to remove the second order trend in gstat?

Dear all, I can remove the linear trend by using ~x+y in gstat, but have no idea to cope with the second order trend. Any help will be appreciated. Cheers, Yong Li ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman

[R-sig-Geo] (no subject)

require(maptools) #read in shape file pts <- readShapePoints(system.file("shapes/baltim.shp", package="maptools")[1]) #define the dimension of grids x.min <- summary(pts)$bbox[1] y.min <- summary(pts)$bbox[2] x.max <- summary(pts)$bbox[3] y.max <- summary(pts)$bbox[4] cellsize <- 20 x.n <- ro

[R-sig-Geo] How to randomly select points in a polygon

I want to overlay points with a polygon, get the points in the polygon and then try to just select a point randomly from them. I have tried using overlay, but it did not give me what I really want. Any help will be appreciated. Yong ___ R-sig-Geo mail

[R-sig-Geo] can't find the function of GE_SpatialGrid

Dear all, The following is the code I try to build grid from polygons, but was hammered by the error in the subject. Any idea helps me. Yong setwd("d:/Rcode/0_Zone") library(maptools) xx <- readShapePoly("Newshape.shp", IDvar="ID", proj4string=CRS("+proj=tmerc

[R-sig-Geo] huge problem with geoR application

Dear All, I have a huge problem with geoR. When I try to build a grid (1000 x 1000) and use krige to do kriging, I am told that R can not handle the vector size of 677mb. I have tried to increase R memory use allocation by -max-mem-size=2500M or memory.limit(size=2500M). It did not work either. C

[R-sig-Geo] ask for help in using geoR

Dear All, When I fit the models to the calculated semivariance using variofit, how can I know the performance of the fitting, say R2? Cheers Yong ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] remote sensing in R

Hi Everybody, Has anybody done or know the remote sensing (e.g. satellite image processing) package in R? Thanks. Yong ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Convert esri shape files to esri ascii grid from R

Hi R spatial experts, Has anybody done the conversion of esri shape files to esri ascii grid based on field values in R? I knew it can be done easily in Arcview using avenue. But if it can be done in R, that will fly my daily work loads. Regards Yong __

[R-sig-Geo] help on "non-square" problem when using write.asciigrid

Dear Friends, Can anybody help me for the subject? When I tried to use write.asciigrid to export a spatialgrid to esri ascii grid, I always failed because of the message "non-square grid is not support" even I knew my grid is 100% squared. Yong ___ R-