Dear Andres, You'll need spTransform() to reproject the layer into WGS84
carto2012x <- unionSpatialPolygons(carto2012, IDs=carto2012$DPA_CANTON) spTransform(carto2012x, CRS("+proj=longlat")) Best regards, ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance Kliniekstraat 25 1070 Anderlecht Belgium To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey 2017-05-10 10:17 GMT+02:00 Andrés Peralta <tiric...@gmail.com>: > Hi to everyone, > > I`m working with the second administrative level cartography from Ecuador > (available in: > http://www.ecuadorencifras.gob.ec//documentos/web-inec/ > Cartografia/2015/Clasificador_Geografico/2012/SHP.zip). > The shape file is called nxcantones.shp. I`ve been having a lot of problems > opening and working with this cartography in R; but i can open it in Q-GIS > and ARCGIS without any problem (I have opened it in both programs and saved > it again). As the cartography has various objects with the same ID - > aparently the same areas but repeated; we had to run the following sintax > in order to have one ID in each area: > > *carto2012 <- readOGR("CANTONES2012/2012CLEAN.shp", "2012CLEAN", > stringsAsFactors=F) * > *proj4string(carto2012) <- CRS("+proj=utm +zone=17 +south +datum=WGS84 > +units=m +no_defs")* > > *carto2012x <- unionSpatialPolygons(carto2012, IDs=carto2012$DPA_CANTON)* > *proj4string(carto2012x) <- CRS("+proj=utm +zone=17 +south +datum=WGS84 > +units=m +no_defs")* > > > *datos <- data.frame(ID=row.names(carto2012x), stringsAsFactors = F)* > *row.names(datos) <- row.names(carto2012x)* > *carto2012x2 <- SpatialPolygonsDataFrame(carto2012x, data=datos)* > > *carto2012 <- carto2012x2 * > > Either using the original (nxcantones) or the modified cartography > (carto2012), I can plot the map in R, but I can´t plot it in R Leaflet. It > just opens the base map. I seems as the projection is lost in the way or > that the cartography has a problem. > > Any ideas? > > -- > > * Andrés Peralta* > Pre-doctoral Researcher > GREDS/EMCONET / ASPB > <https://www.upf.edu/greds-emconet/en/> <http://aspb.cat/> > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo