#Dieses Script schneidet Subsets (hier: Region Draatal) aus NOAA NDVI (!!!) Images (Terascan TDF Format) und #exportiert dieses dann inkl. der Geoinformationen ("AMGLES") ins ASCII Format ("EXPASC"). Mit dem "sed"-Befehl werden #alle "No Data Pixel" im TDF File gesucht und mit "-9" deklariert. #Anschliessend kann man diesen File z.B. mit dem IDL Programm "TDF2ENVI_float_ndvi.pro" in ENVI importieren und #weiterverarbeiten. #! /bin/csh -f cd /usr/data/maroc/all/ foreach file (noaa-1*) subset \ include_vars='ndvi_land' \ 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 ndvi_land'\ all_same_dim='yes'\ list_dims='no'\ $file.sub \ $file.ascii echo removing null data points sed '/\*\*\*/s//-9/g' $file.ascii >$file.idl #mv $file.sub /usr/data/maroc/subset_ndvi_tdf/ rm $file.ascii mv $file.idl /usr/data/maroc/idl_export/envi_maroc_sub_ndvi/ end