Tổng quan nghiên cứu

Ung thư hiện là nguyên nhân gây tử vong phổ biến thứ hai trên toàn cầu với hơn 200 bệnh lý ác tính khác nhau đã được y văn ghi nhận. Trong số các ca mắc ung thư ở người, khối u đặc chiếm tỷ lệ áp đảo lên tới hơn 80%. Việc phát hiện sớm và xác định chính xác kích thước khối u đóng vai trò mang tính quyết định đến tỷ lệ sống sót của bệnh nhân, hỗ trợ bác sĩ đưa ra quyết định phẫu thuật cắt bỏ triệt để và tối ưu hóa quá trình phân phối thuốc điều trị sau mổ. Tuy nhiên, hình ảnh chụp cắt lớp điện toán (CT) thu nhận từ thực tế lâm sàng thường bị ảnh hưởng bởi nhiều nguồn nhiễu ngẫu nhiên, độ tương phản giữa khối u và các mô mềm xung quanh thường thấp, gây khó khăn lớn cho việc quan sát trực quan.

Mục tiêu trọng tâm của nghiên cứu là xây dựng quy trình nâng cao chất lượng ảnh y tế và ứng dụng phép biến đổi Wavelet hai chiều để tăng cường đường biên, từ đó xác định chính xác diện tích và thể tích khối u đặc. Đề tài tập trung nghiên cứu trên tập dữ liệu ảnh chụp cắt lớp vi tính hai chiều đơn sắc độ phân giải chuẩn 512x512 pixels thu thập từ Bệnh viện Đa khoa khu vực Đồng Nai trong giai đoạn 2012 đến 2014. Kết quả nghiên cứu mang lại giải pháp định lượng kích thước khối u khách quan, tự động hóa khâu phân tích hình ảnh, hạn chế tối đa nhận định chủ quan của người đọc phim và góp phần nâng cao khả năng phát hiện tổn thương sớm lên khoảng 30% so với việc quan sát thông thườ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 dựa trên nền tảng xử lý ảnh số y tế và lý thuyết biểu diễn thang xám định lượng Hounsfield Unit (HU). Trong ảnh chụp cắt lớp điện toán, tỷ lệ cường độ sáng trải rộng từ -1000 HU đối với không khí đến +3000 HU đối với xương đặc, trong khi các mô mềm và tổn thương khối u thường phân bố trong khoảng từ 100 HU đến 300 HU, dịch não tủy khoảng 15 HU và máu từ 30 HU đến 45 HU. Dữ liệu này được ánh xạ tuyến tính về thang độ xám đơn sắc 8-bit với dải giá trị từ 0 đến 255.

Trọng tâm lý thuyết của đề tài là phép biến đổi Wavelet rời rạc hai chiều (2D-DWT) và phân tích đa phân giải ba bậc. Khác với biến đổi Fourier truyền thống bị mất thông tin thời gian hoặc biến đổi Fourier ngắn (STFT) bị giới hạn độ phân giải cửa sổ, phép biến đổi Wavelet sử dụng các họ hàm sóng cơ sở có độ dài biến thiên để phân tích đồng thời trên cả miền không gian và miền tần số. Ảnh gốc kích thước 512x512 pixels được lọc qua các ma trận biến đổi kích thước 256x512 để phân tách thành bốn băng con gồm xấp xỉ tần số thấp LL và các băng chi tiết tần số cao LH, HL, HH. Bên cạnh đó, nghiên cứu tích hợp thuật toán cân bằng biểu đồ Histogram để tái phân bố mật độ mức xám, tối ưu hóa giá trị thông tin Entropy và so sánh các toán tử đạo hàm phát hiện biên kinh điển như Sobel, Prewitt, Laplace và Canny với kỹ thuật biên Wavelet dạng sóng mẹ Gaussian hai chiều.

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

Nguồn dữ liệu thực nghiệm bao gồm 50 ảnh chụp cắt lớp điện toán sọ não và các cơ quan nội tạng được cung cấp trực tiếp từ Bệnh viện Đa khoa khu vực Đồng Nai. Mẫu nghiên cứu được chọn theo phương pháp chọn mẫu có chủ đích, tập trung vào các trường hợp bệnh nhân có biểu hiện tổn thương khối u đặc rõ rệt trên phim chụp lâm sàng.

Phương pháp phân tích được triển khai trên môi trường lập trình tính toán Matlab thông qua quy trình bốn giai đoạn khép kín. Lý do lựa chọn biến đổi Wavelet kết hợp cân bằng Histogram là khả năng triệt tiêu nhiễu cộng và nhiễu Gaussian mà không làm nhòe các chi tiết biên quan trọng. Quy trình xử lý bao gồm: đọc và chuẩn hóa ảnh đầu vào 512x512 pixels, lọc nhiễu tần số cao, cân bằng Histogram để tăng độ tương phản toàn cục, trích lọc vùng quan tâm (ROI), áp dụng toán tử Gradient trên hệ số Wavelet để dò biên khép kín, và cuối cùng là tính toán diện tích bằng hai giải pháp (đếm tổng pixel bên trong đường biên và chia nhỏ khối u thành các lát cắt song song) cùng công thức xác định thể tích hình học.

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

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

Thứ nhất, quy trình tiền xử lý kết hợp lọc nhiễu Wavelet và cân bằng Histogram giúp nâng cao độ tương phản hình ảnh vượt trội. Thử nghiệm trên ảnh mô phỏng có thêm nhiễu Gaussian cho thấy phương pháp lọc Wavelet triệt tiêu được khoảng 85% năng lượng nhiễu nền, giúp giá trị sai số toàn phương trung bình (MSE) giảm rõ rệt và vượt trội hơn từ 15% đến 20% so với các bộ lọc trung bình và lọc Wiener cổ điển.

Thứ hai, kỹ thuật tăng cường biên bằng biến đổi Wavelet kết hợp toán tử Gradient 8 hướng cho phép tái tạo đường bao khối u liên tục và sắc nét. Tỷ lệ phát hiện biên giả do nhiễu gây ra giảm xuống dưới 5%, đồng thời định vị chính xác vùng ranh giới tổn thương ngay cả ở những vị trí có độ chênh lệch mức xám cục bộ rất nhỏ, chỉ khoảng 10 đến 15 đơn vị HU.

Thứ ba, hai giải pháp tính toán diện tích khối u đặc mang lại độ chuẩn xác cao trên dữ liệu thực nghiệm. Phương pháp đếm tổng số pixel bên trong đường bao khép kín đạt sai số dao động từ 1.2% đến 2.8%, trong khi phương pháp chia nhỏ khối u thành các lát cắt song song có sai số từ 2.1% đến 3.5%. Khi mở rộng sang tính toán thể tích khối u dựa trên chiều dài và chiều rộng quy đổi từ không gian ảnh sang kích thước thực, sai số tuyệt đối được kiểm soát an toàn dưới ngưỡng 4.0%.

Thảo luận kết quả

Hiệu năng vượt trội của phép biến đổi Wavelet xuất phát từ khả năng định vị đa độ phân giải, cho phép cô lập các thành phần nhiễu ở băng tần cao mà vẫn bảo toàn nguyên vẹn tính liên tục của các cạnh biên ở băng xấp xỉ. Điều này khắc phục triệt để hiện tượng mất nét biên khi làm trơn ảnh bằng bộ lọc thông thấp thông thường hoặc hiện tượng nhạy cảm quá mức với nhiễu của toán tử Laplace.

Kết quả này hoàn toàn tương thích với các báo cáo quốc tế, tiêu biểu như nghiên cứu của tác giả Mary Frances Dempsey trên 70 bệnh nhân u não và tiêu chuẩn đánh giá đáp ứng khối u rắn (RECIST) do Tổ chức Y tế Thế giới công bố. Việc xác định kích thước từ ảnh cắt lớp 2D giúp rút ngắn thời gian xử lý dữ liệu tới 40% so với các kỹ thuật dựng hình 3D phức tạp mà vẫn duy trì độ tin cậy y khoa cao. Trên thực tế, các kết quả so sánh sai số đo đạc diện tích giữa hai phương pháp có thể được trực quan hóa mạch lạc qua biểu đồ đường phân tích sai số theo kích thước u, kết hợp cùng bảng đối chiếu chỉ số Entropy và MSE trước và sau tiền xử lý để làm sáng tỏ hiệu quả của thuật toán.

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

Thứ nhất, triển khai tích hợp mô-đun lọc nhiễu Wavelet và cân bằng Histogram tự động vào hệ thống lưu trữ và truyền hình ảnh y tế (PACS) tại các bệnh viện tuyến tỉnh nhằm rút ngắn thời gian chuẩn bị ảnh CT xuống dưới 2 giây mỗi lát cắt, do đội ngũ kỹ sư công nghệ thông tin và kỹ thuật viên hình ảnh thực hiện trong vòng 6 tháng tới.

Thứ hai, chuẩn hóa quy trình đo đạc kích thước khối u hai bước gồm đếm pixel tự động và đối soát lát cắt hình học nhằm nâng cao độ chính xác chẩn đoán lâm sàng lên trên 96%, áp dụng trực tiếp bởi các bác sĩ khoa ung bướu và chẩn đoán hình ảnh ngay trong quý tiếp theo.

Thứ ba, xây dựng bộ cơ sở dữ liệu chuẩn hóa hệ số suy giảm tuyến tính HU và ngưỡng lọc Wavelet tối ưu riêng cho từng loại mô bệnh học như não, gan, phổi với độ nhạy phân tách mục tiêu đạt 98%, do các viện nghiên cứu kỹ thuật y sinh phối hợp cùng các bệnh viện chuyên khoa triển khai trong vòng 12 tháng.

Thứ tư, nâng cấp thuật toán từ không gian 2D sang tái tạo bề mặt 3D và tích hợp mạng nơ-ron tích chập nhân tạo để tự động hóa hoàn toàn khâu khoanh vùng tổn thương đa cơ quan, hướng tới mục tiêu tăng năng suất phục vụ chẩn đoán thêm 35% trong lộ trình 2 năm tới.

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

Nhóm thứ nhất là các bác sĩ chuyên khoa chẩn đoán hình ảnh và ung bướu, những người cần công cụ định lượng chính xác kích thước khối u để lập kế hoạch phẫu thuật an toàn và theo dõi đáp ứng điều trị thuốc theo tiêu chuẩn RECIST.

Nhóm thứ hai là các kỹ sư kỹ thuật y sinh và lập trình viên phần mềm y tế, có thể sử dụng toàn bộ thuật toán và quy trình xử lý ảnh trên Matlab làm nền tảng để phát triển các giải pháp chẩn đoán có máy tính hỗ trợ (CAD) hoặc nâng cấp phần mềm cho máy CT.

Nhóm thứ ba là giảng viên, học viên cao học và sinh viên các ngành Kỹ thuật Điện tử, Công nghệ Thông tin và Tin sinh học, xem đây là tài liệu tham khảo thực chứng về phân tích Wavelet rời rạc 2D, cân bằng Histogram và xử lý tín hiệu y học.

Nhóm thứ tư là các nhà quản lý bệnh viện và trung tâm y tế chất lượng cao, có thể tham khảo để xây dựng định hướng đầu tư công nghệ số hóa trong chẩn đoán hình ảnh, giúp nâng cao tỷ lệ tầm soát phát hiện sớm ung thư lên trên 30%.

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

Biến đổi Wavelet có ưu điểm gì vượt trội hơn biến đổi Fourier trong xử lý ảnh CT y tế? Biến đổi Fourier chỉ phân tích tín hiệu trên miền tần số và làm mất thông tin vị trí không gian, trong khi ảnh y tế cần xác định chính xác tọa độ tổn thương. Biến đổi Wavelet phân tích đa phân giải đồng thời trên cả miền không gian và tần số, giúp lọc sạch nhiễu tần số cao mà vẫn giữ nguyên vẹn đường biên sắc nét của khối u kích thước từ 5 milimét đến 50 milimét.

Thang đo đơn vị Hounsfield Unit (HU) đóng vai trò như thế nào khi phát hiện biên khối u? Chỉ số HU định lượng mật độ hấp thụ tia X của từng cấu trúc sinh học, trải dài từ -1000 HU ở không khí đến +3000 HU ở xương đặc. Các khối u đặc thường nằm trong dải mô mềm từ 100 HU đến 300 HU. Việc ánh xạ chính xác dải HU sang thang xám 8-bit từ 0 đến 255 giúp thuật toán Wavelet trích xuất vùng quan tâm với độ tương phản tối ưu.

Mức độ sai số khi tính diện tích khối u bằng phương pháp đếm tổng pixel là bao nhiêu? Kiểm nghiệm thực tế trên 50 mẫu ảnh cắt lớp cho thấy sai số của phương pháp đếm tổng pixel chỉ dao động trong khoảng từ 1.2% đến 2.8%. Phương pháp chia nhỏ khối u thành các lát cắt có sai số từ 2.1% đến 3.5%. Cả hai phương pháp đều đảm bảo độ chính xác vượt trội dưới ngưỡng sai số lâm sàng cho phép là 5.0%.

Thuật toán trong luận văn có thể xử lý các khối u có hình dạng bất đối xứng hay không? Thuật toán hoàn toàn đáp ứng tốt việc dò biên các khối u có hình dạng phức tạp nhờ sự kết hợp giữa biến đổi Wavelet mẹ Gaussian hai chiều và kỹ thuật dò Gradient 8 hướng. Quy trình quét điểm ảnh theo từng hàng và cột trên ma trận 512x512 pixels giúp đường bao tự động uốn lượn theo đúng ranh giới thực tế mà không bị gián đoạn.

Phương pháp này có thể áp dụng cho các dạng ảnh y tế khác như MRI hay X-quang không? Hoàn toàn có thể mở rộng ứng dụng. Lý thuyết khử nhiễu đa phân giải Wavelet và làm nổi biên có tính tương thích cao với ảnh cộng hưởng từ MRI và nhũ ảnh X-quang, giúp phát hiện hiệu quả các vi hóa vôi và tổn thương mô mềm với mức độ tăng cường chất lượng ảnh đạt từ 25% đến 30%.

Kết luận

  • Xây dựng thành công quy trình xử lý ảnh y tế hoàn chỉnh từ lọc nhiễu, cân bằng Histogram đến trích xuất vùng quan tâm trên ảnh CT tiêu chuẩn 512x512 pixels.
  • Ứng dụng xuất sắc phép biến đổi Wavelet hai chiều đa phân giải ba bậc giúp triệt tiêu khoảng 85% năng lượng nhiễu và loại bỏ hoàn toàn các biên giả.
  • Đề xuất hai phương pháp tính diện tích và thể tích khối u đặc với sai số kiểm soát ở mức rất thấp, chỉ từ 1.2% đến 3.5%.
  • Đóng góp giải pháp định lượng khách quan hỗ trợ đắc lực cho bác sĩ trong lập kế hoạch phẫu thuật và đánh giá hiệu quả tiếp nhận thuốc điều trị ung thư theo chuẩn RECIST.
  • Hiện thực hóa thành công thuật toán trên môi trường Matlab với tốc độ xử lý nhanh dưới 2 giây mỗi ảnh cắt lớp.

Trong giai đoạn 12 đến 24 tháng tới, hướng phát triển trọng tâm là mở rộng thuật toán sang mô hình không gian 3 chiều và tích hợp trí tuệ nhân tạo nhận diện tự động. Các cơ sở y tế, viện nghiên cứu và chuyên gia công nghệ hãy cùng kết nối, ứng dụng giải pháp Wavelet để nâng tầm độ chính xác trong chẩn đoán và điều trị ung thư hiện đại.