[R-sig-Geo] projection

[Ben Madin]

Gerard,

Until someone else replies, XXXXXXXX might be something like :

+proj=lcc +lat_1=37.8 +lat_0=37.8 +lon_0=34.9134 +k_0=0.99987742 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs

but it might also not be... I'm particularly wary of the x_0 and y_0 values. 

You could also try the proj list?

cheers

Ben



On 18/05/2010, at 20:08 , Heuvelink, Gerard wrote:

> Dear list,
> 
> I hope someone can help me with the following.
> 
> I have an SRTM DEM of Turkey in the "ED 1950 Lambert Conformal Conic"
> projection system.
> 
> However, I do not know how this is referred to in R. For example, I
> might load the dem as:
> 
> DEM = readGDAL("srtm_5km_ascii.txt")
> 
> Next I will want to let know what projection it is by issuing:
> 
> proj4string(DEM) = XXXXXXXXXXX
> 
> What do I write for XXXXXXXX?
> 
> Thanks, Gerard
> 
> 
> Gerard B.M. Heuvelink
> Environmental Sciences Group
> Wageningen University and Research Centre
> P.O. Box 47
> 6700 AA Wageningen
> The Netherlands
> 
> tel +31 317 486538 Mon Tue
> tel +31 317 482716 Wed Thu Fri
> 
> email gerard.heuvelink at wur.nl <mailto:gerard.heuvelink at wur.nl> 
> http://www.lad.wur.nl/UK/ <http://www.lad.wur.nl/UK/> 
> http://www.alterra.wur.nl/NL/onderzoek/Werkveld+Bodem/TBGE/Gerard+Heuvel
> ink/
> <http://www.alterra.wur.nl/NL/onderzoek/Werkveld+Bodem/TBGE/Gerard+Heuve
> link/> 
> 
> 
> 
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo