Bài giảng Xử lý ảnh số - Chương 7: Xử lý ảnh dựa trên Morphological - Ngô Quốc Việt
Giới thiệu Binary Morphology
2. Một số khái niệm về tập hợp
3. Các phép toán cơ bản
Phát triển và bào mòn (Dilation and erosion)
Mở và đóng (Opening and closing)
Biến đổi trúng hay trật (Hit-or-miss transformation)
4. Ứng dụng morphology vào các bài toán xử lý ảnh
Dò biên
Trích các thành phần liên thông (Extraction of Connected
Components)
Tìm bao lồi (Convex Hull)
Tô vùng (Region Filling)
Làm mảnh (Thinning)
Làm dầy (Thickening)
Lọc xương/Lấy cạnh trung tâm (Skeleton)
2. Một số khái niệm về tập hợp
3. Các phép toán cơ bản
Phát triển và bào mòn (Dilation and erosion)
Mở và đóng (Opening and closing)
Biến đổi trúng hay trật (Hit-or-miss transformation)
4. Ứng dụng morphology vào các bài toán xử lý ảnh
Dò biên
Trích các thành phần liên thông (Extraction of Connected
Components)
Tìm bao lồi (Convex Hull)
Tô vùng (Region Filling)
Làm mảnh (Thinning)
Làm dầy (Thickening)
Lọc xương/Lấy cạnh trung tâm (Skeleton)
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Xử lý ảnh số - Chương 7: Xử lý ảnh dựa trên Morphological - Ngô Quốc Việt", để tải tài liệu gốc về máy hãy click vào nút Download ở trên.
File đính kèm:
- bai_giang_xu_ly_anh_so_chuong_7_xu_ly_anh_dua_tren_morpholog.pdf
Nội dung text: Bài giảng Xử lý ảnh số - Chương 7: Xử lý ảnh dựa trên Morphological - Ngô Quốc Việt
- Cho Z là SE và S lưới vuông các vị trí pixlel chứa tập {(r,c), (-r,-c) | (r,c) supp(Z)}. Thì Z(a,b) Z( a, b), (a,b) sup p(Z) Là phản chiếu của structuring element. Z là Z quay 180º quanh gốc của nó. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 11
- Phép toán Dilation đối tượng A theo đối tượng Z (thường là SE) được định nghĩa bởi D(A,Z) AZ a z | a A, z B Đối tượng A tịnh tiến theo từng vectơ (vị trí) trong Z. Kết quả sau khi tịnh tiến tạo ra các đối tượng A1,A2, . An (n là số phần tử trong Z). D(A,Z) là hợp các Ai. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 12
- Một số định nghĩa khác về dilation. Tập tất cả vị trí pixel p, trong I Z p SP (Z p) I . ảnh mà phần giao của Ž+p với I khác rỗng. I Z (Z p). Hội những copies của SE, mỗi bản copy p supp I là dịch chuyển đến vị trí pixel foreground trong ảnh. I Z (I p). Hội những copies của ảnh, mỗi bản là dịch p supp Z chuyển đến vị trí pixel foreground của SE. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 13
- A = {(2,1),(3,1),(4,1),(3,2)} Z = {(0,0),(0,1)} AZ A1 = {(2,1),(3,1),(4,1),(3,2)} A2 = A tịnh tiến bởi vector (0,1) của Z A1 = A tịnh tiến bởi vector (0,0) của Z A2 = {(2,2),(3,2),(4,2),(3,3)} D(A,B) = A1 U A2 = {(2,1),(3,1),(4,1),(3,2), (2,2),(4,2),(3,3)} Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 14
- Sử dụng SE 3x3 – Z8 có dạng cho phép mở rộng đều về cả 8 hướng. 1 1 1 Kết quả dilate trên ảnh mẫu như sau 1 1 1 1 1 1 Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 15
- SE được áp lên mọi pixel của ảnh. Tâm của SE được kết hợp với từng pixel, toàn bộ SE được áp cho pixel đang xét theo cách thay thế pixel đó bằng SE. Khái niệp “áp” là hoạt động “công logic nhị phân giữa các giá trị 0, 1”. SE Ảnh nguồn Đã giãn Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 16
- Tham khảo mã nguồn openCV_Dilate Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 17
- Phép toán Erosion đối tượng A theo đối tượng Z được định nghĩa bởi: E(A,Z) {x | x z A,z Z} Với vector a bất kỳ , khi tịnh tiến a theo từng vectơ trong Z. Kết quả sau khi tịnh tiến a là các vectơ mới: a1,a2, an (n là số lượng phần tử trong Z). Nếu tất cả các vector a1,a2, an đều thuộc A a sẽ thuộc tập E(A,Z). Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 18
- Một số định nghĩa khác Tập tất cả các pixel p, I Z p S Z p I p trong ảnh sao cho Z+p chứa trong I. Phần giao của các copies của reflected. I Z Z p . SE, mỗi cái dịch chuyển đến vị trí pixel p supp I thuộc foreground của ảnh. I Z I p . Phần giao của các copies của mỗi bản p supp Z sao - là dịch chuyển đến vị trí pixel thuộc goreground của refl. SE. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 19
- A = {(1,1),(2,1),(3,1),(4,1),(1,2),(3,2),(1,3),(2,3)} Z= {(1,0),(0,1)} Lấy Vector (1,1) tịnh tiến theo các vector trong Z a1 = (1,1) + (1,0) = (2,1) thuộc A a2 = (1,1) + (0,1) = (1,2) thuộc A (1,1) thuộc E(A,B). Lấy Vecto (3,1) tịnh tiến theo các vector trong Z a1 = (3,1) + (1,0) = (4,1) thuộc A a2 = (3,1) + (0,1) = (3,2) thuộc A (3,1) thuộc E(A,Z). Lấy Vecto (2,2) tịnh tiến theo các vector trong Z a1 = (2,2) + (1,0) = (3,2) thuộc A E(A,Z) = {(1,1),(3,1),(2,2)} a2 = (2,2) + (0,1) = (2,3) thuộc A (2,2) thuộc E(A,Z). Các vectơ khác khi tịnh tiến sẽ có vectơ không thuộc A. Không thuộc E(A,Z) 20
- Xét pixel của ảnh, nếu mọi phần tử của SE trùng với phần ảnh tương ứng, thì logical disjunction (OR operation) được thực hiện giữa tâm của SE với pixel tương ứng để tạo ra pixel trong ảnh output . SE Ảnh nguồn Eroded Các đối tượng nhỏ hơn SE sẽ bị xóa, các đối tượng nối với nhau bởi đường mảnh sẽ tách rời và kích cỡ đối tượng sẽ giảm Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 21
- Tham khảo mã nguồn openCV_Erode Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 22
- Dilate: nhằm tăng kích thước, bịt các lỗ hổng . Có tính giao hoán A(A, B) A B B A D(B, A) . Có tính kết hợp (A B) C A(B C) Erode : co kích thước, mở rộng khoảng hở . Không có tính giao hoán . Không có tính kết hợp Dilate và erode có tính đối ngẫu Có thể cài đặt erode từ dilate và ngược lại c c c c I Z I Z và I Z I Z . Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 23
- Hai phép toán phát triển và bào mòn thường được sử dụng cùng nhau. Từ nhu cầu đó, người ta kết hợp 2 phép toán này để tạo nên một số phép toán có mức độ quan trọng cao hơn: phép mở và phép đóng. Phép mở của tập A bởi phần tử cấu trúc B được định nghĩa như sau: Phép đóng của tập A bởi phần tử cấu trúc B : Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 24
- Opening Closing Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 25
- Phép mở Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 26
- Phép đóng Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 27
- Kết quả của phép mở và phép đóng : Kết quả Phép mở Kết quả Phép đóng Ứng dụng Phép mở có thể sử dụng để loại bỏ các cầu nối, các cành hoặc phần nhô ra của ảnh. Phép đóng có thể sử dụng để lấp đầy các lỗ hổng, các khe hở nhỏ. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 28
- “L” shaped SE Foreground: white pixels Cross-hatched pixels are O marks origin Background: black pixels indeterminate. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 29
- Một số tính chất của phép mở: Tính phản mở rộng : Tính đơn điệu tăng : Tính dừng : Một số tính chất của phép đóng : Tính mở rộng : Tính đơn điệu tăng : Tính dừng : Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 30
- Tham khảo mã nguồn openCV_OpenClose Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 31
- Hit-and-miss được dùng để xác định patterns của foreground và background pixels trong ảnh Xét tập A và một tập hợp B gồm 2 phần tử (B1,B2) với B1 ∩ B2 = Ø, định nghĩa phép biến đổi trúng-trật với công thức như sau (công thức 1) : Một định nghĩa khác về phép biến đổi trúng-trật được xác định bởi (công thức 2): Đụng B1 nhưng không đụng B2. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 32
- Xét một ví dụ ứng dụng của phép biến đổi trúng- trật: Xác định vị trí các điểm góc trên bên phải của tập A Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 33
- Theo công thức 1, kết quả biến đổi là Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 34
- Theo công thức 2, kết quả là Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 35
- Top hat được xác định ℎ 푡 = − ∘ với B là SE Bottom-hat được xác định bởi ℎ 푡 = ∙ − với B là SE Sinh viên làm thực nghiệm với ảnh sau và cho nhận xét (đọc tài liệu Gonzalez từ 695-700) Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 36
- binary image 8-connected SE 4-conn inside bdry original erosion by square difference Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 37
- Gọi biên của đối tượng A là ∂A, lúc đó công thức để tính ∂A như sau: ∂A = A – (A Θ B) B A Θ B ∂A A Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 38
- Ảnh nhị phân 4-connected SE 8-conn inside bdry. original erosion by plus Ảnh hiệu - Difference Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 39
- binary image 4-connected SE 8-conn outside bdry. original dilation by plus difference Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 40
- binary image 8-connected SE 4-conn outside bdry original dilation by square difference Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 41
- binary image erosion by square is 8-conn inside bdry is original in erosion by plus in 4-conn inside bdry Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 42
- binary image dilation by plus is 8-conn outside bdry is original in dilation by square in 4-conn outside bdry Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 43
- inside boundaries are disjoint from outside boundaries 8-bdry/4-bdry/orig all 4 boundaries orig/8-bdry/4-bdry Hình ảnh minh họa trích từ nguồn Vanderbilt University Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 44
- Đa giác A được gọi là đa giác lồi (convex) nếu một đường thẳng nối 2 điểm bất kì thuộc A phải nằm trong A Bao lồi (convex hull) H của một tập bất kì S là một đa giác lồi nhỏ nhất chứa tất cả các điểm thuộc S Tập được tạo thành bởi H-S được gọi là tập thiếu của S Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 45
- Cho Bi, i=1,2,3,4 biểu diễn 4 thành phần cấu trúc Quy trình bao gồm việc thực hiện các phương trình X i (X Bi ) A, i 1,2,3,4; k 1,2,3, k k 1 X i A 0 i i Đặt D X conv , khi bước thực hiện phương trình i i trên đạt được X k X k 1 Bao lồi của A là 4 i C(A) D i 1 Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 46
- 1 2 3 4 Đặt X 0 X 0 X 0 X 0 A Thực hiện với B1 . Thực hiện biến đổi Hit or Miss tập A với B1 cho đến khi không biến đổi được nữa, ta được tập đích Xn. 1 . Hợp Xn với A, gọi là D Thực hiện tương tự với B2, B3, và B4; thu được D2, D3 và D4 Hợp của D1, D2, D3, D4 là bao lồi tìm theo thuật toán này Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 47
- B1 B2 B3 B4 x Don’t care x x x x x x x x x x Background x x x x x x x x x x Foreground A Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 48
- x x x x x Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 49
- 50 x x x x x TS. Ngô Quốc Việt Quốc TS.Ngô - Bài giảng Xử lý ảnh Xử giảngBài
- 51 Việt Quốc Ngô TS. - ảnh lý Xử giảng Bài x x x x x
- Bài Bài giảng Xử lý ảnh - TS. TS. Ngô Quốc Việt x x x x x 52
- D1 D3 4 C(A) Di i 1 D2 D4 Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 53
- Tham khảo mã nguồn openCV_ConvexHull Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 54
- Làm mỏng một tập A theo một phần tử cấu trúc B được xác định: Hit-or-miss transformation Có thể dùng nhiều SE để làm mảnh A Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 55
- Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 56
- Quá trình làm dày A theo phần tử cấu trúc B được xác định Hit-or-miss transformation Quá trình làm dày cũng được xác định qua một chuỗi phần tử cấu trúc. SE của thickening tương tự như SE của thinning nhưng hoán vị vai trò của zero và 1 Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 57
- Tham khảo mã nguồn openCV_ConvexHull Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 58
- Đặt Skel(I,r) là tập các pixel trong I sao cho nếu p Skel(I,r) thì Dp(r), là maximal disk bán kính r trong I. Nghĩa là, Skel(I,r) là quỹ đạo của các tâm của maximal disks bán kính r trong I. Thì S Skel(I,r) r 0 Xương của I là hội của tập tất cả tâm của maximal disk. Maximal disk Dp(r) tiếp xúc với biên của I ít nhất tại hai điểm. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 59
- Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 60
- Ảnh gốc skeleton thô (đỏ) Tỉa nhánh và connected Là tập những tâm skeleton của maximal disks. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 61
- non maximal “disks” maximal disks (red) non max & max disks “disks” are squares non max disks (blue) over skeleton Maximal disk tại pixel p là đĩa lớn nhất fg có chứa p. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 62
- non maximal “disks” maximal disks (red) non max & max disks “disks” là các square non max disks (blue) over skeleton Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 63
- • • • • spurious and disconnected result Note • SE n n n n = = = = = has 3 2 1 0 : : : : Z that is SE SE SE SE 8 = points = = = = the Sq( Sq( Sq( 1 = pixel Sq( 3 7 5 ) ) ) . 3 ) union of all 4 skeleton original Bài Bài giảng Xử lý ảnh to the right top - middle open above erode n=0 w/ Z (n=1) - TS. TS. Ngô Quốc Việt Skel( I, 0 ) top - middle open above erode n=1 w/ Z (n=1) Skel( I, 1 ) top - middle open above erode n=2 w/ Z (n=1) Skel( I, 2 ) top - middle open above erode n=3 w/ Z (n=1) Skel( I, 3 ) 64
- Skeleton thô Xem các pixel giả là Tỉa nhánh Có các pixel “giả” Thành phần liên thông Tốt hơn Ít hơn 3 pix. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 65
- 2 components 2 other components Giao của pruned skel. dilated. of pruned skel dilated. Các thành phần dilated. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 66
- Skeleton thô Skeleton tỉa nhánh reconnected skeleton Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 67
- Một trong những cách biểu diển Skeleton của một tập hợp là biểu diển qua openning và erosion. Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 68
- Quá trình rút xương có thể được biểu diển qua các các thao tác xói mòn (erosion) và mở (opening). Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 69
- Tham khảo mã nguồn openCV_Skeleton Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 70
- Bài giảng Xử lý ảnh - TS. Ngô Quốc Việt 71