Re: [R-sig-Geo] Question re: Intamap package / Interpolate method: Inputs NOT both Spatial Points Data Frame?

n attribute to the Monarch Station Locations File: 'newcol' cbind(monarchstaswith...@data,newcol) } # Write the output file writePointsShape(monarchStasWithGDD,outFileName) } } #### On 10/26/2010 1:09 PM, Edzer Pebesma wrote: Are you

[R-sig-Geo] Question re: Intamap package / Interpolate method: Inputs NOT both Spatial Points Data Frame?

Hello: I note a conflict in the documentation for Intamap/interpolate(), need to know how to resolve: I have two SpatalPointsDataFrame objects: growingDegreeObservations <- readShapePoints("OhioGDDStationsGP.shp") # computed surface phenomena, for individual lat/lon points OrgStaLocs <- r

Re: [R-sig-Geo] reading KML files in R

___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805

Re: [R-sig-Geo] reading KML files in R

-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533

[R-sig-Geo] sp package/ spatial points operation works in R <= 2.8, fails for R >= 2.9 - resend with file download link

s. Question: What has changed in R 2.9 and later? Is this technique deprecated, and is there a better way of creating such a point list? Thanks, Rick R < Attached zip file contains the script and sample data set> -- Rick Reeves Scientific Programmer/Analyst and Data

[R-sig-Geo] Alternate methods: Polygon Algebra / Polygon Overlay with R Spatial objects..

at operate on the Spatial classes? If not, are there any alternatives to the PBSmapping methods? Thanks, Rick R -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 80

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

binary file, not a complete line. Say I want to read values 3, 7, and 10, process them, write the result, and then read values 4, 5, and 17. How do I direct readBin to those specific values? Thanks. You need to seek help. Or help(seek). Yes really. Barry -- Rick Reeves Sci

[R-sig-Geo] maptools/kmlPolygon example: two SPDFs, different results

data set and let me know where I might be mistaken! Thanks, Rick Reeves -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533 ___ R-sig-G

[R-sig-Geo] Preparing KML files for GoogleMaps/Earth: Why no kmlPoint() in maptools package?

n't built the rgdal library on their machines. Or, is there another way? Thanks, Rick R -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533 ___

Re: [R-sig-Geo] matching points to a raster

Elevation") Hope this helps, pls feel free to send more information! Rick Reeves Struve, Juliane wrote: Dear list, I am trying to incoporate spatial information (for example distance) into a movement model by reading an ascii -b converted raster and matching the raster with points of a

[R-sig-Geo] SP Spatial List operation - error in R 2.9 not in R 2.8 (download example)

ot; "y" .. .. ..$ : chr [1:2] "min" "max" ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots .. .. ..@ projargs: chr NA Warning message: In data.row.names(row.names, rowsi, i) : some row.names duplicated: 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37

Re: [R-sig-Geo] Cost distance analysis

v.net.path might work for you. There is a good data interface between R and GRASS: see the spgrass6 package. Rick Reeves Huizhu "Chrissy" Li wrote: Dear All, I am working on a project which involves shortest path, cost distance analysis. I could not find any package in R that c

Re: [R-sig-Geo] (no subject)

___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533

Re: [R-sig-Geo] Smallest distance from one point feature to a pointlayer

https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533 ___ R-sig-Geo mailing list R-sig-Geo@st

Re: [R-sig-Geo] Feature Extraction - object based image analysis

ryrawlins.googlepages.com<http://barryrawlins.googlepages.com/> -- Rick Reeves Scientific Programmer/Analyst and Data Manager National Center for Ecological Analysis and Synthesis UC Santa Barbara www.nceas.ucsb.edu 805 892 2533 ___ R-sig-Geo mailing li

Re: [R-sig-Geo] stratified random sampling

/SampleVectorPolysRastGrid.html You could modify this example to randomly select a point inside each polygon from the point grid (or a similar modification). Hope this helps.. Rick Reeves Heuvelink, Gerard wrote: Dear list, I use readShapePoly to import a shape file. The file has multiple polygons (actually, a

Re: [R-sig-Geo] world map data

Dave: Here is one place to obtain the Digital Chart of the World, which should be more accurate: http://www.maproom.psu.edu/dcw/ (Penn State University DCW Server) Get the shape files, use R maptools routines (e.g., ReadShapePoly) to read the files. Hope this helps, Rick Reeves Dave

Re: [R-sig-Geo] Subsetting a spatial grid

emory limitations. http://www.nceas.ucsb.edu/scicomp/GISSeminar/UseCases/SampleVectorPolygonsRasterGrid/SampleVectorPolysRastGrid.html I might suggest the use of GRASS for this task. Rick Reeves Dylan Beaudette wrote: On Wednesday 10 December 2008, Rob Robinson wrote: Help - please! :-) I have w

Re: [R-sig-Geo] identying shapefile polygons containing long-lat points

Hi: Check out the FindPolys method in the PBSMapping package. here is a code snippet from a larger example: (I will be extracting a simpler use case for distribution soon, but this will give you the general idea HucPolys = readShapePoly("AllHucsMerged") DamPoints = readShapePoints("Part1D

[R-sig-Geo] Solution in R for finding closest points:

Mary: Have a look at this tutorial: http://nceas.ucsb.edu/scicomp/GISSeminar/UseCases/AssignClosestPointsR/AssignClosestPointsR.html Uses R, which has a good set of geospatial analysis functions. This web site has several other examples! Hope this helps, Rick Reeves [EMAIL PROTECTED

Re: [R-sig-Geo] Shortest path around obstacles (OT)

the two systems) You might check out the GIS links on the attached page to get a bird's-eye view of GRASS: http://nceas.ucsb.edu/scicomp/SoftwareLinks.html You would have to be creative in your use of coordinate systems, but you could solve this problem using GRASS GIS. Regards, Rick R

Re: [R-sig-Geo] SpatialPolygonsDataFrame visualization

nd. http://nceas.ucsb.edu/scicomp/GISSeminar/UseCases/MapProdWithRGraphics/OneMapProdWithRGraphics.html Regards, Rick Reeves Mathieu Basille wrote: Hello everyone, I'm not very familiar with visualization of spatial objects in R, especially with the possibilities of the spplot functio

[R-sig-Geo] Method for: Creating contour maps of geocoded point shapefile attributes

classes. Thanks for any insight into this! Rick Reeves begin:vcard fn:Rick Reeves n:Reeves;Rick org:UC Santa Barbara;National Center For Ecological Analysis and Synthesis adr:Suite 300;;735 State Street;Santa Barbara;CA;93101;USA email;internet:[EMAIL PROTECTED] title:Rick Reeves, Scientific

Re: [R-sig-Geo] mapping introduction

map) Best Regards, Rick Reeves Tom Sgouros wrote: Hello all: I was referred to this list when I asked a question on the r-help list about mapping. Unfortunately, I seem to be a little late to the party, and a couple of weeks of monitoring the traffic have left me no more enlightened than I

[R-sig-Geo] Correction to that example...

Debs (and list) Please forgive the oversignt: that post contained two errors: pcrds <- cbind(tritest.vp[[8]][,1],tritest.vp[[8]][,2] ) replaced by pcrds <- cbind(tritest.vp[[i]][,1],tritest.vp[[i]][,2] ) and this line should be removed: points(crds) Thanks, Rick R here is a correct

Re: [R-sig-Geo] Can voronoi polygons calculated in R be converted to ESRI shapefile?

) library(maptools) writePolyShape(SPDF1, "hex") } Debarchaa Ghosh wrote: Hi Rick, Thanks for your help. I'll be looking forward to your script. Thanks, Debs. Rick Reeves wrote: Hello Debs: Yes, this can be done by 'promoting' the polygons generated by the trip

Re: [R-sig-Geo] Plotting point labels along with points using spplot sp.layout list

;> # display list entries for first two points >> # >> pt1 = list("sp.text",c(LatLongs[1,1],LatLongs[1,2]),ss[1],col="white",pos=2) >> pt2 = list("sp.text",c(LatLongs[2,1],LatLongs[2,2]),ss[2],col="white",pos=2) >> # >> #

Re: [R-sig-Geo] center of gravity

since it has the most straightforward input data requirements. Hope this helps, Rick Reeves -- Rick Reeves Scientific Programmer / Analyst National Center for Ecological Analysis and Synthesis UC Santa Barbara [EMAIL PROTECTED] www.nceas.ucsb.edu 805 892 2533 daiconnie wrote: Hi, everyo

Re: [R-sig-Geo] Units of centroid point coordinates generated byget.Pcent and calcCentroid?

proj4string is not NA, then you can use spTransform: d.ll <- spTransform(d, CRS("+proj=longlat")) Otherwise, you can project that matrix of coordinates with (if "m" is a matrix of your polygon centroids) x.ll <- project(m, "insert non-NA p4string here", inv = T) Do

[R-sig-Geo] Units of centroid point coordinates generated by get.Pcent and calcCentroid?

translate them to Lat/Long coordinates? Thanks, Rick Reeves -- Rick Reeves Scientific Programmer / Analyst National Center for Ecological Analysis and Synthesis UC Santa Barbara [EMAIL PROTECTED] www.nceas.ucsb.edu 805 892 2533 begin:vcard fn:Rick Reeves n:Reeves;Rick org:UC Santa Barba

[R-sig-Geo] Issue reading ESRI polygon shapefiles with PBSMapping 'importShapefile'

s shapefile using the PBSmapping tools. Has anyone else encountered and/or repaired this? Plan B would be to convert the Map object TriCounty into a PolySet or PolyData used by PBSmapping package, but so far I have not worked this out. Any suggestions? Thanks, Rick Reeves NCEAS / UC Santa Barbara

[R-sig-Geo] Processing HDF4 data files with R / GeospatialData Abstraction Library (GDAL)

d the issue of processing multiple-band HDF4 files with R / rgdal? Cheers, -- Rick Reeves Scientific Programmer / Analyst National Center for Ecological Analysis and Synthesis UC Santa Barbara [EMAIL PROTECTED] www.nceas.ucsb.edu 805 892 2533 begin:vcard fn:Rick Reeves n:Reeves;Rick org:UC Sa

[R-sig-Geo] simple yet durable error when trying to build GDAL 1.4.0 with HDF4 support

4.2r1-linux, but neither hdf4 or mfhdf lib found." Google search reveals that several others have encountered this issue in the past, but no solution is posted. Has anyone encountered and perhaps resolved this? Thanks in advance, Rick Reeves -- Rick Reeves Scientific Programmer / Analyst

Re: [R-sig-Geo] What is the latest re: reading HDF files with R?

the R HDF5 library routine hdf5load(), which seems to read the new .hdf5 format but NOT the older .hdf format (I tried, unsuccessfully) Question is: Has anyone developed a technique to read HDF files with R? Thanks in advance, Rick R --Rick Reeves Scientific Programmer / Analyst National C

[R-sig-Geo] What is the latest re: reading HDF files with R?

library routine hdf5load(), which seems to read the new .hdf5 format but NOT the older .hdf format (I tried, unsuccessfully) Question is: Has anyone developed a technique to read HDF files with R? Thanks in advance, Rick R -- Rick Reeves Scientific Programmer / Analyst National Center for

Re: [R-sig-Geo] How to get the code of function .spFindCG in package sp?

active R users and I am working to post a series of geospatial data processing tutorials for our residents. Stay tuned for more as soon as I can produce them Regards, Rick Reeves Scientific Programmer National Center for Ecological Analysis and Synthesis/ UC Santa Barbara Wang, Hao wrote: Hi

[R-sig-Geo] PBSMapping package: what are the default units for calcArea() results

t;, # SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] Please see attached jpg for the 'big picture' Many thanks for any suggestions! Regards, Rick Reeves UCSB / NCEAS -- Rick Reeves Scie

[R-sig-Geo] TLA

nyone explored (and returned alive from) the upper limits of R Matrix and SparseM problem sizes? Thanks in advance for any insights! Rick Reeves -- Rick Reeves Scientific Programmer Analyst National Center for Ecological Analysis and Synthesis (NCEAS) University of California, Santa Barbara [

Re: [R-sig-Geo] confused by coordinates

Mr. Zhang: There are many good summaries on the Internet: Google 'spatial coordinate systems' Here is a good example... http://www.environment.sa.gov.au/mapland/sicom/sicom/tp_scs.html Best of luck! Rick Reeves, UC Santa Barbara NCEAS www.nceas.ucsb.edu zhijie zhang wrote: >

[R-sig-Geo] A question re: tripack / voronoi cell area units

o the voronoi.mosaic() method. For example, if the coordinates are expressed in decimal degrees, then the cell side lengths are also expressed in decimal degrees and the cell areas would be in units of decimal degrees ** 2. Can someone affirm or correct this assumption? Thanks, Rick R -- Rick R

[R-sig-Geo] Writing SpatialPolygonsDataFrame object to an ESRI Shape File format? `

both a SpatialPolygon and SpatialPolygonDataFrame objects to use as input... Thanks, RR -- Rick Reeves Scientific ProgrammerAnalyst National Center for Ecological Analysis and Synthesis (NCEAS) University of California, Santa Barbara [EMAIL PROTECTED] 805 892

[R-sig-Geo] Current location of the Srings object (or its replacement)?

rings". Even loading the spmaptools package does not resolve this problem. Any suggestions on resolving this issue? Thanks, Rick Reeves Now, R reports that # --start --- # Srsl <- vector(mode="list", length=length(AllSps)) # # Srings object com

[R-sig-Geo] R function to match closest points to second point set?

tance. Thanks, Rick R -- Rick Reeves Scientific ProgrammerAnalyst National Center for Ecological Analysis and Synthesis (NCEAS) University of California, Santa Barbara [EMAIL PROTECTED] 805 892 2533 ___ R-sig-Geo mailing list R-sig-Geo@stat.math.et

[R-sig-Geo] Just ckecking: voronoi mosaic.cells[[i]]$area[] units?

? Thanks! Rick Reeves Scientific Programmer / Quantitative Analyst National Center for Ecological Analysis and Synthesis University of California, Santa Barbara 805 892 2534 [EMAIL PROTECTED] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https