[R-sig-Geo] wkbMultiPoint and rgdal

Jonathan Greenberg greenberg at ucdavis.edu
Tue Jul 26 23:13:57 CEST 2011


R-sig-geo'ers:

I have an ESRI shapefile which has multi-point geometry, and I'm
trying to get it into R for use with sp() and raster().  Using:

points=readOGR(dsn=getwd(),layer="tnf_sub_2_TreeDig_pts_110725_v02")
returns:

Error in readOGR(dsn = getwd(), layer = "tnf_sub_2_TreeDig_pts_110725_v02",  :
  Incompatible geometry: 4
In addition: There were 50 or more warnings (use warnings() to see the first 50)

the warnings all say:
...
50: In readOGR(dsn = getwd(), layer = "tnf_sub_2_TreeDig_pts_110725_v02",  ... :
  eType not chosen

What I really want is to, in the end, have a SpatialPointsDataFrame
where each multi-point feature is converted to individual points that
have the attribute of what multi-point feature it came from.  What is
the best way to do this?  I'm attaching the zipped shapefile to this
email (hopefully it won't get bounced).

--j

-- 
Jonathan A. Greenberg, PhD
Assistant Project Scientist
Center for Spatial Technologies and Remote Sensing (CSTARS)
Department of Land, Air and Water Resources
University of California, Davis
One Shields Avenue
Davis, CA 95616
Phone: 415-763-5476
AIM: jgrn307, MSN: jgrn307 at hotmail.com, Gchat: jgrn307
-------------- next part --------------
A non-text attachment was scrubbed...
Name: tnf_sub_2_TreeDig_pts_110725_v02.zip
Type: application/zip
Size: 24005 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110726/8444e48e/attachment.zip>


More information about the R-sig-Geo mailing list