会计学1用MATLAB统计图像直方图基本概念灰度直方图:数字图像中各灰度级与其出现的概率的统计关系。可以
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示为且满足第1页/共10页基本概念数字图像的矩阵存储格式第2页/共10页编程思想1、读入图像,cameraman.tif,并显示2、获取图像空间坐标,灰度范围3、统计各个灰度的像素个数4、绘制直方图第3页/共10页绘图: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)第4页/共10页字母颜色标点线型y黄色·点线m粉红○圈线c亮蓝××线r大红++字线g绿色-实线b蓝色*星形线w白色:虚线k黑色-·(--)点划线第5页/共10页(二)阈值平均法为克服邻域平均使图像变模糊的缺点,可以采用加门限的方法来减少这种模糊。具体计算公式是:其中的门限T通常选择为,表示图像的均方差。但实际应用中,门限T要利用经验值和多次试验来获得。这种方法对抑制椒盐噪声比较有效,同时也能较好地保护仅有微小变化差的目标物细节。第6页/共10页非线性平滑处理阈值平均法实例>>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显示第7页/共10页4.2图像的直方图修正计算累积直方图:计算变换后的灰度值:确定灰度变换关系:统计变换后各灰度级的像素个数计算变换后图像的直方图:计算原始直方图:0.110.240.210.250.198455983852102078675,6,773,4625130167776653151.000.980.950.890.810.650.440.1940.020.030.060.080.160.210.250.193801302453336508521020786统计原图像各灰度级像素个数ni276543210列出图像灰度级(i或j)1计算结果计算方法或公式步骤第8页/共10页上机一熟悉MatLab环境及基本操作作业1、图像的基本操作读图像文件及显示图像2、绘制图像的直方图及直方图均衡变换自己编程求解图像的直方图3、自己退化一幅图像,并用维纳滤波复原1、无噪声2、有噪声第9页/共10页