proj: lambert to lat-lon.

Agustin Lobo alobo at ija.csic.es
Sun Jan 14 07:00:00 EST 1996


***************************************************
Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Marti Franques s/n
08028 Barcelona (Spain)

tel. 34 3 330 2716
fax. 34 3 411 0012

alobo at ija.csic.es
alobo at eno.princeton.edu
***************************************************
Delivery-Date: 12 Jan 96 19:59 +0200
From: Agustin Lobo <alobo at ija.csic.es>
Sender: grass-lists-owner at moon.cecer.army.mil.
Reply-To: grassu-list at moon.cecer.army.mil.
Precedence: Bulk
To: gie at charon.er.usgs.gov
Message-ID: <501*alobo at ija.csic.es>
Subject: problem with proj

Thanks for making proj available. I have a problem
and wonder if you could help me.

I need to convert lmbert cartesian coordinates to lat-lon geog.
coordiantes in the Iberian Peninsula. The annotation for my 
(digitized) map is:

Datum: European 195o Postdam.
Ellipsoid: International
Origin:
long: -3d41'15''
lat:40d
secant ("automecoicos" in the original) paralels:
42d49'18.9420'', 37d10'41.0580''
False East: 600 Km
False North: 600 Km

I issue the proj command:
proj -I -v +ellps=intl +proj=lcc +x_0=600000 +y_0=600000 +lat_1="37d10'41.0580"
 +lat_2="42d49'18.9420" 
#Lambert Conformal Conic
#       Conic, Sph&Ell
#       lat_1= and lat_2= or lat_0
# +ellps=intl +proj=lcc +x_0=600000 +y_0=600000 +lat_1=37d10'41.0580
# +lat_2=42d49'18.9420
#--- following specified but NOT used
# +ellps=clrk66 +lat_1=33 +lat_2=45

and get absurd values as:
1110000 860000
3d44'51.55"E    1d49'28.534"N

(1110000 860000 should be like about 6d10' 42d11)

What am I doing wrong?

Thanks

Agus

***************************************************
Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Marti Franques s/n
08028 Barcelona (Spain)

tel. 34 3 330 2716
fax. 34 3 411 0012

alobo at ija.csic.es
alobo at eno.princeton.edu
***************************************************





More information about the grass-user mailing list