[GRASS-user] Landsat image tiffs

[Hamish]

> > Dylan Beaudette wrote:
> Zenon Panoussis wrote:

> > I have had some similar problems wrt. exporting landsat data from
> > GRASS ->  mapsever. here is my approach:
> 
> > #export tiled landsat-shade tiles:
> > for x in `seq 2 1 49`; do g.region rast=ca_landsat_shade_$x;
> > r.out.tiff -t  in=ca_landsat_shade_$x out=ca_landsat_shade_$x.tif;
> > done
> 
> > # re-create geotiffs with gdal_translate (fix bug in tiffs created
> > # with 
> > r.out.tiff)
> > for x in `ls *.tif`; do gdal_translate -a_srs '+proj=aea +x_0=0.0
> > +y_0=0.0  +lon_0=-96.0 +lat_0=40.0 +lat_1=20.0 +lat_2=60.0
> > +datum=NAD83' $x ../new/$x ;  done 
> 
> Hmm. You had to set the projection manually?
..
> No projection. Then perhaps my suspicion that I'm using the wrong
> data type is wrong, it could be the missing projection that's causing
> mapserver to ignore the image.
> 
> > terse i know, but hopefully helpful.
> 
> It is already :) I only need to figure the values of my -a_srs so
> I can test it. And mapserver output your ../new/$x without problems?


hint to get the current mapset's PROJ.4 parameters:

gdal_translate -a_srs "`g.proj -jf`"



Hamish