#!/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/p016r035_7x20000331.ETM-EarthSat/
# requires GRASS 6.3 for recent r.support

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

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

METASTRING="`grep \"SPACECRAFT_ID\|SENSOR_ID\|ACQUISITION_DATE\|WRS_PATH\|WRS_ROW\|SUN_AZIMUTH\|SUN_ELEVATION\|QA_PERCENT_MISSING_DATA\|CLOUD_COVER\|CPF_FILE_NAME\|LMAX_BAND\|LMIN_BAND\|QCALMAX_BAND\|QCALMIN_BAND\" *.met| grep -v '^$' | tr -s '\t' ' ' | sed 's+ ++g' | sed 's+"++g' | tr '\n' ',' | sed 's+,$++g'`"
YEAR=`grep ACQUISITION_DATE *.met| cut -d'=' -f2 | sed 's+ ++g' | cut -d'-' -f1`

#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=lsat7_${YEAR}.$i --o
 r.support map=lsat7_${YEAR}.$i history="$METASTRING"
 echo "Imported as 'lsat7_${YEAR}.$i'"
done

#temp
FILES="${PREFIX1}k${PREFIX2}nn61_${SUFFIX}.tif ${PREFIX1}k${PREFIX2}nn62_${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=lsat7_${YEAR}.$i --o
 r.support map=lsat7_${YEAR}.$i history="$METASTRING"
 echo "Imported as 'lsat7_${YEAR}.$i'"
done

#pan
FILES="${PREFIX1}p${PREFIX2}nn80_${SUFFIX}.tif"

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