#!/bin/sh

#for GLSC Maryland LANDSAT data in gzip (tif.gz)
#MN 10/2003
#MN 2/2004
#MN $Date: 2004-04-26 18:24:01 +0200 (Mon, 26 Apr 2004) $

#Make WGS84/WGS84 -> Gauss-Boaga ${SUFFIX}, cut out region of interest

MYPATH=p016
ROW=r035
MYDATE=20020524
UTM=z17

PREFIX1=${MYPATH}${ROW}_7
PREFIX2=${MYDATE}_${UTM}_

g.region n=228500 s=215000 w=630000 e=645000 res=14.25 -pa
eval `g.region -g`
SUFFIX=nc_spm_wake

############## nothing to change below
#viz/IR
FILES="${PREFIX1}t${PREFIX2}nn10.tif.gz  ${PREFIX1}t${PREFIX2}nn20.tif.gz
 ${PREFIX1}t${PREFIX2}nn30.tif.gz ${PREFIX1}t${PREFIX2}nn40.tif.gz
 ${PREFIX1}t${PREFIX2}nn50.tif.gz ${PREFIX1}t${PREFIX2}nn70.tif.gz"

for k in $FILES ; do
 i=`echo $k | sed 's+'${PREFIX1}'t'${PREFIX2}'++g' | sed 's+.tif.gz++g'`
 echo "Processing $k..."
 gunzip $k
 gdalwarp -order 2 -t_srs "`g.proj -wf`" -te $w $s $e $n -tr 28.5 28.5 ${PREFIX1}t${PREFIX2}$i.tif \
   ${PREFIX1}t${PREFIX2}$i\_${SUFFIX}.tif
 gzip ${PREFIX1}t${PREFIX2}$i.tif
done

#temp
FILES="${PREFIX1}k${PREFIX2}nn61.tif.gz ${PREFIX1}k${PREFIX2}nn62.tif.gz"

for k in $FILES ; do
 i=`echo $k | sed 's+'${PREFIX1}'k'${PREFIX2}'++g' | sed 's+.tif.gz++g'`
 echo "Processing $k..."
 gunzip $k
 gdalwarp -order 2 -t_srs "`g.proj -wf`" -te $w $s $e $n -tr 57.0 57.0 ${PREFIX1}k${PREFIX2}$i.tif \
   ${PREFIX1}k${PREFIX2}$i\_${SUFFIX}.tif
 gzip ${PREFIX1}k${PREFIX2}$i.tif
done

#pan
FILES="${PREFIX1}p${PREFIX2}nn80.tif.gz"

for k in $FILES ; do
 i=`echo $k | sed 's+'${PREFIX1}'p'${PREFIX2}'++g' | sed 's+.tif.gz++g'`
 echo "Processing $k..."
 gunzip $k
 gdalwarp -order 2 -t_srs "`g.proj -wf`" -te $w $s $e $n -tr 14.25 14.25 ${PREFIX1}p${PREFIX2}$i.tif \
   ${PREFIX1}p${PREFIX2}$i\_${SUFFIX}.tif
 gzip ${PREFIX1}p${PREFIX2}$i.tif
done