Loosing information when projecting and writing a shapefile with rgdal
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