数字图像处理几何变换讲稿.ppt

《数字图像处理几何变换讲稿.ppt》由会员分享，可在线阅读，更多相关《数字图像处理几何变换讲稿.ppt（87页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于数字图像处理几何变换关于数字图像处理几何变换第一页，讲稿共八十七页哦本章主要内容和基本要求数字图像的基本代数运算数字图像的几何运算第二页，讲稿共八十七页哦5.1 数字图像处理基本运算的分类一、图像处理算法的分类：单幅图像单幅图像多幅图像单幅图像多幅图像或单幅图像数值/符号（图像分析）二、根据基本运算分类点运算邻域运算第三页，讲稿共八十七页哦5.2 点运算一、数字图像的代数运算1、基本代数运算基础：设图像其运算可以分为图像自身的运算和图像间的运算，从根本上而言可以看成对图像位置的运算或对该位置上灰度级的运算。第四页，讲稿共八十七页哦基本代数运算：图像象素几何位置不变图像灰度级的加、减、乘、除

2、等运算。几何运算：图像象素位置变换后，经过计算，确定该象素灰度的运算。例如图像f和g：第五页，讲稿共八十七页哦 代数运算就是二者点对点的加、减、乘、除运算。对f图所有灰度加c级灰度可作如下处理：J=f+c第六页，讲稿共八十七页哦给整幅图像增加灰度级会使图像亮度得到提高整体偏亮，给个别象素加灰度值可以使目标景物突出，例如：第七页，讲稿共八十七页哦给一幅图像乘以一个系数，可以改变图像的灰度级范围如：原来是20200的灰度级范围，乘以2则有40400的灰度级范围第八页，讲稿共八十七页哦第九页，讲稿共八十七页哦第十页，讲稿共八十七页哦两幅连续的图像运算可以有基本运算组合而成如：Kafbg：第十一页，讲

3、稿共八十七页哦二、代数运算的用途可以对同一场景的多幅图像求平均，以降低加性零值平稳噪声。两幅图像叠加达到二次曝光的效果。减去图中不需要的加性噪声。计算两幅场景的变化。检测轮廓，确定轮廓现。第十二页，讲稿共八十七页哦5.3 图像的几何运算几何运算也叫几何变换：包括图象的平移，旋转，镜象变换，转置，放缩等。一平移(translation)将图像沿水平或垂直方向移动位置后，获得新的图像的变换方式。例如：图像坐标中象素点（x0,y0)移动到（x1,y1)位置灰度不变。第十三页，讲稿共八十七页哦第十四页，讲稿共八十七页哦（x0,y0)点，经过平移（tx,ty），变为点（x1,y1）两点之间的关系 为：x

4、1=x0+tx;y1=y0+ty.以矩阵的形式表示为:第十五页，讲稿共八十七页哦我们更关心的是它的逆变换：因为：我们想知道的是平移后的图象中每个像素的颜色。x0=x1-tx,y0=y1-ty第十六页，讲稿共八十七页哦例如：新图中左上角点的RGB值是多少？很显然，该点是原图中的某一点经过平移后得到的，这两点的颜色肯定是一样的，所以只要知道了原图那点的RGB值即可。那么到底新图中的左上角点对应原图中的哪一点呢？将左上角点的坐标（0,0）代入公式(2)，得到x0=-tx；y0=-ty；所以新图中的(0,0)点的颜色和原图中（-tx,-ty）的一样。设用 表示新图，用 表示原图。第十七页，讲稿共八十七

5、页哦在新图的坐标系中是左上角的象素，原图中的象素在新坐标系中的位置变成了f（tx,-ty)，新的坐标系中g(0,0)f(-tx,-ty)。例如：第十八页，讲稿共八十七页哦平移后出现的问题1、如果新图中有一点（x1,y1），按照公式(2)得到的(x0,y0)不在原图中该怎么办？通常的做法是，把该点的RGB值统一设成(0,0,0)或者(255,255,255)。2、平移后的图象是否要放大？一种做法是不放大，移出的部分被截断，如下图所示，图2为原图，图3为移动后的图。这种处理，文件大小不会改变。第十九页，讲稿共八十七页哦平移前的图第二十页，讲稿共八十七页哦平移后的图像第二十一页，讲稿共八十七页哦还有

6、一种做法是：将图象放大，使得能够显示下所有部分 移动后被放大图像第二十二页，讲稿共八十七页哦这种处理，文件大小要改变。设原图的宽和高分别是w1,h1则新图的宽和高变为w1+|tx|和h1+|ty|，加绝对值符号是因为tx,ty有可能为负.下面给出Translation的源代码。算法的思想是先将所有区域填成白色，然后找平移后显示区域的左上角点(x0,y0)和右下角点(x1,y1)。分几种情况：第二十三页，讲稿共八十七页哦先看x方向(width指图象的宽度)l1.tx-width很显然，图象完全移出了屏幕，不用做任何处理2-widthtx0，如图5所示图像向左上方移动图象区域的x范围从0到widt

7、h-|tx|,对应原图的范围从|tx|到width第二十四页，讲稿共八十七页哦图5第二十五页，讲稿共八十七页哦30txwidth,如图6所示,图像向右下方平移0txwidth,0tyheight的情况图象区域的x范围从tx到width,对应原图的范围从0到width-tx 图6第二十六页，讲稿共八十七页哦4txwidth,图象完全移出了屏幕，不用做任何处理。y方向是对应的（方向是对应的（height表示图象的高度表示图象的高度1ty-height，图象完全移出了屏幕，不用做任何处理2-heightty0，图象区域的y范围从0到height-|ty|,对应原图的范围从|ty|到height30t

8、yheight，图象区域的y范围从ty到height,对应原图的范围从0到height-ty4tyheight，图象完全移出了屏幕，不用做任何处理）第二十七页，讲稿共八十七页哦 这种做法利用了位图存储的连续性，即同一行的像素在内存中是相邻的。利用memcpy函数，从(x0,y0)点开始，一次可以拷贝一整行（宽度为x1-x0），然后将内存指针移到(x0,y0+1)处，拷贝下一行，这样拷贝(y1-y0)行就完成了全部操作，避免了一个一个像素的计算，提高了效率。第二十八页，讲稿共八十七页哦二、旋转(rotation)图像旋转是指图像以某一点作为轴转过一定的角度。通常的做法是以图象的中心为圆心旋转。例

9、如：图7旋转30度（顺时针方向）后的图象如图8所示.可以看出，旋转后图象变大了。第二十九页，讲稿共八十七页哦图8旋转后的图图7原图第三十页，讲稿共八十七页哦若要使旋转后保持原尺寸图像就要进行裁切有些内容就会损失。我们采用第一种做法 第三十一页，讲稿共八十七页哦图像旋转变换矩阵 先来看一下，将一个点顺时针旋转a角后的坐标变换公式，如图10所示，r为该点到原点的距离，在旋转过程中，r保持不变。b为r与x轴之间的角度。x0=r*cos(b)；y0=r*sin(b)图10第三十二页，讲稿共八十七页哦旋转前：x0=r*cos(b)；y0=r*sin(b)旋转a角度后：第三十三页，讲稿共八十七页哦以矩阵的

10、形式表示 第三十四页，讲稿共八十七页哦其逆矩阵的形式表示为第三十五页，讲稿共八十七页哦上面的公式中，坐标系I是以图象的中心为原点，向右为x轴正方向，向上为y轴正方向。它和以图象左上角点，向右为x轴正方向，向下为y轴正方向的坐标系II之间的转换关系如何呢？两种坐标系间的转换关系 第三十六页，讲稿共八十七页哦逆变换为：第三十七页，讲稿共八十七页哦按照上面的公式，我们可以把变换分成3步：第一：将坐标系II变成I；第二：将该点顺时针旋转a角；第三：将坐标系I变回II，这样，我们就得到了变换矩阵，是上面三个矩阵的级联。第三十八页，讲稿共八十七页哦上面公式中Wold表示原图的宽,Hold表示原图 的高,W

11、new新图的宽,Hnew新图的高，可以看出 Wnew=max(|x4-x1|,|x3-x2|);Hnew=max(|y4-y1|,|y3-y2|).第三十九页，讲稿共八十七页哦可以根据以上公式，求出对应的原图中的点，得到它的灰度，如果超出原图范围则填成白色。要注意的是，由于有浮点运算，计算出来点的坐标可能不是整数，采用取整处理，即找最接近的点，这样会带来一些误差（图象可能会出现锯齿），更精确的方法是采用插值。插值算法后面详细介绍。第四十页，讲稿共八十七页哦三、镜象(mirror)镜象分水平镜象和垂直镜象两种 原图第四十一页，讲稿共八十七页哦水平镜像第四十二页，讲稿共八十七页哦垂直镜像第四十三页

12、，讲稿共八十七页哦镜象的变换矩阵很简单，设原图宽为w，高为h，变换后，图的宽和高不变.水平镜象为：垂直镜象为 第四十四页，讲稿共八十七页哦四、转置(transpose)转置是指将x,y坐标对换 原图转制后的图第四十五页，讲稿共八十七页哦要注意的是，转置和旋转90度是有区别的，不信你可以试试，你怎么转，也转不出图15度来。另外，转置后图的宽高对换了。转置的变换矩阵很简单:第四十六页，讲稿共八十七页哦五、缩放（zoom）图像比例缩放是指将给定的图像在x轴方向按比例缩放fx倍，沿y轴方向比例缩放fy。当fx fy时，称为图像的全比例缩放。当fx fy时，图像象素间的相对位置会发生畸变。设原图像中的点

13、P(x0,，y0),比例缩放后在新的图像中对应点P(x,y)，P(x0,，y0)和 P(x,y)之间的对应关系如图所示：第四十七页，讲稿共八十七页哦第四十八页，讲稿共八十七页哦比例缩放前后两点之间的关系用矩阵表示为：第四十九页，讲稿共八十七页哦其逆矩阵为：第五十页，讲稿共八十七页哦插值算法（Interpolation）比例缩放所产生的图像中有可能在原图中找不到相应的象素点,我们采用的做法是找与之最临近的点，我们采用的做法是找与之最临近的点。即利用邻域的像素来估计新的像素值 的处理方法就是插值算法。常用的插值算法有：最邻近插值（Nearest Neighbour Interpolation）。线

14、形插值(Linear Interpolation)。双线形插值（Bilinear Intrepolation）曲线插值（Curvilinear Interpolation）第五十一页，讲稿共八十七页哦（1）最临近插值图像变换后找不到相应的象素点时，直接赋值为和它最相近的的象素值的插值算法。例如：在图像的缩放中，将一幅66的图像按比例缩小，fxfy1/2,图像的长和宽均缩小到原来的一半，成为一幅33的图像。原图像的象素从（0，0）到（5，5），新图像的象素从（0，0）到（2，2）。如下图：第五十二页，讲稿共八十七页哦66的原图f新图中的灰度值利用缩放公式计算得到第五十三页，讲稿共八十七页哦X0=

15、0/0.5=0 y0=0 X1=1/0.5=2 y1=2 X2=2/0.5=4 y2=4 g(0,0)=f(0,0),g(0,1)=f(0,2),g(0,2)=f(0,4),g(1,0)=f(2,0),g(1,1)=f(2,2),g(1,2)=f(2,4)g(2,0)=f(4,0),g(2,1)=f(4,2),g(2,2)=f4,4)33的新图g66的原图f第五十四页，讲稿共八十七页哦图像缩小，可以在原图像中挑选信息，缩小之后的图与原图像比较丢失了很多信息。图像放大处理中，需要对放大后多出来的空格填入适当的象素值，采用最临近插值法举例如下：44的以原图f第五十五页，讲稿共八十七页哦长和宽各放大

16、2倍，利用公式得到88的新图像。新图像中各象素位置X0=0 Y0=0 X1=0.5 Y1=0.5 X2=1 Y2=1 X31.5 Y3=1.5 X4=2 Y4=2 X5=2.5 Y5=2.5 X6=3 Y6=3 X7=3.5 Y7=3.5g(0,0)=f(0,0)g(0,1)=f(0,0.5)g(0,2)=f(0,1)g(0,3)=f(0,1.5)g(0,4)=f(0,2)g(0,5)=f(0,2.5)g(0,6)=f(0,3)g(0,7)=f(0,3.5)g(1,0)=f(0.5,0)g(1,1)=f(0.5,0.5)g(1,2)=f(0.5,1)g(1,3)f(0.5,1.5)g(1,4)

17、=f(0.5,2)g(1,5)=f(0.5,2.5)g(1,6)=f(0.5，3）g(1,7)=f(0.5,3.5)g(2,0)=f(1,0)g(2,1)=f(1,0.5)第五十六页，讲稿共八十七页哦g(2,2)=f(1,1）g(2,3)=f(1,1.5)g(2,4)=f(1,2)g(2,5)=f(1,2.5)g(2,6)=f(2,3）g(2,7)=f(1,3.5)g(3,0)=f(1.5,0)g(3,1)=f(1.5,0.5)g(3,2)=f(f(1.5,1)g(3,3)=f(1.5,1.5)g(3,4)=f(1.5,2)g(3,5)=f(1.5,2.5）g(3,6)=f(1.5,3)g(3

18、,7)=f(1.5,3.5)g(4,0)=f(2,0)g(4,1)f(2,0.5)g(4,2)=f(2,1)g(4,3)=f(2,1.5）g(4,4)=f(2,2)g(4,5)=f(2,2.5)g(4,6)=f(2,3)g(4,7)=f(2,3.5)g(5,0)=f(2.5,0)g(5,1)f(2.5,0.5)g(5,2)=f（2.5,1)g(5,3)=f(2.5,1.5）g(5,4)=f(2.5,2)g(5,5)=f(2.5,2.5)g(5,6)=f(2.5,3)g(5,7)=f(2.5,3.5)g(6,0)=f(3,0)g(6,1)f(3,0.5)g(6,2)=f（3,1)g(6,3)=f

19、(3,1.5）g(6,4)=f(3,2)g(6,5)=f(3,2.5)g(6,6)=f(3,3)g(6,7)=f(3,3.5)g(7,0)=f(3.5,0)g(7,1)f(3.5,0.5)g(7,2)=f(3.5,1)g(7,3)=f(3.5,1.5）g(7,4)=f(3.5,2)g(7,5)=f(3.5,2.5)g(7,6)=f(3.5,3)g(7,7)=f(3.5,3.5)第五十七页，讲稿共八十七页哦放大后88的新图第五十八页，讲稿共八十七页哦(2)线形插值线形插值是使用原图中两个值来构造所求坐标处的值。例如：如图所示，如果已经知道了两点x0,x2处的函数值f(x0),f(x2)，现在要求

20、x1处的函数值f(x1)。我们假设函数是线性的，利用几何知识可以知道:第五十九页，讲稿共八十七页哦第六十页，讲稿共八十七页哦双线性插值在图象处理中需要将线形插值扩展到二维的情况，即采用双线形插值（Bilinear Intrepolation）例如：已知四点 f(x0,y0),f(x2,y2)，f(x2,y0),f(x2,y2)的灰度，要求g（x1,y1)点的灰度，可以先在水平方向上由f(x0,y0),f(x2,y2)线形插值求出g(x0,y1)，f(x2,y0),f(x2,y2)线形插值求出g(x2,y1)，然后在垂直方向上由g(x0,y1),g(x2,y1)线形插值求出g(x1,y1)第六十

21、一页，讲稿共八十七页哦第六十二页，讲稿共八十七页哦双线性插值可以通过三次运算的到。处理效果比较满意。这种处理具有低通滤波性质，使高频分量受到损失，图像轮廓模糊。另外，图像处理中还使用双曲线插值，及三次内插法解决图像缩放及图像旋转后的灰度填充。第六十三页，讲稿共八十七页哦图像几何变换综合训练1、已知33的图像如图所示A、将此图以图形中心为轴顺钟向旋转30度B、将此图放大2倍C、对此图作水平镜像第六十四页，讲稿共八十七页哦1、首先将此图平移至以图像中心为原点的坐标系中第六十五页，讲稿共八十七页哦新坐标系中任一点与原坐标系的点的对应关系如下X1=X-2Y1=Y-2由新图像找原图像中的对应点用拟矩阵第

22、六十六页，讲稿共八十七页哦新图像用g(u,y),表示，原图像用f(x,y),表示有：u0=x2 v0=y2 u(-1）=x1 v(-1)=y1 u1=x3 v1=y3g(-1,-1)=f(1,1)g(-1,0)=f(1,2)g(-1,1)=f(1,3)g(0,-1)=f(2,1)g(0,0)=f(2,2)g(0,1)=f(2,3)g(1,-1)=f(3,1)g(1,0)=f(3,2)g(1,1)=f(3,3)平移后的图像：第六十七页，讲稿共八十七页哦第六十八页，讲稿共八十七页哦第二、将图像以坐标原点为轴顺钟向转过30度第六十九页，讲稿共八十七页哦计算旋转后的图像大小宽度:w=w0cos(a)+

23、H0sin(a)高度:H=H0cos(a)+w0sin(a)因此，旋转后的图像：第七十页，讲稿共八十七页哦旋转后的图像与未旋转图之间的对应关系如下：旋转后的图像坐标用x2,y2表示，未旋转的图像坐标用x1,y1表示。第七十一页，讲稿共八十七页哦第七十二页，讲稿共八十七页哦由新图像找原图像中的对应点用逆矩阵第七十三页，讲稿共八十七页哦新图像用 p(s,t,),表示，原图像用 g(u,v),表示有：S(-2)=u(-2.732)t(-2)=v(-2.732)s(-1)=u(-1.366)t(-1)=v(-1.366)S(0)=u(0)t(0)=v(0)s(1)=u(1.366)t(1)=v(1.3

24、66)s(2)=u(2.732)t(2)=v(2.732)p(-2,-2)=g(-2,-2)p(-2,-1)=g(-2,-1)p(-2,0)g(-2,1)p(-2,1)=g(-2,1)p(-2,2)=g(-2,2)p(-1,-2)=g(-1,-2)p(-1,-1)=g(-1,-1)p(-1,0)g(-1,0)p(-1,1)=g(-1,1)p(-1,2)=g(-1,2)第七十四页，讲稿共八十七页哦P(0,-2)=g(0,-2)p(0,-1)=g(0,-1)p(0,0)g(0,0)p(0,1)=g(0,1)p(0,2)=g(0,2)p(1,-2)=g(1,-2)p(1,-1)=g(1,-1)p(1

25、,0)g(1,1)p(1,1)=g(1,1)p(1,2)=g(1,2)P(2,-2)=g(0,-2)p(2,-1)=g(0,-1)p(2,0)g(0,0)p(2,1)=g(0,1)p(2,2)=g(0,2)经过旋转的图像如下第七十五页，讲稿共八十七页哦用白色填充找不到的象素点得道的图像A，用深色填充找不到的象素点得道的图像BAB第七十六页，讲稿共八十七页哦2、将此图放大2倍M*N=6*6 用g(u,v)表示放大后的图像，f(x,Y)表示原图 u(0)=X(0)v(0)=y(0）u(1)=x(0.5)v(1)=y(0.5)u(2)=x(1)v(2)=y(1)u(3)X(1.5)v3=Y(1.5)

26、u(4)=X(2)v(4)=Y(2)u(5)=X(2.5)v(5)=Y(2.5)第七十七页，讲稿共八十七页哦g(0,0)=f(0,0)g(0,1)=f(0,0.5)g(0,2)=f(0,1)g(0,3)=f(0,1.5)g(0,4)=f(0,2)g(0,5)=f(0,2.5)g(1,0)=f(0.5,0)g(1,1)=f(0.5,0.5)g(1,2)=f(0.5,1)g(1,3)f(0.5,1.5)g(1,4)=f(0.5,2)g(1,5)=f(0.5,2.5)g(2,0)=f(1,0)g(2,1)=f(1,0.5)g(2,2)=f(1,1）g(2,3)=f(1,1.5)g(2,4)=f(1,

27、2)g(2,5)=f(1,2.5)g(3,0)=f(1.5,0)g(3,1)=f(1.5,0.5)g(3,2)=f(f(1.5,1)g(3,3)=f(1.5,1.5)g(3,4)=f(1.5,2)g(3,5)=f(1.5,2.5)g(4,0)=f(2,0)g(4,1)f(2,0.5)g(4,2)=f(2,1)g(4,3)=f(2,1.5）g(4,4)=f(2,2)g(4,5)=f(2,2.5)g(5,0)=f(2.5,0)g(5,1)f(2.5,0.5)g(5,2)=f(2.5,1)g(5,3)=f(2.5,1.5）g(5,4)=f(2.5,2)g(5,5)=f(2.5,2.5)第七十八页，讲

28、稿共八十七页哦第七十九页，讲稿共八十七页哦3、对此图作水平镜像在坐标系中以图像宽度为轴对图像作水平镜像。设原图像中的点为P0(x0,y0),镜像后的点为p(x,y),图像宽度为w。则有：X=w-x0Y=y0矩阵形式第八十页，讲稿共八十七页哦第八十一页，讲稿共八十七页哦2、已知边长为1cm的44的图像f（x，y）如下：A、将它沿p（x，Y）（60，120）与x轴成30度的直线上放大120，写出几何变换公式第八十二页，讲稿共八十七页哦B、分别用最临近插值和双线性插值求f(2.3,3.2)的值A、首先对其平移，平移后的图像中心为p(60,120).平移后的图像用f(x1,y1)表示，原图像用g(x,y).X1=X-60Y1=Y-120矩阵为第八十三页，讲稿共八十七页哦逆运算为：逆矩阵为：第八十四页，讲稿共八十七页哦第二、将图像以坐标原点为轴顺钟向转过30度宽度:w=w0cos(a)+sin(a)高度:H=w0cos(a)+sin(a)因此，旋转后的图像：第八十五页，讲稿共八十七页哦第八十六页，讲稿共八十七页哦感感谢谢大大家家观观看看第八十七页，讲稿共八十七页哦