Tổng quan nghiên cứu

Trong kỹ thuật chẩn đoán hình ảnh y tế hiện đại, ung thư và các khối u mô mềm đang là nguyên nhân gây tử vong hàng đầu với hàng triệu ca mắc mới mỗi năm trên toàn cầu. Các kỹ thuật truyền thống như chụp cắt lớp vi tính hay chụp cộng hưởng từ mang lại hình ảnh chi tiết nhưng lại tồn tại nhược điểm về chi phí cao và nguy cơ phơi nhiễm bức xạ ion hóa. Trong khi đó, siêu âm 2D chế độ B-mode thông thường có vận tốc truyền sóng trung bình trong mô mềm khoảng 1540 m/s với mức phát năng lượng an toàn từ 1 MW đến 10 MW. Tuy nhiên, hình ảnh B-mode phản hồi âm chỉ mang tính chất định tính, chất lượng ảnh phụ thuộc lớn vào tay nghề bác sĩ và rất khó phát hiện các tổn thương có kích thước nhỏ hơn bước sóng âm.

Vấn đề cốt lõi đặt ra là làm thế nào để xây dựng một hệ thống siêu âm cắt lớp tán xạ ngược định lượng vừa đảm bảo độ tương phản cao, vừa tối ưu hóa thời gian tính toán để ứng dụng trong môi trường lâm sàng thực tế. Mục tiêu cụ thể của công trình là nghiên cứu, cải tiến thuật toán lặp vi phân Born bằng phương pháp đặt ngưỡng áp suất tán xạ, giúp loại bỏ các phương trình thu thập dữ liệu dư thừa, từ đó tăng tốc độ hội tụ và giảm độ trễ tái tạo ảnh.

Nghiên cứu được triển khai thực nghiệm mô phỏng tại Trường Đại học Công nghệ, Đại học Quốc gia Hà Nội vào năm 2016. Luận văn tập trung khảo sát mô hình sóng siêu âm tần số 0,46 MHz tác động lên đối tượng hình trụ có đường kính 10 mm với độ chênh lệch vận tốc truyền sóng 2% trong môi trường đồng nhất. Kết quả nghiên cứu mang ý nghĩa thực tiễn lớn khi rút ngắn thời gian tái tạo ảnh trung bình từ 110,7 giây xuống còn 104,4 giây, đạt mức tiết kiệm thời gian 5,7% trên lưới phân giải nhỏ và mở ra tiềm năng tăng tốc gấp nhiều lần đối với các ma trận ảnh y sinh kích thước lớn.

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 dựa trên nền tảng lý thuyết truyền sóng âm toàn phần và bài toán tán xạ ngược phi tuyến trong môi trường không đồng nhất. Trọng tâm lý thuyết là hai mô hình lặp kinh điển: phương pháp lặp Born và phương pháp lặp vi phân Born. Trong đó, phương pháp lặp vi phân Born được đánh giá vượt trội nhờ khả năng cập nhật lại hàm Green nền sau mỗi bước lặp, mang lại tốc độ hội tụ phi tuyến nhanh hơn đáng kể so với phương pháp lặp Born thông thường.

Khung lý thuyết vận hành dựa trên các khái niệm toán học và vật lý cốt lõi:

  1. Hàm mục tiêu biểu diễn sự phân bố độ tương phản âm học giữa môi trường đối tượng và môi trường nền có vận tốc truyền sóng 1480 m/s trong nước.
  2. Vùng quan tâm được chia lưới không gian để rời rạc hóa phương trình tích phân sóng thông qua phương pháp mô-men.
  3. Áp suất tán xạ thu nhận tại các đầu dò sau khi sóng tới tương tác với vật thể.
  4. Toán tử nghịch đảo ma trận giả Moore - Penrose kết hợp chuẩn hóa Tikhonov nhằm giải quyết hệ phương trình tuyến tính quá xác định khi số lượng điểm đo lớn hơn nhiều so với số lượng biến ảnh.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng từ cấu hình hệ thống đo siêu âm tán xạ vòng tròn khép kín 360 độ. Cấu hình phần cứng ảo bao gồm 9 vị trí máy phát sóng siêu âm đặt cách tâm 104 mm, kết hợp với mảng 121 máy thu đặt cách tâm 91 mm quét trên cung góc từ âm 60 độ đến dương 60 độ với bước dịch chuyển 1 độ. Tổng số phương trình đo thu thập được trong một chu kỳ quét là 1089 điểm dữ liệu áp suất tán xạ.

Vùng diện tích quan tâm được rời rạc hóa với cỡ mẫu không gian gồm 196 pixel trên ma trận vuông 14 nhân 14 ô. Phương pháp chọn mẫu dữ liệu áp dụng kỹ thuật quét đa góc đối xứng nhằm triệt tiêu hiện tượng lệch pha tín hiệu. Phương pháp phân tích dữ liệu được lựa chọn là thuật toán lặp vi phân Born tích hợp bộ lọc quyết định ngưỡng áp suất tán xạ tuyệt đối. Lý do lựa chọn giải pháp này là vì trong hệ phương trình đo có kích thước 1089 hàng và 196 cột, phần lớn các giá trị áp suất tán xạ biên độ thấp chỉ đóng góp thành phần nhiễu mà không mang thông tin cấu trúc. Quá trình mô phỏng được thực hiện qua 10 vòng lặp liên tiếp trong môi trường không nhiễu và môi trường có bổ sung nhiễu Gauss 5% để đánh giá toàn diện tính ổn định của hệ thống.

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

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

Quá trình mô phỏng và phân tích định lượng đã mang lại 4 phát hiện quan trọng:

  1. Việc áp dụng ngưỡng áp suất tán xạ giúp cắt giảm trực tiếp kích thước ma trận đo mà không làm suy giảm chất lượng khôi phục đối tượng. Giá trị ngưỡng tối ưu được xác định là A bằng 7 nhân 10 mũ âm 4, tại đó hệ thống loại bỏ được hàng trăm phương trình đo kém ý nghĩa.
  2. Thời gian tái tạo ảnh trung bình của thuật toán cải tiến giảm từ 110,7 giây trong trường hợp không đặt ngưỡng xuống còn 104,4 giây khi đặt ngưỡng tối ưu, tương ứng với hiệu suất cải thiện tốc độ đạt 5,7%.
  3. Trong kịch bản kiểm thử có thêm 5% nhiễu Gauss với hơn 10 lần chạy độc lập, thuật toán lặp vi phân Born có ngưỡng vẫn duy trì sai số tái tạo tương đương trường hợp không đặt ngưỡng, khôi phục chính xác biên dạng hình trụ đường kính 10 mm ở vòng lặp thứ 10.
  4. Khi tiếp tục tăng giá trị ngưỡng lên mức 10 mũ âm 3, sai số tái tạo ảnh tăng vọt trên 25% do ma trận bị khuyết thiếu các thành phần tán xạ quan trọng, chứng minh ngưỡng giới hạn có biên độ tối ưu xác định.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp đạt hiệu quả cao là việc triệt tiêu sự dư thừa thông tin trong bài toán ngược. Khi số lượng máy đo phát sinh 1089 phép đo cho 196 ẩn số, các đầu dò ở vị trí góc xa chỉ ghi nhận tín hiệu tán xạ cực tiểu gần bằng không. Việc duy trì các hàng dữ liệu này trong ma trận tính toán nghịch đảo Moore - Penrose làm gia tăng chi phí xử lý bộ nhớ theo hàm bậc ba mà không gia tăng độ phân giải hình ảnh.

Dữ liệu mô phỏng được thể hiện rõ nét thông qua bảng so sánh kích thước ma trận suy giảm tương ứng với từng mức ngưỡng từ 10 mũ âm 4 đến 10 mũ âm 3 và biểu đồ đường biểu diễn sự hội tụ của sai số tương đối qua 10 bước lặp. Hình ảnh không gian 3 chiều của hàm mục tiêu sau lần lặp thứ 10 cho thấy độ tương phản đạt mức 2% đồng nhất, tiệm cận hoàn hảo với mô hình gốc lý tưởng. So với các nghiên cứu lặp Born cổ điển có số vòng lặp kéo dài từ 20 đến 30 bước, phương pháp đề xuất hội tụ ngay từ bước lặp thứ 4, khẳng định tính khả thi vượt trội trong xử lý dữ liệu lớn.

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

Nhằm hoàn thiện kỹ thuật và đưa giải pháp vào ứng dụng thực tiễn trong ngành thiết bị y tế, 4 khuyến nghị trọng tâm được đề xuất:

  1. Nâng cấp thuật toán xác định ngưỡng từ cơ chế tĩnh sang cơ chế thích nghi động theo từng vòng lặp, hướng tới mục tiêu giảm thêm 15% thời gian tính toán tổng thể trong vòng 6 tháng tới, do các nhóm nghiên cứu phần mềm xử lý tín hiệu chủ trì thực hiện.
  2. Mở rộng quy mô lưới phân giải không gian từ kích thước 14 nhân 14 pixel lên 64 nhân 64 pixel và 128 nhân 128 pixel nhằm nâng cao độ mịn hình ảnh đạt mức dưới 0,5 mm, với lộ trình thực nghiệm 12 tháng do các nghiên cứu sinh chuyên ngành kỹ thuật điện tử đảm nhận.
  3. Tích hợp công nghệ tính toán song song trên phần cứng xử lý đồ họa chuyên dụng và tối ưu hóa hàm toán học khôi phục Gradient phi tuyến, đặt mục tiêu đạt tốc độ xử lý thời gian thực dưới 3 giây cho mỗi lát cắt trong vòng 18 tháng.
  4. Triển khai thử nghiệm lâm sàng trên mô hình vật thể mô phỏng mô người thực tế với dải tần phát từ 1 MHz đến 5 MHz, do các viện trang thiết bị y tế phối hợp cùng các bệnh viện chuyên khoa thực hiện trong giai đoạn 24 tháng.

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

  1. Học viên cao học và nghiên cứu sinh ngành Điện tử - Viễn thông, Kỹ thuật Y sinh: Tài liệu cung cấp cơ sở toán học chi tiết về bài toán tán xạ ngược, phương pháp lặp vi phân Born và kỹ thuật tối ưu hóa ma trận đo trong xử lý ảnh y sinh.
  2. Kỹ sư phát triển phần mềm thiết bị chẩn đoán hình ảnh: Tiếp cận mã nguồn mô phỏng chi tiết bằng phần mềm toán học, áp dụng trực tiếp thuật toán đặt ngưỡng để giải quyết bài toán nghẽn cổ chai tốc độ xử lý dữ liệu mảng đầu dò.
  3. Chuyên gia thiết kế phần cứng siêu âm: Tham khảo cấu hình mảng thu phát 9 nguồn phát và 121 cảm biến thu để tối ưu hóa vị trí đặt đầu dò trên các hệ thống chụp cắt lớp siêu âm thực nghiệm.
  4. Bác sĩ và nhà nghiên cứu ung bướu: Nắm bắt nguyên lý định lượng mô của siêu âm cắt lớp, hỗ trợ đánh giá tiềm năng ứng dụng công nghệ mới trong phát hiện sớm khối u tuyến vú kích thước nhỏ dưới 5 mm.

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

Phương pháp lặp vi phân Born cải tiến có điểm gì khác biệt so với kỹ thuật truyền thống?

Thuật toán truyền thống giữ nguyên toàn bộ 1089 phương trình đo từ 121 cảm biến thu, dẫn đến khối lượng tính toán ma trận rất lớn. Phương pháp cải tiến bổ sung khâu so sánh ngưỡng biên độ áp suất tán xạ, tự động loại bỏ các tín hiệu dưới mức 7 nhân 10 mũ âm 4, giúp giảm kích thước ma trận và tăng tốc độ xử lý 5,7%.

Vì sao giá trị ngưỡng A bằng 7 nhân 10 mũ âm 4 được chọn làm điểm tối ưu?

Qua thực nghiệm kiểm thử trên dải ngưỡng từ 10 mũ âm 4 đến 10 mũ âm 3, mức 7 nhân 10 mũ âm 4 mang lại thời gian tính toán ngắn nhất đạt 104,4 giây và đảm bảo sai số tái tạo ảnh nhỏ nhất, trong khi các mức ngưỡng cao hơn 10 mũ âm 3 làm mất mát dữ liệu và gây sai lệch cấu trúc đối tượng.

Thuật toán có duy trì được độ chính xác khi hệ thống bị ảnh hưởng bởi nhiễu môi trường không?

Thuật toán hoạt động rất ổn định trong môi trường thực nghiệm có 5% nhiễu Gauss. Sau 10 lần lặp kiểm thử độc lập, hình dạng vật thể đường kính 10 mm vẫn được tái tạo rõ nét, tương đương với trường hợp không đặt ngưỡng, chứng minh khả năng lọc nhiễu tốt của bộ nghịch đảo ma trận giả.

Hệ thống thu phát tín hiệu siêu âm trong mô phỏng được cấu hình như thế nào?

Mô hình thực nghiệm ảo bố trí 9 vị trí máy phát sóng siêu âm đặt cách tâm 104 mm cách đều nhau 45 độ, phối hợp cùng mảng 121 máy thu đặt cách tâm 91 mm quét từ góc âm 60 độ đến dương 60 độ với bước dịch 1 độ, đảm bảo thu nhận trọn vẹn dữ liệu tán xạ 360 độ quanh vật thể.

Mức tiết kiệm thời gian 5,7% trên mô hình nhỏ có ý nghĩa như thế nào trong thực tế?

Mức 5,7% đạt được trên ma trận kích thước khiêm tốn 14 nhân 14 ô. Trong các hệ thống chẩn đoán y tế thực tế với độ phân giải hàng trăm nghìn pixel, khối lượng tính toán ma trận tăng theo cấp số nhân, do đó tỷ lệ cắt giảm phương trình đo thừa sẽ giúp tiết kiệm hàng chục phút xử lý cho mỗi ca chụp.

Kết luận

  • Luận văn đã giải quyết thành công bài toán tối ưu hóa tốc độ tái tạo ảnh siêu âm cắt lớp bằng cách tích hợp điều kiện quyết định ngưỡng áp suất tán xạ vào thuật toán lặp vi phân Born.
  • Hệ thống giảm thiểu hiệu quả các phép toán dư thừa, rút ngắn thời gian xử lý 5,7% từ 110,7 giây xuống 104,4 giây trên mô hình đối tượng 10 mm.
  • Chất lượng hình ảnh khôi phục được bảo toàn độ chính xác cao ngay cả trong điều kiện chịu tác động của 5% nhiễu Gauss ngẫu nhiên.
  • Toàn bộ thuật toán xử lý dữ liệu và cấu hình hệ thu phát 9 nguồn phát kết hợp 121 cảm biến thu đã được chuẩn hóa và kiểm chứng trên công cụ tính toán ma trận chuyên dụng.
  • Nghiên cứu đặt nền móng vững chắc cho lộ trình nâng cấp hệ thống lên dải tần 5 MHz và tích hợp vi xử lý tính toán song song trong giai đoạn tiếp theo.

Các viện nghiên cứu và đơn vị công nghệ y tế quan tâm có thể khai thác khung thuật toán này để tiếp tục phát triển các dòng máy chụp siêu âm cắt lớp thế hệ mới phục vụ cộng đồng.