[Proj] polar stereographic coord

[Eric Miller]

You can't project an ascii grid with proj.  It only works with coordinates.  Additionally, the "proj" command only does forward or inverse projections.  You'd need to use "cs2cs" to convert between two coordinate systems.  However, I think gdalwarp might be able to handle projecting the data.

>>> On 8/27/2008 at 9:20 AM, Adriana Huerta-Casas <ahuerta at awi-bremerhaven.de>
wrote:

> 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
>