图像几何变换程序设计-图像处理综合训练

图像⼏何变换程序设计-图像处理综合训练
*******************
实践教案
*******************
兰州理⼯⼤学
计算机与通信学院
计算机图象处理课程设计
题⽬：图像⼏何变换程序设计
⽬录
摘要 (1)
⼀、前⾔ (2)
⼆、算法分析与描述 (2)
三、详细设计过程 (4)
四、调试过程中出现地问题及相应解决办法 (8)
五、程序运⾏截图及其说明 (9)
六、简单操作⼿册 (13)
设计总结 (14)
参考⽂献 (14)
致谢 (15)
附录 (15)
摘要
本次课程设计我地题⽬是图像⼏何变换和图像分割程序设计，计算机图像处理地实现主要以数学模型为基础，通过建⽴合适地算法来实现具体地图像处理，⼏何变换是最常见地图像处理⼿段，通过对变形地图像进⾏⼏何校正，可以得出准确地图像.常⽤地⼏何变换功能包括图像地平移、图像地旋转、图像地缩放、图像地错切等.作为数字图像处理地⼀个重要部分，⼀般⽤MATLAB编程⼯具设计⼀个完整地应⽤程序，实现相应地图像⼏何变换功能.
图像分割是⼀种⾮常重要地图像处理技术，它不仅得到了⼴泛地重视和研究，也在实际中得到了⼤量地应⽤.本⽂主要介绍了计算机图像处理中地图像分割程序地设计与实现.
关键词：图像⼏何变换、平移、旋转、放⼤缩⼩
⼀、前⾔
图形图像处理地应⽤领域涉及⼈类⽣活和⼯作地各个⽅⾯，它是从60年代以来随计算机地技术和VLSI地发展⽽产⽣、发展和不断成熟起来地⼀个新技术领域，理论上和实际应⽤上都并取得了巨⼤地成就.数字图像处理与模拟图像处理地根本不同在于，它不会因图像地存储、传输或复制等⼀系列变换操作⽽导致图像质量地退化，所以图形图像地处理在我们地⽣活中⼜很重要地作⽤.
在对图像地研究和应⽤中，⼈们往往只对图像中地某些部分感兴趣.这些部分通常称为⽬标或前景，它们⼀般对应图像中特定地、具体独特性质地区域.为了辨识和分析⽬标，需要将它们分别提取出来，在此基础上才有可能对⽬标进⼀步利⽤..
⼆、算法分析与描述
2.1图像⼏何变换
2.1.1图像地平移变换
图像地平移变换就是将图像中地像素点按照要求地量进⾏垂直、⽔平移动.图像地⽔平处理，只是改变了原有景物在画⾯上地位置，⽽图像地内容不发⽣变化.
初始坐标为(x0,y0)地点经过平移(tx,ty)(以向右，向下为正⽅向)后，坐标变为(x1,y1).这两点之间地关系是：
x1=x0+tx
y1=y0+ty
使⽤矩阵地形式来表达如下：
x1 y1 1 = x0 y0 1
值得注意地是，⼀个数字图像（灰度图）是以⼀个矩阵来描述地，因此，如果不扩⼤存放处理后地矩阵地⼤⼩，则会出现图像地部分内容移出画⾯地情况.
2.1.2图像地旋转变换
图像地旋转是指以图像中地某⼀点为原点以逆时针或顺时针⽅向翻转⼀定地⾓度，其翻转公式如下：
x1=x0cosa-y0sina；
y1=-x0sina+y0cosa；
⽤矩阵表⽰为：
x1 y1 1 = x0 y0 1
其中，（x1,y1）是原图像地像素点地坐标；（x0,y0）是对应像素点经过翻转变换后地图像地像素点地坐标.
我们旋转所在地坐标系和图像显⽰时对应地Windows屏幕坐标系是不⼀样地，这⾥xoy 为旋转坐标系，x'o'y'为屏幕坐标系.
实际上我们可以分为三步进⾏整个旋转变换：
1.将坐标系x'o'y'变成xoy；
2.将该点顺时针旋转a⾓；
3.将坐标系xoy变回x'o'y'
将上⾯三步变换进⾏合成得到三个矩阵地级联矩阵；
(x0，y0)和(x1,y1)都是x‘o’y‘坐标系中地点；
2.1.3图像地放⼤以及缩⼩
如果⼀幅图像要放⼤k1*k2倍，就是将图像中每个像素复制到k1*k2个像素所构成地⼦块中，这些⼦块再按原来地排列顺序进⾏排列，就可以实现图像地⽅地⽅法.图像缩⼩地⽅法跟图像放⼤地⽅法雷同.
三、详细设计过程
3.1图像⼏何变换
3.1.1图像平移：
1.图像平移地基本原理
图像平移是将⼀副图像中所有地点都按照指定地平移量在⽔平，垂直⽅向移动，平移后地图像与源图像相同.平移后地图像上地每⼀点都可以在原图像中找到对应地点.
X = X0 + X
Y = Y0 + Y
利⽤其次坐标，变换前后图像上地点F0（X0，Y0）和P （X,Y ）之间地关系可以⽤如下地矩阵变换表⽰为：
X0 1 0 X X Y = 0 1 Y Y 1 0 0 1 1
图像放缩地基本理论
图像⽐例缩放是指将给定地图像在X 轴⽅向按⽐例缩放FX 倍，在Y 轴⽅向上按⽐例缩放FY 倍从⽽获得⼀副新地图像.⽐例缩放前后两点F0（X0，
Y0），P （X ，Y ）之间地关系⽤矩形形式可以表⽰为：
X FX 0 0 X0 X=FxX0 Y = 0 Fy 0 y0
1 0 0 0 1 Y=FyY0
三个⼊⼝参数：原图像1，缩放倍数r,插值⽅式mode
最近邻双线性
插值插值
双三次插值3.1.2图像地旋转
3.1.3图像地放⼤
放⼤倍数为1.25倍
D=imread('wa.jpg')。

F=imresize(D,1.25)。

figure,imshow(F),title('放⼤后地图像')。

截图：
四、调试过程中出现地问题及相应解决办
法
通过MA TLAB界⾯运⾏程序进⾏图像⼏何变换时，图像有可能变换不了.可能是所需要变换地图⽚为灰度图像或者⼆值图像.需要先将RGB图像转换⼀下.
五、程序运⾏截图及其说明图像⼏何变换
图像平移：
图5.1 图象低通滤波
图5.2图像地旋转
图5.3图像地放⼤
图5.4 图像地⽔平镜像
图5.5 图像地垂直镜像
六、简单操作⼿册
打开MA TLAB软件，单击File下拉菜单.选择NEW或者Open创建⼀个新地GUI或者打开⼀个已有地GUI（后缀名为.fig地⽂件）.然后根据提⽰进⾏操作.
打开已经创建好地GUI程序.然后进⾏所需要地操作.⽐如我需要进⾏图像地⽔平翻转
设计总结
我地综合训练题⽬是图像⼏何变换程序设计，在这两周地图形图像处理综合训练过程当中，收获很多.在平常地课程学习过程当中，书上地知识感觉不是太难，只要⾃⼰认真听讲，下课后再仔细回顾回顾就能很快掌握，⾃以为什么都明⽩了，但在课设过程中却发现，要将知识实际运⽤起来却是⼗分困难地，所以在今后地学习过程中我⼀定要加强实际运⽤地能⼒训练，更好地将学来地知识转化到解决实际问题之中.通过此次综合训练，使我计算机图像处理地⼯程及其实现有了更深地理解，获得了很多宝贵地经验.特别是怎么样通过理论与实践相结合，把理论知识应⽤到我们做地程序实践过程上去.学会了很多关于计算机图像处理地经验和技巧，更重要地是，知道应如何在困难重重中⼀步⼀步细⼼地发现问题，解决问题,并在遇到问题时同学之间相互帮助中，深刻体会到了团体合作精神地重要性.
参考⽂献
[1] 朱虹.计算机图象处理基础[M]. 科学出版社, 2005
[2] R C.Gonzalez, R E.Woods著,阮秋琦,阮宇智等译.计算机图象处理(第2版).北京:电⼦⼯业出版社,2003
[3] K.R.Castleman. 计算机图象处理.北京:电⼦⼯业出版社,2002
[4] 章毓晋.图像处理与分析-图像⼯程（上册）,清华⼤学,2001
[5] 何斌等编著.Visual C++计算机图象处理.⼈民邮电出版社,2002
[6] 张宏林编著.Visual C++计算机图象模式识别技术及⼯程实践.⼈民邮电出版社,2003.
[7] 黄维通.Visual C++⾯向对象与可视化程序设计.清华⼤学出版社,2003
[8] R C.Gonzalez, R E.Woods, S L. Eddins著,阮秋琦,阮宇智等译.计算机图象处理(MATLAB版).北京:电⼦⼯业出版社,2005
致谢
⾸先我要感谢我地⽼师，在设计过程中给了我很多宝贵意见和丰富地资料并指引我阅读相关地资料和书籍，使我在不熟悉地领域中仍能迅速掌握新地技术.同时感谢我地同学，在设计中遇到问题时和我积极地讨论，并提出了很多建设性地意见，没有他们地帮助，只靠我⾃⼰⼀定很难完成这次任务.在今后地学习和⼯作中，我⼀定要发扬团队精神，齐⼼协⼒完成相应地设计⼯作.
附录
实现图像平移地⽂件地代码：
figure,imshow(F,[])。

title('平移后图形')。

figure,imshow(D)。

function pushbutton6_Callback(hObject, eventdata, handles)
global T
axes(handles.axes2)。

%定义图像显⽰位置
I=handles.img。

%选定操作图
T=getimage。

%获取图像信息
prompt={'输⼊参数:'}。

defans={'0'}。

p=inputdlg(prompt,'input',1,defans)。

%创建数据输⼊窗⼝
p1=str2num(p{1})。

J = double(I)。

%定义为双精度型
H = size(I)。

%定义尺⼨
I_moveresult = zeros(H) %平移后背景变成了⿊⾊地
x_move =p1。

y_move =p1。

I_moveresult(x_move+1:H(1),y_move+1:H(2),1:H(3))= J(1:H(1)-x_move,1:H(2)-y_move,1:H(3))。

imshow(mat2gray(I_moveresult))。

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white')。

end
global T
axes(handles.axes2)。

%定义图像显⽰位置。

I=handles.img。

%选定操作图
T=getimage。

%获取图像信息
prompt={'输⼊参数:'}。

defans={'0'}。

%定义数据放⼤缩⼩
p=inputdlg(prompt,'input',1,defans)。

%创建数据输⼊窗⼝
p1=str2num(p{1})。

f=imresize(I,p1,'nearest')。

%选⽤插值⽅法最近邻法
figure,%新建窗⼝
imshow(f)。

%显⽰结果
title('处理后')。

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white')。

End.
实现图像旋转地⽂件地代码：
function pushbutton8_Callback(hObject, eventdata, handles)
global T
axes(handles.axes2)。

%定义图像显⽰位置。

T=getimage。

%获取图像信息
prompt={'旋转⾓度:'}。

defans={'0'}。

p=inputdlg(prompt,'input',1,defans)。

%创建数据输⼊窗⼝
p1=str2num(p{1})。

%把字符串转换为数值
f=imrotate(handles.img,p1,'nearest')。

%旋转并定义插值⽅法
imshow(f)。

%显⽰结果
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white')。

End
实现图像放⼤地⽂件地代码：
global T
axes(handles.axes2)。

%定义图像显⽰位置。

I=handles.img。

%选定操作图
T=getimage。

%获取图像信息
prompt={'输⼊参数:'}。

defans={'0'}。

%定义数据放⼤缩⼩
p=inputdlg(prompt,'input',1,defans)。

%创建数据输⼊窗⼝
p1=str2num(p{1})。

f=imresize(I,p1,'nearest')。

%选⽤插值⽅法最近邻法
figure,%新建窗⼝
imshow(f)。

%显⽰结果
title('处理后')。

if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white')。

end
global T
axes(handles.axes2)。

%定义图像显⽰位置。

T=getimage。

%获取图像信息
prompt={'旋转⾓度:'}。

defans={'0'}。

p=inputdlg(prompt,'input',1,defans)。

%创建数据输⼊窗⼝
p1=str2num(p{1})。

%把字符串转换为数值
f=imrotate(handles.img,p1,'nearest')。

%旋转并定义插值⽅法
imshow(f)。

%显⽰结果
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white')。

end
实现图像镜像地⽂件地代码：
axes(handles.axes2)。

%定义图像显⽰位置。

A1=handles.img。

%选定操作图
A1=double(A1)。

%定义为双精度型
H=size(A1)。

%定义尺⼨
A3(1:H(1),1:H(2),1:H(3))=A1(1:H(1),H(2):-1:1,1:H(3))。

%⽔平镜像
imshow(uint8(A3))。

%显⽰结果
axes(handles.axes2)。

%定义图像显⽰位置。

A1=handles.img。

%选定操作图
A1=double(A1)。

%定义为双精度型
H=size(A1)。

%定义尺⼨
A2(1:H(1),1:H(2),1:H(3))=A1(H(1):-1:1,1:H(2),1:H(3))。

%垂直镜像
imshow(uint8(A2))。

%显⽰结果。