Tổng quan nghiên cứu

Trong hệ thống y tế hiện đại, khoảng 75% các quyết định can thiệp lâm sàng và phác đồ điều trị phụ thuộc trực tiếp vào kết quả chẩn đoán hình ảnh. Tuy nhiên, sự xuất hiện tất yếu của các dạng nhiễu ngẫu nhiên trong quá trình thu nhận tín hiệu thường làm suy giảm từ 18% đến 30% độ tương phản thực tế giữa các cấu trúc mô mềm và tổn thương bệnh lý. Các thiết bị chẩn đoán quen thuộc như máy chụp X-quang, máy siêu âm (US), hệ thống chụp cắt lớp vi tính (CT) và máy chụp cộng hưởng từ (MRI) đều phải đối mặt với nhiều loại nhiễu phức tạp: từ nhiễu lượng tử photon, dao động nhiệt điện tử mang phân bố Gauss, đến nhiễu đốm hạt trong siêu âm và nhiễu Rician trên tín hiệu từ trường.

Thực trạng này đặt ra yêu cầu cấp thiết về việc phát triển các thuật toán tiền xử lý số học có khả năng loại bỏ tạp âm nhưng không làm suy hao các chi tiết cấu trúc vi mô. Được hoàn thành vào tháng 6 năm 2023 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, luận văn thạc sĩ khoa học máy tính của tác giả Nông Quốc Doãn đã đi sâu giải quyết bài toán nâng cao chất lượng ảnh y học thông qua phương trình khuếch tán phi tuyến phức (Complex Nonlinear Diffusion PDE).

Mục tiêu cụ thể của công trình là xây dựng thuật toán sai phân hữu hạn giải phương trình đạo hàm riêng trong miền phức, nhằm thực hiện song song hai nhiệm vụ: khử triệt để các nguồn nhiễu ngẫu nhiên và bảo tồn nguyên vẹn các đường biên dốc giải phẫu. Nghiên cứu giới hạn phạm vi thực nghiệm trên 4 phương thức ảnh y tế phổ biến trong giai đoạn 2022-2023. Đóng góp của luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc, hướng tới việc tăng chỉ số PSNR thêm khoảng 3 dB đến 6 dB và hạ thấp sai số bình phương trung bình MSE xuống dưới 10%, tạo điều kiện cho các bác sĩ phát hiện sớm các khối u kích thước nhỏ và giảm thiểu đáng kể tỷ lệ chẩn đoán sai lệch.

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

Khung lý thuyết áp dụng

Nền tảng của nghiên cứu dựa trên lý thuyết phương trình đạo hàm riêng parabolic mô tả hiện tượng truyền nhiệt và định luật khuếch tán Fick trong không gian đa chiều. Nhận thấy bộ lọc khuếch tán tuyến tính với hệ số hằng số thường làm mờ cả vùng đồng nhất lẫn ranh giới cơ quan, luận văn đã kế thừa mô hình khuếch tán dị hướng thực của Perona và Malik đề xuất năm 1990, kết hợp tenso khuếch tán tăng cường biên (EED) của Weickert.

Điểm đột phá của đề tài là việc ứng dụng lý thuyết khuếch tán phức do nhà khoa học Guy Gilboa khởi xướng, biến đổi hệ số khuếch tán thực thành số phức dạng $c = r \cdot e^{i\theta}$. Không gian nghiên cứu tập trung vào 3 khái niệm toán học then chốt:

  1. Cơ chế phân tách tín hiệu phức: Phần thực của phương trình đóng vai trò như một bộ lọc làm trơn mịn các vùng đồng nhất để triệt tiêu nhiễu Gauss và nhiễu hạt, trong khi phần ảo tỉ lệ thuận với đạo hàm bậc hai được làm mịn, đóng vai trò như một toán tử Laplacian nhận diện đường bao.
  2. Biên kiểu dốc (ramp edge): Dạng ranh giới mô bệnh học biến thiên độ xám từ từ, vốn rất dễ bị gãy khúc nếu áp dụng đạo hàm bậc một truyền thống.
  3. Hiện tượng bậc thang (staircasing artifact): Tác động tiêu cực thường thấy ở các mô hình khuếch tán phi tuyến cổ điển khi xử lý vùng chuyển tiếp cường độ sáng, được loại bỏ triệt để nhờ hàm dừng phụ thuộc phần ảo.

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

Nghiên cứu sử dụng tập dữ liệu thực nghiệm gồm 120 ảnh y tế chuẩn thuộc 4 nhóm phương thức chẩn đoán: 30 ảnh X-quang xương khớp, 30 ảnh siêu âm mô mềm vùng bụng, 30 lát cắt CT sọ não và ngực, cùng 30 bản ghi MRI thần kinh. Phương pháp chọn mẫu phân tầng có chủ đích được áp dụng nhằm thu thập đầy đủ các dạng mô có độ tương phản thấp và bị tác động bởi các loại nhiễu đặc thù (nhiễu đốm trên siêu âm và nhiễu Rician trên MRI).

Về phương pháp phân tích, tác giả lựa chọn sơ đồ sai phân hữu hạn hiện thời gian thuận - không gian tâm FTCS (Forward-Time Central-Space) trên lưới rời rạc 2 chiều kết hợp điều kiện biên Neumann. Lý do lựa chọn sơ đồ này là tính trực quan, dễ lập trình và cho phép kiểm soát chặt chẽ điều kiện hội tụ số học. Để thuật toán ổn định tuyệt đối trong miền thời gian thực, bước nhảy thời gian được khống chế theo bất đẳng thức $\Delta t \le 0.25 \Delta x^2$ với bước lưới không gian chuẩn hóa và góc pha $\theta = \pi/30$ (khoảng 0.1047 rad). Toàn bộ chương trình được triển khai trên môi trường tính toán MATLAB trong khung thời gian nghiên cứu 12 tháng, từ tháng 6 năm 2022 đến tháng 6 năm 2023, đánh giá thông qua 3 thước đo định lượng khách quan: PSNR, SNR và MSE so với các bộ lọc kinh điển như lọc trung vị 2D và lọc Gauss.

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

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

Quá trình chạy thử nghiệm thuật toán khuếch tán phi tuyến phức trên tập dữ liệu ảnh y khoa mang lại 4 phát hiện quan trọng:

  • Khử triệt để nhiễu đốm trên ảnh siêu âm: So với bộ lọc trung vị kích thước cửa sổ 3x3 pixel, mô hình khuếch tán phức giúp tăng chỉ số PSNR từ mức cơ sở 22.15 dB lên 27.68 dB, tương đương mức cải thiện chất lượng tín hiệu đạt 24.96%. Các bọt nhiễu tán xạ âm thanh bị triệt tiêu mà không làm mất đi các nếp gấp mô mềm.
  • Nâng cao tỷ số tín hiệu trên nhiễu trên ảnh CT và X-quang: Trên các lát cắt CT lồng ngực chịu tác động của nhiễu Gauss nhiệt điện tử, tỷ số tín hiệu trên nhiễu SNR tăng mạnh từ 13.00 dB lên 19.82 dB, tức tăng trưởng hơn 52.46%. Vùng nhu mô phổi và các rãnh xương sườn được làm rõ nét rõ rệt.
  • Khắc phục hoàn toàn hiệu ứng bậc thang trên ảnh MRI: Khác với phương pháp Perona-Malik truyền thống làm giả lập các mảng màu phẳng nhân tạo, thuật toán khuếch tán phức duy trì đường chuyển tiếp tự nhiên giữa chất xám và chất trắng của não bộ, đạt hệ số tương đồng cấu trúc trên 93.5%.
  • Tối ưu hóa chu kỳ lặp và tham số góc pha: Thuật toán đạt trạng thái cân bằng khử nhiễu tối ưu chỉ sau 10 đến 15 bước lặp tính toán với giá trị góc pha nhỏ $\theta = \pi/30$, giúp sai số bình phương trung bình MSE giảm sâu 48.72% so với dữ liệu gốc ban đầu.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương trình khuếch tán phi tuyến phức vượt trội hơn các mô hình truyền thống nằm ở việc tích hợp giá trị ảo trong quá trình vi phân. Khi góc pha $\theta$ đủ nhỏ, phương trình thực chất phân rã thành hai luồng tính toán song song nhưng tương hỗ: phần thực lọc sạch các xung nhiễu ngẫu nhiên, trong khi phần ảo tự động phát hiện vị trí các điểm uốn của biên dốc để lập tức triệt tiêu hệ số dẫn truyền, ngăn chặn việc làm trơn lấn qua biên giới mô tế bào.

Trong khi bộ lọc Gauss làm suy giảm khoảng 35% độ sắc nét ở đường viền giải phẫu và bộ lọc trung vị làm biến dạng các cấu trúc mạch máu nhỏ dưới 2 mm, phương pháp khuếch tán phức bảo toàn trọn vẹn hình thái học của bệnh phẩm.

Về mặt hiển thị trực quan, kết quả nghiên cứu có thể được trình bày sinh động thông qua biểu đồ đường thể hiện tốc độ suy giảm đơn điệu của hàm mất mát MSE qua 20 chu kỳ lặp, kết hợp cùng bảng so sánh đối chuẩn 4 chỉ số (PSNR, SNR, MSE, thời gian xử lý) giữa 5 giải pháp lọc khác nhau. Ngoài ra, việc sử dụng biểu đồ cột biểu diễn mức tăng trưởng phần trăm của PSNR trên từng nhóm bệnh phẩm (siêu âm đạt 24.96%, CT đạt 31.2%, MRI đạt 28.4%) sẽ chứng minh thuyết phục khả năng thích ứng linh hoạt của mô hình với đa dạng nguồn nhiễu y sinh.

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

Dựa trên những luận cứ khoa học và kết quả thực nghiệm đã đạt được, luận văn đưa ra 4 giải pháp ứng dụng và định hướng phát triển có tính thực tiễn cao:

  1. Tích hợp module thuật toán khuếch tán phi tuyến phức vào các hệ thống lưu trữ và truyền hình ảnh y tế PACS tại các cơ sở khám chữa bệnh trong giai đoạn 2024-2025. Chủ thể thực hiện là các kỹ sư công nghệ thông tin y tế, với mục tiêu rút ngắn thời gian lọc tiền xử lý xuống dưới 150 mili-giây cho mỗi lát cắt ảnh kích thước 512x512 pixel.
  2. Chuẩn hóa bộ tham số điều khiển ngưỡng khuếch tán $k$ và hệ số pha $\theta$ tương thích với từng dòng máy siêu âm và máy quét CT chuyên dụng. Nhóm kỹ sư vật lý y sinh cần phối hợp cùng các bác sĩ chẩn đoán hình ảnh hoàn thành bộ tiêu chuẩn cấu hình trước quý 4 năm 2025, hướng đến mục tiêu nâng cao độ tương phản mô mềm tối thiểu 20%.
  3. Nghiên cứu tối ưu hóa mã nguồn thuật toán trên phần cứng xử lý song song GPU sử dụng kiến trúc CUDA hoặc OpenCL. Viện nghiên cứu và các phòng thí nghiệm đại học cần triển khai nhiệm vụ này trong khung thời gian 18 tháng, nhằm đẩy tốc độ khử nhiễu đạt trên 30 khung hình mỗi giây, đáp ứng yêu cầu xử lý hình ảnh thời gian thực trong các ca phẫu thuật nội soi can thiệp.
  4. Xây dựng bộ quy chuẩn đánh giá chất lượng ảnh y học sau khử nhiễu cấp cơ sở, áp dụng đồng thời chỉ số toán học khách quan (PSNR, MSE) và đánh giá mù đôi của hội đồng chuyên môn y khoa. Cơ quan quản lý trang thiết bị y tế cần ban hành văn bản hướng dẫn kiểm định định kỳ vào cuối năm 2025, kiểm soát sai số biến dạng hình học dưới ngưỡng 1.2%.

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

Luận văn thạc sĩ là tài liệu nghiên cứu chuyên sâu, mang lại giá trị thiết thực cho 4 nhóm đối tượng cụ thể:

  • Bác sĩ chuyên khoa chẩn đoán hình ảnh và kỹ thuật viên y tế: Giúp thấu hiểu bản chất vật lý của các dạng nhiễu trên phim X-quang, ảnh chụp CT, siêu âm và MRI. Nắm vững phương pháp này giúp đội ngũ y tế khai thác tối đa độ nét của các vi tổn thương đường kính dưới 3 mm, hỗ trợ chẩn đoán chính xác giai đoạn sớm của bệnh lý và giảm thiểu hơn 15% nguy cơ chỉ định can thiệp xâm lấn không cần thiết.
  • Học viên cao học và nghiên cứu sinh ngành Khoa học máy tính, Kỹ thuật y sinh: Cung cấp khung lý thuyết hoàn chỉnh về phương trình đạo hàm riêng ứng dụng trong thị giác máy tính. Đây là tài liệu tham khảo giá trị với hệ thống hơn 40 công thức toán học rời rạc hóa và mã nguồn mẫu trên MATLAB, rút ngắn từ 3 đến 6 tháng thực nghiệm cho các đề tài xử lý ảnh nâng cao.
  • Đội ngũ kỹ sư phát triển phần mềm thiết bị y tế: Cung cấp giải pháp kỹ thuật tối ưu để nhúng vào bộ vi xử lý của các máy quét thế hệ mới, giúp giảm tải bộ nhớ đệm khoảng 25% so với các mạng nơ-ron học sâu phức tạp nhưng vẫn bảo toàn đường biên tối ưu.
  • Giảng viên và chuyên gia nghiên cứu toán ứng dụng: Sử dụng làm tài liệu giảng dạy chuyên đề mô hình hóa hiện tượng vật lý trong không gian số và chuyển giao công nghệ cho các dự án hợp tác y tế - công nghệ thông tin.

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

Tại sao phương trình khuếch tán phức lại vượt trội hơn khuếch tán thực trong việc bảo toàn biên?

Trong khuếch tán thực, hệ số khuếch tán phụ thuộc vào gradient bậc một, dễ dẫn đến hiện tượng bậc thang làm gãy khúc các đoạn biên dốc tự nhiên. Ngược lại, khuếch tán phức khai thác phần ảo tương đương đạo hàm bậc hai được làm trơn. Nhờ đó, thuật toán nhận diện chính xác các điểm uốn ranh giới mà không làm chậm quá trình khuếch tán tại sườn dốc, giúp tăng chỉ số PSNR thêm khoảng 3 dB đến 5 dB so với mô hình thực.

Thuật toán có làm biến dạng các chi tiết giải phẫu nhỏ trong ảnh y khoa không?

Không. Với thiết kế hàm dừng khuếch tán giảm theo độ lớn phần ảo và số chu kỳ lặp được kiểm soát chặt chẽ từ 10 đến 15 bước, thuật toán tự động đóng băng quá trình làm mịn khi tiếp cận các cấu trúc biên sắc nét. Dữ liệu thực nghiệm trên 120 ảnh y tế cho thấy sai số biến dạng hình học luôn duy trì dưới 1.5%, bảo tồn nguyên vẹn các mạch máu và vi tổn thương có kích thước nhỏ từ 1 đến 2 mm.

Góc pha theta ảnh hưởng như thế nào đến độ ổn định của thuật toán?

Góc pha $\theta$ quyết định sự tương tác giữa phần thực và phần ảo trong hệ số khuếch tán phức $c = e^{i\theta}$. Nghiên cứu đã chứng minh thuật toán chỉ đạt trạng thái ổn định và hội tụ toán học cao khi $\theta$ tiến dần về 0 (thực nghiệm chọn $\theta = \pi/30$). Nếu tăng $\theta$ tiệm cận giá trị $\pi/2$, bước thời gian bị thu hẹp đáng kể, gây mất ổn định số học và phát sinh lỗi tính toán nghiêm trọng.

Phương pháp này có xử lý được ảnh chụp màu y tế hay chỉ áp dụng cho ảnh thang độ xám?

Thuật toán ban đầu được thiết kế chuyên biệt cho ảnh thang độ xám 2D, vốn là định dạng tiêu chuẩn của hơn 90% dữ liệu X-quang, CT và MRI. Tuy nhiên, mô hình hoàn toàn có thể mở rộng cho ảnh màu (như ảnh nội soi tiêu hóa) bằng cách phân tách không gian màu thành các kênh độc lập, hoặc xử lý trên kênh độ sáng mà vẫn giữ nguyên vẹn thông tin màu sắc gốc.

Thời gian tính toán của mô hình khuếch tán phức trên máy tính thông thường là bao lâu?

Trên cấu hình máy tính cá nhân tiêu chuẩn sử dụng bộ vi xử lý đa nhân thông dụng, thuật toán giải sai phân FTCS cho một ảnh kích thước 512x512 pixel mất khoảng 0.25 đến 0.45 giây cho 15 vòng lặp. Khi được chuyển đổi sang ngôn ngữ C++ hoặc tăng tốc phần cứng thông qua chip đồ họa GPU, thời gian xử lý có thể giảm xuống dưới 30 mili-giây, hoàn toàn đáp ứng yêu cầu lâm sàng trực tiếp.

Kết luận

  • Luận văn thạc sĩ khoa học máy tính năm 2023 của tác giả Nông Quốc Doãn đã xây dựng thành công giải pháp tăng cường chất lượng ảnh y học dựa trên nền tảng phương trình khuếch tán phi tuyến phức.
  • Công trình giải quyết triệt để sự đánh đổi nan giải giữa nhiệm vụ triệt tiêu tạp âm ngẫu nhiên và yêu cầu bảo toàn đường biên dốc giải phẫu trên 4 phương thức ảnh: X-quang, siêu âm, CT và MRI.
  • Kết quả thực nghiệm chứng minh tính ưu việt vượt trội của mô hình với chỉ số PSNR tăng khoảng 24.96% trên ảnh siêu âm và tỷ số SNR tăng trên 52% đối với ảnh chụp cắt lớp vi tính.
  • Kế hoạch phát triển tiếp theo tập trung vào việc lập trình tối ưu hóa song song trên GPU và mở rộng thuật toán sang mô hình dữ liệu thể tích khối 3D trong vòng 12 đến 24 tháng tới.
  • Để nâng cao độ chính xác trong công tác sàng lọc bệnh học, các đơn vị y tế và chuyên gia công nghệ hãy áp dụng thử nghiệm mô hình khuếch tán phức vào quy trình tiền xử lý ảnh y sinh ngay hôm nay.