TWI234640B - Contactless three dimensional measure method - Google Patents

Contactless three dimensional measure method Download PDF

Info

Publication number
TWI234640B
TWI234640B TW91125314A TW91125314A TWI234640B TW I234640 B TWI234640 B TW I234640B TW 91125314 A TW91125314 A TW 91125314A TW 91125314 A TW91125314 A TW 91125314A TW I234640 B TWI234640 B TW I234640B
Authority
TW
Taiwan
Prior art keywords
image
coordinates
dimensional space
space
measurement
Prior art date
Application number
TW91125314A
Other languages
English (en)
Inventor
Ming Jang
Gau-Huei Lin
Wei-Je Jang
Original Assignee
Ming Jang
Gau-Huei Lin
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ming Jang, Gau-Huei Lin filed Critical Ming Jang
Priority to TW91125314A priority Critical patent/TWI234640B/zh
Application granted granted Critical
Publication of TWI234640B publication Critical patent/TWI234640B/zh

Links

Description

1234640 案號91125314 年 月 曰 修正 五、發明說明(1) 本發明非接觸式三維空間量測方法,係指一種應用於 三維座標之量測技術,可由數位影像直接求得空間中待測 物體之三維座標。
快速精準的三次元量測,無論在產品設計、工業製造 、產品檢驗以及商業多媒體資訊方面都有極大的需求。一 般而言,物體三次元外型的量測技術上可分為接觸式與非 接觸式兩大類。接觸式量測主要是利用機械式探針,裝置 於三次元量床(CMM )上,藉由平台的移動裝置及探針接 觸物體表面的逐點方式,來量測工件物之外型,其缺點為 量測速度慢及容易對軟質的物體表面造成破壞,而非接觸 式量測的方法則可以彌補此項之缺憾。 在非接觸式三維量測技術中,光學量測技術是最重要 的方式之一。光學非接觸式量測技術方法可分為主動式及 被動式兩類,其主要的差異在被動式光學非接觸式量測不 提供特定光源,而主動式光學非接觸式量測則需要,立體 視覺法為被動式量測中最具代表的方式,其主要的關鍵在 於量測校正技術,良好的量測校正才能確保量測的精度。
立體視覺校正技術一般最常用的是三維空間透視投影 的針孔(p i η - h ο 1 e )模型,此p i η - h ο 1 e數學模型只利用單 純的幾何線性轉換關係,並沒有考慮鏡頭的非線性像差因 素,故只適合於長焦距鏡頭CCD攝影機模型或較粗略的量 測上。另一種考慮軸對稱像差的CCD攝影機模型乃由Tsai [R. Y. Tsai, ” A versatile camera calibration for high-accuracy 3D machine vision metrology using
第4頁 1234640 _案號91125314_年月日__ 五、發明說明(2) off-the-shelf TV cameras and lenses", IEEE Journal of Robotics and Automation, pp. 323-344, 1987 ; R. Y. Tsai, and R. K. Lenz, ” A new technique for fully autonomous and efficient 3D robotics hand/ey e calibration 丨丨,IEEE Trans. Robot.
Automat. , vol· 5,no. 3,pp. 3 45 - 3 5 8,Jun. 1 9 8 9. 】所提出,可求得C C D攝影機的特性參數如影像中心、旋 轉角等線性係數以及球面像差等非線性係數。其後陸續經 由M· Robert’s 及A.J· Naftel [Μ· Roberts and A. J. Naitel,M A genetic algorithm approach to camera calibration in 3D machine vision", IEE Colloquium on Genetic Algorithms in Image Processing and
Vision’, pp· 12/1-12/5, 1994.】修改Tsai 的模型,因 其參數計算仍太複雜,故不適合快速影像量測應用。 除了以上針對攝影機校正所提到的空間透視投影轉換 及軸對稱像差關係的轉換外,A · N i 1 s s ο η及P. Η ο 1 m b e r g 【A· Nilsson and P· Holmberg, "Combining a state 2D vision camera and a π ultrasonic range detector for 3D position estimation", IEEE Instrumentation and Measurement Technology Conference, no. 18-20, p p · 4 5 3 - 4 5 8, M a y 1 9 9 3 ·】提出二維像差影像與二維線性 轉換影像的多項式函數關係式,省略求得c c D攝影機的特 性參數’大量簡化了 T s a i模型參數的計算;其後μ · c h a n g 及W.C. Tai【Μ· Chang and W. C. Tai,π 360-deg
1234640 案號 91125314_± 修正 五、發明說明(3)
profile noncontact measurement using a neura network" , Optical Engineering, v o1 . 34, no. pp. 3572-3576, 1995 ; W. C. Tai and M. Chang,
Noncontact published in Optical Engineering, vo1. 3 5, η o . 9, 1 9 9 6】更直接提出二維平面空間與二維影像 座標的映射函數轉換,其函數轉換的過程中,並不需要考 慮C C D攝影機參數的特性,簡化了校正的程序。然而,此 二維 礎, 故在 方法 特徵 徵記 係式 映射 座標 空間 係式 ,不 升量 式三 影像座 並以機 二維量 有鑑於 ,藉由 記號與 號座標 ,任何 函數關 值。 由是, 量測方 的運算 僅有效 測的速 為達到 維空間 .....十山門比,丨不W里閃丞 構位移疊合之方式求得三維空間影像之座標, 測與二維校正方面亦有空間疊合之誤差。 Ϊ ^ ΐ,係提供一種非接觸式三維空間量測 1β π糸列已知座標的二維空間校正網格 其對映成像於雙CCD攝影機二維影像平^ ’即可建構出真實空間與雙CCD影像介 、 在CCD影像空間的 :j的關 係式的運算1到其在真實^間 J只際位置 t赉:之主要目的,即在於提供 法2雙CCD攝影機與三維空接觸式三維 、上找到其在真實空間中的實際、射函數關 产i ^ =尺寸量測時所須之龐大運曾置座標值 度與精度。 更可提 ί 3二的,本發明是這樣實現的:一 里測方法,係包含: 種非接觸
1234640 案號91125314 年 月 曰 修正 4 五、發明說明(4) 一影像之擷取,乃由至少雙部影像感測元件、雷射光度控 制器與三軸移動平台所構成,利用雷射光度控制器產生 一特徵記號於量測物上,並配合三軸移動平台與影像感 測元件,可獲得特徵記號之空間實際座標與影像平面座 標;
一三維空間座標與影像平面座標之運算,藉由一函數係數 值可進行二維平面座標與三維空間座標之轉換;與 一電腦輔助處理,負責處理影像資料與進行數值資料之運 算,故,先透過影像擷取系統獲得特徵記號之實際座標 與影像平面座標,可求出三維空間映射函數之函數係數 值,再利用該函數係數值與量測物之影像平面座標,即 可獲得量測物之三維空間座標。 為使 貴審查委員進一步了解本發明之結構特徵及功 效,茲藉由下述具體之實施例,並配合所附之圖式,對本 發明做一詳細之說明,說明如后: 量測系統校正原理
在影像量測系統中,「影像量測空間」是指在空間中 的物體位置座標,經CCD攝影機鏡頭投影在CCD攝影機影像 感測元件上,可清楚成像所構成的空間集合。而「三維空 間影像量測」是指在一影像量測空間中的位置點,藉由其 成像於C C D感測元件的影像位置運算得到其實際的空間座 標值。請參考第1 a圖所示,為本發明雙攝影機量測架構 之示意圖,在X座標1 0 0、Y座標2 00與Z座標3 0 0 所組成之三維空間1 〇 〇 〇中,設定一個虛擬矩形量測空
第7頁 1234640 案號9Π25314 年 月 曰 修正 五、發明說明(5) 間3 0 ,在前方放置CCD1攝影機1 〇以及CCD2攝影機2 0 ,使量測空間可完全投影至攝影機影像平面;如第1 b圖 所示,係為攝影機之影像示意圖,當進行三維空間影像量 測時,先將特徵記號4 0放置於該矩形量測空間3 0中, 藉由此特徵記號4 0對映在影像平面的位置,分別為: CCD1攝影機1 〇可得左側影像圖1 2 ,CCD2攝影機2 0可 得右侧影像圖2 2 ,藉此即可運算求出其相對於攝影機空 間之實際位置座標。
若要由影像的位置運算求其相對於攝影機空間中之位 置,量測系統需先經過校正的程序,主要是利用已知空間 中點座標與其在CCD攝影機之影像點座標的對映函數關係 ,來建立系統的校正參數,並據以為後續量測的參數。 以下就線性透視投影轉換的推演來探討高次多項式映 射函數在三維空間影像量測校正應用。 【雙攝影機線性透視投影轉換】 當空間點座標(X,Y,Z )投影至兩個影像平面時,設其 透視投影轉換為線性關係,其多項式可表示為: + α2 ^Υ-\-α3 *Ζ (la) y{ =bO +¾ ^ X + b2 ^Y + h3
x2=c[i+cl^ Χ + ο2^ΥΛ-ο3^Ζ (lb)
夕2 = + 半Z + *^3 * Z
第8頁 1234640 _案號 91125314_年月日_ 五、發明說明(6) 為式型 -\ej矩換轉 2 ^2 Γ . 1I-J ζ 7 r 0 12 3 ώ ώ ο ,1 2 3 ο λι 2 ά ώ ά ά ο fl 2 3 2 或 C I I 5-- 3 U 72 x2 rL II I I rL 中其
I-1 z r X 1X - L II I J B c 0 12 3 Ja Γα ώ Ja 0 12 3 έο άι ά2 03 A ¾¾ 若已知空間中至少4個不共平面的點座標以及所相對 應的4組投影影像座標,則可求出投影轉換矩陣[斗 m—v] (4) 若已知空間中多於4個不共平面點座標以及所相對應 的投影影像座標則可用最小平方誤差法來求平均解
第9頁 1234640 案號 91125314 年月曰 修正 五、發明說明(7) (5) [[ΒΥΐΒ^ΒηΛ (6) 另一方面亦可求出兩組影像點至空間座標的轉換[小 U] [C]—1 (7)【多攝影機線性透視投影轉換】當空間點座標(X,Υ,Ζ)投影至多個影像平面時,設其透視 投影轉換為線性,其多項式可表示為: (8a) (8b) (8c)
j^2 = cIq ^ Jf d 2^ y Z (8d)
— Sq + + ^ Z (8e) Λ = /0+Λ*π + /2*7+/3*Ζ (8f)
第10頁 1234640 « 91125314 , , a 修正 五、發明說明(8) 轉換矩陣型式為 k 3^2 λ] = [! γ ζ] 或 U'卜[S ][ C'] 其中
B ] = [l X T Z b - - 知 已 若 4 的 應 投 (¾ c1。 α{ bx c{ b2 C2 a3 b3 c3 (10) 9 /0Λ/2Λ 0 12 3 e"^e 0 12 3 Ju Ju Ja 對 相 所 及 以陣 標矩 座換 點轉 的影 面投 平出 共求 不可 _-^個則 /0/4/34 , 94 票 至彳 ^04^2^3座中、 ^0^1¾ 4 像 /yy io 61¾¾空影 α。αι A A -影 __- 4 Η於可 C 多剧 丨中, 間標 空座 知像 已 若影 影 投 的 應 對 相 所 丨及 1 以 (1標 座 點 面 平 1共 1 不 個
S 最 用 均 平 求 來 法 差 誤 方 平
第11頁 1234640 案號91125314 年 月 曰 修正 五、發明說明(9) [s]rU] = [s]r[s][c] (X 2) [C ] - Ib]t[b ]]_1[ b]t[a] (13) 以此類推亦可求出多組影像點至空間座標的轉換如下 U" ] = [s ][cff ] ][c]r=[s ][c][c]r (1 4)
[B]=[A][cV[[c][cY}~1 (1 5) 上述雙攝影機線性透視投影轉換與多攝影機線性透視 投影轉換之方程式乃習知三維空間量測技術之原理,然而 其藉由線性推導的方式並無法解決攝影鏡頭的像差以及其 他因素造成投影關係為非線性之情況。 【雙攝影機高次多項式三維空間映射函數】 當考慮攝影鏡頭的像差以及其他因素造成投影關係為 非線性時,本發明提出以高次多項式三維空間映射函數轉 換來解決此問題。
三維空間影像量測與攝影機影像點座標的對映關係假 設以高次多項式三維空間映射函數轉換來做分析時,其轉 換與逆轉換之間並不像線性轉換般容易計算,所以分別用 高次多項式最小平方誤差擬合的方式來進行。 當由空間中的點座標轉換至影像座標時,以下列的函
第12頁 1234640 _案號 91125314_年月日_ifi 五、發明說明(10) 數來做轉換 n n-k n-k-j (16a) λ=0 ;=0 Ϊ=0 η n-k n-k-j y^G^XJ^Z)^ Σ ΣΒΦχ1γ]Ζ" λ=0 j=0 i=0 η n-k n-k-j x2=f2(x?f?z) = 2 Σ YJcl]kxiY^k (16b) λ=0 7=0 z=0 n n-k n-k-j y^G^XJ.Z)^ Σ
λ=0 7=0 ϊ==0 其中(X,Y,Z )為空間點座標,(a,y〇 為攝影機CCD1的影像座 標,(W2) 為攝影機C C D 2的影像座標,及為空間中的點 座標轉換至攝影機CCD 1影像座標的映射函數係數,及 為空間中的點座標轉換至攝影機CCD2影像座標的映射 函數係數,η為多項式函數的階數。 若由雙攝影機影像座標轉換至空間中的點座標時,以 下列的函數來做轉換
第13頁 1234640 案號91125314 年 曰 修正 五、發明說明(11) n n-l n-l-k n-l—k-j
^ = /(^1^1^2^2) = 2 Σ Σ Σ aijkixi y/x2kyJ /=0 λ=0 ;=0 Ϊ=0 n n-l n-l-k n-l-k-j ^ = g = Σ Σ Σ Σ bijkixiy\ x2kyi /=0 λ=0 >0 Ϊ=0 (17) n n-l n-l-k n-l-k-j z = /2 (xpj1?x2^2) = 2 Σ Σ Σ cijkixi yi x2kyi /=0 k=[\ j=0 i= 0 其中^^及·為影像座標轉換至空間中的點座標 的映射函數係數。 【多攝影機高次多項式三維空間映射函數】 當由空間中的點座標轉換至多攝影機影像座標時,以 下列的函數來做轉換 n n-k n-k-j x^Fx{XJ.Z)^ Σ Σ4^Ζ" λ=0 7=0 ϊ=0 η n-k n-k-j y^G.iXJ^Z)^ Σ ΣΒ^ΧΎ;Ζ" λ=0 ;=0 ϊ=0 η rt-k n-k-j x2^F2{XJ^)^ Σ ZC^T;Z" (18) λ=0 彳=0 ϊ=0 η n-k n-k-j ^2 = σ2(χ5Γ5ζ) = Σ Σ Σ°^χίγ^ λ=0 j=0 ΐ=0
第14頁 1234640 案號91125314 年 月 曰 修正 五、發明說明(12) n n-k n-k-j x^Fn{xj,z)-z Σ ΣΝΦχίγ'ζ" 允=〇 J=0 ϊ=0 n n-k n-k-j λ = Α(χ,7,ζ) = Σ Σ λ=0 7=0 Σ=0 其中(Χ«,Λ)為攝影機CCDn的影像座標,及ΜΪ^為空 間中的點座標轉換至攝影機CCDn影像座標的映射函數係 數。
若由多攝影機影像座標轉換至空間中的點座標時,以 下列的函數來做轉換 (19) X = f(xl,yl,x2,y2,...jc^yJ Z-h {x^y^x2,y2,...xn,yn) 【空間映射函數影像量測解析度】
三維空間投影在攝影機影像平面時,其投影影像的大 小、位置會因不同的光學設計而不同,所以必須先進行影 像量測解析度分析,這也是光學設計時的重要考量因素之 一。在本發明中三維空間影像量測解析度的定義為三維空 間投影在單位影像像素點的距離,單位為m m / p i X e 1。
第15頁 案號 91125314 1234640 修正 五、發明說明(13) 影像量測解析度又可細分 空間解析度,其計算的方式如;·,,各軸向解析度以及 X軸向解析度的計算係先在十 定點(X + dX,Y,Z)及(X〜dX γ 疋義的量測空間中取特 換至影像座標的映射函數,’可^將,7人空間中點座標轉 座標值,則X轴向解析度為 f X ’ y )及(X 2,y 2 )的影像 X resolution 二^/ U:, )2 + Ui)2 (2 0) 同理’將(X,Y + dY,Z)及(X,Y-dY,z)代入空間中點座 至影像座標的映射函數,可得(X 3,y 3 )及(X 4,y 4 )的偾: 標值’將(X,Y,Z + dZ)及(X,Y,Z-dZ)代入空間中點座f = f 至影像座標的映射函數,可得(X5,y5)及(X6,y6)的f 標值,則Y軸向解析度及Z軸向解析度分別為 〜m度 V(X4 -¾)2 + (凡-少3 )2 (21) y resohdtion 二 2 * 沿] Z resolution = yl(x6-x5f +(y6-y5f (22)
第16頁 1234640 案號91125314 年 月 曰 修正 五、發明說明(14) 而空間影像量測解析度則定為 resolution二 4X —resolution2 +Y —resolution2 +Z — resolution2 因為空間影像量測是透視投影的關係,其解析度不是 均一的,所以要將量測空間中的位置點逐一取點來分析。 三維空間影像量測實例
在本發明中,用雙攝影機來進行量測及三轴移動平台 來建立三維影像量測空間,並使用雙線雷射來形成十字特 徵記號,操作流程參考第2圖所示,係為非接觸式三維空 間量測流程示意圖,涵蓋了: 一影像之擷取,乃由至少雙部影像感測元件、雷射光度控 制器與三軸移動平台所構成,利用雷射光度控制器產生 一特徵記號於量測物上,並配合三軸移動平台與影像感 測元件,可獲得特徵記號之空間實際座標與影像平面座 標;
一三維空間座標與影像平面座標之運算,藉由一函數係數 值可進行二維平面座標與三維空間座標之轉換;與 一電腦輔助處理,負責處理影像資料與進行數值資料之運 算,故,先透過影像擷取系統獲得特徵記號之實際座標 與影像平面座標,可求出三維空間映射函數之函數係數 值,再利用該函數係數值與量測物之影像平面座標,即
第17頁 1234640 案號 91125314 曰 修正 五、發明說明(15) 可獲得量測物之三維空間座標。 首先形成特徵記號5 0 0後,分別經由三轴移動平台 之定位5 1 0與雙CCD攝影機之取像5 2 0,可獲得多組^ 三維空間座標(X,Y,Z ) ’與對應之雙c C D影像平面^座標 (α,λ)、(〜λ),代入三維空間映射函數校正5 3 0 ( 1 7 ) 式,可求得該式中的各項函數係數值5 4 〇 ,完成校正雙 攝影機高次多項式三維空間映射函數5 5 0 ,此時,(1 7 )式中的各項函數係數值則為已知。量測物6 Q Q則是 經雙CCD攝影機之取像5 2 0取得雙CCD影像平面標(Wl)、 〇2,乃)後’代入雙攝影南次多項式三維空間映射函 數5 5 0 ,即可求得三維空間座標(X,γ,z ) 6 1 0 。 【雙攝影機量測設置】 本發明中所使用的影像感測元件係為電荷耦合器 (Charge-Coup 1ed Device, CCD ) 〇 量測空間範圍通常與 視角與景深有關,當焦距 變形量也越大,當光圈越 ’所以要根據所需要的量 鏡頭,如第3圖所示,使 影機CCD1 1 〇 與CCD2 2 0 雙攝影機架設距離a為2 4 0 使用的鏡頭為AVENIR公司 光圈最大為F1. 2,水平視 攝影機鏡頭焦距及光圈所形成的 越短則視角越大,但相對的影像 小則景深越深,但進光量也變少 測空間範圍來選用合適的攝影機 用敏通公司的M7V561中解析度攝 ,有效像素為573(H)*505(V), mm,與基準線夾角b為50 ° ,所 所有的SSE0812,其焦距為8mm, 角c約為3 4 ° ,垂直視角約為
第18頁 1234640 _案號 91125314_年月日__ 五、發明說明(16) 3 0° ,調整光圈及像距為F 8 · 0,使景深範圍約為1 0 0 mm以 上,如此建立的影像量測空間e,可涵蓋所設定的X軸 8 0 m m * Y軸3 0 m m以及Z軸4 8 m m的特定量測空間f,再者,雷射 二極體1 5架設於CCD1 1 〇與CCD2 2 0中間。 【特徵記號設置】 特徵記號設置的目的是為了量測系統的校正,其型式 可依不同需求來設計,主要的特點要在影像中能容易辨識 以及能精確的計算其影像位置,如第4圖所示,可為下列 數種形式:圓形、環型、方形及十字形等,本發明中選用 十字形為空間特徵記號。 【三軸移動平台設置】 在建構量測系統的空間映射函數時,需要多個三維空 間特徵記號與影像點座標的相對應關係,所以使用精密的 三軸移動平台移動特徵記號及雙攝影機探頭,利用移動平 台精密的位移來得到特徵記號空間位置與其成像於雙攝影 機探頭之像素點的相對位置變化量。 【映射函數校正實驗】 首先將十字特徵記號放置於三維量測空間中相對零點 位置,其由攝影機CCD1及攝影機CCD2擷取的影像如第5圖 所示。設定X軸移動平台行程8 0 m m,Y軸行程3 0 m m,Z軸行 程4 8 m m,當移動間距X軸1 0 m m、Y軸5 m m、Z軸8 m m時,則X軸 移動9次,Y軸移動7次,Z軸移動7次,共計4 4 1次,可得到 4 4 1組三維空間十字特徵記號位置投影至雙攝影機影像點 座標的原始資料,如第6圖所示,係為將攝影機所得到的
第19頁 1234640 _案號 91125314_年月日__ 五、發明說明(17) 全部影像點資料同時顯示在同一平面的結果,其中,第6 a圖為左CCD影像;第6 b圖為右CCD影像。若以網格直線連 結時,如第7圖所示,第7a圖為左CCD影像;第7b圖為 右C C D影像。 【特徵記號影像位置重複性分析】
將空間網格特徵記號影像擷取程序執行8次,共計可 得7 0 5 6 ( 8 * 4 4 1 * 2 )個影像資料點,將各空間點的影像點位 置求其平均值。X軸影像位置與γ軸影像位置誤差圖及影像 位置誤差圖分別如第8圖、第9圖與第1 0圖所示。其X 軸影像最大誤差正值為0 . 7 3 9 4 p i X e 1 ,負值為-1 . 4 3 2 2 p i X e 1,標準差為0 · 0 7 9 8。Y軸影像最大誤差正值為 0.0480pixel ,負值為-0.2111 pixel ,標準差為 0.0480。 影像位置距離誤差為0 · 1. 4 4 1 7 p i X e 1 ,標準差為0 . 0 5 7 3。 由此可知X軸影像的重覆度略差於Y軸,可能的原因是 影像在本系統仍為類比的傳輸,在水平同步傳輸的準確有 誤差,而Y軸可視為數位傳輸,所以重覆度較好。以特徵 記號影像重覆度而言,在三個標準差的考量下,X軸可達 0 . 2 3 9像素點,而Y軸達0 · 1 7 2像素點,所以影像系統可達 1 / 4〜1 / 6次像素點的解析能力。 【映射函數擬合誤差分析】
將平均後的空間點座標與影像點位置座標代入映射函 數(1 7 )式,以最小平方法擬合資料,求得影像座標轉換 至空間中的點座標的映射函數係數。在本實驗中選用4階 擬合函數,可得到7 2項係數,表1列出前3 5項的映射函數
第20頁 1234640 _案號91125314_年月曰 修正_ 五、發明說明(18) 係數,因影像的位置輸入值範圍在X軸為0〜6 3 9,y軸為 0〜4 7 9,所以在高次項時雖係數值很小,但仍對空間的誤 差值有影響。 計算擬合誤差時,將影像座標代入映射函數(1 7 )式 中求得空間中的點座標,並計算其與原始空間座標的誤 差。X軸擬合結果如第1 1圖所示,最大擬合正誤差為
0.070 mm,最大擬合負誤差為-0.061mm,標準差為0.019。 Y軸擬合結果如第1 2圖所示,最大擬合正誤差為〇. 〇45mm ,最大擬合負誤差為-0. 0 4 5mm,標準差為0. 014。Z軸擬合 結果如第1 3圖,最大擬合正誤差為〇.〇36mm,標準差為 0 . 0 1 0。空間位置擬合結果如第1 4圖所示,最大擬合誤 差為0 . 0 7 6 m m,標準差為0 . 0 1 2。由此可知,系統擬合誤差 在三個標準差時在0.04mm以下。
第21頁 1234640 _案號91125314_年月日_ 五、發明說明(19) 表1左右攝影機影像點映射至空間點(χπζ)多項式擬合結果 項次 函式 X擬合係數 Y擬合係數 Ζ擬合係數 1 1 79.906695665944 14.982505437418 47.944017399549 2 0.000233597822 0.000043799591 0.000140158693 3 y. -0.002287320965 -0.000428872680 -0.001372392580 4 0.000233043251 0.000043695609 0.000139825951 5 乃 0.003153670719 0.00059131 3259 0.001892202433 6 2 ^1 -0.000000065389 -0.000000012261 -0.000000039234 7 0.000004351014 〇.〇〇〇〇〇〇815815 0.000002610609 8 又丨“2 -0.000000103231 -0.000000019356 -0.000000061938 9 -0.000006805409 -0.000001276014 -0.000004083245 10 y.2 -0.000035131039 -0.000006587070 -0.000021078624 11 儿* & 0.000001910260 〇.〇〇〇〇〇〇358174 0.000001146156 12 7 丨 *y2 0.000081828337 0.000015342813 0.000049097002 13 2 -0.000000039094 -0.000000007330 -0.000000023457 14 尤2*y2 -0.000004396477 -0.000000824339 -0.000002637886 15 2 y2 -0.000048578686 -0.000009108504 -0.000029147212 16 3 -0.000000000119 -0.000000000022 -0.000000000071 17 xi2*y. -0.000000001394 -0.000000000261 -0.000000000836 18 i “2 -0.000000000394 -0.000000000074 -0.000000000236 19 〇.〇〇〇〇〇〇〇〇2868 〇.〇〇〇〇〇〇000538 〇.〇〇〇〇〇〇〇〇1721 _ιιι 第22頁 1234640 _案號91125314_年月日 修正 五、發明說明(20) 項次 函式 x擬合係數 Y擬合係數 Ζ擬合係數 20 x' * y? 〇.〇〇〇〇〇〇〇47658 〇.〇〇〇〇〇〇〇〇8936 〇.〇〇〇〇〇〇〇28595 21 X丨木3^丨* 〇.〇〇〇〇〇〇000805 〇.〇〇〇〇〇〇〇〇〇151 〇.〇〇〇〇〇〇〇〇〇483 22 U * h -0.000000117608 -0.000000022051 -0.000000070565 23 * Χ2 -0.000000000434 -0.000000000081 -0.000000000260 24 又丨木^ 〇.〇〇〇〇〇〇〇〇2086 〇.〇〇〇〇〇〇〇〇〇391 〇.〇〇〇〇〇〇〇〇1252 25 〇.〇〇〇〇〇〇〇75802 〇.〇〇〇〇〇〇〇14213 〇.〇〇〇〇〇〇〇4 5481 26 3 yi -0.000000150259 -0.000000028174 -0.000000090155 27 yi2 〇.〇〇〇〇〇〇023087 〇.〇〇〇〇〇〇〇〇4329 〇.〇〇〇〇〇〇013852 28 2 v yi * 〇.〇〇〇〇〇〇553057 〇.〇〇〇〇〇〇103698 〇.〇〇〇〇〇〇331834 29 v 2 〇.〇〇〇〇〇〇〇〇2145 〇.〇〇〇〇〇〇〇〇〇402 〇.〇〇〇〇〇〇001287 30 yi -0.000000057977 -0.000000010871 -0.000000034786 31 y. * y^2 -0.000000661368 -0.000000124006 -0.000000396821 32 3 -0.000000000158 -0.000000000030 -0.000000000095 33 2 v & -0.000000000715 -0.000000000134 -0.000000000429 34 〇.〇〇〇〇〇〇〇40862 〇.〇〇〇〇〇〇〇〇7662 〇.〇〇〇〇〇〇〇24517 35 3 〇.〇〇〇〇〇〇258261 〇.〇〇〇〇〇〇〇4 84 24 〇.〇〇〇〇〇〇154956 【映射函數影像量測解析度分析】 將平均後的空間點座標與其影像點位置座標代入映射 函數(1 6 )式,求得空間中的點座標轉換至影像座標的映 射函數係數,再依照【空間映射函數影像量測解析度】之 定義計算出空間中的影像量測解析度,其結果如第1 5圖 及第1 6圖所示,X軸依遠近的不同,約由0.17mm/pixel 至0.43mm/pixel,而Y軸及Z軸約在0.2mm/pixel以内。整
第23頁 1234640 案號 91125314 Λ_η 曰 修正 五、發明說明(21) 體而言,空間解析度在0.3mm/pixel至0.5mm/pixel間。第 1 7圖及第1 8圖是以3D立體圖來表示解析度與量測空間 的關係,其中,第1 7圖為CCD 1影像點映射至空間點=^ 解析度立體圖;第1 8圖為CCD 2影像點映射至空間點空間 解析度立體圖。 由此可清楚看出解析度與量測空間距(:(:^攝影機 土 近有決定性的影響。越近則解析度越好,反之,則較逐 所以選擇量測範圍時亦決定了解析度的範圍,並且 二 時要盡量靠近CCD,才能得到較好的結果。其像素里測 分布於0· 4mm/Pixel 至 〇. 6 mm/Pixei 間。 厅度
Rn綜合以上結果可知,在量測系統所設定的 〇mm*30mm*48mm的量測空間中,利 特徵記號,其所得到影像位置重覆】十字 下,可達到1/4次像素點的重覆产, 一^車差的情況 度分析中,可知道因量測处 又旦 里測空間的解析 3的影響,其空間解析度至?5雙攝影機遠 〇·ΐ25_間,而映射函 啊度為〇·〇75ππη至 〇. 〇4mm以下,合乎李鞭0决差在三個標準差時為 ,升量= = 需;,整體而言,在本系 ί : J ::大小、提高CCD像素點數目、及切以的數做法有:慎選 J以次能t等方式。 輪 發先几件與控制電路元件所ί度j:器’主要由 〃 兔光元件係為
第24頁 1234640 _案號91125314_年月日__ 五、發明說明(22) 定式或旋轉式的雷射二極體,而所投射之特徵記號,其結 構光源可為點光源、線光源或面光源其中之一。在本發明 中,係以簡單之實驗設計,乃利用固定式的雷射二極體產 生十字特徵記號,於實際應用上,可利用線光源或面光源 作為特徵記號之結構光源,搭配旋轉式雷射二極體,進行 量測物之非接觸式三維空間量測,其系統乃不超脫本發明 所提之方法與裝置。
承前所述,本發明乃提出一種非接觸式三維空間量測 方法,透過雙C CD攝影機與三維空間映射函數關係式的運 算,可以找到其在真實空間中的實際位置座標值,不僅有 效減少三維尺寸量測時所須之龐大運算,更可提升量測的 速度與精度堪稱具創作性與進步性,符合發明專利之法定 要件,爰依法提出發明專利申請。 雖本發明以一較佳實施例揭露如上,但並非用以限定 本發明實施之範圍。任何熟習此項技藝者,在不脫離本發 明之精神和範圍内,當可作些許之更動與潤飾,即凡依本 發明所做的均等變化與修飾,應為本發明專利範圍所涵蓋 ,其界定應以申請專利範圍為準。
第25頁 1234640 _案號9Π25314_年月曰 修正_ 圖式簡單說明 第1 a圖係為雙攝影機量測架構示意圖。 第1 b圖係為攝影機之影像示意圖。 第2圖係為非接觸式三維空間量測流程示意圖。 第3圖係為雙攝影機設計示意圖。 第4圖係為各種特徵記號形式示意圖。 第5圖係為特徵記號擷取影像示意圖。 第6圖係為特徵記號投影至雙攝影機點之原始座標資料示 意圖:(a )左CCD影像,(b )右CCD影像。 第7圖係為特徵記號在雙攝影機點座標資料之網格示意圖 :(a )左CCD影像,(b )右CCD影像。 第8圖係為特徵記號X轴影像位置誤差圖。 第9圖係為特徵記號Y軸影像位置誤差圖。 第1 0圖係為影像位置誤差圖。 第1 1圖係為特徵記號X軸映射誤差圖。 第1 2圖係為特徵記號Y軸映射誤差圖。 第1 3圖係為特徵記號Z軸映射誤差圖。 第1 4圖係為空間位置擬合誤差圖。 第1 5圖係為CCD1在量測空間的解析度。 第1 6圖係為CCD2在量測空間的解析度。 第1 7圖係為CCD 1影像點映射至空間點(X,Y,Z )空間解析 度(mm/pixel)立體圖。 第1 8圖係為CCD2影像點映射至空間點(X,Υ,Z)空間解析 度(mm/pixel)立體圖。
第26頁 1234640 案號 91125314 曰 修正 圖式簡單說明 圖 號 簡 單 說 明 a • 雙 攝 影 機 架 設 距 離 b 夾 角 C 水 平 視 角 e • · •影像 量 測 空 間 f • 特 定 量 測 空 間 1 0 · • CCD1 攝 影 機 1 2 • 左 側 影 像 圖 1 5 · • •雷 射 二 極 體 2 0 CCD2 攝 影 機 2 2 · • •右 側 影 像 圖 3 0 矩 形 量 測 空 間 4 0 · • · · 特 徵 記 號 1 0 0 • • X座標 2 0 0 • · · • Y座標 3 0 0 5 0 0 5 10 5 2 0 5 3 0 5 4 0 5 5 0 6 0 0 6 10 10 0 0 • · Ζ座標 特徵記號 三軸移動平台之定位 雙CCD攝影機之取像 三維空間映射函數校正 求得該式中的各項函數係數值 高次多項式三維空間映射函數 量測物 三維空間座標(X,Υ,Ζ ) 三維空間
第27頁

Claims (1)

1234640 案號91125314 年月日 修正 六、申請專利範圍 1 · 一種非接觸式三維空間量測方法,係包含: 一影像之擷取,乃由至少雙部影像感測元件、雷射光 度控制器與三軸移動平台所構成,利用雷射光度控 制器產生一特徵記號於量測物上,並配合三轴移動 平台與影像感測元件,可獲得特徵記號之空間實際 座標與影像平面座標; 一三維空間座標與影像平面座標之運算,藉由一三維 空間映射函數係數值可進行二維平面座標與三維空 間座標之轉換;與
一電腦輔助處理,負責處理影像資料與進行數值資料 之運算,故,先透過影像擷取系統獲得特徵記號之 實際座標與影像平面座標,將三維空間影像量測與 雙或多影像感測元件影像點座標的對應關係以三維 空間映射函數轉換,可求出三維空間映射函數之函 數係數值,再利用該函數係數值與量測物之影像平 面座標,即可獲得量測物之三維空間座標。 2 ·如申請專利範圍第1項所述之非接觸式三維空間量測 方法,其中,影像感測元件係為電荷耦合器(Charge -Coupled Device, CCD ) 〇 3 ·如申請專利範圍第1項所述之非接觸式三維空間量測
方法,其中,雷射光度控制器主要由發光元件與控制 電路元件所構成。 4 ·如申請專利範圍第3項所述之非接觸式三維空間量測 方法,其中,發光元件係為雷射二極體。
第28頁 1234640 __案號 91125314_年月日__ 六、申請專利範圍 5 ·如申請專利範圍第3項所述之非接觸式三維空間量測 方法,其中,發光元件所投射的特徵記號,其結構光 源係為點光源、線光源或面光源其中之一。 6 ·如申請專利範圍第3項所述之非接觸式三維空間量測 方法,其中,發光元件係為固定式裝置。 7 ·如申請專利範圍第3項所述之非接觸式三維空間量測 方法,其中,發光元件係為旋轉式裝置。
第29頁
TW91125314A 2002-10-25 2002-10-25 Contactless three dimensional measure method TWI234640B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW91125314A TWI234640B (en) 2002-10-25 2002-10-25 Contactless three dimensional measure method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW91125314A TWI234640B (en) 2002-10-25 2002-10-25 Contactless three dimensional measure method

Publications (1)

Publication Number Publication Date
TWI234640B true TWI234640B (en) 2005-06-21

Family

ID=36597888

Family Applications (1)

Application Number Title Priority Date Filing Date
TW91125314A TWI234640B (en) 2002-10-25 2002-10-25 Contactless three dimensional measure method

Country Status (1)

Country Link
TW (1) TWI234640B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI393853B (zh) * 2009-12-29 2013-04-21 Metal Ind Res & Dev Ct 三維空間座標量測裝置及其方法
TWI484138B (zh) * 2013-04-02 2015-05-11 Track laser scanning equipment and method
TWI662694B (zh) * 2017-12-20 2019-06-11 緯創資通股份有限公司 三維影像攝取方法及系統

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI393853B (zh) * 2009-12-29 2013-04-21 Metal Ind Res & Dev Ct 三維空間座標量測裝置及其方法
TWI484138B (zh) * 2013-04-02 2015-05-11 Track laser scanning equipment and method
TWI662694B (zh) * 2017-12-20 2019-06-11 緯創資通股份有限公司 三維影像攝取方法及系統
US10542247B2 (en) 2017-12-20 2020-01-21 Wistron Corporation 3D image capture method and system

Similar Documents

Publication Publication Date Title
CN111536902B (zh) 一种基于双棋盘格的振镜扫描系统标定方法
CN108120392B (zh) 气液两相流中气泡三维测量系统及方法
CN104315995B (zh) 基于虚拟多立方体标准目标的tof深度相机三维坐标标定装置和方法
JP6537237B2 (ja) 情報処理装置および方法
CN105716527B (zh) 激光焊缝跟踪传感器标定方法
CN110378969B (zh) 一种基于3d几何约束的汇聚式双目相机标定方法
CN104316083B (zh) 一种虚拟多球体球心定位的tof深度相机三维坐标标定装置和方法
CN109373912B (zh) 一种基于双目视觉的非接触式六自由度位移测量方法
WO2018201677A1 (zh) 基于光束平差的远心镜头三维成像系统的标定方法及装置
TW201812700A (zh) 多維度尺寸量測系統及其方法
Zhou et al. A novel laser vision sensor for omnidirectional 3D measurement
CN109272555B (zh) 一种rgb-d相机的外部参数获得及标定方法
Fernandez et al. Planar-based camera-projector calibration
CN102980526A (zh) 用黑白相机获取彩色图像的三维扫描仪及其扫描方法
WO2020199439A1 (zh) 基于单双目混合测量的三维点云计算方法
JP2015106287A (ja) キャリブレーション装置及び方法
JP2010276433A (ja) 撮像装置、画像処理装置及び距離計測装置
CN102881040A (zh) 一种数码相机移动拍摄三维重建方法
CN113804128A (zh) 一种双轴承孔同轴度误差视觉测量装置及测量方法
Xu et al. 3-D reconstruction of binocular vision using distance objective generated from two pairs of skew projection lines
CN112489109A (zh) 一种三维成像系统方法、装置及三维成像系统
TWI234640B (en) Contactless three dimensional measure method
Im et al. A solution for camera occlusion using a repaired pattern from a projector
Chang et al. Non-contact scanning measurement utilizing a space mapping method
CN114993207B (zh) 基于双目测量系统的三维重建方法

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees