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