#!/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_7x20020524.ETM-EarthSat-Orthorectified
# 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
r.in.gdal p016r035_7t20000331_z17_nc_spm_wake.tif out=lsat7_${YEAR} --o
g.rename lsat7_2000.1,lsat7_2000.10 --o
g.rename lsat7_2000.2,lsat7_2000.20 --o
g.rename lsat7_2000.3,lsat7_2000.30 --o
g.rename lsat7_2000.4,lsat7_2000.40 --o
g.rename lsat7_2000.5,lsat7_2000.50 --o
g.rename lsat7_2000.6,lsat7_2000.70 --o

#temp
r.in.gdal p016r035_7k20000331_z17_nc_spm_wake.tif out=lsat7_${YEAR} --o
g.rename lsat7_2000,lsat7_2000.61 --o

#pan
r.in.gdal p016r035_7p20000331_z17_nc_spm_wake.tif out=lsat7_${YEAR} --o
g.rename lsat7_2000,lsat7_2000.80 --o

#viz/IR
for k in 10 20 30 40 50 70 ; do
 echo "Processing $k..."
 r.support map=lsat7_${YEAR}.$k history="$METASTRING"
 echo "Imported as 'lsat7_${YEAR}.$k'"
done

#temp
for k in 61 ; do
 echo "Processing $k..."
 r.support map=lsat7_${YEAR}.$k history="$METASTRING"
 echo "Imported as 'lsat7_${YEAR}.$k'"
done

#pan
for k in 80 ; do
 echo "Processing $k..."
 r.support map=lsat7_${YEAR}.$k history="$METASTRING"
 echo "Imported as 'lsat7_${YEAR}.$k'"
done