[R-sig-Geo] Transform decimal degrees units in km - ReadOGR & owin{spatstat}

Aurelie C.Godin godina at dal.ca
Tue Apr 3 20:38:05 CEST 2012


Dear list,

I am trying to figure out how I can transform my unit of length from decimal
degree to square km when I import the following simple ArcGIS polygon shape. 
Ultimately, I would like to have this information available when I assign
this spatial polygon to my window in spatstat [i.e., Window area = xx square
unit
Unit of length: xx km].

div0A <- readOGR(dsn="Archive", layer="Projections")
> summary(div0A)
Object of class SpatialPolygonsDataFrame
Coordinates:
     min       max
x -80.00 -57.63430
y  66.25  78.16666
Is projected: FALSE 
proj4string :
[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
Data attributes:
       Id      NPAzimutha         UTM20            UTM19           AlberEA      
 Min.   :0   Min.   :348877   Min.   :349232   Min.   :349162   Min.  
:348656  
 1st Qu.:0   1st Qu.:348877   1st Qu.:349232   1st Qu.:349162   1st
Qu.:348656  
 Median :0   Median :348877   Median :349232   Median :349162   Median
:348656  
 Mean   :0   Mean   :348877   Mean   :349232   Mean   :349162   Mean  
:348656  
 3rd Qu.:0   3rd Qu.:348877   3rd Qu.:349232   3rd Qu.:349162   3rd
Qu.:348656  
 Max.   :0   Max.   :348877   Max.   :349232   Max.   :349162   Max.  
:348656 

library(spatstat)
newW2<-as(div0A,"owin")
> newW2
window: polygonal boundary
enclosing rectangle: [-80, -57.6343] x [66.25, 78.16666] units 

Thank you very much in advance,
Aurelie

-----
Aurelie Cosandey-Godin
Ph.D. Student, Department of Biology, Dalhousie University
Industrial Graduate Fellow, WWF-Canada
Email: godina at dal.ca | Web: wormlab.biology.dal.ca

--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Transform-decimal-degrees-units-in-km-ReadOGR-owin-spatstat-tp7434090p7434090.html
Sent from the R-sig-geo mailing list archive at Nabble.com.



More information about the R-sig-Geo mailing list