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

[Michael Barton]

Zbigniew,

I updated the r.in.aster script to read Level 1A ASTER imagery. I have 
a set and have tested it. It seems to work fine. I include a copy here. 
I am also posting it on my website in a new GRASS script exchange that 
I am going to launch (with some trepidation) in a followup posting. 
Thanks for your help in testing this.

Michael Barton

-------------- next part --------------
A non-text attachment was scrubbed...
Name: r.in.aster.tar.gz
Type: application/x-gzip
Size: 1985 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20040601/6cf6e5f7/r.in.aster.tar.gz
-------------- next part --------------


______________________________
Michael Barton, Professor & Curator
School of Human Origins, Cultures, & Societies
Arizona State University
Tempe, AZ  85287-2402
USA

voice: 480-965-6262; fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton
On May 28, 2004, at 2:23 PM, Zbigniew Perski wrote:

> Dear Paul,
>
> As I did, you are trying to georeference subdataset 1. This is NOT  
> image band 1. This is the metadata and GCP's for image bands 1-3 VNIR.
>
> So, now I understand my frustrations!
>
>  To actually georeference image band 1, you need to indicate 
> subdataset  2. The relationship between images and subdatasets is 
> listed at the  bottom of the output to gdal-info. In brief, it is
>
>>
>> band 1=subdataset 2, 2=3, 3n=4, 3b=5, 4=8, 5=9, 6=10, 7=11, 8=12, 
>> 9=13,  10=16, 11=17, 12=18, 13=19, 14=20. Subdatasets 1, 5, 6, 7, 14, 
>> and 15  are reserved for metadata and GCP's.
>
> OK.
>
>>
>>
>> If Level 1A can be georeferenced, I can easily make the change to the 
>>  script to let you select them via a flag.
>
> This I think will be helpfull for many people since L1A data are more 
> frequent in ASTER database (as I remember 2 years ago L1B was realy 
> rare.).
>
>> Also, I suspect that MODIS  images work much the same way as they 
>> come from the same organizations  and use the same hdf format. If so 
>> and someone could send me the output  of gdal-info for a MODIS image,
>
> For MODIS not, but for ASTER 1A I can send you some outputs in monday 
> if you wisch. I Have not MODIS data in this moment.
>
>> I could try to either include MODIS in  an expanded r.in.aster or 
>> make a companion r.in.modis.
>>
> It will be great too I think!
>
>
> Many thanks
>
> Zbigniew
>
>
>