Thực nghiệm thu ảnh 3D của chi tiết cơ khí

n Ngu nồ chi u sángế q p q 1 p q p P P P

Hình 2-18 Biểu diễn chùm tia chiếu tới vật thể dưới dạng mặt phẳng

Mặt phẳng P được biểu diễn dưới dạng 1 điểm gốc q và 2 vec tơ1 và2 Với bất kỳ điểm p thuộc mặt phẳng P đều thỏa mãn phương trình sau:

� = {� = � + (( + !! ∶ (�! ∈ ℝ} (2-50) Phương trình (2-50) không phải là cách biểu diễn duy nhất, điểm q có thể thay thế

bằng một điểm khác nằm trên mặt phẳng và hai véc tơ1,2 có thể thay thế bằng bất kỳ tổ hợp tuyến tính khác của1,2

Trong trường hợp khác, người ta thường biểu diễn mặt phẳng dưới dạng một điểm q nằm trên mặt phẳng và một véc tơ pháp tuyến nℝ):

� = {�: �/(� − �) = 0} (2-51)

Cũng tương tự như phương trình (2-50) phương trình (2-51) cũng không phải là các biểu diễn duy nhất, điểm q có thể thể thay thế bằng bất kỳ điểm nào khác trong mặt phẳng P và véc tơ pháp tuyến n bằng bất kỳ bội số vô hướng �� 0

Để chuyển đổi giữa phương trình (2-50) và phương trình (2-51) chúng ta có thể lấy � =( × ! là tích của 2 véc tơ cơ sở1 và2, và đây là 2 véc tơ cần tìm Trong thực tế, chỉ cần tìm một véc tơ1 trực giao với n, véc tơ2 = n1

Như vậy, một đường thẳng L cũng có thể được biểu diễn bằng giao của 2 mặt phẳng theo công thức sau:

� = {�: �(/(� − �) = �!/(� − �) = 0} (2-52) Trong đó, 2 véc tơ pháp tuyến �1 và �2 độc lập tuyến tính, khi đó ta xác định được điểm chung của 2 mặt phẳng thay vì 2 điểm khác nhau Như vậy phương trình (2-52) được coi là phương trình giao điểm giữa mặt phẳng của hệ chiếu sáng tại các điểm giao p trên bề mặt vật

Thông thường các hệ chiếu sáng được sử dụng trong quá trình quét Laser và ánh sáng cấu trúc chứa các đường hoặc điểm có thể xác định được Hệ chiếu sáng được đơn giản hóa theo mô hình lỗ nhỏ, tức là một đường chiếu tạo ra một mặt phẳng ánh sáng; mặt phẳng này là duy nhất chứa đường thẳng trên mặt phẳng hình ảnh và tâm chiếu, với mỗi điểm chiếu được tạo ra bởi một tia sáng (tia này tạo bởi tâm của hệ chiếu sáng và

điểm chiếu trên vật) và cũng là một điểm duy nhất, giao điểm của một mặt phẳng ánh sáng với vật thể thành nhiều đoạn cong được chiếu sáng (hình 2-18) là tập bởi các điểm được chiếu Các điểm này được hệ máy ảnh thu được và xác định được tia ảnh duy nhất từ 1 điểm được chiếu sáng và tâm của hệ máy ảnh Hiện tại, chúng ta giả định rằng vị trí và hướng của máy chiếu và máy ảnh được biết đến theo hệ tọa độ toàn cầu Theo giả thiết này, phương trình của mặt phẳng và tia chiếu, cũng như phương trình của tia máy ảnh tương ứng với các điểm được chiếu sáng, được xác định bằng các tham số có thể đo được Từ các phép đo này, vị trí của các điểm được chiếu sáng có thể được khôi phục bằng cách giao nhau giữa các mặt phẳng hoặc tia sáng với tia hệ thu ảnh tương ứng với các điểm được chiếu sáng Thông qua các quá trình như vậy, sự không rõ ràng về độ sâu do phép chiếu mô hình lỗ nhỏ có thể được loại bỏ, cho phép khôi phục mô hình bề mặt 3D

Gọi �G là tâm của hệ chiếu sáng, P là mặt phẳng tạo bởi �G và các tia chiếu sáng tới vật, n là véc tơ pháp tuyến của mặt phẳng P Gọi L là đường thẳng đi qua tâm �< của hệ thu ảnh và điểm p là giao điểm giữa mặt phẳng chiếu sáng và vật (theo hình 2-19) Đường thẳng L được biểu diễn ở dạng tham số:

� = {� = �< + �� ∶ � ∈ ℝ} Và mặt phẳng P được biểu diễn dưới dạng:

� = x�: �/{� − �G | = 0y

(2-53)

(2-54)

Hình 2-19 Tam giác đạc giao điểm giữa mặt phẳng và đường thẳng (nguồn: [79])

Trong một số trường hợp đặc biệt, đường thẳng và mặt phẳng song song nên không có điểm p giao nhau, hay nói cách khác véc tơ � và n trực giao �/�=0 khi này cũng xảy ra trường hợp đường thẳng L nằm trong mặt phẳng P Do vậy, trong bài toán này, chúng ta chỉ xét trong trường hợp véc tơ � và n không trực giao thỏa mãn điều kiện: �/� 0 khi đó giao của đường thẳng L và mặt phẳng P luôn có một giao điểm p Vì điểm p nằm trên đường thẳng L nên ta có thể viết � = �< + ��, và cần xác định giá trị của � Mặt khác, điểm p cũng nằm trên mặt phẳng cho nên giá trị � cần thỏa mãn phương trình tuyến tính:

�� !� − ��" = �� !�� + �� − ��" = 0 Và giá trị � được xác định rút ra từ công thức (2-55)

(2-55)

� = �/ { �G − � < | �/�

(2-56)

Như đã nói ở trên, công thức (2-57) chỉ đúng trong điều kiện đường thẳng L và mặt phẳng P không song song hay nói cách khác �/� 0 Trong thực thế sẽ có một tập hợp các giao điểm p trên bề mặt vật, chúng ta cần xác định các điểm p đó trong điều kiện hẹp, nghĩa là p là giao điểm của 2 đường thẳng giao nhau trong không gian 3D theo hình 2-20

Hình 2-20 Tam giác đạc bởi giao điểm giữa 2 đường thẳng (nguồn: [79])

Giá sử ta có đường thẳng L1 chứa điểm q1 là gốc của hệ chiếu sáng và p nằm trên vật, đường thẳng L2 chứa điểm q2 là gốc của hệ thu ảnh và điểm p Như vậy biểu diễn 2 đường thẳng L1 và L2 dưới dạng:

�( = {� = �( + �(�( ∶ �( ∈ ℝ}

�! = {� = �! + �!�! ∶ �! ∈ ℝ}

(2-57)

(2-58) Xét trong một số trường hợp, véc tơ �( và véc tơ �! có thể là phụ thuộc tuyến tính

(nếu một trong hai véc tơ là bội số vô hướng của véc tơ còn lại) hoặc độc lập tuyến tính Hai véc tơ �( và véc tơ �! phụ thuộc tuyến tính trong trường hợp hai đường thẳng song song, trong trường hợp này chúng không cắt nhau và không trùng nhau Trong trường hợp khác nếu chúng độc lập tuyến tính thì véc tơ �( và véc tơ �! có thể cắt nhau hoặc không

Xét trong trong trường hợp chúng cắt nhau thì giao điểm p đó là giao điểm duy nhất thỏa mãn điều kiện cần và đủ để khi véc tơ �( và véc tơ �! là độc lập tuyến tính và các giá trị vô hướng �( �à �! sao cho:

�( + �(�( = �! + �!�! (2-59)

Hoặc tương tự véc tơ �( − �! là phụ thuộc tuyến tính vào véc tơ �( và véc tơ �! Xét trường hợp hai đường thẳng không chắt nhau, ta có xác định giao điểm gần đúng là “điểm gần nhất với 2 đường thẳng” Hay nói cách khác, cho dù hai đường thẳng cắt nhau hay không thì vẫn có thể xác định được giao điểm gần đúng là điểm p thỏa mãn điều kiện tổng bình phương khoảng cách nhỏ nhất đến cả 2 đường thẳng:

�(�, �(, �! ) = ‖�( + �(�( − �‖! + ‖�! + �!�! − �‖! (2-60) Để xác định giá trị của p, hàm �(�, �(, �! ) được coi là một hàm không âm, 5 biến số bao gồm: Giá trị tọa độ 3 chiều của điểm p và hai biến số vô hướng �( và �!

Hình 2-21 Xác định giao điểm gần đúng của 2 đường thẳng (nguồn: [79])

Trước hết cần xác định hai biến số �( và �! của hàm không âm bậc 2 Giả sử điểm �( = �( + �(�( là điểm nằm trên đường thẳng �( và điểm �! = �! + �!�! nằm trên đường thẳng �! Trung điểm của đoạn thẳng nối giữa �( và �! là �(!, ta cố phương trình sau:

1 1 2 2

(2-61) Điều kiện để cho (�, �(, �! ) của hàm � là nhỏ nhất khi xét đạo hàm riêng của �

theo 5 giá trị là nhỏ nhất, theo đó ta xét đạo hàm riêng của � (công thức 2-60) theo 3 biến của tọa độ điểm p:

�

Φ

�� = (� − �() + (� − �!) = 0

(2-62)

Giá trị p là nhỏ nhất nếu nó là trở thành trung điểm �(! của đoạn thẳng nối giữa �( và �!; Như vậy cần xác định bình phương tối thiểu khoảng cách từ một điểm �( trên đoạn thẳng �( và một điểm �! trên đoạn thẳng �! Trong trường hợp này xét tìm nhỏ nhất của phương trình không âm 2 biến:

(

Ψ �(, �!) = 2�(�(!, �(, �!) = (‖ �! + �!�!) − (�( + �(�()‖! (2-63) Từ công thức (2-22) ta xác định đạo hàm riêng của từ biến số �(, �! nhỏ nhất tiến dần tới 0: � Ψ ��( � Ψ ��! = �(/(�(�( − �!�! + �( − �!) = �(‖�(‖! − �!�(/�! + �(/(�( − �!) = �!/(�!�! − �(�( + �! − �() = �!‖�!‖! − �!�!/�( + �!/(�! − �() (2-64) (2-65)

Hai phương trình tuyến tính (2-64), (2-65) biến �(, �! có thể được thu gọn và biểu diễn dưới dạng ma trận như sau:

‖�(‖! −�!/�( −�(/�! ‖�!‖! � � / (�! − �() £ (2-66)

Phương trình (2-25) phụ thuộc vào tuyến tính độc lập của 2 véc tơ �( và véc tơ �!, ma trận 2x2 có thể được thu gọn lại:

�! −�!/�(

−�(/�!

‖�!‖! 8(

£ (2-67)

Phương trình (2-67) có thể viết dưới dạng:

� �! 1 ‖�‖! ‖�(‖!‖�!‖! − (�(/�!)! �!�( �(/�! ‖�!‖! � /(�! − �() �!(�! − �() £ (2-68) £ ¢ (/ ¢ /! ¢ ( £ = �!(�! − �() � /(�! − �() ¢ (/ ¢ ( £ = ¢ � ‖�(‖! £ ¢ ( £ = ¢ (/ �(! = �( + (�! − �() = �! + (�( − �!)

Theo công thức (2-68) giao điểm gần đúng p có thể xác định được sau khi tìm được giá trị của biến số �( ℎ�ặ� �!

2 4 2 Mô hình lỗ nhỏ của máy ảnh

Đối với một hệ thống đo sử dụng phương pháp ánh sáng cấu trúc, máy ảnh có vai trò thu nhận hình ảnh các tia sáng phản xạ trên bề mặt vật từ hệ chiếu sáng Để thể hiện quá trình tạo ảnh vật thể của một hệ thu ảnh người ta sử dụng mô hình lỗ nhỏ, mô hình này bao gồm một mặt phẳng ảnh I và một tâm điểm Oc, ta có thể biểu diễn mô hình như sau:

Hình 2-22 Mô hình hệ thu ảnh lỗ nhỏ (nguồn: [79])

Tuy nhiên, với mô hình như trên hình 2-22, ảnh của vật bị đảo ngược, do đó, người ta chuyển sang một mô hình tương đương để dễ dàng hơn trong quá trình tính toán Trên hình 2-23, đổi vị trí giữa mặt phẳng ảnh và mặt phẳng lỗ nhỏ

Hình 2-23 Mô hình hệ thu ảnh lỗ nhỏ chuyển đổi (nguồn: [79])

Xét mô hình lỗ nhỏ (hình 2-22) với hệ tọa độ máy ảnh gồm tâm của phép chiếu q và hệ trực chuẩn cơ sở {�(, �!, �)} và mặt phẳng ảnh tạo bởi 2 véc tơ cơ bản tại tiêu cự f =1 Với bất kỳ điểm p trong không gian 3 chiều của hệ tọa độ thực có tọa độ (�(, �!, �)) đều có mối liên hệ với tọa độ của máy ảnh Hình chiếu của điểm p lên mặt phẳng ảnh là điểm ảnh 2 chiều u được xác định bởi thông số �( và �), ta có thể viết tọa độ điểm u dưới dạng véc tơ 3 chiều � = (�(, �!, 1), giả sử ta có là một tham số vô hướng không âm, mối liên hệ giữa điểm p thực và điểm ảnh u có thể được biểu diễn như sau:

�(

¥�!¦ = � §�!¨ �) 1

(2-69) Gọi o là tâm chiếu của gốc tọa độ của hệ tọa độ thực có tọa độ (0, 0, 0)/, điểm q và véc tơ �( và �! được xác định theo:

1 [�(|�!|�] = §0 0 0 1 0 0 0¨ 1 (2-70) Với mỗi hệ thu ảnh được gắn với một hệ tọa độ ảnh, một điểm p trong hệ tọa độ

thực được biểu diễn dưới dạng véc tơ �K = (�K(, �K!, �K))/, trong hệ tọa độ ảnh cũng có dạng véc tơ �9 = (�9(, �9!, �9))/, mối quan hệ giữa 2 véc tơ �K và �9 được biến đổi bằng một ma trận dịch chuyển � ∈ ℝ) và một ma trận quay � ∈ ℝ)4):

�9 = ��K + � (2-71)

Xét trong hệ tọa độ của máy ảnh, mối quan hệ giữa tọa độ điểm trong không gian 3 chiều và tọa độ ảnh 2 chiều của phép chiếu bằng mô hình lỗ nhỏ với công thức (2-69) với �� = �9 , mối liên hệ với tọa độ thực với tọa độ ảnh được thể hiện bằng công thức:

�� = ��K + � (2-72)

Như vậy, các tham số (R, T) được gọi là các tham số ngoại của máy ảnh biểu diễn vị trí và hướng của máy ảnh so với hệ tọa độ thực

Giả sử đơn vị đo độ dài trên mặt phẳng ảnh bằng với hệ tọa độ thực, khoảng cách từ tâm hình chiếu đến mặt phẳng ảnh bằng một đơn vị độ dài và gốc của hệ tọa độ ảnh có tọa độ �( = 0 và �! = 0 Tuy nhiên trong thực tế, trường hợp này không xảy ra vì đơn vị đo trên mặt phẳng ảnh là điểm ảnh, còn với hệ toạ độ thực là đơn vị chiều dài (mét, milimet hoặc inch,…), khoảng cách từ tâm chiếu đến mặt phẳng ảnh có thể thay đổi tùy theo từng máy ảnh và theo quy định hệ tọa độ của ảnh là nằm trên cùng bên trái Ngoài ra, mặt phẳng ảnh có thể bị nghiêng so với mặt mô hình mặt phẳng ảnh lý tưởng, do vậy để thể hiện sự không đồng nhất giữa mô hình lý tưởng và mô hình thực tế, người ta sử dụng một ma trận � ∈ ℝ)4) trong phương trình chiếu để mô tả các thông số nội tại của một hệ máy ảnh:

Ma trận K có dạng: �� = �(��K + �) (2-73) � �( � = § 0 0 � �L � �! 0 �( �!¨ 1 (2-74) �(

Trong đó f là chiều dài tiêu cự (được xác định là khoảng cách giữa tâm chiếu đến mặt phẳng ảnh) Các tham số �(, �! là tham số thể hiện tỷ lệ tọa độ thứ nhất và tọa độ thứ 2 do các hệ máy ảnh có kích thước điểm ảnh là hình chữ nhật Tham số �L để bù cho mặt phẳng ảnh bị nghiêng (�(, �!)/ là tọa độ ảnh của giao điểm của đường thẳng đứng trong tọa độ máy ảnh với mặt phẳng ảnh, điểm này gọi là tâm ảnh hoặc điểm chính Các tham số của ma trận K không phụ thuộc vào các vị trí đặt máy ảnh so với vật mà nó thể hiện các đặc tính vật lý liên quan đến thiết kế cơ học và quang học của máy ảnh Nếu không có sự điều chỉnh khoảng cách giữa các thành phần linh kiên trong máy ảnh, thì ma trận K không thay đổi, do vậy chúng ta chỉ cần xác định thông qua quá trình hiệu chuẩn ban đầu và được sử dụng trong suốt quá trình tính toán sau này mà không cần phải thực hiện lại Như vậy sau khi xác định được ma trận K thì ta cũng đã có được phương trình (2-73) thể hiện mối quan hệ của một điểm 3 chiều trong tọa độ thực với tọa độ ảnh 2 chiều

Với mỗi điểm ảnh có tọa độ � = (�(, �!, 1)/ theo hình 2-20 và công thức (2-53) � = {� = �< + �� ∶ � ∈ ℝ} ta luôn xác định được một đường thẳng duy nhất nối với tâm chiếu, để tính toán tọa độ 3D của một điểm ảnh là cần xác các tham số trong hệ phương trình này Như ta đã biết, với mỗi một điểm trong hệ tọa độ thực luôn nằm trên một đường thẳng nối dài từ tâm chiếu đi qua mặt phẳng ảnh, do đó để thể hiện tọa độ của một điểm thực đó từ công thức (2-32) ta có thể viết lại:

�K = (−�/�) + �(�/�) (2-75)

Với �8( = �/ Từ công thức (2-53) và công thức (2-75) ta rút ra kết luận đường thẳng đi qua một điểm q có tọa độ trong hệ tọa độ thực là �K = −�/� nếu q là tâm chiếu và véc tơ � trong hệ toạ độ thực là �K = �/�

Hình 2-24 Mặt phẳng ảnh được xác định bởi đường thẳng trên ảnh và tâm chiếu

Một đường thẳng trên ảnh được biểu diễn bằng 8 tham số trên tọa độ ảnh qua công thức sau:

� = � ∶ �/� = �(�( + �!�! + �) = 0

52

Trong đó � = (�(, �!, �))/ là một véc tơ với �( ≠ 0 hoặc �! ≠ 0, khi sử dụng hệ chiếu sáng để chiếu các vân lên bề mặt vật quét sẽ bao gồm các đường theo chiều dọc và chiều ngang, do đó ta có thể biểu diễn đường thẳng nằm ngang dưới dạng:

�M = {�: �/� = �! − � = 0} (2-77) Trong đó � là tọa độ thứ 2 của một điểm trên đường thẳng Trong trường hợp này ta có thể coi � = (0, 1, −� )/ Tương tự như vậy đối với đường thẳng nằm đứng:

�N = {�: �/� = �( − � = 0} (2-78)

Trong đó � là tọa độ thứ 1 của một điểm trên đường thẳng Trong trường hợp này ta có thể coi � = (1, 0, −� )/ Từ đó ta có thể xác định được 1 mặt phẳng duy nhất tạo bởi đường thẳng L và gốc chiếu Với mỗi điểm ảnh với tọa độ ảnh u trên đường thẳng L, đường thẳng sẽ luôn bao gồm một điểm thuộc mặt phẳng P và tâm chiếu Giả sử p là