Problems With Cartography - Help

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|

Problems With Cartography - Help

Andrés Peralta
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
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: Problems With Cartography - Help

Thierry Onkelinx
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 <[hidden email]>:

> 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
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: Problems With Cartography - Help

Andrés Peralta
Thank you Thierry, it worked perfectly.

Best regards,

Andrés

On Wed, May 10, 2017 at 10:25 AM, Thierry Onkelinx <[hidden email]
> wrote:

> 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 <[hidden email]>:
>
>> Hi to everyone,
>>
>> I`m working with the second administrative level cartography from Ecuador
>> (available in:
>> http://www.ecuadorencifras.gob.ec//documentos/web-inec/Carto
>> grafia/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
>> [hidden email]
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>


--

*             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
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo