0001 function [fourier_coeff_2d, sorted_values_prof, x, y, template, mean_pdf_profile, no_values_x, no_values_y] = load_2d_data_vertical(number)
0002
0003
0004 switch number
0005 case 1
0006
0007
0008 load('/home/bengt/CIWSIR/Dataset_gen/CloudSAT_data/test_2d23d.mat')
0009 template=IWC;
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023 otherwise
0024
0025
0026 load data_cloud_prof.mat
0027 template = double(lwc);
0028 end
0029
0030 [no_values_y, no_values_x] = size(template);
0031 x = 1:no_values_x;
0032 y = 1:no_values_y;
0033 mean_pdf_profile = squeeze(mean(template, 2));
0034
0035
0036 template = remove_average_profile(template, mean_pdf_profile);
0037 sorted_values_prof = sort(template, 2);
0038 sorted_values_prof = sorted_values_prof';
0039 total_variance_pdf = std(sorted_values_prof(:)).^2;
0040
0041
0042 fourier_coeff_2d = abs(ifft2( template ));
0043 power = fourier_coeff_2d.^2;
0044 total_variance_spec = sum(sum(power));
0045 power = power * total_variance_pdf / total_variance_spec;
0046 fourier_coeff_2d = sqrt(power');