#Dieses Script prozessiert im Verzeichnis /yukon/ftp/NOAA/  die via fpt
#zugesandten NOAA FILES, d.h. sie werden navigiert und georeferenziert, um so
#eine  erste Qualitätsanalyse der Files z.B. in TVIS durchführen zu können.
#! /bin/csh -fx

cd /yukon/ftp/NOAA/
foreach file (NSS.LHRR.NL.*)

echo hrptin

hrptin \
        header_size='16384' \
        hrpt_format='klmlacp' \
        hrpt_types='avhrr' \
        on_pass_disk='n' \
        tape_device=$file \
        pass_year='2002' \
        channels='1 2 3 4 5' \
        fix_missing='yes' \
        sharpen='y y y y y' \
        calibrate='yes' \
        temp_units='celsius' \
        byte_output='no' \
        full_report='no' \
        solar_corr='no' \
        geo_correct='y' \
        delta_line='1' \
        delta_sample='1' \
        use_master='yes' \
        master_file='/wolga/TERA_SCRIPT/MOROCCO/Master_marokko_utmwgs84.mas' \
$file \
$file.avhrr


navbox2 \
	variable= \
	coast_file=dcw.coast \
	box_height=50 \
	box_width=50 \
	min_cross_dev=2.5 \
	min_feature=15 \
	max_regions=4 \
	min_region_perc=10 \
	min_coast_pts=40 \
	min_sat_elev=0 \
	min_good_boxes=1 \
$file.avhrr \
navboxes

nav2 \
        expr_vars='avhrr_ch1 avhrr_ch2' \
        expression='x1' \
        box_file='navboxes' \
        coast_file='dcw.coast' \
        max_coast_pts='10000' \
        max_line_shift='15' \
        max_samp_shift='15' \
        min_good_boxes='1' \
        fix_attitude='yes' \
        reset_sat='yes' \
        use_time='yes' \
        method='correlate' \
        min_corr='60' \
        interp_corr='yes' \
        min_feature='4' \
        fixed_thresh='no' \
        stdev_units='0' \
        num_bins='256' \
        estimate_range='yes' \
$file.avhrr

nav2 \
        expr_vars='avhrr_ch1 avhrr_ch2' \
        expression='x1' \
        box_file='navboxes' \
        coast_file='dcw.coast' \
        max_line_shift='15' \
        max_samp_shift='15' \
        min_good_boxes='1' \
        fix_attitude='yes' \
        reset_sat='yes' \
        use_time='yes' \
        method='edge' \
        min_coast_pts='40' \
	min_edge_diff='1' \
	min_edge_perc='90' \
$file.avhrr

echo fastreg2

fastreg2 \
        master_file='/wolga/TERA_SCRIPT/MOROCCO/Master_marokko_utmwgs84.mas' \
        master_var= \
        include_vars= \
        brute_force='no' \
        poly_size='100' \
        interpolate='nn' \
        off_protect='no' \
        whole_input='no' \
	add_to_output='no' \
$file.avhrr \
$file.geo

set file1=`echo $file.geo | cut -c14-18`
cp $file.geo ${file1}_utmwgs84

set satellite=`getinfo attr_name=satellite ${file1}_utmwgs84` 
cp ${file1}_utmwgs84 ${file1}_${satellite}_utmwgs84

set projection_name=`getinfo attr_name=projection_name ${file1}_utmwgs84`
mv ${file1}_${satellite}_utmwgs84 {$file1}_${satellite}_${projection_name}


rm ${file1}_utmwgs84
rm $file.avhrr
rm $file.geo
rm navboxes
end