图像信息处理之图像旋转和阈值分割

1、图像信息处理上机实验报告图像旋转、阈值分割1、实验内容：（1）写出以图像中心点作为支点旋转任意角度，图像仍保持完整的旋转程序；（2）对Fig33图，采用灰度直方图峰谷方法分离出米粒，要求编写程序自动求直方图峰谷对应的分割阈值;(3)练习PS中把jimoxiong变换到jimoxiong2，或者反向恢复，并证明你的操作是成功的，写出基本步骤。2、具体实验一、图像旋转：实验原理：（1） ：计算公式：i=icos-jsin j=icos+jcos(i,j)为原图像素点坐标，（i,j）是对应像素点经旋转变换后图像像素点坐标；（2） ：对计算所得坐标值取整；（3） 对取整后坐标值范围进行画布扩大；（4）

2、 对图像出现的空穴进行实验程序：function r=imrotate(I,theta)I = imread(L:图像信息处理17周上机Fig33.BMP);row,col= size(I);theta=(30*pi/180); %旋转30度 cosa=cos(theta); sina=sin(theta); %原图四个顶点坐标 srcx1=-col*0.5;srcy1=row*0.5; srcx2=col*0.5;srcy2=row*0.5; srcx3=-col*0.5;srcy3=-row*0.5; srcx4=col*0.5;srcy4=-row*0.5; %图像旋转后四个顶点坐标 d

3、stx1=cosa*srcx1-sina*srcy1; dsty1=sina*srcx1+cosa*srcy1; dstx2=cosa*srcx2-sina*srcy2; dsty2=sina*srcx2+cosa*srcy2; dstx3=cosa*srcx3-sina*srcy3; dsty3=sina*srcx3+cosa*srcy3; dstx4=cosa*srcx4-sina*srcy4; dsty4=sina*srcx4+cosa*srcy4; %计算旋转后图像宽度和高度 %旋转后图像高为对应顶点纵坐标之差较大值 h=max(abs(dsty4-dsty1),abs(dsty2-d

4、sty3)+0.5; %旋转后图像宽为对应顶点纵坐标之差较大值 w=max(abs(dstx4-dstx1),abs(dstx2-dstx3)+0.5; %对新图高宽进行四舍五入 h=floor(h); w=floor(w); r=zeros(h,w); f1=-w*0.5*cosa-h*0.5*sina+0.5*col; f2=w*0.5*sina-h*0.5*cosa+0.5*row; for x=1:w for y=1:h x0=floor(x*cosa+y*sina+f1); y0=floor(-x*sina+y+cosa+f2);%对图像进行填充 if x00&x00&y0=T S1

5、=S1+double(Im(i,j); %大于阈域值图像点灰度值累加n1=n1+1; %大于阈域值图像点个数累加else S0=S0+double(Im(i,j); %小于阈域值图像点灰度值累加n0=n0+1; %小于阀域值图像点个数累加 end end end T0=S0/n0;T1=S1/n1; if abs(T-(T0+T1)/2)0.1 %迭代至前后两次阀域值相差几乎为0时停止迭代。break; else T=(T0+T1)/2; %在阈值T下，迭代阈值的计算过程end end i1=im2bw(Im,T/255); % 图像在最佳阈值下二值化 subplot(1,3,3),imsho

6、w(i1) title(结果)zd = 255zx = 0T = 75.7379实验结果：三、ps原图：实验结果：步骤：编辑新建画布选择合适尺寸，并将将背景设为白色 打开原图将图层解锁调整原图大小点击移动将原图拖动到新建图层中 水平翻转自由变换调整旋转到合适位置 合并图层并导出3、总结与分析：实验中由于自己对软件的不熟悉，使用时就总会遇到一些很细节性的问题，比如全角、半角符号、函数名、使用规则、存放路径等等很细小的问题，往往会在这上面花费很多很多时间；另一个问题便是自己的逻辑思维能力了，拿到题目，用脑子能想到怎么做，但是具体要编程转化成语言时就很难了，自己往往纠结半天也出不来结果，总要借助辅导书才能完成。Phtoshop部分自己做得时候往往也都是辅导书和网络一起辅助进行，这个过程，自己确实也熟悉了它的运行环境