UTM to LatLon conversion in libraries?

UTM to LatLon conversion in libraries?

Post by iavar.. » Sun, 07 Jun 1998 04:00:00



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

 
 
 

UTM to LatLon conversion in libraries?

Post by William Capeha » Thu, 18 Jun 1998 04:00:00




>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://www.veryComputer.com/.*ia.edu/~eww6n/math/MapProjection.html

that can help with the UTM conversions

================================================================

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://www.veryComputer.com/~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)
;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;