Tổng quan nghiên cứu

Theo thống kê của Tổ chức Y tế Thế giới, các khối u hệ thần kinh trung ương tuy chỉ chiếm khoảng 2% tổng số các ca ung thư trên toàn cầu nhưng lại nằm trong nhóm bệnh lý có tỷ lệ tử vong cao hàng đầu. Đáng chú ý, có khoảng 33,3% trong tổng số các ca u não được phân loại là u ác tính ở cấp độ III và cấp độ IV, đe dọa trực tiếp đến tính mạng bệnh nhân nếu không được can thiệp kịp thời. Kỹ thuật chụp cộng hưởng từ (MRI) đa kênh là tiêu chuẩn vàng trong chẩn đoán lâm sàng nhờ khả năng hiển thị tương phản vượt trội giữa các cấu trúc mô mềm. Tuy nhiên, việc phân đoạn thủ công từng lát cắt thể tích 3D tiêu tốn từ 30 đến 60 phút mỗi ca chụp, đồng thời phụ thuộc lớn vào kinh nghiệm chủ quan của bác sĩ chẩn đoán hình ảnh.

Luận văn tập trung giải quyết bài toán tự động hóa hoàn toàn quy trình phân đoạn và bóc tách các thành phần bất thường trên ảnh MRI não người. Mục tiêu then chốt của nghiên cứu là xây dựng mô hình phân lớp thông minh có khả năng nhận diện chính xác toàn bộ khối u cũng như bóc tách 4 cấu trúc vi mô phức tạp: nhân hoại tử, phần u mở rộng, phần u không mở rộng và vùng phù nề. Đề tài được triển khai thực nghiệm toàn diện trong giai đoạn từ tháng 1 năm 2016 đến tháng 12 năm 2016 tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM, sử dụng bộ dữ liệu chuẩn quốc tế gồm 160 ca chụp MRI thể tích não bộ. Kết quả nghiên cứu mang lại giá trị khoa học và ứng dụng vượt bậc khi nâng hệ số tương đồng không gian DICE lên mức 85,1% cho toàn bộ khối u và 73,0% cho vùng khối u mở rộng, hỗ trợ đắc lực cho công tác lập kế hoạch phẫu thuật thần kinh và xạ trị chuẩn xác.

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 xây dựng trên nền tảng lý thuyết học máy có giám sát với thuật toán Rừng ngẫu nhiên (Random Forests), kết hợp mô hình tối ưu hóa năng lượng trường ngẫu nhiên Markov (Markov Random Fields) và kỹ thuật Cắt đồ thị (Graph-Cut). Thuật toán Rừng ngẫu nhiên sử dụng tập hợp 64 cây quyết định được huấn luyện độc lập và song song, giúp cân bằng phương sai và giảm thiểu hiện tượng quá khớp (overfitting) trên tập dữ liệu đa chiều.

Các khái niệm cốt lõi được áp dụng trong luận văn bao gồm: chuỗi xung cộng hưởng từ đa kênh (T1, T1-gad, T2, T2-FLAIR), đơn vị thể tích ảnh (voxel 3D), phân đoạn ngữ nghĩa ảnh y tế, và hệ số tương đồng DICE. Nhằm khắc phục nhược điểm thiếu thông tin tương quan không gian của các mô hình truyền thống, đề tài tích hợp lý thuyết kết cấu thông qua ma trận đồng xuất hiện mức xám (Gray-Level Co-occurrence Matrix - GLCM). Ma trận này được thiết lập trong không gian 3D trên cửa sổ lân cận kích thước 21x21x21 voxel với 13 hướng không gian độc lập, lượng hóa ảnh về 64 mức xám phi tuyến tính theo hàm logarit tự nhiên. Từ đó, 4 chỉ số kết cấu then chốt gồm độ tương phản (Contrast), năng lượng (Energy), độ đồng chất (Homogeneity) và độ tương quan (Correlation) được trích xuất nhằm mô tả chính xác ranh giới mô bệnh lý.

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

Nguồn dữ liệu thực nghiệm được trích xuất từ cơ sở dữ liệu chuẩn hóa của cuộc thi phân đoạn khối u não MICCAI BRATS. Cỡ mẫu nghiên cứu gồm 160 tập ảnh chụp 3D MRI đa kênh của bệnh nhân u não, đã được các chuyên gia y tế đầu ngành chú giải và gán nhãn thủ công (ground truth). Phương pháp chọn mẫu áp dụng cơ chế phân chia ngẫu nhiên độc lập 50-50, bao gồm 80 mẫu thể tích dùng cho quá trình huấn luyện mô hình và 80 mẫu hoàn toàn tách biệt dùng cho quá trình kiểm thử đánh giá.

Phương pháp phân tích được thực hiện theo quy trình 2 giai đoạn bằng ngôn ngữ Python với thư viện scikit-learn. Giai đoạn 1 sử dụng 18 đặc trưng cơ bản (cường độ sáng 4 kênh, bộ lọc Gaussian, Sobel, Laplacian, thống kê cục bộ 5x5x5 và 3 tọa độ không gian tương đối so với điểm sáng cực đại) kết hợp ngưỡng xác suất 0,6 để tách toàn bộ khối u. Giai đoạn 2 sử dụng 70 đặc trưng (bổ sung 52 đặc trưng từ ma trận đồng xuất hiện) để phân loại 4 thành phần u chi tiết. Toàn bộ thuật toán được thực thi trên hạ tầng máy chủ điện toán đám mây Amazon AWS cấu hình 4 nhân CPU, 32GB RAM và 30GB bộ nhớ trong thời gian 12 tháng. Việc lựa chọn Random Forests kết hợp Graph-Cut mở rộng giúp bảo đảm tính liên tục của đường biên mô học và tối ưu hóa độ chính xác toàn cục.

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

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

Quá trình kiểm thử thực nghiệm trên 80 tập dữ liệu ảnh 3D MRI độc lập đã đem lại 4 phát hiện khoa học quan trọng:

  • Phát hiện 1 (Phân đoạn toàn bộ khối u): Mô hình đề xuất kết hợp tọa độ tương đối và ma trận đồng xuất hiện đạt hệ số DICE 85,1%, cải thiện 3,1% so với phương pháp tham chiếu nền tảng của Chawla năm 2015 (đạt 82,0%). Riêng phương pháp tích hợp tọa độ đơn lẻ đạt 84,0% và phương pháp chỉ dùng ma trận đồng xuất hiện đạt 84,7%.
  • Phát hiện 2 (Bóc tách nhân khối u hoại tử): Độ chính xác phân đoạn vùng nhân hoại tử đạt chỉ số DICE 70,0%, tăng trưởng 3,0% so với mức 67,0% của mô hình cơ sở. Các đặc trưng kết cấu GLCM đóng vai trò quyết định trong việc nhận diện vùng hoại tử trung tâm.
  • Phát hiện 3 (Nhận diện phần khối u mở rộng): Kết quả phân vùng cấu trúc u mở rộng đạt chỉ số DICE 73,0%, vượt trội 5,0% so với phương pháp gốc (đạt 68,0%) và tăng 3,0% so với mô hình chỉ sử dụng đặc trưng tọa độ (đạt 70,0%).
  • Phát hiện 4 (Hiệu quả thời gian tính toán): Thời gian thực thi mô hình chỉ dùng tọa độ là 0,8 giờ, trong khi mô hình kết hợp đầy đủ đòi hỏi 3,8 giờ trên cùng cấu hình máy chủ. Việc chỉ kích hoạt 52 đặc trưng GLCM ở bước 2 đã tiết kiệm hơn 75% tổng thời gian tính toán so với việc áp dụng trên toàn bộ thể tích sọ não.

Thảo luận kết quả

Sự gia tăng vượt trội về độ chính xác phân đoạn bắt nguồn từ việc mô hình hóa thành công cấu trúc giải phẫu học của khối u não. Trong tự nhiên, khối u não thường có xu hướng phân tầng hình thái từ trong ra ngoài: nhân hoại tử ở trung tâm, tiếp đến là phần khối u mở rộng, phần u không mở rộng và ngoài cùng là vùng phù nề. Việc bổ sung 3 tọa độ không gian tương đối so với điểm có cường độ sáng cực đại đã cung cấp thông tin định vị khoảng cách hình học, giúp thuật toán phân biệt chính xác lõi u và vỏ u. Đồng thời, 52 đặc trưng ma trận đồng xuất hiện GLCM 3D giúp nắm bắt chi tiết sự biến thiên kết cấu cục bộ giữa các voxel lân cận, loại bỏ hiện tượng phân lớp nhầm lẫn giữa vùng phù nề và dịch não tủy.

Dữ liệu thực nghiệm của nghiên cứu được hệ thống hóa rõ ràng thông qua bảng đối chiếu đa tiêu chí và biểu đồ cột so sánh tương quan giữa 4 cấu hình thuật toán. Đồ thị cho thấy rõ xu hướng tăng trưởng tuyến tính của hệ số DICE trên cả 3 cấu trúc vùng (toàn bộ khối u, nhân u và u mở rộng) khi lần lượt tích hợp các lớp đặc trưng mới. So với các giải pháp truyền thống như Support Vector Machine (thường nhạy cảm với dữ liệu mất cân bằng) hay Fuzzy C-Means (dễ bị nhiễu biên), việc kết hợp Random Forests với bước hậu xử lý lọc hình thái 3D (phần tử cầu kích thước 13x13x13 voxel, loại bỏ khối liên thông dưới 3000 điểm) và thuật toán Graph-Cut đã triệt tiêu hoàn toàn các đốm nhiễu ngẫu nhiên, tạo ra đường biên u mềm mại và có độ tin cậy y khoa cao.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp kỹ thuật và ứng dụng thực tiễn mang tính khả thi cao:

  • Tối ưu hóa thuật toán trích xuất đặc trưng bằng tính toán song song: Bộ phận R&D tại các viện nghiên cứu công nghệ y tế cần triển khai tái cấu trúc mã nguồn tính toán ma trận đồng xuất hiện GLCM 13 hướng trên nền tảng xử lý đồ họa GPU thông qua CUDA. Mục tiêu cụ thể là rút ngắn thời gian xử lý toàn bộ ca chụp từ 3,8 giờ xuống dưới 0,5 giờ trong lộ trình 6 tháng tới.
  • Tích hợp kiến trúc học sâu đa tầng: Đội ngũ kỹ sư khoa học dữ liệu nên nghiên cứu kết hợp mô hình phân tầng Random Forests hiện tại với mạng nơ-ron tích chập 3D U-Net trong giai đoạn 12 tháng tiếp theo. Đích ngắm là nâng cao chỉ số DICE của vùng nhân hoại tử từ mức 70,0% hiện tại lên vượt ngưỡng 78,0%.
  • Hoàn thiện và chuẩn hóa phần mềm giao diện tương tác lâm sàng: Nhóm phát triển phần mềm y tế cần tích hợp sâu module phân đoạn tự động vào hệ thống ITK-SNAP thông qua giao thức truyền thông mạng máy chủ - máy trạm. Thời gian phản hồi xử lý ảnh cần được tối ưu dưới 15 giây mỗi ca bệnh nhằm đáp ứng trực tiếp nhu cầu hội chẩn của bác sĩ trong vòng 9 tháng.
  • Triển khai thử nghiệm lâm sàng đa trung tâm: Các bệnh viện chuyên khoa ngoại thần kinh và ung bướu cần phối hợp thực hiện thử nghiệm diện rộng trên cỡ mẫu tối thiểu 500 ca bệnh mới trong vòng 18 tháng. Hoạt động này nhằm chuẩn hóa ngưỡng xác suất phân lớp 0,6 trên nhiều dòng máy quét MRI khác nhau như 1.5 Tesla và 3.0 Tesla.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật chuyên sâu cho 4 nhóm đối tượng chính:

  • Bác sĩ Chẩn đoán hình ảnh và Phẫu thuật thần kinh: Cung cấp cơ chế phân tích tự động 4 vùng mô bệnh lý trên các chuỗi xung MRI đa kênh, giúp bác sĩ tiết kiệm khoảng 70% thời gian đo đạc thủ công và nâng cao độ an toàn khi xác định ranh giới bóc tách khối u.
  • Học viên cao học và Nghiên cứu sinh ngành Khoa học máy tính, Trí tuệ nhân tạo: Cung cấp tài liệu tham khảo chi tiết về kỹ thuật trích xuất véc-tơ đặc trưng 70 chiều, phương pháp kết hợp thuật toán phân loại Rừng ngẫu nhiên với tối ưu hóa năng lượng trường ngẫu nhiên Markov (MRF/Graph-Cut) trên dữ liệu ảnh 3D.
  • Kỹ sư phát triển phần mềm y tế và Công nghệ y tế (HealthTech): Tham khảo kiến trúc tích hợp hệ thống giữa thư viện mã nguồn mở ITK/ITK-SNAP với máy chủ điện toán đám mây Amazon AWS, phục vụ thiết kế các hệ thống lưu trữ và truyền hình ảnh y tế (PACS) thế hệ mới.
  • Nhà quản lý bệnh viện và Chuyên gia chuyển đổi số y tế: Cung cấp luận cứ kỹ thuật và kinh tế để đánh giá hiệu quả đầu tư trang thiết bị máy chủ xử lý phân tích hình ảnh (tiêu chuẩn 4 nhân CPU, 32GB RAM trở lên), thúc đẩy tiến trình số hóa bệnh viện thông minh.

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

  • Tại sao nghiên cứu lại chọn thuật toán Rừng ngẫu nhiên thay vì mô hình SVM hay K-means?
    Thuật toán Rừng ngẫu nhiên với 64 cây quyết định sở hữu khả năng xử lý xuất sắc tập dữ liệu đa chiều gồm 70 đặc trưng và chống nhiễu vượt trội. Trong khi K-means có độ chính xác thấp và SVM tốn nhiều tài nguyên khi phân loại đa lớp, Rừng ngẫu nhiên cho phép tính toán xác suất phân lớp tin cậy và hỗ trợ xử lý song song hiệu quả.

  • Đặc trưng ma trận đồng xuất hiện (GLCM) đóng góp cụ thể như thế nào vào kết quả phân đoạn?
    GLCM lượng hóa kết cấu không gian 3D qua 13 hướng xung quanh mỗi voxel trong vùng lân cận 21x21x21 điểm ảnh. Nhờ 52 đặc trưng kết cấu bổ sung này, hệ số DICE của vùng u mở rộng đã tăng mạnh 5,0% (từ 68,0% lên 73,0%), giúp phân định ranh giới phức tạp giữa mô u sống và mô lành.

  • Làm thế nào để mô hình kiểm soát được thời gian tính toán khi xử lý dữ liệu thể tích 3D lớn?
    Nghiên cứu áp dụng quy trình phân đoạn 2 giai đoạn thông minh. Giai đoạn 1 chỉ sử dụng 18 đặc trưng cơ bản với thời gian xử lý 0,8 giờ để định vị khối u. Giai đoạn 2 mới tính toán 52 đặc trưng GLCM trên vùng voxel thuộc khối u, giúp tiết kiệm hơn 75% khối lượng tính toán so với cách làm thông thường.

  • Vai trò của phương pháp hậu xử lý Markov Random Fields và Graph-Cut là gì?
    Các thuật toán học máy điểm ảnh độc lập thường tạo ra các voxel nhiễu rời rạc. Kỹ thuật Markov Random Fields kết hợp thuật toán Graph-Cut tối ưu hóa hàm năng lượng không gian, tự động loại bỏ các khối vụn dưới 3000 voxel và làm trơn đường bao ranh giới u một cách tự nhiên.

  • Công cụ trực quan hóa tích hợp ITK-SNAP hoạt động theo cơ chế nào?
    Công cụ được xây dựng dưới dạng module mở rộng kết nối trực tiếp với giao diện ITK-SNAP qua địa chỉ IP máy chủ. Người dùng chỉ cần tải dữ liệu MRI 3D, nhấn lệnh phân đoạn để hệ thống gửi dữ liệu lên máy chủ xử lý và nhận về kết quả hiển thị 4 màu sắc trực quan theo từng phân vùng u.

Kết luận

  • Luận văn đã xây dựng thành công quy trình tự động hóa hoàn toàn 2 giai đoạn phân đoạn khối u não trên ảnh 3D MRI bằng thuật toán Rừng ngẫu nhiên 64 cây quyết định.
  • Đề xuất đột phá hai nhóm đặc trưng mới gồm tọa độ không gian tương đối và 52 chỉ số kết cấu ma trận đồng xuất hiện GLCM 3D trên 13 hướng không gian.
  • Nâng cao vượt bậc độ chính xác phân vùng với hệ số tương đồng DICE đạt 85,1% cho toàn bộ khối u (tăng 3,1%) và 73,0% cho vùng khối u mở rộng (tăng 5,0%).
  • Hiện thực hóa hoàn chỉnh công cụ phần mềm trực quan trên nền tảng ITK-SNAP kết nối điện toán đám mây, phục vụ đắc lực cho công tác chẩn đoán hình ảnh.
  • Thiết lập nền tảng khoa học vững chắc để phát triển các hệ thống trợ lý AI y tế thế hệ mới kết hợp mô hình học sâu và tính toán song song GPU trong lộ trình 12 đến 18 tháng tới. Quý độc giả, các nhà nghiên cứu và kỹ sư y tế quan tâm có thể khai thác trực tiếp khung phương pháp luận của luận văn để mở rộng nghiên cứu sang các dạng bệnh lý tổn thương thần kinh phức tạp khác.