dimage_gaussiantest.m

来自「包括相邻帧差分结果的高斯性和独立性检验实验结果演示程序。」· M 代码 · 共 84 行

clear;
clc;

%读入连续四帧图像
im1=double(imread('img1.bmp'));
im2=double(imread('img2.bmp'));
im3=double(imread('img3.bmp'));
im4=double(imread('img4.bmp'));

[ny,nx]=size(im1);
ax=subplot(4,4,1);
colormap(gray(256));
image(im1);
title('原图像');

ax=subplot(4,4,5);
image(im2);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);
    
ax=subplot(4,4,9);
image(im3);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);
    
ax=subplot(4,4,13);
image(im4);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);


ax=subplot(4,4,6);
image(im1);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);
title('估计图像');

ax=subplot(4,4,10);
image(im2);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);

ax=subplot(4,4,14);
image(im3);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);

ax=subplot(4,4,7);

tmp1=(im2-im1);
yn=tmp1((ny/2):(ny/2+99),(nx/2):(nx/2+99));


image(tmp1);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);   
title('残留噪声');

ax=subplot(4,4,8);
[tstar,meanval]=kendallrankcorr(tmp1(120:140,160:170));
title ('残留噪声经验概率密度函数')

ax=subplot(4,4,11);

tmp2=im3-im2;
yn=tmp2((ny/2):(ny/2+99),(nx/2):(nx/2+99));

image(tmp2);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);   
ax=subplot(4,4,12);
[tstar,meanval]=kendallrankcorr(tmp2(120:140,160:170));

    
ax=subplot(4,4,15);
tmp3=im4-im3;
yn=tmp3((ny/2):(ny/2+99),(nx/2):(nx/2+99));

image(tmp3);
set(ax,'Xtick',[],'Xticklabel',[],...
        'Ytick',[],'Yticklabel',[]);   
ax=subplot(4,4,16);
[tstar,meanval]=kendallrankcorr(tmp3(120:140,160:170));