[R-sig-Geo] Loosing information when projecting and writing a shapefile with rgdal

Barry Rowlingson b.rowlingson at lancaster.ac.uk
Sat Oct 23 00:32:34 CEST 2010


On Fri, Oct 22, 2010 at 8:01 PM, Manuel Spínola <mspinola10 at gmail.com> wrote:

> cob.transf = spTransform(cob, CRS("+proj=tmerc +lon_0=-84 +lat_0=0
> +x_0=500000 +k=0.9999 +datum=WGS84")) # to CRTM05

 Is it just a typo? You called it 'cob.transf' above, but below you've
written 'corr.transf':

> writeOGR(corr.transf,
> dsn="C:/ProyectoRespacial/CBSS/Shapefiles_proyectados",
> layer="coberturaCRTM05", driver="ESRI Shapefile")

?

--
blog: http://geospaced.blogspot.com/
web: http://www.maths.lancs.ac.uk/~rowlings
web: http://www.rowlingson.com/
twitter: http://twitter.com/geospacedman
pics: http://www.flickr.com/photos/spacedman



More information about the R-sig-Geo mailing list