load('convex_subpixel_ptmean_03_SNR_30_sigma_1_Nb_2_Train')
wavelength=wavelength/1000;%为啥除1000呢
ZeroV = zeros([1 1000])
indices = find(GT==1)
for i=1:1000
z =(E_t-X(1:211,indices(i))).*(E_t-X(1:211,indices(i)));
dist = sqrt(sum(z));
ZeroV(i)=dist
end
[n,im]=min(ZeroV);
cols=im;
pix_col=indices(cols)
for i=1:1000
ZeroV(i)=acos(dot(E_t,X(1:211,indices(i)))/(norm(E_t)*norm(X(1:211,indices(i)))))*180/pi
end
[n,im]=min(ZeroV);
cols=im; %cols 为最小值所在列数
pix_col=indices(cols)