% demo_mnt_pitkin: 3D visualization of DEM.

% Set I/O paths and files
homedir=Interpolation_MNT;
input_dir=[homedir filesep 'data' filesep];
output_dir=[input_dir 'results' filesep];
mnt_file='pitkin';
occl_file=[mnt_file '.occl7'];


% Read Images
I = double(imread([input_dir mnt_file '.tif']));
Ioccl = double(imread([input_dir occl_file  '.tif']));


% Select subimage to display    
dX = 70:190; dY = 150:350;


mask=Ioccl==0;
figure(1); view3Dmask(0.5*I(dX,dY),mask(dX,dY));
           title('Original with occlusion');
