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.heuvel...@wur.nl <mailto:gerard.heuvel...@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@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo