[GRASSLIST:3533] Re: script for importing Terra-ASTER data

[Zbigniew Perski]

I found it. The problem was caused by Level 1A data.

r.in.aster is tryng to read the data as 1B level:
HDF4_SDS:ASTER_L1B:"/home/pers/Projekty/Svalbard/Ny_alesund/ASTER/pg-PR1A0000-2000122501_029_009":2
when I tried to use gdalwarp stand alone:

gdalwarp -t_srs '+proj=latlong +a=6378137 +rf=298.257223563 +no_defs'
HDF4_SDS:ASTER_L1A:"pg-PR1A0000-2000122501_029_009":1
pg-PR1A0000-2000122501_029_009_2.tif
ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates for
HDF4_SDS:ASTER_L1A:pg-PR1A0000-2000122501_029_009:1.
There is no affine transformation and no GCPs.


In fact the georeferencing should be possible because the data conaitns
131 GCP-s:
Coordinate System is `'
GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.2572235629972,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree (supplier to define representation)",0.01745329251994328],AUTHORITY["EPSG","4326"]]
GCP[  0]: Id=1, Info=
          (0.5,-119.5) -> (13.7399,79.2685,0)
GCP[  1]: Id=2, Info=
          (410.5,-119.5) -> (13.9667,79.2295,0)
......


Regards

Zbigniew




W liście z pią, 28-05-2004, godz. 07:00, Michael Barton pisze: 
> Paul,
> 
> Thanks for pointing this out.
> 
> Michael
> ____________________
> C. Michael Barton, Professor
> School of Human Origins, Cultures, & Societies
> PO Box 872402
> Arizona State University
> Tempe, AZ  85287-2402
> USA
> 
> Phone: 480-965-6262
> Fax: 480-965-7671
> www: <www.public.asu.edu/~cmbarton>
> On May 27, 2004, at 9:16 PM, Paul Kelly wrote:
> 
> > On Thu, 27 May 2004, Michael Barton wrote:
> >
> >> The script should work in either 5.3 or 5.7, but the GUI interface 
> >> only comes up in 5.7.
> >
> > It uses g.proj, which is only available in 5.7.
> >
> > Paul
> >
>