[R-sig-Geo] zscore for krige.cv with IDW

2009-12-09 Thread Tobin Cara
Hello, I am trying to compare the zscores of using various kriging methods with inverse distance weighting. When I run crossval_idw - krige.cv(Prec00hr~1, DataCoordhr, nfold=2) I get all NA's for the zscore. As the residual is a valid number, I am assuming that the standard deviation is not

[R-sig-Geo] transforming a 3d SpatialGridDataFrame into a mesh of some kind

2009-12-09 Thread Mark Connolly
I have a 3d SpatialGridDataFrame I created using krige (IDW). I am using rgl facilities to create an interactive plot of the points, but I'd really like to have something more volumetric than points. My copy of the truly wonderful Applied Spatial Data Analysis with R (shameless pander) is

Re: [R-sig-Geo] convert KML file into shape file or CSV file

2009-12-09 Thread Barry Rowlingson
On Wed, Dec 9, 2009 at 1:22 AM, rusers.sh rusers...@gmail.com wrote: Hi all,  Does anybody know how to import a KML file in R and them save it as a shape or CSV file. Take the following KML file as an example, http://www.nature.com/nature/googleearth/avianflu1.kml  I have tried several

[R-sig-Geo] Rshhs() import problem

2009-12-09 Thread Maik Renner
Dear R-Geographers, I am trying to use the GSHHS shoreline data for annotating my maps. I used Rgshhs() from the package maptools to import it into R as a Spatial object. However it seems that there is a problem due to the new version of gshhs_2.0.zip. In the help gshhs_1.11.zip is referred.

[R-sig-Geo] reading select values from a binary file

2009-12-09 Thread Jonathan Thayn
Hopefully there is a simple solution to this problem. I process really large satellite images in R and, in order to not offend the memory-limitation demi-gods, I open each row of the imagery one at a time, process it, write the results, and then open the next line. I use readBin() to open each

Re: [R-sig-Geo] reading select values from a binary file

2009-12-09 Thread Barry Rowlingson
On Wed, Dec 9, 2009 at 4:58 PM, Jonathan Thayn jth...@ilstu.edu wrote: Hopefully there is a simple solution to this problem. I process really large satellite images in R and, in order to not offend the memory-limitation demi-gods, I open each row of the imagery one at a time, process it,

Re: [R-sig-Geo] transforming a 3d SpatialGridDataFrame into a mesh of some kind

2009-12-09 Thread Robert J. Hijmans
Mark, If spgdf is a SpatialGridDataFrame object # install.packages(raster, repos=http://R-Forge.R-project.org;) library(raster) plot3D(raster(spgdf)) Or you can have a look at the code: getMethod(plot3D, 'RasterLayer') Robert On Wed, Dec 9, 2009 at 5:10 AM, Mark Connolly mark_conno...@acm.org

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

2009-12-09 Thread Zia Ahmed
Dear Garcia, Thank you so much! Before running my data, I have tested your code using meuse data. It works fine. But it produce following output like below from V1 to V3103. I do not understand what does the values in rows mean (1 to 4). Here is the code that I used. I have questions How do

Re: [R-sig-Geo] transforming a 3d SpatialGridDataFrame into a mesh of some kind

2009-12-09 Thread Mark Connolly
I _think_ raster may be solving a different problem. I am looking for a voxel-type raster. Let me know if I am wrong. plot3d(raster(kriged_sp)) Error in as.double(x) : cannot coerce type 'S4' to vector of type 'double' Let me describe the SpatialGridDataFrame a bit better. The coordinates

Re: [R-sig-Geo] reading select values from a binary file

2009-12-09 Thread Barry Rowlingson
On Wed, Dec 9, 2009 at 6:20 PM, rick reeves ree...@nceas.ucsb.edu wrote: Hello All: Faced with a similar challenge, and NOT wanting to resort to writing a C language function employing fseek(), I just used two readBin calls: One to read (and implicitly discard) data up to the spot I actually

Re: [R-sig-Geo] Rshhs() import problem

2009-12-09 Thread Roger Bivand
On Wed, 9 Dec 2009, Maik Renner wrote: Dear R-Geographers, I am trying to use the GSHHS shoreline data for annotating my maps. I used Rgshhs() from the package maptools to import it into R as a Spatial object. However it seems that there is a problem due to the new version of gshhs_2.0.zip. In

Re: [R-sig-Geo] Spatial R Courses

2009-12-09 Thread Tomislav Hengl
David Rossiter runs this excellent distance education course: http://www.itc.nl/Pub/Study/Courses/C10-AES-DE-02 *Next one starts in 6 weeks! We will run a 5-day trainig course in June (which is more affordable, but shorter): 27 June - 4 July 2010, GEOSTAT 2010 Summer School for PhD

Re: [R-sig-Geo] reading select values from a binary file

2009-12-09 Thread Michael Sumner
If you can read the whole file at once it might be better to push the results into a matrix for indexing in R. This could also be useful for reading chunks of the file, and processing in parts. seek is fine but if you are moving around for three values each read it will be slow. If the

Re: [R-sig-Geo] convert KML file into shape file or CSV file

2009-12-09 Thread rusers.sh
It is the problem of the file. I should have thought of this because i have tried several tools. Thanks. Based on your replies, i think it is easy to convert the files between KML and SHP if the files are correct and have simple structure. We only need readOGR and writeOGR in rgdal package to

Re: [R-sig-Geo] convert KML file into shape file or CSV file

2009-12-09 Thread Alex Mandel
Almost, in my sleepiness I left out an important part from that pseudo code. You need to store the data in a variable/object and then write it out. importeddata - readOGR(pathtofile,layer.kml) writeOGR(importeddata ,pathtooutput,driver=ESRI Shapefile,layer=output.shp) Make sure to read the help

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

2009-12-09 Thread Jaime R. Garcia M.
On Wed, 2009-12-09 at 11:46 -0500, Zia Ahmed wrote: Thank you so much! Before running my data, I have tested your code using meuse data. It works fine. But it produce following output like below from V1 to V3103. I do not understand what does the values in rows mean (1 to 4). Here is the

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

2009-12-09 Thread Zia Ahmed
Thank you so much! Before running my data, I have tested your code using meuse data. It works fine. But it produce following output like below from V1 to V3103. I do not understand what does the values in rows mean (1 to 4). Here is the code that I used. I have questions How do I plot

[R-sig-Geo] creating cluster with between-points arbitrary distances

2009-12-09 Thread Alexandre VILLERS
Security, version of virus signature database 4674 (20091209) __ The message was checked by ESET Mail Security. http://www.eset.com ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] creating cluster with between-points arbitrary distances

2009-12-09 Thread Dan Putler
Hi Alexandre, Greetings from yesterday evening. I think I know what you want, and I've been working on the same problem in a very different substantive context. To make sure I understand the question, you are basically stating that if the gap between nearest birds of the same species exceeds

Re: [R-sig-Geo] creating cluster with between-points arbitrary distances

2009-12-09 Thread Rainer M Krug
On Thu, Dec 10, 2009 at 9:23 AM, Dan Putler dan.put...@sauder.ubc.cawrote: Hi Alexandre, Greetings from yesterday evening. I think I know what you want, and I've been working on the same problem in a very different substantive context. To make sure I understand the question, you are

Re: [R-sig-Geo] creating cluster with between-points arbitrary distances

2009-12-09 Thread Alexandre VILLERS
Niort Phone +33 (0)5 49 09 96 13 Fax +33 (0)5 49 09 65 26 __ Information from ESET Mail Security, version of virus signature database 4674 (20091209) __ The message was checked by ESET Mail Security. http://www.eset.com ___ R-sig