Tổng quan nghiên cứu

Sự phát triển của công nghệ thông tin trong y học hiện đại đã tạo ra bước đột phá với các thiết bị chẩn đoán hình ảnh tiên tiến như máy chụp cắt lớp vi tính (CT-Scanner) và máy chụp cộng hưởng từ (MRI). Tại các cơ sở y tế công lập như Bệnh viện A Thái Nguyên – nơi tiếp nhận và điều trị nội trú cho gần 1.000 bệnh nhân mỗi ngày – nhu cầu chẩn đoán hình ảnh chính xác đóng vai trò quyết định đến hiệu quả điều trị. Tuy nhiên, phần lớn các hệ thống máy tính đi kèm máy chụp CT thế hệ cũ chỉ cung cấp chuỗi ảnh cắt lớp 2D rời rạc, gây khó khăn cho bác sĩ trong việc hình dung toàn diện không gian giải phẫu và các tổn thương phức tạp.

Trong khi các phần mềm thương mại ngoại nhập chuyên dụng như 3D-Doctor, eFilm, hay Vitrea2 có chi phí bản quyền lên tới hàng chục nghìn USD (thường vượt quá 20.000 USD cho mỗi trạm làm việc PACS), các bệnh viện tuyến tỉnh và địa phương rất khó tiếp cận do ngân sách hạn chế. Bài toán đặt ra là cần nghiên cứu và làm chủ kỹ thuật xây dựng mô hình 3D từ chuỗi ảnh chụp cắt lớp y tế trên các dòng máy tính phổ thông.

Mục tiêu cụ thể của nghiên cứu là xây dựng quy trình thuật toán hoàn chỉnh từ khâu đọc dữ liệu ảnh chuẩn DICOM, tiền xử lý khử nhiễu, phân đoạn tách biên cho đến tái tạo bề mặt 3D dưới dạng lưới đa giác tam giác. Phạm vi nghiên cứu tập trung vào các bộ dữ liệu ảnh cắt lớp sọ não và chấn thương hàm mặt thu thập tại Bệnh viện A Thái Nguyên trong giai đoạn 2015–2016. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm 100% chi phí phần mềm thương mại, đồng thời rút ngắn thời gian hội chẩn trước mổ khoảng 45% và nâng cao độ chính xác nhận diện thương tổn xương lên trên 95%.

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 ứng dụng nền tảng lý thuyết tái tạo Radon (1917) và nguyên lý đo hệ số suy giảm tia X của Cormack và Hounsfield (1967) – công trình đạt giải Nobel Y học năm 1979. Nguyên lý chụp cắt lớp vi tính dựa trên khả năng hấp thụ tia X khác nhau giữa các mô, cho phép máy tính lượng hóa độ cản quang để phân biệt các cấu trúc cơ thể có độ chênh lệch tỷ trọng chỉ từ 0,5%.

Mô hình nghiên cứu áp dụng kỹ thuật biểu diễn bề mặt đối tượng 3 chiều (Boundary Representation) thông qua cấu trúc dữ liệu đồ họa bao gồm 3 thành phần cốt lõi:

  • Tập đỉnh (VertexList): Tập hợp các vector tọa độ 3 chiều $(x, y, z)$ xác định vị trí các điểm ranh giới của đối tượng trong không gian.
  • Tập mặt (FaceList): Tập hợp các chỉ số liên kết 3 đỉnh liền kề tạo thành lưới các mặt tam giác (triangle mesh), tối ưu hóa tốc độ kết xuất bề mặt.
  • Tập UV (UVs): Tập hợp các vector 2 chiều quy định tọa độ trải chất liệu (texture mapping) lên bề mặt mô hình.

Các khái niệm chính được chuẩn hóa gồm: định dạng ảnh y tế DICOM (Digital Imaging and Communications in Medicine) với các trường dữ liệu quan trọng như Bit Allocated (16 bit), Bit Stored (12-16 bit), High Bit; hệ thống lưu trữ và truyền thông hình ảnh PACS; và phương pháp nội suy láng giềng tự nhiên (Natural Neighbor Interpolation - NNI) để làm mịn bề mặt giữa các lát cắt liên tiếp.

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

Nghiên cứu sử dụng nguồn dữ liệu thực tế gồm hơn 150 lát cắt CT định dạng DICOM thu thập từ các hồ sơ bệnh án điều trị nội trú về chấn thương sọ não và hàm mặt tại Bệnh viện A Thái Nguyên. Phương pháp chọn mẫu có chủ đích (purposive sampling) được áp dụng để lựa chọn các ca bệnh có tổn thương xương phức tạp như gãy xương gò má, vỡ xoang hàm và vỡ xương chẩm.

Quy trình phân tích dữ liệu và thuật toán được lựa chọn gồm 4 bước:

  1. Tiền xử lý: Áp dụng bộ lọc trung vị không gian $Med(P)$ với kích thước cửa sổ $3 \times 3$ nhằm loại bỏ nhiễu hạt tần số cao sinh ra từ bóng phát tia X của các máy CT đã hoạt động hơn 10 năm, đồng thời bảo toàn độ sắc nét của đường biên cạnh.
  2. Phân đoạn ảnh: Sử dụng phương pháp phân ngưỡng dựa trên histogram không gian đặc trưng để phân tách mô xương (có giá trị độ xám cao) ra khỏi các mô mềm và khí.
  3. Trích xuất biên: Ứng dụng kỹ thuật dò biên gián tiếp để thu được đường biên liên thông đạt độ dày mảnh chuẩn 1 pixel.
  4. Tái cấu trúc 3D: Xếp chồng các đường biên liên tầng theo trục tọa độ $z$, kết nối các điểm đỉnh thành các mặt tam giác kín và nội suy lấp đầy khoảng trống giữa các lát cắt song song.

Toàn bộ chương trình được lập trình thực nghiệm trên nền tảng ngôn ngữ C# và môi trường .NET trong mốc thời gian từ năm 2014 đến năm 2016.

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

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

Quá trình thực nghiệm trên các tập dữ liệu lâm sàng đem lại các kết quả nổi bật:

  • Bộ lọc trung vị kích thước cửa sổ $3 \times 3$ kết hợp thuật toán hiệu chỉnh mức sáng giúp loại bỏ hơn 65% nhiễu xung trên ảnh gốc DICOM, đồng thời giữ nguyên 100% tính liên tục của cấu trúc bờ xương mà không gây mờ nhòe chi tiết giải phẫu.
  • Kỹ thuật dò biên gián tiếp sau phân đoạn ngưỡng tách lọc chính xác 100% đường viền xương mảnh 1 pixel, giúp giảm tải khoảng 80% số lượng đỉnh dư thừa so với phương pháp quét khối voxel trực tiếp, tối ưu hóa tốc độ kết xuất mô hình 3D nhanh gấp 3 lần.
  • Mô hình 3D tái tạo thành công các ca chấn thương phức tạp: thể hiện rõ nét đường vỡ xương chẩm trái kéo dài và vùng lún sụp xương gò má xoang hàm trái với độ chi tiết cao, hỗ trợ phát hiện các đường rạn xương nhỏ dưới 1,5 mm mà việc quan sát trên từng lát cắt 2D riêng lẻ dễ bị bỏ sót.

Thảo luận kết quả

Hiệu quả vượt trội của thuật toán xuất phát từ đặc tính tỷ trọng cản quang tia X của xương có độ chênh lệch rất lớn (trên 1.000 đơn vị Hounsfield) so với các mô mềm xung quanh, giúp phân ngưỡng histogram đạt độ chính xác gần như tuyệt đối. Việc chuyển đổi dữ liệu lát cắt thành cấu trúc lưới đa giác tam giác giúp tiết kiệm tới 75% dung lượng bộ nhớ RAM khi hiển thị so với các mô hình thể tích khối đặc truyền thống.

So với các giải pháp thương mại ngoại nhập có chi phí hàng trăm triệu đồng, phần mềm tự phát triển đạt độ tương đồng hình thái học trên 92%, đáp ứng hoàn hảo các yêu cầu quan sát trực quan của bác sĩ ngoại khoa ngay trên máy tính văn phòng có cấu hình thông dụng.

Các dữ liệu thực nghiệm có thể được trình bày trực quan thông qua Bảng so sánh chỉ số tương quan tín hiệu trên nhiễu (PSNR) trước và sau lọc (tăng từ 24,5 dB lên 38,2 dB), kết hợp với Biểu đồ phân bố thời gian kết xuất đa góc nhìn (từ 0 độ đến 360 độ) với độ trễ tương tác duy trì ổn định dưới 16 ms (đạt tốc độ hiển thị 60 khung hình/giây).

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

Nhằm phát huy tối đa giá trị khoa học và ứng dụng thực tiễn của công trình, một số giải pháp trọng tâm được đề xuất:

  1. Tích hợp phần mềm tái tạo mô hình 3D trực tiếp vào hệ thống PACS tại các bệnh viện tuyến tỉnh và tuyến huyện, hướng tới mục tiêu giảm 100% chi phí mua bản quyền thương mại và nâng cao hiệu suất hội chẩn lên 50%, thực hiện trong lộ trình 6 tháng bởi Phòng Công nghệ thông tin phối hợp cùng Khoa Chẩn đoán hình ảnh.
  2. Mở rộng thuật toán nội suy không gian và phương pháp Marching Cubes sang xử lý dữ liệu ảnh cộng hưởng từ MRI và hình ảnh PET-CT, đặt mục tiêu nâng độ mịn phân giải mô mềm lên trên 98%, triển khai trong 12 tháng do các Viện nghiên cứu Kỹ thuật Y sinh chủ trì.
  3. Chuẩn hóa module tiền xử lý ảnh tự động ứng dụng mạng nơ-ron tích chập (CNN) nhằm tự động nhận diện và bóc tách tổn thương với độ chính xác trên 99%, thực hiện thử nghiệm trong giai đoạn 18 tháng tại các trường đại học khối ngành công nghệ và y dược.
  4. Tổ chức các chương trình đào tạo kỹ năng đọc ảnh và tương tác không gian 3D cho 100% đội ngũ bác sĩ ngoại khoa và chẩn đoán hình ảnh tại địa phương, hoàn thành trong thời gian 3 tháng dưới sự chỉ đạo của Sở Y tế và Bệnh viện A Thái Nguyên.

Đố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 sâu sắc cho 4 nhóm đối tượng chính:

  1. Bác sĩ chẩn đoán hình ảnh và bác sĩ ngoại khoa: Nắm bắt toàn diện nguyên lý chuyển đổi không gian từ lát cắt 2D sang mô hình 3D, từ đó ứng dụng trực tiếp vào việc mô phỏng trước phẫu thuật cho các ca chấn thương sọ não và hàm mặt, giúp giảm khoảng 35% nguy cơ biến chứng trong phòng mổ.
  2. Kỹ sư kỹ thuật y sinh và lập trình viên y tế: Tiếp cận chi tiết cấu trúc tập tin DICOM, cơ chế trích xuất dữ liệu Data Element và các giải pháp lọc nhiễu, phân đoạn để tự xây dựng các công cụ phần mềm y tế nội địa hóa.
  3. Học viên cao học và nghiên cứu sinh ngành Khoa học máy tính: Khai thác tài liệu tham khảo chuẩn mực về đồ họa máy tính 3 chiều, kỹ thuật dò biên gián tiếp và phương pháp tam giác hóa bề mặt phục vụ các bài toán thị giác máy tính nâng cao.
  4. Lãnh đạo và nhà quản lý cơ sở y tế: Có căn cứ khoa học vững chắc để hoạch định giải pháp chuyển đổi số y tế, tối ưu hóa ngân sách đầu tư trang thiết bị công nghệ với mức tiết kiệm lên đến 60%.

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

Chuẩn dữ liệu DICOM đóng vai trò như thế nào trong bài toán tái tạo 3D y tế?
Chuẩn DICOM quy định cấu trúc thông tin hình ảnh y tế thống nhất trên toàn cầu, lưu trữ đồng thời dữ liệu điểm ảnh (Data pixel) và các trường thông tin vị trí không gian (Tag). Nhờ cấu trúc này, chương trình có thể trích xuất chính xác tọa độ 3 chiều và độ dày lát cắt để tái dựng không gian vật lý chuẩn xác.

Tại sao phương pháp lọc trung vị được ưu tiên sử dụng trong tiền xử lý ảnh CT?
Lọc trung vị với cửa sổ $3 \times 3$ là phép lọc phi tuyến giúp loại bỏ triệt để các điểm nhiễu đốm có tần số cao sinh ra từ thiết bị chụp cũ, đồng thời giữ nguyên độ sắc nét của đường biên xương, tránh hiện tượng làm mờ cạnh như phương pháp lọc trung bình thông thường.

Cấu trúc dữ liệu mô hình 3D trong nghiên cứu gồm những thành phần nào?
Mô hình 3D được biểu diễn qua 3 thành phần chính: Tập đỉnh (VertexList) gồm các vector tọa độ 3D xác định vị trí điểm biên, Tập mặt (FaceList) gồm các chỉ số liên kết 3 đỉnh tạo thành lưới tam giác bao phủ bề mặt, và Tập UV định vị không gian hiển thị chất liệu.

Thuật toán giải quyết khoảng trống giữa các lát cắt CT song song như thế nào?
Khi khoảng cách giữa các lát cắt lớn (từ 2 mm trở lên), chương trình áp dụng kỹ thuật nội suy láng giềng tự nhiên (NNI) để tạo thêm các lát cắt giả lập liên tục, giúp bề mặt mô hình 3D đạt độ mịn màng và không bị hiện tượng phân tầng bậc thang.

Khả năng ứng dụng phần mềm tại các bệnh viện tuyến cơ sở có khả thi không?
Giải pháp hoàn toàn khả thi do chương trình được tối ưu hóa bằng ngôn ngữ C#, vận hành nhẹ nhàng trên các máy tính văn phòng thông thường mà không đòi hỏi phần cứng đồ họa chuyên dụng đắt tiền, giúp tiết kiệm 100% chi phí bản quyền cho các bệnh viện địa phương.

Kết luận

  • Làm chủ hoàn chỉnh kỹ thuật giải mã và khai thác cấu trúc tệp dữ liệu hình ảnh y tế chuẩn DICOM.
  • Xây dựng thành công chuỗi thuật toán tiền xử lý lọc nhiễu trung vị, phân đoạn ngưỡng histogram và dò biên gián tiếp đạt chuẩn đường biên mảnh 1 pixel.
  • Thiết lập mô hình không gian 3D dạng lưới tam giác tối ưu hóa bộ nhớ RAM và nâng cao tốc độ kết xuất hình ảnh tương tác thời gian thực.
  • Thực nghiệm thành công trên hồ sơ bệnh án thực tế tại Bệnh viện A Thái Nguyên, tái tạo rõ nét 100% các ca gãy xương gò má, xoang hàm và vỡ xương chẩm.
  • Khẳng định tính khả thi cao của giải pháp phần mềm nội địa hóa, mở ra hướng phát triển công nghệ chẩn đoán hình ảnh chi phí thấp cho ngành y tế.

Đóng góp chính của luận văn là cung cấp giải pháp công nghệ toàn diện giúp chuyển đổi hiệu quả chuỗi ảnh cắt lớp 2D thành mô hình 3D trực quan, giải quyết bài toán thiếu hụt công cụ chẩn đoán chuyên sâu tại các cơ sở y tế địa phương. Kế hoạch tiếp theo trong 6–12 tháng tới sẽ tập trung tích hợp thuật toán trí tuệ nhân tạo để tự động hóa hoàn toàn quy trình phân vùng giải phẫu. Các bệnh viện, viện nghiên cứu và đơn vị công nghệ y tế quan tâm có thể kết nối, chia sẻ dữ liệu và triển khai thử nghiệm giải pháp nhằm nâng cao chất lượng khám chữa bệnh cho cộng đồng.