用MATLAB统计图像直方图 10页

用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=[111;111;111];>>h=h/9;>>K=conv2(J,h);>>M=double(J);>>fori=1:242forj=1:308L(i,j)=J(i,j);if(abs(K(i,j)-J(i,j))>100)L(i,j)=K(i,j);endendend;>>subplot(2,2,1),imshow(I);subplot(2,2,2),imshow(J);>>subplot(2,2,3),imshow(K,[]);subplot(2,2,4),imshow(L,[]);IJLK显示 4.2图像的直方图修正计算累积直方图：计算变换后的灰度值：确定灰度变换关系：统计变换后各灰度级的像素个数计算变换后图像的直方图：计算原始直方图：0.110.240.210.250.198455983852102078675,6,773,4625130167776653151.000.980.950.890.810.650.440.1940.020.030.060.080.160.210.250.193801302453336508521020786统计原图像各灰度级像素个数ni276543210列出图像灰度级（i或j）1计算结果计算方法或公式步骤 上机一熟悉MatLab环境及基本操作作业1、图像的基本操作读图像文件及显示图像2、绘制图像的直方图及直方图均衡变换自己编程求解图像的直方图3、自己退化一幅图像，并用维纳滤波复原1、无噪声2、有噪声