[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