Tổng quan nghiên cứu

Theo thống kê từ ngành y tế, khoảng 90% người dân Việt Nam mắc các bệnh lý liên quan đến răng miệng, trong đó 75% trường hợp chịu ảnh hưởng bởi sâu răng vĩnh viễn và hơn 60% mắc các tổn thương cấu trúc cứng như mòn răng, vỡ thân răng hoặc viêm tiêu xương ổ răng. Trong bối cảnh men răng chứa đến 96% hàm lượng chất vô cơ và ngà răng chiếm khoảng 75%, việc phát hiện sớm các tổn thương ngầm đóng vai trò sống còn trong việc bảo tồn răng thật. Để phục vụ công tác chẩn đoán và điều trị, ảnh X-quang nha khoa (bao gồm ảnh cắn cánh bitewing, ảnh quanh chóp periapical và ảnh toàn cảnh panoramic) là công cụ chỉ định phổ biến nhất nhờ chi phí hợp lý và khả năng ghi nhận cấu trúc mô học sâu.

Tuy nhiên, bài toán phân đoạn ảnh X-quang nha khoa tự động đối mặt với nhiều rào cản kỹ thuật lớn do độ tương phản cục bộ thấp, nhiễu tia X và sự chồng lấn mức xám phức tạp giữa ba vùng: mô răng (độ sáng cao nhất), xương ổ răng - nướu (mức xám trung bình) và vùng nền (mức xám tối nhất). Các thuật toán học máy không giám sát truyền thống như Fuzzy C-Means (FCM) thường dễ rơi vào điểm cực tiểu cục bộ và rất nhạy cảm với dữ liệu nhiễu. Luận văn thạc sĩ chuyên ngành Khoa học máy tính của tác giả Lê Thị Mai Hương, được hoàn thành năm 2017 tại Trường Đại học Công nghệ Thông tin và Truyền thông – Đại học Thái Nguyên dưới sự hướng dẫn khoa học của TS. Nguyễn Đình Dũng, đã tập trung giải quyết triệt để hạn chế này. Mục tiêu cốt lõi của công trình là nghiên cứu các giải thuật phân cụm bán giám sát mờ hiện đại và đề xuất mô hình lai ghép đa tầng Otsu - FCM - eSFCM, giúp tự động hóa quá trình phân tách vùng răng với độ chính xác phân đoạn đạt trên 92%, hỗ trợ giảm thiểu 40% thời gian phân tích hình ảnh của bác sĩ lâm sàng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng vững chắc của lý thuyết Khám phá tri thức trong cơ sở dữ liệu (KDD) 5 giai đoạn: Trích chọn dữ liệu, Tiền xử lý, Chuyển đổi, Khai phá dữ liệu và Đánh giá tri thức. Luận văn vận dụng sâu các mô hình toán học và lý thuyết học máy chuyên sâu sau:

  • Lý thuyết tập mờ và Thuật toán phân cụm mờ FCM (Bezdek, 1981): Khắc phục nhược điểm của phân cụm rõ (crisp clustering) khi đường biên giữa các vùng mô sinh học không tách biệt hoàn toàn. Mỗi điểm ảnh $x_k$ được gán một độ thuộc $u_{kj} \in [0, 1]$ đối với từng cụm $j$ thông qua việc tối thiểu hóa hàm mục tiêu có trọng số mờ hóa $m$ (thông thường chọn $m = 2$).
  • Thuật toán phân cụm nửa giám sát chuẩn SSSFC và eSFCM: Mô hình eSFCM kết hợp thành phần hiệu chỉnh Entropy vào hàm mục tiêu cùng ma trận khoảng cách Mahalanobis nhằm tối ưu hóa quá trình phân bổ độ thuộc, tận dụng triệt để thông tin bổ trợ định hướng từ chuyên gia hoặc từ bước tiền phân cụm.
  • Thuật toán phân ngưỡng Otsu (1979): Kỹ thuật nhị phân hóa dựa trên việc tối đa hóa phương sai giữa các lớp (inter-class variance), cho phép bóc tách tự động vùng nền tối với tốc độ xử lý nhanh và độ ổn định cao.

Các khái niệm then chốt trong nghiên cứu bao gồm: Không gian metric khoảng cách (Euclide, Manhattan, Minkowski), Ma trận phân hoạch mờ $U$, Tập tâm cụm $V$, Ràng buộc Must-link và Cannot-link, cùng các chỉ số hợp lệ cụm chuyên dụng.

Phương pháp nghiên cứu

  • Nguồn dữ liệu và Cỡ mẫu: Nghiên cứu sử dụng tập dữ liệu thực nghiệm gồm 120 ảnh X-quang nha khoa số hóa đa định dạng (ảnh toàn cảnh panoramic kích thước trung bình $1500 \times 800$ pixel, ảnh cắn cánh và ảnh quanh chóp $600 \times 400$ pixel) được thu thập từ các cơ sở khám chữa bệnh nha khoa và bộ dữ liệu y khoa chuẩn.
  • Phương pháp chọn mẫu: Luận văn áp dụng phương pháp chọn mẫu phân tầng có chủ đích (purposive stratified sampling). Dữ liệu được chia thành các nhóm đại diện cho từng hình thái giải phẫu: nhóm cấu trúc răng chuẩn (40 ảnh), nhóm răng chèn ép/mất răng (35 ảnh), nhóm có can thiệp vật liệu trám kim loại gây nhiễu cản tia (25 ảnh) và nhóm bệnh lý viêm quanh chóp/tiêu xương (20 ảnh).
  • Lý do lựa chọn phương pháp phân tích: Mô hình kết hợp lai ghép Otsu - FCM - eSFCM được lựa chọn vì thuật toán Otsu giúp triệt tiêu 100% không gian tính toán của vùng nền đen ngoài miệng; thuật toán FCM đóng vai trò khởi tạo thông tin bổ trợ sơ bộ không giám sát; và thuật toán eSFCM tối ưu hóa biên độ thuộc mờ, triệt tiêu nhiễu cục bộ và tránh rơi vào bẫy cực trị địa phương.
  • Quy trình và Timeline nghiên cứu: Công trình được thực hiện trong thời gian 12 tháng (từ tháng 6/2016 đến tháng 6/2017). Ứng dụng thực nghiệm được xây dựng hoàn chỉnh trên môi trường MATLAB, tích hợp đầy đủ quy trình từ đọc ảnh, tiền xử lý, phân cụm lai ghép đến xuất kết quả đánh giá định lượng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  • Hiệu quả khử nền vượt trội: Việc áp dụng thuật toán tách ngưỡng Otsu ở giai đoạn tiền xử lý đã loại bỏ chính xác 98.5% diện tích vùng nền đen và mô mềm bên ngoài cung răng, giúp giảm 35% dung lượng dữ liệu cần tính toán ma trận độ thuộc ở các bước sau.
  • Nâng cao chất lượng phân tách cụm: Mô hình lai ghép nửa giám sát mờ eSFCM với thông tin bổ trợ từ ma trận $U_{FCM}$ đạt chỉ số PBM (Pakhira-Bandyopadhyay-Maulik) trung bình là 3.65, tăng 28.3% so với thuật toán phân cụm mờ FCM truyền thống (chỉ số PBM chỉ đạt 2.84).
  • Tối ưu hóa tốc độ hội tụ: Thuật toán eSFCM có bổ trợ thông tin rút ngắn số vòng lặp hội tụ từ mức trung bình 48 chu kỳ ở FCM xuống còn 19 đến 22 chu kỳ lặp với cùng ngưỡng dừng $\varepsilon = 10^{-5}$, giúp tiết kiệm khoảng 42% thời gian xử lý CPU trên mỗi bức ảnh.
  • Độ chính xác nhận diện biên răng bệnh lý: Trên các mẫu ảnh phức tạp có răng nêm chèn hoặc chứa vật liệu trám răng nhân tạo, thuật toán đề xuất giữ vững độ tách biệt biên thân răng và xương hàm với tỷ lệ lỗi gán nhãn pixel giảm xuống dưới 6.8%, khắc phục hiện tượng nhòe đường nối men - xê măng mà các phương pháp truyền thống thường mắc phải.

Thảo luận kết quả

Sự vượt trội của mô hình lai ghép bắt nguồn từ cơ chế kiểm soát bất định đa tầng: Thuật toán Otsu đóng vai trò "ngắt kênh" đối với các pixel nền có cường độ xám tiệm cận 0, giúp FCM không bị phân tán tâm cụm vào các vùng vô nghĩa. Tiếp đó, việc loại bỏ các giá trị hàm thuộc cực tiểu trong ma trận $U_{FCM}$ trước khi nạp vào eSFCM đã tạo ra một ma trận tri thức bổ trợ sắc nét, giúp thuật toán eSFCM định hướng hàm mục tiêu hội tụ nhanh chóng vào cấu trúc giải phẫu thực của thân và chân răng.

+-------------------------------------------------------------------------+
|                  SO SANH HIEN THI CHI SO DANH GIA                       |
|                                                                         |
|  Chi so DB (Davies-Bouldin) - Cang nho cang tot:                       |
|  FCM truyen thong:   [##############                    ] 1.48          |
|  SSSFC ban giam sat: [##########                        ] 1.05          |
|  Mo hinh Lai ghep:   [######                            ] 0.68 (-54.0%) |
|                                                                         |
|  Chi so PBM - Cang lon cang tot:                                        |
|  FCM truyen thong:   [############                      ] 2.84          |
|  SSSFC ban giam sat: [###############                   ] 3.18          |
|  Mo hinh Lai ghep:   [##################                ] 3.65 (+28.3%) |
+-------------------------------------------------------------------------+

Dữ liệu phân nghiệm được minh họa rõ nét qua Bảng so sánh định lượng các chỉ số hình học cụm và Biểu đồ đường biểu diễn sự suy giảm của hàm mục tiêu theo từng bước lặp. Kết quả cho thấy đồ thị hàm mất mát của mô hình lai ghép dốc đứng và đạt trạng thái bình ổn chỉ sau 15 bước, trong khi đồ thị của FCM và SSSFC dao động kéo dài đến bước thứ 45. Điều này khẳng định thuật toán lai ghép không chỉ cải thiện chất lượng hình ảnh đầu ra mà còn mang tính khả thi cao trong triển khai thực tế.

Đề xuất và khuyến nghị

  • Nâng cấp giải thuật xử lý trên không gian 3 chiều: Đội ngũ kỹ sư nghiên cứu cần mở rộng mô hình phân cụm nửa giám sát mờ sang cấu trúc voxel để xử lý dữ liệu ảnh chụp cắt lớp vi tính hình nón (Cone Beam CT), hướng tới mục tiêu nâng độ chính xác tái tạo bề mặt 3D lên 96% trong lộ trình 12 tháng tới.
  • Tích hợp bộ lọc thích nghi khử nhiễu tiền xử lý: Nhóm phát triển phần mềm y tế cần thiết lập module lọc nhiễu Wavelet đa mức hoặc lọc khuếch tán bất đẳng hướng (Anisotropic Diffusion) nhằm triệt tiêu thêm 15% nhiễu tán xạ tia X trước khi thực hiện tách ngưỡng Otsu, hoàn thành trong quý 2 năm tiếp theo.
  • Chuẩn hóa kho dữ liệu mở về ảnh X-quang nha khoa: Các bệnh viện răng hàm mặt trung ương cùng các trường đại học chuyên ngành công nghệ thông tin cần phối hợp xây dựng kho dữ liệu 5.000 ảnh X-quang có nhãn chuyên gia chuẩn hóa trong thời gian 18 tháng để phục vụ huấn luyện và kiểm thử benchmark.
  • Tối ưu hóa mã nguồn và nhúng vào hệ thống phần mềm phòng khám: Doanh nghiệp công nghệ y tế (MedTech) cần tái cấu trúc thuật toán sang ngôn ngữ C++ hoặc Python Tensor, mục tiêu đạt tốc độ phân đoạn thời gian thực dưới 0.8 giây/ảnh trên các thiết bị máy tính nha khoa tiêu chuẩn trong giai đoạn 24 tháng.

Đối tượng nên tham khảo luận văn

  • Học viên cao học và Nghiên cứu sinh ngành Khoa học máy tính: Tiếp cận khung lý thuyết toán học chuyên sâu về tập mờ, phương pháp giải tối ưu Lagrange và mô hình tích hợp ràng buộc entropy vào hàm mục tiêu trong học máy nửa giám sát.
  • Kỹ sư thị giác máy tính và Xử lý ảnh y tế: Khai thác quy trình kiến trúc thực tế, sơ đồ usecase, sơ đồ tuần tự và toàn bộ mã nguồn MATLAB mẫu được đính kèm ở phụ lục luận văn để áp dụng vào các dự án phân đoạn ảnh y khoa.
  • Bác sĩ nha khoa và Chuyên gia chẩn đoán hình ảnh: Nắm bắt cơ chế nhận dạng tự động cấu trúc răng, men răng và xương ổ răng của các thuật toán trí tuệ nhân tạo, từ đó nâng cao độ tin cậy khi sử dụng các hệ thống hỗ trợ ra quyết định lâm sàng.
  • Nhà phát triển sản phẩm công nghệ y tế (MedTech): Sử dụng các kết quả thực nghiệm và bảng đối sánh hiệu năng của luận văn làm căn cứ kỹ thuật để thiết kế tính năng tiền phân vùng tự động trong các phần mềm quản lý và chẩn đoán nha khoa.

Câu hỏi thường gặp

Thuật toán phân cụm mờ FCM truyền thống gặp rào cản gì khi phân đoạn ảnh X-quang nha khoa? FCM truyền thống xem tất cả các pixel đều có vai trò ngang nhau và chỉ dựa vào khoảng cách Euclidean đến tâm cụm. Do ảnh X-quang nha khoa có độ tương phản kém và nhiễu hạt cao, thuật toán FCM thường bị rơi vào điểm cực tiểu địa phương, khiến đường biên phân đoạn răng bị mờ và nhầm lẫn giữa vùng ngà răng với xương ổ răng có mức xám xấp xỉ nhau.

Mô hình lai ghép Otsu - FCM - eSFCM hoạt động theo nguyên lý nào? Mô hình hoạt động qua ba chặng tuần tự: Đầu tiên, thuật toán tách ngưỡng Otsu loại bỏ hoàn toàn vùng nền đen để cô lập vùng giải phẫu chính. Tiếp theo, thuật toán FCM phân cụm sơ bộ để tạo ra ma trận độ thuộc ban đầu. Cuối cùng, thuật toán eSFCM sử dụng ma trận này làm tri thức bổ trợ kết hợp với hiệu chỉnh entropy để phân cụm chính xác các đường biên răng.

Thông tin bổ trợ trong thuật toán bán giám sát mờ eSFCM được trích xuất như thế nào? Thông tin bổ trợ không đòi hỏi chuyên gia phải gán nhãn thủ công từng pixel. Thay vào đó, nó được tự động trích xuất trực tiếp từ ma trận độ thuộc $U_{FCM}$ của bước phân cụm FCM trước đó, sau khi đã áp dụng quy tắc lọc bỏ các giá trị hàm thuộc nhỏ nhất tại mỗi điểm ảnh nhằm loại trừ các yếu tố nhiễu không xác định.

Thuật toán lai ghép xử lý các ca lâm sàng có chất hàn trám kim loại hoặc răng mọc lệch như thế nào? Vật liệu trám kim loại thường cản quang mạnh, tạo ra các vệt sáng chói lóa trên ảnh X-quang. Nhờ hàm khoảng cách Mahalanobis và cơ chế điều chuẩn entropy trong eSFCM, thuật toán tự động cân chỉnh ma trận hiệp phương sai của từng cụm, hạn chế hiện tượng biến dạng tâm cụm và giúp duy trì độ bao tách thân răng chính xác với sai số dưới 7%.

Làm thế nào để tái sử dụng mã nguồn MATLAB trong luận văn cho các bài toán phân đoạn ảnh khác? Toàn bộ mã nguồn MATLAB trong phụ lục luận văn được thiết kế theo cấu trúc module hóa rõ ràng, bao gồm hàm đọc ảnh, hàm tiền xử lý Otsu, hàm tính toán độ thuộc mờ và hàm đo lường chỉ số DB/PBM. Người lập trình có thể dễ dàng thay đổi số lượng cụm $C$, hệ số mờ $m$ hoặc thay thế dữ liệu đầu vào bằng ảnh X-quang phổi hoặc ảnh MRI não.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở toán học của lý thuyết tập mờ, thuật toán FCM và các biến thể phân cụm bán giám sát mờ hiện đại như SSSFC và eSFCM.
  • Đề xuất thành công mô hình lai ghép ba giai đoạn Otsu - FCM - eSFCM, giải quyết triệt để bài toán phân đoạn vùng răng trên ảnh X-quang nha khoa phức tạp.
  • Đạt hiệu quả vượt trội trong thực nghiệm với chỉ số PBM tăng 28.3%, chỉ số DB giảm xuống 0.68 và rút ngắn trên 40% số chu kỳ lặp so với các thuật toán độc lập.
  • Xây dựng hoàn chỉnh ứng dụng thử nghiệm trên MATLAB với giao diện tương tác đồ họa thân thiện, hỗ trợ đầy đủ các tính năng nạp ảnh, phân đoạn và đánh giá định lượng.
  • Công trình đóng góp nền tảng kỹ thuật vững chắc cho việc phát triển các hệ sinh thái phần mềm nha khoa thông minh hỗ trợ bác sĩ chẩn đoán và lập kế hoạch điều trị tự động.

Trong giai đoạn 12 đến 24 tháng tới, hướng phát triển kế tiếp cần tập trung vào việc mở rộng không gian thuật toán sang ảnh CT 3 chiều và tích hợp kiến trúc mạng nơ-ron tích chập sâu. Các nhà nghiên cứu và kỹ sư quan tâm được khuyến khích kế thừa khung giải thuật lai ghép này để tiếp tục phát triển các công cụ hỗ trợ y tế thiết thực.