Xử lý ảnh - Lọc ảnh ppsx

1 Đ4. Kỹ thuật lọc ảnh Quá trình số hoá ảnh thực có thể gây nhiễu cho ảnh số. Nhiễu trong ảnh số đa dạng và phức tạp. Trong phần này ta sẽ hiểu nhiễu là những phần tử ảnh mà giá trị của nó trội so với các phần tử xung quanh. Xét theo tần số, các phần tử nhiễu có tần số cao so với các điểm xung quanh. Tuy nhiên, có thể thấy rõ là rất khó phân biệt điểm đặc biệt của ảnh với nhiễu. Ví dụ, vết nhăn trên quần áo và vết gập của ảnh có trớc khi số hoá. Cơ sở để xử lý điểm ảnh bị nhiễu là dùng thông tin của các điểm ảnh lân cận. Thay thế giá trị điểm ảnh cần xử lý bằng tổ hợp giá trị các điểm ảnh lân cận. Việc thay thế dựa trên giả định là các điểm ảnh lân cận có giá trị gần giống với điểm ảnh cần xét. Về kỹ thuật, ngời ta sử dụng một mặt nạ di chuyển khắp ảnh gốc. ở mỗi vị trí của ảnh, tính toán theo các điểm trong mặt nạ và thay vào giá trị cho điểm ở tâm mặt nạ. Theo các loại mặt nạ khác nhau mà có các cách tính khác nhau, tổ hợp giá trị các điểm lân cận điểm đợc xét 1. Lọc tuyến tính Công thức tính giá trị mới cho điểm ảnh tại (x,y): G(x,y) = W 1 M 1 + W 2 M 2 + W 3 M 3 + W 4 M 4 + W 5 M 5 + W 6 M 6 + W 7 M 7 + W 8 M 8 +W 9 M 9 trong đó W k , M k là các ma trận 98 7 654 321 WWW WWW WWW W 98 7 654 321 MMM MM MMM MM Đặt mặt nạ sao cho tâm của mặt nạ, M 5 trùng với điểm (x,y), khi này M 1 trùng với điểm (x-1, y-1), M 9 trùng với điểm (x+1,y+1). Khi này W 5 = f(x,y); W 1 = f(x-1,y-1), )1,1()1,()1,1( ),1(),(),1( )1,1()1,()1,1( 98 7 654 321 yxfyxfyxf yxfyxfyxf yxfyxfyxf WWW WWW WWW W Một số mặt nạ thờng đợc sử dụng: 2 111 111 111 9 1 1 H 111 121 111 10 1 2 H 121 242 121 16 1 3 H Kỹ thuật thực hiện: 7 7 7 7 7 7 7 7 7 6 Tính g(6,4) theo H1: Các giá trị trong cửa sổ khi đặt tâm cửa sổ vào điểm (6,4), hàng 4, cột 6: 7 7 7 7 7 1 7 7 7 G(4,6) = [(7+7+7+7+1+1+7+7+7)/9] = [51/9 ] = [5.66 ] = 6 6 6 6 6 7 7 7 7 7 7 5 5 5 6 7 7 7 7 7 7 5 5 5 6 7 1 1 7 7 7 4 4 5 6 7 7 7 7 7 7 4 4 5 6 6 7 7 7 6 6 3 3 4 5 6 7 7 7 6 6 4 4 5 6 7 7 0 7 6 6 4 4 5 6 7 7 7 7 6 6 4 4 5 6 7 7 7 7 6 6 ảnh gốc I 7 7 7 7 7 7 7 7 7 6 7 7 7 7 7 7 7 7 7 7 6 6 6 6 7 7 7 7 7 7 6 6 6 6 7 7 7 7 7 7 5 5 5 6 7 7 7 7 7 7 5 5 6 6 7 6 6 6 7 7 5 5 5 6 7 7 1 7 7 7 5 5 5 6 7 6 6 6 7 7 4 4 5 6 7 7 7 7 7 7 4 5 5 6 7 6 6 6 7 7 4 4 5 6 6 7 7 7 6 6 4 4 5 6 6 7 7 7 7 6 3 3 4 5 6 7 7 7 6 6 4 4 5 6 6 6 6 6 6 6 4 4 5 6 7 7 0 7 6 6 4 4 5 6 6 6 6 6 6 6 4 4 5 6 7 7 7 7 6 6 4 4 5 6 7 6 6 6 6 6 4 4 5 6 7 7 7 7 6 6 4 4 5 6 7 7 7 7 6 6 ảnh gốc I ảnh I qua mặt nạ H1 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 6 6 6 6 6 7 7 7 7 7 7 6 6 6 6 7 7 7 7 7 7 5 5 6 6 7 6 6 6 7 7 5 5 5 6 7 7 6 7 7 7 5 5 5 6 7 6 6 6 7 7 5 5 5 6 7 6 6 6 7 7 4 5 5 6 7 6 6 6 7 7 4 4 5 6 7 7 6 7 7 7 4 4 5 6 6 7 7 7 7 6 4 4 5 6 6 7 7 7 6 6 4 4 5 6 6 6 6 6 6 6 4 4 5 5 6 6 6 6 6 6 4 4 5 6 7 6 6 6 6 6 4 4 5 6 7 6 5 6 6 6 4 4 5 6 7 6 6 6 6 6 4 4 5 6 7 7 6 6 6 6 3 4 4 5 6 7 7 7 7 6 6 4 4 5 6 7 7 7 7 6 6 ảnh I qua mặt nạ H2 ảnh I qua mặt nạ H3 2. Lọc phi tuyến Kỹ thuật lọc phi tuyến coi một điểm ảnh kết quả không phải là tổ hợp tuyến tính của các điểm lân cận. Bộ lọc phi tuyến thờng dùng là lọc trung vị. A) B) Hình A) ảnh gốc; B) ảnh qua lọc H 1 Hình C) ảnh lọc qua H 2 D) ảnh lọc qua H 3 4 a. Lọc trung vị Khái niệm Phần tử b của dãy có 2m+1 phần tử x 1 , x 2 , , x m , x m+1 , x m+2 , , x 2m , x 2m+1 (1) đợc gọi là trung vị nếu có m phần tử lớn hơn hoặc bằng b và m phần tử nhỏ hơn hoặc bằng b. Ví dụ 1: Dãy {1, 2, 9, 5, 4} có trung vị bằng 4, vì có 2 phần tử { 1, 2 , 9, 5, 4} bé hơn và 2 phần tử { 1, 2 , 9, 5 , 4} lớn hơn. Ví dụ 2: Dãy {15, 17,18, 16, 78, 17, 17, 15, 20} có trung vị bằng 17, vì có 4 phần tử (viết nhỏ hơn và gạch chân) bé hơn hoặc bằng { 15 , 17 ,18, 16 , 78, 17, 17, 15 , 20} và 4 phần tử (đợc đánh dấu bằng màu đỏ) lớn hơn hoặc bằng: { 15 , 17 , 18 , 16 , 78 , 17 , 17, 15 , 20 } Thuật toán: Di chuyển một cửa sổ 3 3 trên mặt phẳng ảnh, từ trái qua phải, từ trên xuống dới. Mỗi lần di chuyển một cột từ trái qua phải, hết một hàng, di chuyển một hàng từ trên xuống dới. Tại mỗi vị trí, xét 9 phần tử ảnh trong cửa sổ, thay thế phần tử ở giữa bằng trung vị của 9 phần tử trong cửa sổ. Ngoài cửa sổ 3 3 có thể sử dụng cửa sổ kích thớc 5 5 hoặc 7 7. Thủ tục lọc có thể đợc thực hiện nhiều lần, cho đến khi thủ tục lọc không làm thay đổi ảnh. Thông thờng, cửa sổ 3 3 hay đợc lựa chọn để sử dụng. Ví dụ 3: Cột x 15 17 18 Hàng y 16 78 17 17 15 20 5 Đặt cửa sổ sao cho tâm trùng với điểm (x,y). Giá trị của phần tử ảnh tại (x,y) là 78. Các phần tử ảnh trong cửa sổ có giá trị lập thành dãy 15, 17, 18, 16, 78, 17, 17, 15, 20 Theo ví dụ trên, trung vị của dãy này là 17. Giá trị mới của phần tử ảnh tại điểm (x,y) sẽ là 17. Cột x 15 17 18 Hàng y 16 17 17 17 15 20 Kỹ thuật thực hiện: a) Sắp xếp thành một dãy X tăng dần hoặc giảm dần các giá trị lấy đợc trong cửa sổ: 15 15 16 17 17 17 18 20 78 1 2 3 4 5 6 7 8 9 b) Đặt lại f(x,y) := X(5) Lu ý rằng ảnh mới qua xử lý khác với ảnh gốc. b. Lọc giả trung vị Với lọc trung vị, số lợng tính toán khá lớn (có thể bằng số mũ của kích thớc cửa sổ lọc). Vì vậy, để khắc phục nhợc điểm này, ngời ta dùng một phơng pháp khác lọc giả trung vị (Pseudo Hình E) ảnh gốc F) ảnh qua lọc trung vị 6 Median Filter). Thí dụ với dãy số: a, b, c, d, e lọc giả trung vị đợc định nghĩa nh sau: Giả trung vị(a,b,c,d,e)= ),,(),,,(),,,( ),,(),,,(),,,( 2 1 edcMaxdcbMaxcbaMaxMIN edcMindcbMincbaMinMAX Thuật toán: 1. Lấy các phần tử trong cửa sổ ra mảng một chiều (L phần tử). 2. Tìm min của lần lợt các chuỗi con rồi lấy max: gọi m1 là giá trị này. 3. Tìm max của lần lợt các chuỗi con rồi lấy min: gọi m2 là giá trị tìm đợc. 4. Gán giá trị điểm đang xét là trung bình cộng của m1 và m2. Lọc giả trung vị có nhiều điểm giống nh lọc trung vị. Dãy lấy ra không cần sắp xếp và giá trị gọi là trung vị lại đợc tính theo trung bình cộng của max của min và min của max. Hai loại mặt nạ hay dùng là mặt nạ vuông và mặt nạ chữ thập. Thực tế, mặt nạ vuông không làm biến dạng ảnh mà lại hiệu quả. Tuy nhiên trong lọc giả trung vị, mặt nạ chữ thập cho kết quả khả quan hơn nhiều. a) Mặt nạ chữ thập b) Mặt nạ vuông 5 x5 Các kỹ thuật lọc trình bày trên là lọc thông thấp. Nó đợc dùng để lọc nhiễu. Ngoài lọc thông thấp, ngời ta còn sử dụng lọc thông cao. Lọc thông cao dùng để làm nổi bật các chi tiết có tần số thấp. Các phần tử có tần số không gian cao sẽ sáng hơn, còn các phần tử có tần số không gian thấp sẽ đen đi. Kỹ thuật lọc thông cao cũng đợc thực hiện nhờ thao tác nhân chập. Các mặt nạ hay đợc dùng nh: 111 191 111 010 151 010 121 252 121 7 Bài 2. Nâng cao chất lợng ảnh. Nâng cao chất lợng ảnh là một bớc quan trọng, tạo tiền đề cho xử lý ảnh. Mục đích chính là làm nổi bật một số đặc tính của ảnh nh thay đổi độ tơng phản, lọc nhiễu, nổi biên, làm trơn ảnh, khuếch đại ảnh, I. Cải thiện ảnh dùng toán tử điểm Toán tử điểm là toán tử không bộ nhớ, ở đó một mức xám u [0,N] đợc ánh xạ sang một mức xám v [0,N]: v=f(u). ánh xạ f tuỳ theo các ứng dụng khác nhau có các dạng khác nhau: Tăng độ tơng phản. Tách nhiễu phân ngỡng. Biến đổi âm bản. Cắt theo mức. Trích chọn bít. Trừ ảnh. Nén dải độ sáng. 1. Tăng độ tơng phản. Độ tơng phản biểu diễn sự thay đổi độ sáng của đối tợng so với nền, hay nói cách khác, độ tơng phản là độ nổi của điểm ảnh hay vùng ảnh so với nền. ảnh số là tập hợp các điểm, mà mỗi điểm có giá trị độ sáng khác nhau. ở đây, độ sáng để mắt ngời dễ cảm nhận ảnh song không phải là quyết định. Thực tế chỉ ra rằng hai đối tợng có cùng độ sáng nhng đặt trên hai nền khác nhau sẽ cho cảm nhận khác nhau. Vì vậy, ta có thể thay đổi độ tơng phản của ảnh sao cho phù hợp. Để điều chỉnh độ tơng phản của ảnh, ta điều chỉnh biên độ của ảnh trên toàn dải hay trên dải có giới hạn bằng cách biến đổi tuyến tính biên độ đầu vào (dùng hàm biến đổi là hàm tuyến tính) hay phi tuyến (hàm mũ hay hàm logarit). Khi dùng hàm tuyến tính các độ dốc , , phải chọn lớn hơn 1 trong miền cần giãn. Các cận a, b có thể chọn khi xem xét lợc đồ xám của ảnh. v b v a a b L u v 8 Lubvbu buavau auu uf b a )( )( 0 )( Chú ý: Nếu giãn độ tơng phải bằng hàm tuyến tính ta có: = = ảnh kết quả trùng với ảnh gốc. , , >1 giãn độ tơng phản. , , <1 co độ tơng phản. Hàm mũ hay dùng giãn độ tơng phản có dạng F = (X[m,n]) p ; p thờng chọn bằng 2 2. Tách nhiễu và phân ngỡng. ubL buau au uf 00 )( Tách nhiễu là trờng hợp đặc biệt của giãn độ tơng phản khi hệ số góc ==0. Tách nhiễu đợc ứng dụng trong trờng hợp khi biết tín hiệu nằm trên khoảng [a, b]. Phân ngỡng là trờng hợp đặc biệt của tách nhiễu khi a=b=const. Trong trờng hợp này ảnh đầu ra là ảnh nhị phân (vì chỉ có hai mức). Phân ngỡng hay dùng trong kỹ thuật in ảnh hai màu vì ảnh gần nhị phân không thể cho ra ảnh nhị phân khi quét ảnh bởi có sự xuất hiện của nhiễu do bộ cảm biến và sự biến đổi của nền. Thídụ nh trờng hợp ảnh vân tay. 3. Biến đổi âm bản. f(u) = L u Biến đổi âm bản rất có ích khi hiện các ảnh Y học. 4. Cắt theo mức. Kỹ thuật này dùng hai phép ánh xạ khác nhau cho trờng hợp có nền và không nền. u v a b L u v 9 Có nền diKhacu buaL uf )( Không nền diKhac buaL uf 0 )( Biến đổi này cho phép phân đoạn một số mức xám từ phần còn lại của ảnh. Nó hữu dụng khi nhiều đặc tính khác nhau của ảnh nằm trên nhiều miền mức xám khác nhau. 5. Trích chọn bít. Mỗi điểm ảnh đợc mã hoá trên B bít. Nếu B=8 có ảnh 2 8 =256 mức xám (ảnh nhị phân ứng với B=1). Trong các bít mã hoá này, ngời ta chia làm hai loại: bít bậc thấp và bít bậc cao. Với bít bậc cao, độ bảo toàn thông tin cao hơn nhiều so với bit bậc thấp. U = k 1 2 B-1 + k 2 2 B-2 + + k B-1 2 + k B Nếu ta muốn trích chọn bit thứ n và hiện chúng, ta dùng biến đổi: diKhac kneuL uf n 0 1 )( 6. Trừ ảnh Trừ ảnh dùng để tách nhiễu khỏi nền. Ta quan sát ảnh ở hai thời điểm khác nhau, so sánh chúng để tìm ra sự khác nhau, sau đó dóng thẳng hai hai ảnh rồi trừ đi và thu đợc ảnh mới. ảnh mới đó chính là sự khác nhau. Kỹ thuật này hay đợc dùng trong dự báo thời tiết hay trong y học. 7. Nén dải độ sáng Đôi khi do dải động của ảnh lớn, việc quan sát không thuận tiện, cần phải thu nhỏ dải độ sáng lại mà ngời ta gọi là nén dải độ sáng. Ngời ta dùng phép biến đổi Loga sau: v(m,n)=c log 10 (+u(m,n)) Với c là hằng số tỷ lệ, rất nhỏ so với u(m,n). Thờng chọn cỡ 10 -3 . II. Cải thiện ảnh dùng toán tử không gian a b u v 10 Cải thiện ảnh là làm cho ảnh có chất lợng tốt hơn theo ý đồ sử dụng. Thờng là ảnh thu nhận có nhiễu cần phải loại bỏ nhiễu hay ảnh không sắc nét bị mờ hoặc cần làm rõ các chi tiết nh biên. Các toán tử không gian dùng trong kỹ thuật tăng cờng ảnh đợc phân theo nhóm theo công dụng: làm trơn nhiễu, nổi biên. Để làm trơn nhiễu hay tách nhiễu ngời ta sử dụng các bộ lọc tuyến tính (lọc trung bình, thông thấp) hay lọc phi tuyến (trung vị, giả trung vị, lọc đồng hình). Do bản chất của nhiễu ứng với tần số cao và cơ sở lý thuyết của bộ lọc là bộ lọc chỉ cho tín hiệu có tần số nào đó thông qua (dải tần bộ lọc). Do vậy để lọc nhiễu ta dùng lọc thông thấp (theo quan điểm tần số không gian) hay lấy tổ hợp tuyến tính để san bằng (lọc trung bình). Để làm nổi cạnh (ứng với tần số cao), ngời ta dùng các bộ lọc thông cao, Laplace. 1. Làm trơn nhiễu bằng lọc tuyến tính, lọc trung bình và lọc dải thông thấp a) Lọc trung bình không gian Mỗi điểm ảnh đợc thay thế bằng lọc trung bình trọng số của các điểm lân cận: wlk lnkmylkanmv ),( ),(),(),( Nếu ta dùng các trong số nh nhau, công thức trở thành. wlk w lnkmy N nmv ),( ),( 1 ),( Với w kl N a 1 ; N w : là số điểm ảnh trong cửa sổ lọc. y(m,n) : ảnh đầu vào. v(m,n) : ảnh đầu ra. w(m,n): cử sổ lọc a(m,n) : trọng số lọc Lọc trung bình có trọng số chính là thực hiện chập ảnh đầu vào với nhân chập H. Nhân chập H trong trờng hợp này có dạng: 111 111 111 9 1 H [...]... hình thực hiện lấy logarit của ảnh quan sát log(X(m,n)) = log(X (m,n))+log((m,n)) Rõ ràng nhiễu nhân có trong ảnh sẽ bị giảm 2 Làm trơn nhiễu bằng lọc phi tuyến: Lọc trung vị Lọc giả trung vị Lọc ngoài a) Lọc trung vị Nhắc lại rằng khái niệm trung vị đã nêu trong bài 1 và được viết: v(m,n)= Trung vi(y(m-k,n-l) với (k,l) W ( 4-8 ) Kỹ thuật này đòi hỏi giá trị các điểm ảnh trong cửa sổ phải xếp theo... tính của lọc trung bình ta có: 11 Y ( m, n ) 1 Nw x Hay: Y (m, n) goc (m k , n l ) [m, n] ( k , l )w 1 Nw x goc (m k , n l ) ( k ,l )w 2 n Nw Như vậy nhiễu cộng trong ảnh giảm đi Nw lần c) Lọc đồng hình Kỹ thuật lọc này hiệu quả đối với nhiễu nhân Thực tế ảnh quan sát được gồm ảnh gốc nhân với một hệ số nhiễu X(m,n) = X (m,n).(m,n) X (m,n) : ảnh gốc X (m,n) : ảnh thu được (m,n) : nhiễu Lọc đồng... biên> Tính chất của lọc trung vị: - Lọc trung vị là phi tuyến vì: 12 Trung vị ((x(m) + y (m))# Trung vị(x(m)) + Trung vị(y(m)) - Hữu ích cho việc loại bỏ các điểm các điểm ảnh hay các hàng mà vẫn bảo toàn độ phân giải - Hiệu quả giảm khi số điểm nhiễu trong cửa sổ lớn hơn hay bằng một nửa số điểm trong cửa sổ Điều này dễ giải thích vì trung vị là (Nw+1)/2 giá trị lớn nhất nếu Nw lẻ Lọc trung vị cho trường...Trong lọc trung bình đôi khi người ta ưu tiên cho các hướng để bảo vệ biên của ảnh khỏi bị mờ đi do làm trơn ảnh Ví dụ: Giả sử đầu vào biểu diễn bởi ma trận I: 4 5 I 6 5 5 7 2 7 1 7 1 7 1 6 1 8 3 7 5 7 1 7 6 1 2 ảnh số thu được bởi lọc trung bình Y=HI có dạng: 26 35 1 Y 36 9 36 24 36 31 19 16 39 46 31 27 43 49 34 27 43 48 34 12 35 33 22 11 Một bộ lọc trung bình không... giảm dần so với giá trị trung vị Kích thước cửa sổ thường được chọn sao cho số điểm ảnh trong cửa sổ là số lẻ Các cửa sổ hay dùng là cửa sổ 3 x3, 5x5, hay 7x7 Thí dụ: Nếu y(m)={2,3,8,4,2} và cửa sổ W= (-1 ,0,1), ảnh kết quả thu được sau lọc trung vị sẽ là v(m)= (2,3,4,4,2) Thực vậy: mỗi lần ta so sánh một dãy 3 điểm ảnh đầu vào với trung vị, không kể điểm biên Do đó: v[0]=2 v[1]= Trung... 1]} 2 4 H là ma trận 2*2 Mỗi điểm ảnh kết quả có giá trị bằng trung bình cộng của nó với 4 điểm lân cận b) Lọc thông thấp: Lọc thông thấp thường sử dụng để làm trơn nhiễu Sử dụng một số nhân chập sau: 0 1 0 1 H 1 2 1 8 0 1 0 1 b 1 b b 2 Hb 2 (b 2) 1 b 1 b 1 Khi b=1 Hb chính là nhân chập H1 (lọc trung bình) Để hiểu rõ hơn, ta viết lại phương trình thu nhân ảnh dưới dạng: Xqs[m,n] = Xgoc[m,n]... trường hợp 2 chiều coi như lọc trung vị tách được theo từng chiều, có nghĩa là người ta tiến hành lọc trung vị cho cột tiếp theo cho hàng b) Lọc ngoài (Outlier Filter) Giả thiết rằng có một mức ngưỡng nào đó cho các mức nhiễu (có thể dựa vào lược đồ xám) Tiến hành so sánh giá trị của một điểm ảnh với trung bình số học 8 lân cận của nó Nếu sự sai lệch này lớn hơn ngưỡng, điểm ảnh này được coi như nhiễu... giá trị của điểm ảnh bằng giá trị trung bình 8 lân cận vừa tính được Các cửa sổ tính toán thường là 3x3 Tuy nhiên cửa sổ có thể mở rộng đến 5x5 hay 7x7 để đảm bảo tính tương quan giữa các điểm ảnh Vấn đề quan trọng là xác định ngưỡng để loại nhiễu mà vẫn không làm mất thông tin 1 Mở và đọc Header của ảnh I Ghi Header vào tệp mới MI 2 Chép nguyên hàng đầu từ I sang MI 3 For i:=2 to BMP.H-1 a) Chép điểm... làm mất thông tin 1 Mở và đọc Header của ảnh I Ghi Header vào tệp mới MI 2 Chép nguyên hàng đầu từ I sang MI 3 For i:=2 to BMP.H-1 a) Chép điểm đầu của hàng từ I sang MI b) For j:=2 toBMP.W-1 Đọc từ tệp: cửa sổ (j-1,i-1,j+1,i+1) vào W Tính G(W); Ghi G(W) vào tệp MI; c) Chép điểm cuối của hàng từ tệp I sang MI 4 Chép nguyên hàng cuối cùng từ tệp I sang tệp MI 13 . phi tuyến thờng dùng là lọc trung vị. A) B) Hình A) ảnh gốc; B) ảnh qua lọc H 1 Hình C) ảnh lọc qua H 2 D) ảnh lọc qua H 3 4 a. Lọc trung vị Khái niệm. N w : là số điểm ảnh trong cửa sổ lọc. y(m,n) : ảnh đầu vào. v(m,n) : ảnh đầu ra. w(m,n): cử sổ lọc a(m,n) : trọng số lọc Lọc trung bình có trọng số chính là thực hiện chập ảnh đầu vào với. f(x,y) := X(5) Lu ý rằng ảnh mới qua xử lý khác với ảnh gốc. b. Lọc giả trung vị Với lọc trung vị, số lợng tính toán khá lớn (có thể bằng số mũ của kích thớc cửa sổ lọc) . Vì vậy, để khắc phục