|
|
3. 标准相关
%%%%%%%%%%%%%%%%%%%%% 标准相关 %%%%%%%%%%%%%%
f1=mean(mean(fenetre.wind_coul(:,:,1)));
f2=mean(mean(fenetre.wind_coul(:,:,2)));
f3=mean(mean(fenetre.wind_coul(:,:,3)));
Sf=sqrt((fenetre.wind_coul(:,:,1)-f1).^2+(fenetre.wind_coul(:,:,2)-f2).^2+(fenetre.wind_coul(:,:,3)-f3).^2);
fchapo1=(fenetre.wind_coul(:,:,1)-f1)./Sf;
fchapo2=(fenetre.wind_coul(:,:,2)-f2)./Sf;
fchapo3=(fenetre.wind_coul(:,:,3)-f3)./Sf;
fichier2='acquisition165.jpg';
x2=imread(fichier2);
%%%%%%%% x2 亮度转型显示 %%%%%%%%%%%%%
%%%%%%%% x3=f(x2) %%%%%%%%%%%%%%%%%%
a=1.2;
b=20;
x3=a.*x2+b;
fenetrei=fenetre;
xx=fenetrei.centre(1)
yy=fenetrei.centre(2);
for i=1:20
fenetrei.centre(1)=fenetrei.centre(1)-1;
fenetrei.centre(2)=yy;
for j=1:20
fenetrei.sommets=[];
fenetrei.centre(2)=fenetrei.centre(2)-1;
[fenetrei]=InitFen(fenetrei,x2);
[fenetrej]=InitFen(fenetrei,x3);
figure(6)
clf;
figure(7)
clf;
%%%%%%%%%%% 计算标准相关 %%%%%%%%%%%%%%
g1=mean(mean(fenetrei.wind_coul(:,:,1)));
g2=mean(mean(fenetrei.wind_coul(:,:,2)));
g3=mean(mean(fenetrei.wind_coul(:,:,3)));
Sg=sqrt((fenetrei.wind_coul(:,:,1)-g1).^2+(fenetrei.wind_coul(:,:,2)-g2).^2+(fenetrei.wind_coul(:,:,3)-g3).^2);
gchapo1=(fenetrei.wind_coul(:,:,1)-g1)./Sg;
gchapo2=(fenetrei.wind_coul(:,:,2)-g2)./Sg;
gchapo3=(fenetrei.wind_coul(:,:,3)-g3)./Sg;
sommefg=(fchapo1-gchapo1).^2+(fchapo2-gchapo2).^2+(fchapo3-gchapo3).^2;
SSD(i,j)=sum(sum(sommefg));
<span lang="EN-US">if SSD(i,j) |
-
8.jpg
(35.55 KB, 下载次数: 4)
-
9.jpg
(59.24 KB, 下载次数: 6)
|