Dear R FAQ Team :p
So I have been able to use the readOGR function, could anyone transfer to its
maintainer a request to add writeOGR in the "see also" section of help?
Anyway, after having done that for my shapefile, I tried to deal with the
points I had to plot on, here is what I tried:
On Mon, 30 Jun 2014, Barry Rowlingson wrote:
Not your fault that readShapeSpatial ignores the projection!
I think Roger recently spoke about fixing these things...
The underlying shapelib code will not be maintained, anmd
maptools::readShape* will be deprecated very soon. Always use
rgdal::
Not your fault that readShapeSpatial ignores the projection!
I think Roger recently spoke about fixing these things...
Barry
On Mon, Jun 30, 2014 at 4:33 PM, Gilles Benjamin Leduc wrote:
> Thanks Barry,
>
> That's it… Sorry if my question was too simple… I'm not a geograph, I'm a
> molecular
Thanks Barry,
That's it… Sorry if my question was too simple… I'm not a geograph, I'm a
molecular biologist… so I'm not used to these packages … I started from the
simple script I had for normal maps… (and it is monday, brain isn't fully
operational!)
Benjamin
Le Lundi 30 Juin 2014 15
On Mon, Jun 30, 2014 at 3:18 PM, Sarah Goslee wrote:
> Hi,
>
> At a guess, your Front dataset doesn't have a CRS set, so
Probably because you've used readShapeSpatial which ignores projection
info from shapefiles...
>> Front<-readShapeSpatial("ISL_adm0.shp")
>> class(Front)
>> [1] "SpatialPolygo
Hi,
At a guess, your Front dataset doesn't have a CRS set, so
spTransform() doesn't know what to transform *from*.
See the examples in ?spTransform for the whole process, including
setting the CRS for your data.
Sarah
On Mon, Jun 30, 2014 at 10:10 AM, Gilles Benjamin Leduc wrote:
> Hi,
>
> I n
Hi,
I need to convert some geographical coordinate it UTM.
First I have a shapefile:
Front<-readShapeSpatial("ISL_adm0.shp")
class(Front)
[1] "SpatialPolygonsDataFrame"
then… I tried this after googling:
DataUTM <- spTransform(Front, CRS("+proj=utm +zone=27 +north +datum=WGS84"))
Erreur dans