matlab 图像二维直方图

clc;
clear all;
close all;
tic
I=imread('corner.png');
% I=double(I);
 f=rgb2gray(I);
 figure,imshow(f),title('原图像');
  figure,imhist(f),title('一维直方图');
[m,n,d]=size(f)
A=f;
for i=1:1:256
    for j=1:1:256
        twodim2(i,j)=0;
    end
end
 J = imnoise(f,'gaussian',0.05);
 figure,imshow(J),title('加躁后的图像');
 figure,imhist(J),title('一维直方图');
H=fspecial('average');
 k=imfilter(J,H);
% k=medfilt2(A)
for j=1:1:n;
    for i=1:1:m;
        m1=J(i,j);
        m2=k(i,j);
        twodim2(m1+1,m2+1)= twodim2(m1+1,m2+1)+1;
        
    end
end
i=1:1:256;
j=1:1:256;
figure;
 mesh(i,j,twodim2(i,j)/(m*n));

title('二维直方图')
toc


rusult:

 <