#!/bin/sh

#import preprocessed GLSC Maryland LANDSAT data 
# along with relevant metadata
# Markus Neteler 3,8/2007

# ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/WRS2/p016/r035/p016r35_5t871014.TM-EarthSat-Orthorectified
# p016r035_5t19871014_z17_nn10_nc_spm_wake.tif

# requires GRASS 6.3 for recent r.support

MYPATH=p016
ROW=r035
MYDATE=19871014
UTM=z17
SUFFIX=nc_spm_wake

############## nothing to change below
PREFIX1=${MYPATH}${ROW}_5
PREFIX2=${MYDATE}_${UTM}_

dos2unix *.met
METASTRING="`grep "IMAGE_ID\|PATH\|ROW \|DATE\|PLATFORM" *.met | tr ' ' '=' | tr '\n' ','| sed 's+,$++g'`"
YEAR="19`grep DATE *.met| cut -d'/' -f3`"

#viz/IR
FILES="${PREFIX1}t${PREFIX2}nn10_${SUFFIX}.tif  ${PREFIX1}t${PREFIX2}nn20_${SUFFIX}.tif
 ${PREFIX1}t${PREFIX2}nn30_${SUFFIX}.tif ${PREFIX1}t${PREFIX2}nn40_${SUFFIX}.tif
 ${PREFIX1}t${PREFIX2}nn50_${SUFFIX}.tif ${PREFIX1}t${PREFIX2}nn70_${SUFFIX}.tif"

for k in $FILES ; do
 i=`echo $k | sed 's+'${PREFIX1}'t'${PREFIX2}nn'++g' | sed 's+_'${SUFFIX}'.tif++g'`
 echo "Processing $k..."
 r.in.gdal $k out=lsat5_${YEAR}.$i --o
 r.support map=lsat5_${YEAR}.$i history="$METASTRING"
 echo "Imported as 'lsat5_${YEAR}.$i'"
done

#temp
FILES="${PREFIX1}k${PREFIX2}nn60_${SUFFIX}.tif"

for k in $FILES ; do
 i=`echo $k | sed 's+'${PREFIX1}'k'${PREFIX2}nn'++g' | sed 's+_'${SUFFIX}'.tif++g'`
 echo "Processing $k..."
 r.in.gdal $k out=lsat5_${YEAR}.$i --o
 r.support map=lsat5_${YEAR}.$i history="$METASTRING"
 echo "Imported as 'lsat5_${YEAR}.$i'"
done