#Dieses Script deklariert alle Pixel im Meeresbreich als "badval", so dass bei Klassifikation der
#Landoberfläche keine Komplikationene auftreten können. NOAA  16 (avhrr_ch3a!!!)
#! /bin/csh -f

cd /usr/data/masked_files/
foreach file (noaa-16*)
emathp \
	expr_vars='avhrr_ch1 avhrr_ch2 avhrr_ch3a avhrr_ch4 avhrr_ch5 latitude longitude' \
	num_exprs='12' \
	y1_expr='x1' \
	y2_expr='x2' \
	y3_expr='x3' \
	y4_expr='x4' \
	y5_expr='x5' \
	y6_expr='x6' \
	y7_expr='x7' \
	y8_expr='land ? y1 : badval' \
	y9_expr='land ? y2 : badval' \
	y10_expr='land ? y3 : badval' \
	y11_expr='land ? y4 : badval' \
	y12_expr='land ? y5 : badval' \
	use_degrees='n' \
	poly_size='100' \
	save_exprs='8,9,10,11,12' \
	var_names='avhrr_ch1 avhrr_ch2 avhrr_ch3a avhrr_ch4 avhrr_ch5' \
	var_units='albedo*100% albedo*100% temp_deg_c temp_deg_c temp_deg_c' \
	var_types='sh sh sh sh sh' \
	scale_offsets='0 0 0 0 0' \
	scale_factors='1 1 1 1 1' \
$file $file.land

rm $file
mv -f $file.land $file

end



#mv $file.land /usr/data/maroc/masked_files/ohne_meer/
#cd /usr/data/maroc/masked_files/ohne_meer/
#mv -f $file.land $file