第8章 幾何運算

1、DIP 用于改變圖像中各物體之間的空間關(guān)用于改變圖像中各物體之間的空間關(guān) 系。系。相當于物體在圖像內(nèi)移動相當于物體在圖像內(nèi)移動。 空間變換要有約束，否則變換后的圖像會空間變換要有約束，否則變換后的圖像會“支支 離破碎離破碎”。 定義：定義： 其中：其中： f(x, y)為輸入圖像，為輸入圖像， g(x, y)為輸出圖像。函數(shù)為輸出圖像。函數(shù) a(x, y)和和b(x, y)唯一描述了空間變換唯一描述了空間變換，若它們是連續(xù)，若它們是連續(xù) 的，則連通關(guān)系將在圖像中得到保持。的，則連通關(guān)系將在圖像中得到保持。 DIP f(x, y)的灰度值僅在整數(shù)位置的灰度值僅在整數(shù)位置(x, y)被定義被定義（

2、原始原始 圖像各像素是位于各離散的整數(shù)坐標點上圖像各像素是位于各離散的整數(shù)坐標點上），但變換），但變換 后的后的。 這就意味著，這就意味著，g(x, y)的灰度值通常會由處于非整的灰度值通常會由處于非整 數(shù)坐標點上的數(shù)坐標點上的f(x, y)的值來決定。的值來決定。 因此，如果把因此，如果把f映射為映射為g，則則f中的一個像素可能中的一個像素可能 會映射到會映射到g中幾個像素之間的位置；反之亦然。中幾個像素之間的位置；反之亦然。 所有像素都正好位于采樣?xùn)鸥竦恼麛?shù)坐標所有像素都正好位于采樣?xùn)鸥竦恼麛?shù)坐標 處處。 DIP 當實現(xiàn)一個幾何運算時，可以想象當實現(xiàn)一個幾何運算時，可以想象將輸入圖像的將輸

3、入圖像的 灰度一個一個像素地轉(zhuǎn)移到輸出圖像中灰度一個一個像素地轉(zhuǎn)移到輸出圖像中，若一個像素，若一個像素 被映射到被映射到4 4個輸出像素之間的位置，則其灰度值就要個輸出像素之間的位置，則其灰度值就要 按照插值算法，在按照插值算法，在4 4個輸出像素之間進行分配。這便個輸出像素之間進行分配。這便 是是。 x y y x f(x, y)g(x, y) x, y非整形非整形x, y整形整形 DIP 還可以有另一種與上述相反的映射過程，即：還可以有另一種與上述相反的映射過程，即：將將 輸出圖像的灰度一個一個像素地反向映射到輸入圖像輸出圖像的灰度一個一個像素地反向映射到輸入圖像 中中，因此一個輸出像素可

4、能會被映射到，因此一個輸出像素可能會被映射到4 4個輸入像素個輸入像素 之間的位置，輸出像素的灰度值就要按照插值算法，之間的位置，輸出像素的灰度值就要按照插值算法， 由由4 4個輸入像素的某種組合來決定。這便是個輸入像素的某種組合來決定。這便是 。 x y y x f(x, y)g(x, y) x, y整形整形x, y非整形非整形 DIP 輸出像素被映射到輸入圖像中的非整數(shù)位置，輸出像素被映射到輸入圖像中的非整數(shù)位置， 即即4 4個像素之間，因此需要插值運算。個像素之間，因此需要插值運算。 輸出像素的灰度值等于離它所映射到的位置最輸出像素的灰度值等于離它所映射到的位置最 近的輸入像素的灰度值近

5、的輸入像素的灰度值。 簡單，但精度不高。簡單，但精度不高。 DIP 函數(shù)函數(shù)f(x, y)在單位正方形頂點的值已知，希望通在單位正方形頂點的值已知，希望通 過插值找到正方形內(nèi)任意點的值過插值找到正方形內(nèi)任意點的值。 顯然，令由雙線性方程顯然，令由雙線性方程來來 定義的一個定義的一個與已知灰度值的與已知灰度值的4個頂點擬合，個頂點擬合， 即可根據(jù)這即可根據(jù)這4個頂點的灰度值求出個頂點的灰度值求出a、b、c、d 這這4個系個系 數(shù)。從而得到用于計算灰度級插值的雙線性方程。數(shù)。從而得到用于計算灰度級插值的雙線性方程。 f(0) f(1) 01 DIP 也可將上三式合并：也可將上三式合并： 由上可見：

6、前由上可見：前3個式子共個式子共3次乘、次乘、6次加次加，而第，而第4個式子個式子 則有則有4次乘、次乘、8次加次加。因此，用前。因此，用前3個式子分步計算要比只個式子分步計算要比只 用第用第4個式子一次計算還要快些。個式子一次計算還要快些。 DIP (x,1) (x,y) (0,y) (x,0) (1,1) (0,1) (0,0) (1,0) x y x y 由圖可見，通過雙線性插值得到的由圖可見，通過雙線性插值得到的與與4個個 頂點頂點通常是不一樣的。通常是不一樣的。 DIP 放大放大263.9% 最臨近插值最臨近插值雙線性插值雙線性插值 DIP 雙線性插雙線性插值具有平滑作用，會使圖像的

7、值具有平滑作用，會使圖像的細節(jié)退化細節(jié)退化， 特別是在放大處理時，這種影響將更為明顯特別是在放大處理時，這種影響將更為明顯。 另外，由雙線性插值得到的雙曲拋物面在鄰域的另外，由雙線性插值得到的雙曲拋物面在鄰域的 邊界處相接吻合，但該處的斜率并不吻合，即邊界處相接吻合，但該處的斜率并不吻合，即導(dǎo)數(shù)在導(dǎo)數(shù)在 邊界處不連續(xù)邊界處不連續(xù)。 因此，高階插值相當于因此，高階插值相當于用一張大的蒙皮恰好蓋住用一張大的蒙皮恰好蓋住 各個灰度樁值各個灰度樁值，而不是象雙線性插值那樣，僅是對每，而不是象雙線性插值那樣，僅是對每 相鄰的相鄰的4 4個樁蓋一張蒙皮。個樁蓋一張蒙皮。 灰度樁值灰度樁值 DIP 根據(jù)前面

8、的公式：根據(jù)前面的公式： 令令，則則 令令，則則 此時，圖像相對于基準點被平移了此時，圖像相對于基準點被平移了 用齊次矩陣表示為用齊次矩陣表示為 22 00 xy 0 0 ( ,)10 ( ,)01 10011 a x yxx b x yyy DIP 令令，則圖像在則圖像在x軸向和軸向和y軸向分軸向分 別放大別放大c和和d倍。（倍。（左上角的圖像原點在圖像縮放時不動左上角的圖像原點在圖像縮放時不動） 令上述令上述，則則，產(chǎn)生關(guān)于產(chǎn)生關(guān)于y軸的鏡像。軸的鏡像。 令上述令上述，則則，產(chǎn)生關(guān)于產(chǎn)生關(guān)于x軸的鏡像。軸的鏡像。 ( ,)1/00 ( ,)01/0 10011 a x ycx b x yd

9、y DIP 令令， 則圖像繞原點順時針轉(zhuǎn)則圖像繞原點順時針轉(zhuǎn) 角。角。 ( ,)cos( )sin( )0 ( ,)sin( )cos( )0 10011 a x yx b x yy 當圖像繞任意點當圖像繞任意點旋轉(zhuǎn)時，可以先將圖像平移旋轉(zhuǎn)時，可以先將圖像平移 到原點，經(jīng)旋轉(zhuǎn)變換后，再移回到到原點，經(jīng)旋轉(zhuǎn)變換后，再移回到位置。即：位置。即： 00 00 ( , )1 0cos( )sin( )01 0 ( , )0 1sin( )cos( )00 1 10 010010 011 a x yxxx b x yyyy DIP 如果圖像僅作如果圖像僅作平移平移或或縮放縮放操作，輸出像素的坐標位操作，

10、輸出像素的坐標位 置就僅置就僅分別依賴于分別依賴于x和和y。可以。可以在這兩個方向分別進行變在這兩個方向分別進行變 換換，效率還可能更高。例如，先在水平方向產(chǎn)生一個中，效率還可能更高。例如，先在水平方向產(chǎn)生一個中 間圖像，然后把間圖像，然后把中間圖像中間圖像作為輸入，再在垂直方向進行作為輸入，再在垂直方向進行 變換，得到最終結(jié)果。變換，得到最終結(jié)果。 等寬等寬 DIP 對于對于旋轉(zhuǎn)變換旋轉(zhuǎn)變換，Catmull和和Smith提出了兩步法：提出了兩步法： 即根據(jù)即根據(jù)導(dǎo)出：導(dǎo)出： 再代入到再代入到，得：，得： 因此，在運算的第一步（只在水平方向）先把因此，在運算的第一步（只在水平方向）先把 與與聯(lián)

11、立（聯(lián)立（即變換圖像即變換圖像 ，其中，其中 仍為原來的值仍為原來的值）。）。 然后再在運算的第二步（只在垂直方向）將然后再在運算的第二步（只在垂直方向）將 式與式與聯(lián)立（聯(lián)立（即即，其中其中 剛剛變化完，不需再變剛剛變化完，不需再變）。）。 ( ,)sin() cos() a x yy x ( ,)sin( ) ( ,) cos( ) a x yy b x y DIP 在上述可分離旋轉(zhuǎn)在上述可分離旋轉(zhuǎn)變換中，要求變換中，要求 值要小，否則當值要小，否則當 時因時因而失靈。但通常情況的旋轉(zhuǎn)角度而失靈。但通常情況的旋轉(zhuǎn)角度 都很小，由于都很小，由于及及的整數(shù)倍旋轉(zhuǎn)可以通過簡單的的整數(shù)倍旋轉(zhuǎn)可以通

12、過簡單的 行列交換而實現(xiàn)，因此，旋轉(zhuǎn)角度可以限制在行列交換而實現(xiàn)，因此，旋轉(zhuǎn)角度可以限制在之之 間。間。 由以上可見，平移、放大、旋轉(zhuǎn)均可以分為兩步由以上可見，平移、放大、旋轉(zhuǎn)均可以分為兩步 的一維運算。的一維運算。 DIP 在許多圖像處理應(yīng)用中，所需的在許多圖像處理應(yīng)用中，所需的空間變換都相空間變換都相 當復(fù)雜當復(fù)雜，無法用前述無法用前述那樣那樣 的簡便數(shù)學(xué)表達式表示的簡便數(shù)學(xué)表達式表示。另外，。另外，所需的空間變換經(jīng)所需的空間變換經(jīng) 常要從對實際圖像的測量中獲得常要從對實際圖像的測量中獲得，因此更希望用這，因此更希望用這 些些測量結(jié)果測量結(jié)果而不是函數(shù)形式來而不是函數(shù)形式來描述幾何變換描述

13、幾何變換。 例如，攝像畫面上部的畸變，需用矩形柵格目例如，攝像畫面上部的畸變，需用矩形柵格目 標對其數(shù)字化，再對不同的區(qū)域用相應(yīng)的空間變換標對其數(shù)字化，再對不同的區(qū)域用相應(yīng)的空間變換 來校正。對同一部有畸變失真的攝像機來說，上述來校正。對同一部有畸變失真的攝像機來說，上述 空間變換可用于對所有從該攝像機輸出的圖像的校空間變換可用于對所有從該攝像機輸出的圖像的校 正。正。 DIP 指定圖像的一系列指定圖像的一系列控制點作空間變換，控制點作空間變換， 其他非控制點的位移由插值來決定其他非控制點的位移由插值來決定。 將圖像分成許多多邊形區(qū)域，并對每個將圖像分成許多多邊形區(qū)域，并對每個 區(qū)域使用雙線性

14、映射。區(qū)域使用雙線性映射。 DIP 找到變換式找到變換式中中 和和的函數(shù)表達式，該表達式通常為的函數(shù)表達式，該表達式通常為多項多項 式式，其參數(shù)選擇能使多項式與控制點及其位移量相，其參數(shù)選擇能使多項式與控制點及其位移量相 吻合。（吻合。（參考文獻參考文獻1中使用了高達中使用了高達5階的多項式卷階的多項式卷 繞函數(shù)繞函數(shù)）。）。 如果控制點數(shù)多于多項式項數(shù)，就必須采用多如果控制點數(shù)多于多項式項數(shù)，就必須采用多 項式擬合來決定多項式系數(shù)。這時要求項式擬合來決定多項式系數(shù)。這時要求空間變換整空間變換整 體上與指定的控制點最佳擬合體上與指定的控制點最佳擬合。但對于每一個控制。但對于每一個控制 點則不一

15、定嚴格吻合。點則不一定嚴格吻合。 DIP DIP 輸入的輸入的控制點形成柵格控制點形成柵格，此柵格，此柵格映射映射為輸出圖像為輸出圖像 中中的的矩形柵格矩形柵格。輸入控制點對應(yīng)矩。輸入控制點對應(yīng)矩 形頂點，輸入多邊形內(nèi)的點則對應(yīng)矩形內(nèi)各點。形頂點，輸入多邊形內(nèi)的點則對應(yīng)矩形內(nèi)各點。 DIP 通常采用通常采用雙線性插值雙線性插值對柵格進行插值，這是因為雙線對柵格進行插值，這是因為雙線 性插值計算簡單，且可以產(chǎn)生一個能保持連續(xù)性和連通性性插值計算簡單，且可以產(chǎn)生一個能保持連續(xù)性和連通性 的光滑映射的光滑映射。 雙線性空間變換的一般表達式為雙線性空間變換的一般表達式為： 可見該變換由可見該變換由這這

16、8個系數(shù)來定義。個系數(shù)來定義。 輸入四邊形的輸入四邊形的4個頂點個頂點輸出四邊形的輸出四邊形的4個頂點個頂點 這樣，從這樣，從到到 的映射，得到的映射，得到4個含個含的方程的方程 從從 到到 的映射，得到的映射，得到4個含個含的方程的方程 根據(jù)上述根據(jù)上述8個方程，即可解出個方程，即可解出這這8 個系數(shù)，個系數(shù)， 。 DIP 還有另一種更方便、高效的定義還有另一種更方便、高效的定義： 上式中的上式中的和和表示像素的位移量，是表示像素的位移量，是 和和 的雙線性函數(shù)。的雙線性函數(shù)。 xa,ya xb,yb xc,yc xd,yd x1,y1x2,y1 x2,y2x1,y2 x,y1 x,y2 x

17、,y dx(x2,y2) dy(x2,y2) 輸入四邊形 輸出四邊形 DIP 由于是雙線性，那么由于是雙線性，那么沿著每一輸出行（或每一沿著每一輸出行（或每一 輸出列）都是線性的輸出列）都是線性的。對每行定義增量。對每行定義增量 x（像素間像素間 隔為隔為1），因此對），因此對dx：，：， 同樣對同樣對dy: (1, )( , )dx xydx x yx ( ,1)( , )dy x ydy x yy 增量增量 x隨著行的不同而變化隨著行的不同而變化，但在輸出矩形頂，但在輸出矩形頂 點的位移值可容易地計算（對應(yīng)點約束），即點的位移值可容易地計算（對應(yīng)點約束），即通過通過 給定點的位移進行插值給

18、定點的位移進行插值來求得。來求得。 僅需兩個加法（僅需兩個加法（dx和和dy各一）。各一）。 DIP 1 11 1 4 44 456 23 233 3 6 665 上圖是定義了上圖是定義了6個相鄰矩形的輸出圖像。其空間個相鄰矩形的輸出圖像。其空間 變換如前所述。外面矩形的雙線性系數(shù)來自其最鄰近變換如前所述。外面矩形的雙線性系數(shù)來自其最鄰近 的內(nèi)部矩形。另外，相鄰矩形的頂點要重合。的內(nèi)部矩形。另外，相鄰矩形的頂點要重合。 通常一個由四邊形到矩形的映射不足以描述所期通常一個由四邊形到矩形的映射不足以描述所期 望的空間變換，這就需要望的空間變換，這就需要 DIP 消除由于攝像機導(dǎo)致的數(shù)字圖像的幾何畸

19、變。消除由于攝像機導(dǎo)致的數(shù)字圖像的幾何畸變。 校正從衛(wèi)星上或飛機側(cè)視雷達上得到的發(fā)生幾何畸校正從衛(wèi)星上或飛機側(cè)視雷達上得到的發(fā)生幾何畸 變的數(shù)字圖像。變的數(shù)字圖像。 DIP 一些圖像系統(tǒng)使用非矩形的一些圖像系統(tǒng)使用非矩形的 像素坐標，用普通顯示器觀察圖像素坐標，用普通顯示器觀察圖 像時需轉(zhuǎn)換為矩形像素坐標。像時需轉(zhuǎn)換為矩形像素坐標。 例：例：Viking Lander宇宙飛船宇宙飛船 使用使用角掃描攝像機角掃描攝像機來數(shù)字化火星來數(shù)字化火星 的全景圖像。采用球面坐標系統(tǒng)，的全景圖像。采用球面坐標系統(tǒng)， 掃描線之間以相等的俯仰角隔開，掃描線之間以相等的俯仰角隔開， 其像素間距代表相等的方位角增其

20、像素間距代表相等的方位角增 量。量。校正時，需要將球面投影到校正時，需要將球面投影到 一個相切的平面上，即投影從球一個相切的平面上，即投影從球 面中心點發(fā)出，并將球面上的點面中心點發(fā)出，并將球面上的點 投影到平面上。投影到平面上。 用矩形系統(tǒng)顯示角掃描圖像用矩形系統(tǒng)顯示角掃描圖像 DIP 對相似的圖像進行配準，以便進行比較。例如，對相似的圖像進行配準，以便進行比較。例如， 利用圖像相減來檢測運動或變化利用圖像相減來檢測運動或變化。如果相減的圖像有。如果相減的圖像有 微小的位移，則相減后的差圖像將會有明顯的偏導(dǎo)數(shù)微小的位移，則相減后的差圖像將會有明顯的偏導(dǎo)數(shù) 成分，這就掩蓋了真正有意義的圖像差別

21、。靜止圖像成分，這就掩蓋了真正有意義的圖像差別。靜止圖像 用固定位置的攝像機拍攝多幅，則各圖像是配準的。用固定位置的攝像機拍攝多幅，則各圖像是配準的。 否則，否則，在圖像相減前要進行配準在圖像相減前要進行配準。 簡單的配準用平移，而旋轉(zhuǎn)或更復(fù)雜的變形需要簡單的配準用平移，而旋轉(zhuǎn)或更復(fù)雜的變形需要 復(fù)雜的幾何運算。復(fù)雜的幾何運算。 對各圖像的細小特征點進行定位，可用于定義控對各圖像的細小特征點進行定位，可用于定義控 制點，而后基于控制點進行圖像的減運算。制點，而后基于控制點進行圖像的減運算。 DIP 將圖像以另一種形式進行表示，以便于解釋。將圖像以另一種形式進行表示，以便于解釋。 例如，通過幾何

22、運算將果蠅染色體的原始圖像拉例如，通過幾何運算將果蠅染色體的原始圖像拉 平（原始圖像由顯微鏡所拍攝的照片拼接而成），其平（原始圖像由顯微鏡所拍攝的照片拼接而成），其 原理是將輸入圖像的每個染色體疊加上一個四邊形的原理是將輸入圖像的每個染色體疊加上一個四邊形的 控制柵格，并使其中的兩條邊平行于染色體的軸線，控制柵格，并使其中的兩條邊平行于染色體的軸線， 因而這些四邊形被映射成輸出圖像中水平排列的矩形因而這些四邊形被映射成輸出圖像中水平排列的矩形 串。串。 DIP 例如，從宇例如，從宇 宙飛船上傳回的宙飛船上傳回的 地球、月球圖片地球、月球圖片 是球面圖片，而是球面圖片，而 拼接成地圖時是拼接成地

23、圖時是 平面的，因此需平面的，因此需 要幾何運算。這要幾何運算。這 是一個是一個“球面球面” “平面平面”的幾何的幾何 變換過程變換過程。 DIP 1、沿某些直線保持、沿某些直線保持不變，則地圖不變，則地圖等距等距 （equidistant）。 2、投影中區(qū)域的、投影中區(qū)域的保持不變，則地圖保持不變，則地圖等值等值 （equivalence）。 3、投影中、投影中保持不變（表面上直線相交的保持不變（表面上直線相交的 夾角與地圖上投影線之間的夾角相同），則地夾角與地圖上投影線之間的夾角相同），則地 圖圖保角保角（conformal）或或同形同形（orthomolphic）。 DIP 正交投影法正

24、交投影法 球極平面投影法球極平面投影法 Mercator投影法投影法 Lambert保角圓錐投影法保角圓錐投影法 各種方法的生成技術(shù)、地圖性質(zhì)等均不同。各種方法的生成技術(shù)、地圖性質(zhì)等均不同。 DIP 在影視作品中常見，表現(xiàn)為在影視作品中常見，表現(xiàn)為 。由于在變換過程中，物體中。由于在變換過程中，物體中 的特征點會從其起始位置平滑地移向終止位置，因此的特征點會從其起始位置平滑地移向終止位置，因此 可產(chǎn)生為生動的效果?？僧a(chǎn)生為生動的效果。 這這不同于影視制作技術(shù)中的不同于影視制作技術(shù)中的技術(shù)技術(shù)。 因為淡入淡出過程中，第一幅畫面漸漸隱去，隨后第因為淡入淡出過程中，第一幅畫面漸漸隱去，隨后第 二幅畫面才淡淡地顯現(xiàn)出來。兩幅畫面沒有融合的過二幅畫面才淡淡地顯現(xiàn)出來。兩幅畫面沒有融合的過 程，因而不會產(chǎn)生逼真的視覺轉(zhuǎn)換。程，因而不會產(chǎn)生逼真的視覺轉(zhuǎn)換。 DIP 通常所說的變形都是指由通常所說的變形都是指