[R-sig-Geo] Conversion geographical coodinate in UTM

[Barry Rowlingson]

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 <gbl1 at hi.is> wrote:
> 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:12 GMT, Barry Rowlingson <b.rowlingson at lancaster.ac.uk> a écrit:
>
>> On Mon, Jun 30, 2014 at 3:18 PM, Sarah Goslee <sarah.goslee at gmail.com> 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] "SpatialPolygonsDataFrame"
>>
>> Use rgdal:
>>
>> require(rgdal)
>> Front = readOGR(".","ISL_adm0")
>>
>> then check Front has a proj4string. summary(Front) should mention it.
>>
>> This is getting to be the R-sig-geo equivalent of R FAQ 7.31..
>>
>> Barry
>
>
>