|
wolffd@0
|
1 function [dis_s12 dis_s13] = conv_euclidean_dist(index_list,vector_list,index_map,w,norm)
|
|
wolffd@0
|
2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
wolffd@0
|
3 % compute the distance vector from 2 vector using convolution euclidean %
|
|
wolffd@0
|
4 % w = convolution window %
|
|
wolffd@0
|
5 % sontran2013 %
|
|
wolffd@0
|
6 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
wolffd@0
|
7 s_num = size(index_list,1);
|
|
wolffd@0
|
8 s_len = size(vector_list,2);
|
|
wolffd@0
|
9
|
|
wolffd@0
|
10 dis_s12 = zeros(s_num,s_len);
|
|
wolffd@0
|
11 dis_s13 = zeros(s_num,s_len);
|
|
wolffd@0
|
12
|
|
wolffd@0
|
13 % Convert list of assigned indices to ordered indices
|
|
wolffd@0
|
14 [checker index_list] = ismember(index_list,index_map);
|
|
wolffd@0
|
15 cv_vec = ones(1,w);
|
|
wolffd@0
|
16 if all(all(checker(:,1:3)))
|
|
wolffd@0
|
17 dis_s12 = sqrt(conv2([1],cv_vec,(vector_list(index_list(:,1),:) - vector_list(index_list(:,2),:)).^2));
|
|
wolffd@0
|
18 dis_s13 = sqrt(conv2([1],cv_vec,(vector_list(index_list(:,1),:) - vector_list(index_list(:,3),:)).^2));
|
|
wolffd@0
|
19 % dis_s12 = dis_s12(:,w:end-w);
|
|
wolffd@0
|
20 % dis_s13 = dis_s13(:,w:end-w);
|
|
wolffd@0
|
21 if norm
|
|
wolffd@0
|
22 dis_s12 = dis_s12./(conv2([1],cv_vec,ones(s_num,s_len)));
|
|
wolffd@0
|
23 dis_s13 = dis_s13./(conv2([1],cv_vec,ones(s_num,s_len)));
|
|
wolffd@0
|
24 end
|
|
wolffd@0
|
25 else
|
|
wolffd@0
|
26 fprintf('Some indices are not mapped\n');
|
|
wolffd@0
|
27 end
|
|
wolffd@0
|
28 end |
