[Proj] polar stereographic coord

Adriana Huerta-Casas ahuerta at awi-bremerhaven.de
Wed Aug 27 12:20:37 EDT 2008


Hi,
                                                                                
I'm new using proj4 and I can't really figure out how to do this properly:
                                                                                
I need to convert data given in polar stereographical coordinates (from
BEDMAP, 71S lat true scale,0E central meridian,5km nominal grid spacing)
to cartesian coordinates. The input file is ascii and contain headers.
I'm after a file x,y,data.
                                                                                
The first lines of the file to read and convert are:
ncols         1200
nrows         1000
xllcorner     -3000000
yllcorner     -2500000
cellsize      5000
NODATA_value  -9999
-9999 -9999 -9999 -999 ...
                                                                                
                                                                                
I'm doing the follwing:
                                                                                
 cat icethic.asc | proj +proj=stere +lon_0=0 +lat_0=-90 +lat_ts=-71
+ellps=WGS84 +datum=WGS84
                                                                                
but I'm wondering how to write a file and read the headers.
                                                                                
Thanks a lot!
adriana


-- 
----------------------------
Dr. Adriana Huerta-Casas
Alfred-Wegener-Institute for Polar 
 and Marine Research
Bussestrasse 24 (Building F-406)
D-27570 Bremerhaven, Germany
Telephone:+49(471)4831-1855





More information about the Proj mailing list