#Dieses Script schneidet Subsets (Region Draatal) aus NOAA 1x  Channel (!!!) Images (Terascan TDF Format) und
#exportiert dieses dann inkl. der Geoinformationen ins ASCII Format.
#Anschliessend kann man diesen File z.B. mit dem IDL Programm "TDF2ENVI_byte_channel.pro" in ENVI importieren und
#weiterverarbeiten. "line" und "sample" sowie Ein- ud Ausgabe im Envi Sript müssen angepasst werden.
#Infos zu "line" und "sample" bekommt man mit dem Befehl "contents" und dem entsprechenden File. (=>"contents filexyz")
#RSRG 29.05.2002 

#! /bin/csh -fx

cd /usr/data/maroc/masked_files/

foreach file (noaa-11*)

subset \
        include_vars='avhrr_ch1 avhrr_ch2 avhrr_ch3 avhrr_ch4 avhrr_ch5' \
        clip_dims='line sample' \
        range_pairs='550 950 350 860' \
        instantiate='yes' \
$file \
$file.sub

angles \
	latitude='yes' \
	longitude='yes' \
	sat_zenith='n' \
	sun_zenith='n' \
	scatter_phase='n' \
	sun_reflection='n' \
	rel_azimuth='n' \
	store_cosines='n' \
	real_output='y' \
	poly_size='100' \
$file.sub

expasc \
	include_vars='latitude longitude avhrr_ch1 avhrr_ch2 avhrr_ch3 avhrr_ch4 avhrr_ch5'\
	all_same_dim='yes'\
	list_dims='no'\
$file.sub \
$file.ascii

echo removing null data points

sed '/\*\*\*/s//0/g' $file.ascii >$file.idl


rm $file.ascii
mv $file.idl /usr/data/maroc/idl_export/envi_maroc_sub_channel/
end