echo on;
% The program WAVELET3D computes and visualizes a 3-D Wavelet
% of the specific data set.
% 
% Data loaded is:    'p71301.dat'- 'p713100.dat'.
% Note: see marked lines in the code to change cases.
%
% Calling sequence-
% wavelet3d
%
% N. E. Huang (NASA GSFC)		14.Feb.2001
% --- wavelet3d.m --- Version  14.Feb.2001 ------------------


clear;						% Start Fresh

load p71301.dat; 				% Load Data
load p71302.dat; 				% Load Data
load p71303.dat; 				% Load Data
load p71304.dat; 				% Load Data
load p71305.dat; 				% Load Data
load p71306.dat; 				% Load Data
load p71307.dat; 				% Load Data
load p71308.dat; 				% Load Data
load p71309.dat; 				% Load Data
load p713010.dat; 				% Load Data
load p713011.dat; 				% Load Data
load p713012.dat; 				% Load Data
load p713013.dat; 				% Load Data
load p713014.dat; 				% Load Data
load p713015.dat; 				% Load Data
load p713016.dat; 				% Load Data
load p713017.dat; 				% Load Data
load p713018.dat; 				% Load Data
load p713019.dat; 				% Load Data
load p713020.dat; 					
load p713021.dat; 	
load p713022.dat; 	
load p713023.dat; 	
load p713024.dat; 	
load p713025.dat; 	
load p713026.dat; 	
load p713027.dat; 	
load p713028.dat; 	
load p713029.dat; 	
load p713030.dat;
load p713031.dat; 
load p713032.dat; 
load p713033.dat; 
load p713033.dat; 
load p713034.dat; 
load p713035.dat; 
load p713036.dat; 
load p713037.dat; 
load p713038.dat; 
load p713039.dat;  	
load p713040.dat; 	
load p713041.dat; 
load p713042.dat; 
load p713043.dat; 
load p713044.dat; 
load p713045.dat; 
load p713046.dat; 
load p713047.dat; 
load p713048.dat; 
load p713049.dat; 
load p713050.dat; 
load p713051.dat; 
load p713052.dat; 
load p713053.dat; 
load p713054.dat; 
load p713055.dat; 
load p713056.dat; 
load p713057.dat; 
load p713058.dat; 
load p713059.dat; 
load p713060.dat;
load p713061.dat; 
load p713062.dat; 
load p713063.dat;  
load p713064.dat; 
load p713065.dat; 
load p713066.dat; 
load p713067.dat; 
load p713068.dat; 
load p713069.dat; 
load p713070.dat; 
load p713071.dat; 
load p713072.dat; 
load p713073.dat; 
load p713074.dat; 
load p713075.dat; 
load p713076.dat; 
load p713077.dat; 
load p713078.dat; 
load p713079.dat; 
load p713080.dat;
load p713081.dat; 
load p713082.dat; 
load p713083.dat; 
load p713084.dat; 
load p713085.dat; 
load p713086.dat; 
load p713087.dat; 
load p713088.dat; 
load p713089.dat; 
load p713090.dat; 
load p713091.dat;  
load p713092.dat;  
load p713093.dat;  
load p713094.dat;  
load p713095.dat;  
load p713096.dat;  
load p713097.dat;  
load p713098.dat;  
load p713099.dat;  
load p7130100.dat;

w1=wt2(p71301, 5.5, 4, 512, 100);	  % Compute the Wavelet for each line
w2=wt2(p71302, 5.5, 4, 512, 100);	
w3=wt2(p71303, 5.5, 4, 512, 100);	
w4=wt2(p71304, 5.5, 4, 512, 100);	
w5=wt2(p71305, 5.5, 4, 512, 100);	
w6=wt2(p71306, 5.5, 4, 512, 100);	
w7=wt2(p71307, 5.5, 4, 512, 100);	
w8=wt2(p71308, 5.5, 4, 512, 100);	
w9=wt2(p71309, 5.5, 4, 512, 100);	
w10=wt2(p713010, 5.5, 4, 512, 100);	
w11=wt2(p713011, 5.5, 4, 512, 100);
w12=wt2(p713012, 5.5, 4, 512, 100);	
w13=wt2(p713013, 5.5, 4, 512, 100);	
w14=wt2(p713014, 5.5, 4, 512, 100);	
w15=wt2(p713015, 5.5, 4, 512, 100);	
w16=wt2(p713016, 5.5, 4, 512, 100);	
w17=wt2(p713017, 5.5, 4, 512, 100);	
w18=wt2(p713018, 5.5, 4, 512, 100);	
w19=wt2(p713019, 5.5, 4, 512, 100);	
w20=wt2(p713020, 5.5, 4, 512, 100);
w21=wt2(p713021, 5.5, 4, 512, 100);		
w22=wt2(p713022, 5.5, 4, 512, 100);
w23=wt2(p713023, 5.5, 4, 512, 100);		
w24=wt2(p713024, 5.5, 4, 512, 100);
w25=wt2(p713025, 5.5, 4, 512, 100);		
w26=wt2(p713026, 5.5, 4, 512, 100);
w27=wt2(p713027, 5.5, 4, 512, 100);		
w28=wt2(p713028, 5.5, 4, 512, 100);
w29=wt2(p713029, 5.5, 4, 512, 100);			
w30=wt2(p713030, 5.5, 4, 512, 100);
w31=wt2(p713031, 5.5, 4, 512, 100);
w32=wt2(p713032, 5.5, 4, 512, 100);
w33=wt2(p713033, 5.5, 4, 512, 100);	
w34=wt2(p713034, 5.5, 4, 512, 100);
w35=wt2(p713035, 5.5, 4, 512, 100);	
w36=wt2(p713036, 5.5, 4, 512, 100);
w37=wt2(p713037, 5.5, 4, 512, 100);	
w38=wt2(p713038, 5.5, 4, 512, 100);
w39=wt2(p713039, 5.5, 4, 512, 100);			
w40=wt2(p713040, 5.5, 4, 512, 100);
w41=wt2(p713041, 5.5, 4, 512, 100);
w42=wt2(p713042, 5.5, 4, 512, 100);
w43=wt2(p713043, 5.5, 4, 512, 100);	
w44=wt2(p713044, 5.5, 4, 512, 100);
w45=wt2(p713045, 5.5, 4, 512, 100);	
w46=wt2(p713046, 5.5, 4, 512, 100);
w47=wt2(p713047, 5.5, 4, 512, 100);	
w48=wt2(p713048, 5.5, 4, 512, 100);
w49=wt2(p713049, 5.5, 4, 512, 100);	
w50=wt2(p713050, 5.5, 4, 512, 100);
w51=wt2(p713051, 5.5, 4, 512, 100);
w52=wt2(p713052, 5.5, 4, 512, 100);
w53=wt2(p713053, 5.5, 4, 512, 100);	
w54=wt2(p713054, 5.5, 4, 512, 100);
w55=wt2(p713055, 5.5, 4, 512, 100);	
w56=wt2(p713056, 5.5, 4, 512, 100);
w57=wt2(p713057, 5.5, 4, 512, 100);	
w58=wt2(p713058, 5.5, 4, 512, 100);
w59=wt2(p713059, 5.5, 4, 512, 100);	
w60=wt2(p713060, 5.5, 4, 512, 100);
w61=wt2(p713061, 5.5, 4, 512, 100);
w62=wt2(p713062, 5.5, 4, 512, 100);
w63=wt2(p713063, 5.5, 4, 512, 100);	
w64=wt2(p713064, 5.5, 4, 512, 100);
w65=wt2(p713065, 5.5, 4, 512, 100);	
w66=wt2(p713066, 5.5, 4, 512, 100);
w67=wt2(p713067, 5.5, 4, 512, 100);	
w68=wt2(p713068, 5.5, 4, 512, 100);
w69=wt2(p713069, 5.5, 4, 512, 100);	
w70=wt2(p713070, 5.5, 4, 512, 100);
w71=wt2(p713071, 5.5, 4, 512, 100);	
w72=wt2(p713072, 5.5, 4, 512, 100);
w73=wt2(p713073, 5.5, 4, 512, 100);
w74=wt2(p713074, 5.5, 4, 512, 100);
w75=wt2(p713075, 5.5, 4, 512, 100);	
w76=wt2(p713076, 5.5, 4, 512, 100);
w77=wt2(p713077, 5.5, 4, 512, 100);	
w78=wt2(p713078, 5.5, 4, 512, 100);
w79=wt2(p713079, 5.5, 4, 512, 100);	
w80=wt2(p713080, 5.5, 4, 512, 100);
w81=wt2(p713081, 5.5, 4, 512, 100);
w82=wt2(p713082, 5.5, 4, 512, 100);
w83=wt2(p713083, 5.5, 4, 512, 100);	
w84=wt2(p713084, 5.5, 4, 512, 100);
w85=wt2(p713085, 5.5, 4, 512, 100);	
w86=wt2(p713086, 5.5, 4, 512, 100);
w87=wt2(p713087, 5.5, 4, 512, 100);	
w88=wt2(p713088, 5.5, 4, 512, 100);
w89=wt2(p713089, 5.5, 4, 512, 100);	
w90=wt2(p713090, 5.5, 4, 512, 100);
w91=wt2(p713091, 5.5, 4, 512, 100);	
w92=wt2(p713092, 5.5, 4, 512, 100);
w93=wt2(p713093, 5.5, 4, 512, 100);
w94=wt2(p713094, 5.5, 4, 512, 100);	
w95=wt2(p713095, 5.5, 4, 512, 100);
w96=wt2(p713096, 5.5, 4, 512, 100);	
w97=wt2(p713097, 5.5, 4, 512, 100);	
w98=wt2(p713098, 5.5, 4, 512, 100);
w99=wt2(p713099, 5.5, 4, 512, 100);					
w100=wt2(p7130100, 5.5,  4, 512, 100);

% Construct the 3-d Wavelt; the axies and the Wavelet volume

wx=linspace(0, 26.54, 170)';

k=linspace( 1/26.54, 512/(26.54*4), 100)';



W=zeros(100,170, 100);			% Make & Assemble Data Volume

W(:, :, 1)=w1;
W(:, :, 2)=w2;
W(:, :, 3)=w3;
W(:, :, 4)=w4;
W(:, :, 5)=w5;
W(:, :, 6)=w6;
W(:, :, 7)=w7;
W(:, :, 8)=w8;
W(:, :, 9)=w9;
W(:, :, 10)=w10;
W(:, :, 11)=w11;
W(:, :, 12)=w12;
W(:, :, 13)=w13;
W(:, :, 14)=w14;
W(:, :, 15)=w15;
W(:, :, 16)=w16;
W(:, :, 17)=w17;
W(:, :, 18)=w18;
W(:, :, 19)=w19;
W(:, :, 20)=w20;
W(:, :, 21)=w21;
W(:, :, 22)=w22;
W(:, :, 23)=w23;
W(:, :, 24)=w24;
W(:, :, 25)=w25;
W(:, :, 26)=w26;
W(:, :, 27)=w27;
W(:, :, 28)=w28;
W(:, :, 29)=w29;
W(:, :, 30)=w30;
W(:, :, 31)=w31;
W(:, :, 32)=w32;
W(:, :, 33)=w33;
W(:, :, 34)=w34;
W(:, :, 35)=w35;
W(:, :, 36)=w36;
W(:, :, 37)=w37;
W(:, :, 38)=w38;
W(:, :, 39)=w39;
W(:, :, 40)=w40;
W(:, :, 41)=w41;
W(:, :, 42)=w42;
W(:, :, 43)=w43;
W(:, :, 44)=w44;
W(:, :, 45)=w45;
W(:, :, 46)=w46;
W(:, :, 47)=w47;
W(:, :, 48)=w48;
W(:, :, 49)=w49;
W(:, :, 50)=w50;
W(:, :, 51)=w51;
W(:, :, 52)=w52;
W(:, :, 53)=w53;
W(:, :, 54)=w54;
W(:, :, 55)=w55;
W(:, :, 56)=w56;
W(:, :, 57)=w57;
W(:, :, 58)=w58;
W(:, :, 59)=w59;
W(:, :, 60)=w60;
W(:, :, 61)=w61;
W(:, :, 62)=w62;
W(:, :, 63)=w63;
W(:, :, 64)=w64;
W(:, :, 65)=w65;
W(:, :, 66)=w66;
W(:, :, 67)=w67;
W(:, :, 68)=w68;
W(:, :, 69)=w69;
W(:, :, 70)=w70;
W(:, :, 71)=w71;
W(:, :, 72)=w72;
W(:, :, 73)=w73;
W(:, :, 74)=w74;
W(:, :, 75)=w75;
W(:, :, 76)=w76;
W(:, :, 77)=w77;
W(:, :, 78)=w78;
W(:, :, 79)=w79;
W(:, :, 80)=w80;
W(:, :, 81)=w81;
W(:, :, 82)=w82;
W(:, :, 83)=w83;
W(:, :, 84)=w84;
W(:, :, 85)=w85;
W(:, :, 86)=w86;
W(:, :, 87)=w87;
W(:, :, 88)=w88;
W(:, :, 89)=w89;
W(:, :, 90)=w90;
W(:, :, 91)=w91;
W(:, :, 92)=w92;
W(:, :, 93)=w93;
W(:, :, 94)=w94;
W(:, :, 95)=w95;
W(:, :, 96)=w96;
W(:, :, 97)=w97;
W(:, :, 98)=w98;
W(:, :, 99)=w99;
W(:, :, 100)=w100;

% --- Averaging Section --- %

W(:, :, 1) = (W(:, :, 1) + W(:, :, 2)) ./ 2.;
W(:, :, 100) = (W(:, :, 99) + W(:, :, 100)) ./ 2.;
W(:, :, 2) = (4. .* W(:, :, 2) + W(:, :, 1) + W(:, :, 3)) ./ 6.;
W(:, :, 99) = (4. .* W(:, :, 99) + W(:, :, 98) + W(:, :, 100)) ./ 6.;

for i = 3:98;
	W(:, :,i) = (6. .* W(:, :, i) + 2. .* (W(:, :, i-1) +  W(:, :, i+1)) ...
	             +  W(:, :, i-2) +  W(:, :, i+2)) ./ 12.;
end;

% --- Save result for further plotting ---- %

save Wvlt W wx k;

tt=linspace(15, 25, 100)';

% --- Change this Value -----> x.xxx <---------- Change this Value -- %

p = patch(isosurface(wx,k,tt,W, 0.02));

% ------------------------------------------------------------------- %

isonormals(wx,k, tt, W, p)
set(p, 'FaceColor', 'cyan', 'EdgeColor', 'none');
daspect([.5 .2 0.25])
campos([-45 10 90]);
lighting phong;
grid off;
camlight left

% --- Change the Title ---> Case   Frames  Line No. Contour Level --- %
title( 'X Distance: 0-26 cm; Y Distance: 15-25 cm; Wavenumber: 0-4.5 1/cm')
toptitle('3D Wavelet Spectrum: NEH-07 F130 L193 to 492 C.02')
% ------------------------------------------------------------------- %

% --- swavelet3d.m Ends Normally --- %