UTM to LatLon conversion in libraries? [message #11955] |
Sat, 06 June 1998 00:00  |
Ferdinando Iavarone
Messages: 10 Registered: October 1997
|
Junior Member |
|
|
I'm looking for a UTM to Lat-Lon conversion function.
I need it to convert a Digital Elevation Model given in UTM coordinates
to Geographic coordinates, so that it can be related with satellite
positions.
Anybody heard of an IDL function implemented somewhere that can help?
best regards,
/NaN do
--
Nando Iavarone
Gruppo Telerilevamento
Dipartimento Interateneo di Fisica
Via Amendola, 173
70126 BARI
ITALY
Tel.: +39 80 5443166
Fax.: +39 80 5443165
e_mail: iavarone@axpba0.ba.infn.it
|
|
|
Re: UTM to LatLon conversion in libraries? [message #12016 is a reply to message #11955] |
Wed, 17 June 1998 00:00  |
wcapehar
Messages: 12 Registered: May 1998
|
Junior Member |
|
|
In article <3578F459.C11E6BCF@axpba0.ba.infn.it>,
<iavarone@axpba0.ba.infn.it> wrote:
> I'm looking for a UTM to Lat-Lon conversion function.
> I need it to convert a Digital Elevation Model given in UTM coordinates
> to Geographic coordinates, so that it can be related with satellite
> positions.
> Anybody heard of an IDL function implemented somewhere that can help?
I have a PVWAVE function that does goes from Lambert Conformal Conic to LL and it
works in IDL. It was easy enough to write since UTM is more straight forward than
LCC you should have a better time. It's hamfisted but it gets the job done.
There are good sites our there esp at
ftp://mapping.usgs.gov/pub/software/current_software/gctpc/ and
http://astsun.astro.virginia.edu/~eww6n/math/MapProjection.h tml
that can help with the UTM conversions
============================================================ ====
Bill Capehart <wjc@ias.sdsmt.edu> Land Surface Processes
Institute of Atmospheric Sciences
South Dakota School of Mines and Technology Ph: +1-605-394-1994
Rapid City, SD 57701-3995 Fax: +1-605-394-6061
=============== http://sky.ias.sdsmt.edu/~wjc ==================
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;
;
; Georeferencing data
;
; World Meteorological Organization Grid #212 (Lam Conf Con)
;
pi4 = !dpi / 4.0;
radius = 6371000.0;
;
dy = 40.63525*1000.;
dx = 40.63525*1000.;
;
centlat = 35.0000 * !dpi / 180.00;
centlon = -95.0000 * !dpi / 180.00;
centralx = 105-1
centraly = 49-1
;
reflat1 = 25.0000 * !dpi / 180.00;
reflat2 = 25.0000 * !dpi / 180.00;
;
ycorner = -(centraly)*dy;
xcorner = -(centralx)*dx;
;
eastings = xcorner + dindgen(nx)*dx
northings = ycorner + dindgen(ny)*dy
;
latitude = dblarr(nx,ny)
longitude = dblarr(nx,ny)
;
; Projection Specific Area Begins
;
; [UTM Conversion info goes here]
;
Projection Specific Area Ends
;
centlat = centlat / !dpi * 180.00;
centlon = centlon / !dpi * 180.00;
;
reflat1 = reflat1 / !dpi * 180.00;
reflat2 = reflat2 / !dpi * 180.00;
;
print, 'Grid 212 Addresses' print, ' LLC:
',latitude(0,0), longitude(0,0), $
northings(0), eastings(0)
print, ' ULC: ',latitude(0,ny-1), longitude(0,ny-1), $
northings(ny-1), eastings(0)
print, ' URC: ',latitude(nx-1,ny-1),longitude(nx-1,ny-1), $
northings(ny-1),eastings(nx-1)
print, ' LRC: ',latitude(nx-1,0), longitude(nx-1,0), $
northings(0), eastings(nx-1)
print, ' '
print, ' CPt: ',latitude(centralx,centraly), $
longitude(centralx,centraly), $
northings(centraly), eastings(centralx)
;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;
|
|
|