用MATLAB统计图像直方图.pptVIP

用MATLAB统计图像直方图 安玉磊 基本概念 灰度直方图：数字图像中各灰度级与其出现的概率的统计关系。可以表示为 且满足 基本概念 数字图像的矩阵存储格式 编程思想 1、读入图像，cameraman.tif，并显示 2、获取图像空间坐标，灰度范围 3、统计各个灰度的像素个数 4、绘制直方图 绘图：plot（x,y) x = -pi:pi/10:pi; y = tan(sin(x)) - sin(tan(x)); plot(x,y,--rs,LineWidth,2,... MarkerEdgeColor,k,... MarkerFaceColor,g,... MarkerSize,10) 字母 颜色 标点 线型 y 黄色 · 点线 m 粉红 ○ 圈线 c 亮蓝 × ×线 r 大红 ＋ ＋字线 g 绿色 － 实线 b 蓝色 * 星形线 w 白色 ： 虚线 k 黑色 －· (--) 点划线 （二）阈值平均法 为克服邻域平均使图像变模糊的缺点，可以采用加门限的方法来减少这种模糊。具体计算公式是： 其中的门限T通常选择为 ， 表示图像的均方差。但实际应用中，门限T要利用经验值和多次试验来获得。这种方法对抑制椒盐噪声比较有效，同时也能较好地保护仅有微小变化差的目标物细节。 非线性平滑处理 阈值平均法 实例 I=imread(eight.tif); J=imnoise(I,salt pepper); h=[1 1 1; 1 1 1;1 1 1]; h=h/9; K=conv2(J,h); M=double(J); for i=1:242 for j=1:308 L(i,j) = J(i,j); if (abs(K(i,j)-J(i,j)) 100) L(i,j) = K(i,j); end end end; subplot(2,2,1),imshow(I); subplot(2,2,2),imshow(J); subplot(2,2,3),imshow(K,[]);subplot(2,2,4),imshow(L,[]); I J L K 显示 4.2 图像的直方图修正 计算累积直方图： 计算变换后的灰度值： 确定灰度变换关系： 统计变换后各灰度级的像素个数 计算变换后图像的直方图： 计算原始直方图： 0.11 0.24 0.21 0.25 0.19 8 455 983 852 1020 786 7 5,6,7?7 3,4?6 2?5 1?3 0?1 6 7 7 7 6 6 5 3 1 5 1.00 0.98 0.95 0.89 0.81 0.65 0.44 0.19 4 0.02 0.03 0.06 0.08 0.16 0.21 0.25 0.19 3 80 130 245 333 650 852 1020 786 统计原图像各灰度级像素个数ni 2 7 6 5 4 3 2 1 0 列出图像灰度级（i或j） 1 计算结果 计算方法或公式 步骤 上机一 熟悉MatLab环境及基本操作 作业 1、图像的基本操作 读图像文件及显示图像 2、绘制图像的直方图及直方图均衡变换 自己编程求解图像的直方图 3、自己退化一幅图像，并用维纳滤波复原 1、无噪声 2、有噪声