Bộ lọc hình thái khác biệt như thế nào so với bộ lọc tuyến tính truyền thống?

Bộ lọc tuyến tính truyền thống dựa trên phép tích chập số học và lọc tín hiệu theo miền tần số phổ, dễ làm nhòe các biên cạnh đối tượng. Ngược lại, bộ lọc hình thái là toán tử phi tuyến dựa trên lý thuyết tập hợp và đại số lưới, lọc bỏ nhiễu và biến đổi tín hiệu theo hình thái hình học mà vẫn bảo toàn biên cạnh sắc nét.

Vì sao hai phép toán mở và đóng lại giữ vai trò trụ cột trong lọc hình thái?

Các phép co và giãn đơn thuần làm thay đổi kích thước của đối tượng (thu nhỏ hoặc phóng to). Phép mở (co trước giãn sau) và phép đóng (giãn trước co sau) là các toán tử lũy đẳng và tăng, cho phép loại bỏ nhiễu sáng hoặc tối mà vẫn duy trì kích thước tổng thể của các cấu trúc chính.

Tính chất lũy đẳng (idempotence) của bộ lọc hình thái có ý nghĩa kỹ thuật gì?

Tính lũy đẳng đảm bảo rằng việc áp dụng bộ lọc lặp lại lần thứ hai không làm thay đổi thêm tín hiệu đã lọc. Điều này giúp bộ lọc loại bỏ dứt điểm các thành phần không mong muốn ngay trong lần quét đầu tiên mà không làm suy hao tín hiệu ở các bước xử lý kế tiếp.

Bộ lọc hình thái: Cơ sở toán học, toán tử và ứng dụng

Bộ lọc hình thái (tiếng Anh: morphological filter) là một lớp các toán tử xử lý tín hiệu và hình ảnh phi tuyến tính thỏa mãn hai tiên đề đại số về tính tăng và tính lũy đẳng trên cấu trúc lưới, trong đó việc trích xuất và biến đổi hình học của tín hiệu được thực hiện thông qua sự tương tác với một hình mẫu định trước gọi là phần tử cấu trúc (structuring element). Trong khi các phép toán sơ cấp như phép co và phép giãn làm thay đổi quy mô hình học, các bộ lọc hình thái thực thụ như phép mở, phép đóng và bộ lọc xen kẽ tuần tự cho phép loại bỏ nhiễu và bảo toàn biên cạnh mà không làm dịch chuyển ranh giới đối tượng.

Bối cảnh ra đời và cơ sở đại số của hình thái học toán học

Hình thái học toán học ban đầu được xây dựng trên lý thuyết tập hợp cho các ảnh nhị phân. Năm 1986, Jean Serra đã công bố một bài báo tổng quan nền tảng trên tạp chí Computer Vision, Graphics, and Image Processing, hệ thống hóa hình thái học toán học từ không gian tập hợp Euclid sang lý thuyết lưới đầy đủ (complete lattices), qua đó mở rộng việc áp dụng các phép toán hình thái từ ảnh nhị phân sang tín hiệu và ảnh mức xám liên tục cũng như rời rạc.

Tiếp nối nền tảng đó, năm 1987, Petros Maragos và Ronald W. Schafer đã công bố công trình nghiên cứu hệ thống gồm hai phần trên tạp chí IEEE Transactions on Acoustics, Speech, and Signal Processing. Trong phần một, các tác giả đã phân tích bản chất lý thuyết tập hợp của các bộ lọc hình thái và thiết lập mối quan hệ tương hỗ với các bộ lọc tuyến tính bất biến theo dịch chuyển. Trong phần hai, công trình làm sáng tỏ liên hệ toán học chặt chẽ giữa bộ lọc hình thái với bộ lọc trung vị, bộ lọc thống kê thứ tự và bộ lọc xếp chồng.

Cũng vào năm 1987, Robert M. Haralick cùng Stanley R. Sternberg và Xinhua Zhuang đã công bố bài báo quan trọng trên IEEE Transactions on Pattern Analysis and Machine Intelligence, chuẩn hóa các thuật toán hình thái học cho thị giác máy tính và phân tích ảnh công nghiệp, tạo tiền đề cho việc ứng dụng rộng rãi các bộ lọc hình thái trong khoa học kỹ thuật hiện đại.

Hai phép toán cơ sở: Phép co và phép giãn

Hệ thống các bộ lọc hình thái được xây dựng từ hai phép toán đại số cơ bản: phép co (erosion) và phép giãn (dilation). Cho tín hiệu đầu vào $f(x)$ và phần tử cấu trúc $b(x)$ xác định trên không gian Euclid:

1. Phép giãn mức xám

Phép giãn của hàm $f$ bởi phần tử cấu trúc $b$ , ký hiệu là $f \oplus b$ , được định nghĩa thông qua phép cộng Minkowski mở rộng và toán tử cận trên:

(f \oplus b)(x) = \sup_{y} \{ f(x - y) + b(y) \}

Trong trường hợp phần tử cấu trúc phẳng có giá trị bằng không trên miền xác định $B$ , phép giãn tương đương với phép tìm giá trị cực đại cục bộ:

(f \oplus b)(x) = \max_{y \in B} \{ f(x - y) \}

Về mặt hình học, phép giãn làm mở rộng các vùng sáng, lấp đầy các thung lũng hoặc lỗ thủng tối có kích thước nhỏ hơn phần tử cấu trúc và kết nối các vùng sáng bị đứt đoạn.

2. Phép co mức xám

Phép co của hàm $f$ bởi phần tử cấu trúc $b$ , ký hiệu là $f \ominus b$ , được định nghĩa thông qua phép trừ Minkowski mở rộng và toán tử cận dưới:

(f \ominus b)(x) = \inf_{y} \{ f(x + y) - b(y) \}

Với phần tử cấu trúc phẳng, phép co trở thành phép tìm giá trị cực tiểu cục bộ:

(f \ominus b)(x) = \min_{y \in B} \{ f(x + y) \}

Phép co làm thu hẹp các vùng sáng, mở rộng các vùng tối và loại bỏ hoàn toàn các cấu trúc sáng cô lập có kích thước nhỏ hơn phần tử cấu trúc.

Hai toán tử lọc cơ bản: Phép mở và phép đóng

Mặc dù phép co và phép giãn loại bỏ được các thành phần nhiễu tương ứng, chúng làm thay đổi kích thước hình học tổng thể của đối tượng và không có tính lũy đẳng, vì vậy chưa phải là bộ lọc hình thái hoàn chỉnh theo nghĩa đại số. Để khắc phục điều này, việc kết hợp liên tiếp hai phép toán cơ sở tạo thành hai toán tử lọc cốt lõi: phép mở (opening) và phép đóng (closing).

1. Phép mở hình thái

Phép mở của tín hiệu $f$ bởi phần tử cấu trúc $b$ , ký hiệu là $f \circ b$ , được định nghĩa là một phép co liên tiếp sau đó bởi một phép giãn cùng phần tử cấu trúc:

f \circ b = (f \ominus b) \oplus b

Về mặt trực quan, phép mở tương đương với việc lăn phần tử cấu trúc phía dưới bề mặt đồ thị tín hiệu. Mọi đỉnh nhọn sáng nhỏ hơn phần tử cấu trúc sẽ bị cắt gọt, trong khi hình dạng tổng thể của các vùng sáng lớn hơn được bảo tồn nguyên vẹn. Nhờ đó, phép mở loại bỏ hiệu quả nhiễu xung sáng dạng muối tiêu mà không làm mờ biên cạnh chuyển mức lớn.

2. Phép đóng hình thái

Phép đóng của tín hiệu $f$ bởi phần tử cấu trúc $b$ , ký hiệu là $f \bullet b$ , được định nghĩa là một phép giãn liên tiếp sau đó bởi một phép co:

f \bullet b = (f \oplus b) \ominus b

Cơ chế hình học của phép đóng tương đương với việc lăn phần tử cấu trúc phía trên bề mặt tín hiệu. Những thung lũng hẹp và khe rãnh tối nhỏ hơn phần tử cấu trúc sẽ được lấp đầy, trong khi độ cao của các đỉnh núi lớn không bị suy giảm. Phép đóng là công cụ chuyên dụng để triệt tiêu nhiễu hạt tối và hàn gắn các vết nứt gãy.

Tính chất đại số tiên đề của bộ lọc hình thái

Theo định nghĩa đại số của Jean Serra và Petros Maragos, một toán tử ánh xạ $\psi$ được công nhận là một bộ lọc hình thái khi và chỉ khi thỏa mãn hai tính chất tiên đề sau đây:

Tính tăng (Increasing property): Nếu $f(x) \le g(x)$ với mọi $x$ , thì $\psi(f)(x) \le \psi(g)(x)$ . Tính chất này bảo đảm rằng bộ lọc luôn tôn trọng mối quan hệ thứ tự sắp xếp của cường độ tín hiệu.
Tính lũy đẳng (Idempotence): Áp dụng bộ lọc hai lần liên tiếp cho kết quả đồng nhất với một lần áp dụng: $\psi(\psi(f)) = \psi(f)$ . Tính lũy đẳng thể hiện trạng thái nghiệm dừng ổn định, nghĩa là tín hiệu sau khi lọc đã đạt trạng thái cân bằng hình học đối với phần tử cấu trúc và không bị suy biến thêm nếu tiếp tục lọc.

Khi phần tử cấu trúc chứa gốc tọa độ, phép mở sở hữu thêm tính phản bành trướng: $f \circ b \le f$ , còn phép đóng sở hữu tính bành trướng: $f \le f \bullet b$ .

Bảng so sánh các toán tử hình thái cơ bản

Toán tử	Ký hiệu toán học	Bản chất biến đổi	Tính bành trướng (khi chứa gốc tọa độ)	Tính lũy đẳng (Bộ lọc đại số)
Phép giãn	Toán tử cộng Minkowski $\oplus$	Cực đại cục bộ	Bành trướng (khi chứa gốc tọa độ)	Không lũy đẳng (Toán tử sơ cấp)
Phép co	Toán tử trừ Minkowski $\ominus$	Cực tiểu cục bộ	Phản bành trướng (khi chứa gốc tọa độ)	Không lũy đẳng (Toán tử sơ cấp)
Phép mở	Toán tử mở $\circ$	Cắt đỉnh sáng cục bộ	Phản bành trướng (với mọi phần tử)	Lũy đẳng (Bộ lọc hình thái)
Phép đóng	Toán tử đóng $\bullet$	Lấp rãnh tối cục bộ	Bành trướng (với mọi phần tử)	Lũy đẳng (Bộ lọc hình thái)

Các cấu trúc lọc nâng cao và biến đổi phái sinh

Từ các toán tử cơ sở, nhiều kiến trúc bộ lọc hình thái phức hợp và toán tử phát hiện biên đã được phát triển:

Bộ lọc xen kẽ tuần tự (Alternating Sequential Filters - ASF): Thay vì lọc trực tiếp bằng một phần tử cấu trúc lớn, ASF áp dụng một chuỗi các phép mở và phép đóng đan xen với phần tử cấu trúc có kích thước tăng dần đều. Cấu trúc này loại bỏ nhiễu từ thang đo mịn đến thô, duy trì hình dạng hình học mượt mà và giảm thiểu hiện tượng dịch chuyển biên.
Biến đổi Top-hat: Phép trừ giữa tín hiệu gốc và tín hiệu đã lọc mở tạo ra biến đổi White Top-hat: $f - (f \circ b)$ , dùng để tách các chi tiết sáng nhỏ trên nền không đồng đều. Ngược lại, biến đổi Black Top-hat: $(f \bullet b) - f$ dùng để trích xuất các thung lũng tối hẹp.
Gradient hình thái (Morphological Gradient): Hiệu số giữa phép giãn và phép co: $G(f) = (f \oplus b) - (f \ominus b)$ . Cần lưu ý rằng gradient hình thái là một toán tử phát hiện biên cạnh phái sinh chứ không phải là bộ lọc hình thái đại số, vì nó không thỏa mãn tính tăng và tính lũy đẳng. Toán tử này tạo ra các đường bao biên khép kín mà không phụ thuộc vào xấp xỉ đạo hàm vi phân tuyến tính.

Ranh giới và hạn chế của bộ lọc hình thái

Mặc dù sở hữu nhiều ưu điểm vượt trội trong bảo toàn biên cạnh, bộ lọc hình thái có những giới hạn kỹ thuật cần được cân nhắc kỹ lưỡng:

Độ nhạy cảm với hình học và kích thước phần tử cấu trúc: Hiệu quả của bộ lọc phụ thuộc quyết định vào sự tương thích hình học giữa phần tử cấu trúc và đặc trưng cần xử lý. Một phần tử cấu trúc có kích thước cố định không thể thích ứng đồng thời với các chi tiết có tỷ lệ biến thiên phức tạp trong cùng một bức ảnh.
Kém hiệu quả trước nhiễu cộng Gauss: Khác với bộ lọc trung bình tuyến tính hoặc bộ lọc Wiener vốn tối ưu hóa khử nhiễu Gauss phân bố đều, toán tử tìm cực trị max/min của hình thái học có thể bị lệch pha hoặc tạo thành các mảng giả tạo nếu tín hiệu bị ô nhiễm bởi nhiễu trắng mật độ dày.
Méo dạng hình học khi phần tử cấu trúc bất đối xứng: Nếu phần tử cấu trúc không có tính đối xứng qua gốc tọa độ hoặc không chứa gốc tọa độ, các phép co giãn sẽ gây ra hiện tượng dịch chuyển vị trí trung tâm của đối tượng xử lý.
Độ phức tạp tính toán với phần tử cấu trúc không phân tách: Khi mở rộng sang không gian ba chiều hoặc sử dụng các phần tử cấu trúc phi phẳng kích thước lớn, chi phí tính toán tăng nhanh trừ khi phần tử cấu trúc có thể phân tách thành tích Minkowski của các phần tử nhỏ hơn.

Tóm lại, bộ lọc hình thái đại diện cho một bước chuyển đổi mô thức từ giải tích tần số Fourier sang phân tích hình học không gian trong xử lý tín hiệu. Hiểu rõ các tiên đề đại số cùng các giới hạn vật lý của chúng là cơ sở để kết hợp linh hoạt giữa lọc hình thái và các kỹ thuật xử lý hình ảnh hiện đại.