[R-sig-Geo] spTransform: datum syntax?
Agustin Lobo
Agustin.Lobo at ija.csic.es
Thu Dec 6 14:46:54 CET 2007
Dear list,
I have to transform the datum of an spatial object.
First, I create the object with:
pressp <- SpatialPointsDataFrame(pres[,c(10,6)], data=pres[,c(1,2,11)],
+ proj4string = CRS("+proj=longlat +ellps=WGS84"))
Then I have to transform to UTM but using the International or Hayford
datum.
The following works:
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31))
but then I keep WGS84.
The following also works:
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31 +datum=NAD27"))
but this is not the datum I want
Then I've tried:
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31 +datum=ED50"))
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31
+datum=International"))
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31 +datum=Int"))
presspUTM <- spTransform(pressp, CRS("+proj=utm +zone=31 +datum=Hayford"))
and many other similar names, also tried with ellps=
The doc refers to the proj.4 doc, but cannot find any list of datums
there, including the pdf docs of proj.4
Does anyone know where I could find a list of supported datums in proj.4
and their correct names for the +datum= argument?
Thanks!
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: Agustin.Lobo at ija.csic.es
http://www.ija.csic.es/gt/obster
More information about the R-sig-Geo
mailing list