Tổng quan luận án

Bệnh ung thư vú là một trong những nguyên nhân gây tử vong hàng đầu ở phụ nữ trên phạm vi toàn cầu. Theo số liệu thống kê được tác giả trích dẫn, mỗi năm thế giới ghi nhận khoảng 1,4 triệu ca mắc mới và 458.000 trường hợp tử vong [1]. Tại Việt Nam, tỷ lệ mắc mới ước tính khoảng 7.000 ca mỗi năm với tỷ lệ tử vong xấp xỉ 35%, đồng thời bệnh có xu hướng trẻ hóa qua các năm. Việc phát hiện sớm khối u khi kích thước còn nhỏ hơn 5 mm mang ý nghĩa quyết định, giúp tỷ lệ điều trị thành công đạt tới 90%. Hiện nay, kỹ thuật chụp X-quang tuyến vú (mammography) phổ biến đối với phụ nữ hậu mãn kinh nhưng bộc lộ hạn chế rõ rệt ở phụ nữ dưới 50 tuổi do mật độ mô vú dày đặc [67]. Trong khi đó, kỹ thuật siêu âm B-mode thông thường chỉ cung cấp thông tin định tính về hàm trở kháng âm, phụ thuộc lớn vào kinh nghiệm chủ quan của bác sĩ chẩn đoán và khó phát hiện tổn thương ở kích thước nhỏ.

Kỹ thuật siêu âm cắt lớp (Ultrasound Tomography - UST) dựa trên lý thuyết tán xạ ngược và phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM) có khả năng định lượng các đặc tính cơ học và tốc độ truyền âm của mô, cho phép phát hiện các khối u có kích thước nhỏ hơn bước sóng tới. Tuy nhiên, rào cản lớn nhất khiến các thiết bị siêu âm cắt lớp tán xạ chưa được thương mại hóa rộng rãi nằm ở độ phức tạp tính toán rất cao, thời gian tái tạo ảnh kéo dài, dung lượng bộ nhớ yêu cầu lớn và độ phân giải ảnh dễ bị ảnh hưởng bởi nhiễu. Bên cạnh đó, việc áp dụng kỹ thuật lấy mẫu nén ngẫu nhiên truyền thống gặp khó khăn khi hiện thực hóa trên phần cứng hệ thống do đòi hỏi các bộ tạo số ngẫu nhiên phức tạp, còn các giải thuật kết hợp đa tần số chưa xác định được thời điểm chuyển tiếp tối ưu giữa các dải tần số. Luận án được thực hiện nhằm giải quyết các khoảng trống kỹ thuật này.

Mục tiêu nghiên cứu của luận án được xác định gồm:

  1. Xây dựng và phát triển các phương pháp, mô hình, giải thuật và công cụ nhằm tăng tốc độ lấy mẫu và nâng cao chất lượng tạo ảnh siêu âm cắt lớp bằng kỹ thuật kết hợp hai tần số (DF-DBIM).
  2. Phát triển kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS) tích hợp vào phương pháp DBIM nhằm giảm thiểu số lượng đầu dò thu nhận, giảm chi phí tính toán và bảo đảm khả năng thực thi trên phần cứng của các hệ thống chụp cắt lớp siêu âm.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Sự thay đổi tốc độ truyền âm khi sóng siêu âm lan truyền qua môi trường không đồng nhất (mô bệnh lý, khối u); kỹ thuật tạo ảnh dựa trên trường tán xạ ngược; phương pháp lặp vi phân Born (DBIM); các thuật toán xử lý tín hiệu kết hợp hai tần số và lấy mẫu nén giả ngẫu nhiên.
  • Phạm vi nghiên cứu: Luận án giới hạn việc khảo sát trên mô hình tạo ảnh lát cắt hai chiều (2D) trong hệ tọa độ trụ; đối tượng tán xạ được mô hình hóa dưới dạng hình trụ tròn đặt tại tâm hệ đo trong môi trường nền đồng nhất (nước); chênh lệch tốc độ truyền sóng giữa khối u và môi trường nền ở mức nhỏ (khoảng 2%); tần số kích thích khảo sát trong dải từ 0,64 MHz đến 2 MHz; dữ liệu thực nghiệm kiểm chứng được thu thập từ nhóm nghiên cứu tại Đại học Illinois (UIUC, Hoa Kỳ).

Tổng quan tài liệu và vị trí của luận án

Tổng quan lịch sử phát triển cho thấy công nghệ tạo ảnh y sinh khởi nguồn từ phát hiện tia X của Wilhelm Roentgen năm 1885, sau đó mở rộng sang các kỹ thuật hiện đại như MRI, CT, PET, SPECT và siêu âm [66], [100]. Trong lĩnh vực âm học, nguyên lý sonar phát triển từ năm 1910 đã đặt nền móng cho kỹ thuật siêu âm B-mode [14]. Để trích xuất thông tin định lượng từ ảnh B-mode, nhiều nhóm nghiên cứu đã áp dụng các tham số kết cấu bậc một và bậc hai [28], [85], [99], mô hình thống kê đường bao [71], [98] hoặc phân tích biến đổi sóng con (wavelet) [11], [37], tuy nhiên các thông số này không độc lập với hệ thống đo.

Đối với kỹ thuật siêu âm định lượng, hai hướng tiếp cận chính được phát triển trong các thập kỷ qua:

  1. Siêu âm cắt lớp dựa trên lý thuyết chiếu tia: Xuất hiện từ cuối thập niên 1970 và đầu thập niên 1980 [30], [41], [42], [70]. Tiêu biểu là hệ thống CURE của Viện Ung thư Karmanos sử dụng phương pháp tia thẳng [60] kết hợp hiệu chỉnh khúc xạ [16]; hệ thống HUTT tại Đại học Nam California phân tích độ suy hao theo tần số [44], [45], [46]; các hệ thống 3D tại Viện Công nghệ Karlsruhe sử dụng mảng đầu dò hình trụ [35], [76]; và nghiên cứu của Zhang cùng cộng sự về thông số phi tuyến từ sóng hài bậc hai [21], [22]. Hạn chế cốt lõi của lý thuyết chiếu tia là chỉ tái tạo chính xác các cấu trúc có kích thước lớn hơn bước sóng từ 2 đến 5 lần [79], [102], không đáp ứng yêu cầu phát hiện tổn thương giai đoạn sớm.
  2. Siêu âm cắt lớp dựa trên lý thuyết tán xạ: Hướng tiếp cận này ước lượng các phân bố âm học như tốc độ truyền âm, suy hao và mật độ [19], [33], [36], [48], [51], [52], [54], [65]. Điển hình cho thiết bị thương mại là hệ thống chụp ảnh vú SoftVue [62] với 2048 biến tử, 512 kênh phát và 512 kênh thu, tần số 3 MHz. Về mặt giải thuật nghịch đảo sóng, hai phép xấp xỉ phổ biến gồm xấp xỉ Born bậc một [5] và xấp xỉ Rytov bậc một [4]. Nghiên cứu của Slaney và cộng sự [57] chỉ ra rằng xấp xỉ Rytov tuy không giới hạn kích thước vùng tán xạ nhưng gặp hiện tượng gập pha (phase wrapping) gây sai số lớn khi nghịch đảo [13], trong khi các giải thuật tháo gập pha [55], [89] chưa đạt hiệu quả cao. Ngược lại, xấp xỉ Born bậc một hoạt động ổn định khi điều kiện lệch pha cực đại thỏa mãn $|\Delta \phi| < \pi$.

Đối với các công trình cải tiến phương pháp lặp vi phân Born (DBIM):

  • Cải tiến thuật toán lặp: Nhóm tác giả tại [94] kết hợp phương pháp lặp Born (BIM) ở các vòng đầu để kháng nhiễu và chuyển sang DBIM ở các vòng sau để tăng tốc hội tụ. Một số nghiên cứu bổ sung việc khôi phục mật độ $\rho$ song song với tốc độ truyền âm [74], [75], [81], [83].
  • Cải tiến cấu hình đo và tối ưu hóa: Sử dụng chuẩn hóa Tikhonov truyền thống [31], thay thế bằng bài toán bình phương tối thiểu chuẩn hóa $l_1$ ($l_1$-LSP) [96], hoặc bố trí đầu dò ngẫu nhiên xung quanh đối tượng [95].
  • Tích hợp các kỹ thuật xử lý tín hiệu: Sử dụng phép nội suy láng giềng gần nhất (Nearest neighbor) từ lưới thô $6 \times 6$ lên lưới tinh $12 \times 12$ để giảm chi phí tính toán ban đầu [68], [93]. Đối với kỹ thuật kết hợp tần số [9], [58], [81], [84], [86], [87], [101], phương pháp đa tần số (MF-DBIM) [7] đòi hỏi lưu trữ lượng dữ liệu quá lớn, trong khi kỹ thuật hai tần số (DF-DBIM) chưa giải quyết được thời điểm chuyển tiếp tối ưu giữa tần số thấp và cao. Đối với kỹ thuật lấy mẫu nén (CS) theo lý thuyết của Candes, Tao [17] và Donoho [18], các công trình [6], [12], [34], [47], [69], [80], [90] và nhóm nghiên cứu tại Đại học Công nghệ, ĐHQGHN [20], [64], [91], [92], [97] đã bước đầu ứng dụng lấy mẫu ngẫu nhiên vào siêu âm cắt lớp nhưng vướng phải trở ngại khi thiết kế mạch phần cứng do tính chất ngẫu nhiên của vị trí đầu dò.

Khoảng trống luận án tập trung giải quyết: Xác định quy tắc và số vòng lặp chuyển giao tối ưu cho thuật toán hai tần số DF-DBIM; đề xuất phương pháp lấy mẫu nén giả ngẫu nhiên (DCS-DBIM) dựa trên chuỗi tất định nhằm giảm số kênh đo, giảm thời gian tính toán mà vẫn đảm bảo tính khả thi khi chế tạo phần cứng.

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

Cơ sở lý thuyết âm học và mô hình toán học

Sóng siêu âm là sóng dao động đàn hồi có tần số $f > 20\text{ kHz}$. Vận tốc truyền sóng biến đổi tùy theo đặc tính cơ học của từng mô sinh học. Sự biến thiên này được thể hiện trong dữ liệu đối sánh của luận án:

Cơ quan / Mô sinh học Tốc độ truyền sóng siêu âm (m/s)
Màng mỡ / Mô mềm 1540
Não bộ 1550
Gan 1565
Thận 1575
Máu 1570
Cơ bắp 1600
Thần kinh thị giác 1615
Thấu kính mắt 1630
Xương 4090
Phổi 630
Nước nền (tại nhiệt độ phòng) 1484

Trường áp suất âm tổng $p(\vec{r})$ khi truyền qua môi trường không đồng nhất thỏa mãn hệ thức: $$\vec{p} = \vec{p}{inc} + \vec{p}{sc}$$ trong đó $\vec{p}{inc}$ là áp suất sóng tới và $\vec{p}{sc}$ là áp suất sóng tán xạ. Phương trình vi phân Helmholtz mô tả sự lan truyền sóng có dạng: $$(\nabla^2 + k_0^2) p_{sc}(\vec{r}) = -k_0^2 O(\vec{r}) p(\vec{r})$$ với $k_0 = 2\pi f / c_0$ là số sóng trong môi trường nền và $k(\vec{r}) = 2\pi f / c(\vec{r})$ là số sóng của môi trường không đồng nhất. Hàm Green trong không gian hai chiều đồng nhất được biểu diễn bằng hàm Hankel loại 1 bậc 0: $$G_0(k_0 |\vec{r} - \vec{r}'|) = \frac{j}{4} H_0^{(1)}(k_0 |\vec{r} - \vec{r}'|)$$ Hàm mục tiêu $O(\vec{r})$ đặc trưng cho sự tương phản vận tốc âm được xác định: $$O(\vec{r}) = \frac{k^2(\vec{r}) - k_0^2}{k_0^2} \approx \frac{c_0^2}{c^2(\vec{r})} - 1$$

Phương trình tích phân liên tục được rời rạc hóa bằng phương pháp Moment (Method of Moments - MoM) sử dụng các hàm cơ bản sinc và hàm delta trên lưới khảo sát kích thước $N \times N$ điểm ảnh:

  • Áp suất tổng tại các điểm bên trong vùng chia lưới: $$\bar{p} = (\bar{I} - \bar{C} \operatorname{diag}(\bar{O}))^{-1} \bar{p}_{inc}$$
  • Áp suất tán xạ đo được tại $N_r$ đầu thu từ $N_t$ đầu phát (kích thước vector $N_t N_r \times 1$): $$\bar{p}_{sc} = \bar{B} \operatorname{diag}(\bar{p}) \bar{O}$$ trong đó $\bar{B}$ và $\bar{C}$ là các ma trận hệ thống chứa các tích phân thành phần của hàm Green.

Theo xấp xỉ Born bậc một, khi giả định $\vec{p}{sc} \ll \vec{p}{inc}$, phương trình tuyến tính hóa giữa độ sai khác áp suất tán xạ $\Delta \bar{p}{sc}$ và lượng cập nhật hàm mục tiêu $\Delta \bar{O}$ được thiết lập: $$\Delta \bar{p}{sc} = \bar{B} \operatorname{diag}(\bar{p}{inc}) \Delta \bar{O} = \bar{M} \Delta \bar{O}$$ Nghiệm của bài toán nghịch đảo được giải thông qua kỹ thuật chuẩn hóa Tikhonov: $$\Delta \bar{O} = \arg\min{\bar{O}} \left( |\Delta \bar{p}{sc} - \bar{M} \Delta \bar{O}|^2 + \gamma |\Delta \bar{O}|^2 \right)$$ Chi phí tính toán và yêu cầu bộ nhớ của thuật toán Tikhonov đạt mức $O(N_t N_r N^2)$. Sóng tới từ máy phát được mô hình hóa bằng chùm Bessel bậc 0: $$p{inc}(\vec{r}_k) = J_0(k_0 |\vec{r}_k - \vec{r}t|) = \sum{m=0}^{\infty} \frac{(-1)^m}{m! , \Gamma(m+1)} \left( \frac{k_0 |\vec{r}_k - \vec{r}_t|}{2} \right)^{2m}$$

Phương pháp nghiên cứu và dữ liệu thực nghiệm

Tác giả kết hợp phương pháp mô hình hóa toán học, mô phỏng số trên máy tính và kiểm chứng bằng dữ liệu đo thực nghiệm. Bộ dữ liệu thực nghiệm được cung cấp bởi nhóm nghiên cứu tại Đại học Illinois (UIUC, Hoa Kỳ) [82], thực hiện trên mô hình đối tượng thực nghiệm gồm một quả bóng cao su có đường kính $0,3\text{ mm} \pm 0,05\text{ mm}$ đặt trong bể nước ở nhiệt độ $18,1^\circ\text{C}$ ($c_0 = 1484\text{ m/s}$). Khoảng cách từ tâm hệ đo đến đầu phát là $104\text{ mm}$, khoảng cách đến đầu thu là $91\text{ mm}$, tần số hoạt động $0,64\text{ MHz}$. Đầu thu ghi nhận áp suất tán xạ khi quét từng bước $10^\circ$ trong phạm vi cung góc $120^\circ$.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương này hệ thống hóa toàn diện cơ sở của các kỹ thuật tạo ảnh y sinh và đi sâu phân tích kỹ thuật siêu âm cắt lớp. Tác giả phân định rõ hai hướng tiếp cận: tạo ảnh dựa trên lý thuyết chiếu tia và tạo ảnh dựa trên lý thuyết tán xạ. Phân tích chỉ ra rằng lý thuyết chiếu tia vấp phải giới hạn vật lý không thể khôi phục các chi tiết nhỏ hơn bước sóng, trong khi kỹ thuật tán xạ ngược giải quyết được bài toán này. Chương 1 tổng kết các nghiên cứu về xấp xỉ Born bậc một và xấp xỉ Rytov bậc một, luận giải nguyên nhân lựa chọn xấp xỉ Born dựa trên điều kiện pha $|\Delta \phi| < \pi$. Đồng thời, tác giả khảo sát các công trình liên quan đến phương pháp DBIM: giải pháp kết hợp BIM-DBIM [94], kỹ thuật nội suy lưới [93], kết hợp đa tần số [7], [81] và lấy mẫu nén ngẫu nhiên [80]. Từ việc chỉ ra độ phức tạp tính toán lớn của DBIM và tính bất khả thi trong chế tạo phần cứng của CS ngẫu nhiên, chương này thiết lập định hướng nghiên cứu tập trung vào hai giải thuật: DF-DBIM (hai tần số) và DCS-DBIM (lấy mẫu nén giả ngẫu nhiên).

Chương 2: Kỹ thuật tạo ảnh siêu âm cắt lớp

Nội dung chương 2 tập trung xây dựng mô hình giải tích và quy trình thuật toán cho phương pháp lặp vi phân Born (DBIM). Quy trình phát - thu được mô tả theo cấu hình vòng tròn khép kín $360^\circ$ quanh đối tượng: lần lượt từng máy phát trong số $N_t$ đầu phát hoạt động, toàn bộ $N_r$ đầu thu ghi nhận tín hiệu tán xạ đồng thời, tạo ra $N_t \times N_r$ phép đo độc lập. Tác giả xây dựng chi tiết thuật toán DBIM (Giải thuật 2.1), trong đó ma trận hệ thống và trường Green được tính toán lại sau mỗi vòng lặp dựa trên hàm mục tiêu ước lượng ở bước trước.

Phần thực nghiệm kiểm chứng tiến hành so sánh dữ liệu tán xạ mô phỏng ở các mức phân giải lưới khác nhau ($N = 16, 32, 64$) với dữ liệu đo thực nghiệm từ UIUC trên đối tượng quả bóng đường kính $0,3\text{ mm}$:

  • Ở độ phân giải $N = 16$, sai lệch giữa mô phỏng và thực nghiệm xuất hiện rõ rệt do lưới quá thô.
  • Ở độ phân giải $N = 32$ và $N = 64$, đường cong áp suất tán xạ mô phỏng bám rất sát dữ liệu thực nghiệm, đặc biệt ở các vùng có biên độ tán xạ lớn ($\bar{p}_{sc} > 0,4$). Khi áp dụng thêm bộ lọc trung bình để làm trơn dữ liệu mô phỏng, độ tương đồng đạt mức cao.
  • Trên cơ sở cân đối giữa độ chính xác và chi phí tính toán ($O(N_t N_r N^2)$), tác giả lựa chọn độ phân giải $N = 32$ làm chuẩn tham số cho các mô phỏng tiếp theo trong luận án.

Chương 3: Thuật toán kết hợp tần số dùng trong tạo ảnh siêu âm cắt lớp

Chương 3 đề xuất giải thuật hai tần số DF-DBIM (Giải thuật 3.1) nhằm giải quyết hiện tượng mất hội tụ của xấp xỉ Born trong môi trường tán xạ mạnh mà không làm tăng quá mức chi phí bộ nhớ như phương pháp đa tần số MF-DBIM. Nguyên lý của giải thuật là sử dụng tần số thấp $f_1$ trong $N_{f1}$ vòng lặp đầu để mở rộng miền hội tụ và thu được cấu trúc bao quát của đối tượng, sau đó chuyển sang tần số cao $f_2$ trong $N_{f2}$ vòng lặp còn lại ($N_{total} = N_{f1} + N_{f2} = 8$ vòng lặp) để tái tạo chi tiết và nâng cao độ phân giải không gian.

Vấn đề cốt lõi được tác giả giải quyết là khảo sát tìm giá trị chuyển giao $N_{f1}$ tối ưu thông qua 4 kịch bản mô phỏng số với số lượng máy phát $N_t$, máy thu $N_r$ và số biến tính toán khác nhau:

Kịch bản Tần số $f_1$ / $f_2$ (MHz) Cấu hình đo ($N_t \times N_r$) Số điểm lưới ($N$) Số phép đo / Số biến Giá trị $N_{f1}$ tối ưu
Kịch bản 1 1,0 / 2,0 Khảo sát theo $N_t, N_r$ 22 Phép đo < Số biến $N_{f1} = 3$
Kịch bản 2 1,0 / 2,0 Khảo sát theo $N_t, N_r$ 22 Phép đo $\approx$ Số biến $N_{f1} = 3$
Kịch bản 3 1,0 / 2,0 Cố định cấu hình 22 Phép đo > Số biến $N_{f1} = 4$
Kịch bản 4 1,0 / 2,0 Cấu hình mở rộng 22 Khảo sát độc lập $N_{f1} = 4$

Kết quả so sánh sai số thặng dư tương đối (RRE) sau 8 vòng lặp khẳng định:

  • Nếu chuyển giao quá sớm ($N_{f1} = 1$ hoặc $2$), thuật toán dễ rơi vào cực tiểu cục bộ do tần số thấp chưa kịp hội tụ ổn định.
  • Nếu chuyển giao quá muộn ($N_{f1} \ge 6$), số vòng lặp ở tần số cao không đủ để làm sắc nét ảnh khôi phục.
  • Tại điểm chuyển tiếp tối ưu ($N_{f1} = 3$ hoặc $4$ tùy theo tỷ lệ số phép đo trên số biến), sai số chuẩn hóa của giải thuật DF-DBIM đạt mức thấp hơn rõ rệt so với phương pháp DBIM đơn tần số truyền thống, đồng thời biên dạng mặt cắt dọc hàm mục tiêu bám sát giá trị lý tưởng.

Chương 4: Thuật toán lấy mẫu nén dùng trong tạo ảnh siêu âm cắt lớp

Chương 4 tập trung xử lý trở ngại về số lượng kênh đo và thời gian tính toán của DBIM bằng kỹ thuật lấy mẫu nén giả ngẫu nhiên (Deterministic Compressed Sensing - DCS). Thay vì bố trí đầu thu ngẫu nhiên dựa trên các bộ phát số ngẫu nhiên phần cứng (HRNG) khó chế tạo, tác giả sử dụng chuỗi ánh xạ logistic tất định để xác định vị trí phân bố của các máy thu xung quanh vòng tròn đo.

Mô hình DCS-DBIM (Giải thuật 4.1) kết hợp cấu hình thu giả ngẫu nhiên với thuật toán giải bài toán tối ưu hóa chuẩn hóa $l_1$ ($l_1$-regularized Least Squares Problem - $l_1$-LSP) nhằm khôi phục tín hiệu mục tiêu thưa. Hiệu năng của giải pháp được kiểm chứng thông qua so sánh giữa cấu hình đo đều truyền thống và cấu hình DCS đề xuất:

Phương pháp Cấu hình đầu dò ($N_t \times N_r$) Tổng số phép đo Tỷ số nén ($r$) Thời gian tính toán sau 8 vòng lặp Mức độ sai số chuẩn hóa (RRE)
DBIM truyền thống $22 \times 22$ 484 1,0 (đo đầy đủ) Chuẩn quy chiếu cao Mức sai số cơ sở
DCS-DBIM đề xuất $16 \times 16$ 256 0,53 Giảm đáng kể Tương đương hoặc thấp hơn DBIM
DBIM truyền thống $30 \times 30$ 900 1,0 (đo đầy đủ) Rất cao Mức sai số cơ sở
DCS-DBIM đề xuất $20 \times 20$ 400 0,44 Giảm hơn 50% Thấp hơn DBIM truyền thống

Kết quả mô phỏng sau 8 vòng lặp cho thấy DCS-DBIM với 256 phép đo đạt chất lượng tái tạo ảnh và sai số tương đương phương pháp DBIM truyền thống sử dụng 484 phép đo. Đặc biệt, khi so sánh trực tiếp giữa DCS-DBIM và CS-DBIM (lấy mẫu ngẫu nhiên thuần túy) ở các tỷ số nén khác nhau, DCS-DBIM cho thấy độ ổn định sai số tương đương nhưng vượt trội về tính khả thi thực tế khi triển khai phần cứng.

Kết quả và những đóng góp mới

Luận án đã đạt được các đóng góp khoa học và kỹ thuật cụ thể sau:

  1. Đóng góp về lý thuyết và thuật toán:

    • Đề xuất hoàn chỉnh quy trình giải thuật DF-DBIM kết hợp hai tần số. Luận án đã giải quyết bài toán xác định thời điểm chuyển giao tần số ($N_{f1}$), chứng minh bằng định lượng rằng việc chuyển đổi từ tần số thấp sang tần số cao ở vòng lặp thứ 3 hoặc thứ 4 (trong chuỗi 8 vòng lặp) giúp tối ưu hóa sự cân bằng giữa tốc độ hội tụ và độ phân giải không gian của ảnh siêu âm cắt lớp.
    • Xây dựng thành công thuật toán DCS-DBIM kết hợp phương pháp lặp vi phân Born với kỹ thuật lấy mẫu nén giả ngẫu nhiên sử dụng chuỗi ánh xạ logistic và chuẩn hóa $l_1$-LSP. Đây là giải pháp thay thế cho lấy mẫu nén ngẫu nhiên, loại bỏ yêu cầu sử dụng mạch tạo số ngẫu nhiên phức tạp trên phần cứng.

  2. Đóng góp về mặt kỹ thuật và thực tiễn:

    • Giảm đáng kể số lượng phần tử đầu dò và số phép đo thu thập dữ liệu (từ 484 phép đo xuống 256 phép đo; từ 900 phép đo xuống 400 phép đo) mà vẫn duy trì hoặc cải thiện chất lượng ảnh khôi phục, giúp tiết kiệm bộ nhớ tính toán và rút ngắn thời gian tạo ảnh.
    • Kiểm chứng thành công mô hình mô phỏng DBIM với dữ liệu thực nghiệm áp suất tán xạ từ quả bóng đường kính $0,3\text{ mm}$ do Đại học Illinois cung cấp, xác lập tham số chia lưới tối ưu ($N = 32$) bảo đảm độ tin cậy khoa học khi áp dụng thuật toán vào các bài toán tái tạo ảnh tiếp theo.

  3. Giải pháp và kiến nghị đề xuất:

    • Luận án kiến nghị áp dụng cấu hình mảng đầu dò giả ngẫu nhiên kết hợp thuật toán DCS-DBIM cho các thiết bị chụp cắt lớp siêu âm thế hệ mới nhằm hạ thấp giá thành phần cứng và hướng tới khả năng xử lý thời gian thực trong sàng lọc lâm sàng ung thư vú.

Hạn chế và hướng nghiên cứu tiếp

Căn cứ trên các nội dung văn bản trình bày, những hạn chế kỹ thuật của luận án bao gồm:

  • Giới hạn chiều không gian: Các thuật toán mới chỉ được thử nghiệm và đánh giá trên các lát cắt hai chiều (2D), chưa được mở rộng thực thi trực tiếp trên mô hình không gian ba chiều (3D) hoàn chỉnh.
  • Môi trường mô phỏng đơn giản hóa: Đối tượng nghiên cứu trong các kịch bản mô phỏng được giả định là hình trụ tròn lý tưởng đặt trong môi trường nền đồng nhất (nước), chưa bao hàm đầy đủ cấu trúc giải phẫu không đồng nhất phức tạp của vú thật (như hệ thống ống dẫn sữa, mô mỡ, mô liên kết).
  • Thuộc tính âm học khảo sát: Luận án mới chỉ tập trung xử lý tham số tốc độ truyền âm ($c$), chưa tích hợp đồng thời việc tái tạo hệ số suy hao âm ($\alpha$) và phân bố mật độ mô ($\rho$) vào cùng một khung thuật toán.
  • Mức độ kiểm chứng thực nghiệm: Dữ liệu thực nghiệm mới dừng lại ở mô hình phantom vật lý đơn giản (quả bóng cao su đường kính $0,3\text{ mm}$), chưa có điều kiện thử nghiệm trên mẫu mô sinh học thật hoặc thử nghiệm tiền lâm sàng trên bệnh nhân.

Hướng nghiên cứu tiếp theo được xác định là mở rộng thuật toán sang mô hình 3D, tích hợp khôi phục đồng thời đa tham số âm học và từng bước chế tạo hệ thống phần cứng thực nghiệm để thử nghiệm trên các phantom sinh học phức tạp.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị chuyên môn đối với:

  • Nghiên cứu sinh, học viên cao học và giảng viên thuộc các chuyên ngành Kỹ thuật Điện tử, Xử lý tín hiệu, Kỹ thuật Viễn thông và Kỹ thuật Y sinh: Cung cấp tài liệu tham khảo chi tiết về cơ sở toán học của bài toán tán xạ ngược, kỹ thuật rời rạc hóa phương trình tích phân bằng phương pháp Moment (MoM) và kỹ thuật lập trình giải thuật DBIM tại Chương 2.
  • Các kỹ sư và nhóm nghiên cứu phát triển thiết bị y tế: Có thể tham khảo phương pháp xác định điểm chuyển tiếp tần số tại Chương 3 để tối ưu hóa quy trình thu nhận tín hiệu đa tần, cũng như ứng dụng sơ đồ bố trí đầu dò giả ngẫu nhiên theo chuỗi logistic tại Chương 4 nhằm giảm bớt số kênh thu phát khi thiết kế hệ thống phần cứng siêu âm cắt lớp.

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

1. Tại sao tác giả lựa chọn kỹ thuật xấp xỉ Born bậc một thay vì xấp xỉ Rytov bậc một?
Trả lời: Mặc dù xấp xỉ Rytov không bị giới hạn bởi kích thước đối tượng tán xạ, nhưng phép nghịch đảo dựa trên xấp xỉ Rytov thường gặp lỗi gập pha (phase wrapping) nghiêm trọng khi ước tính pha tổng từ dữ liệu đo, dẫn đến sai số đột ngột. Trong khi đó, xấp xỉ Born bậc một hoạt động ổn định và tin cậy khi điều kiện độ lệch pha cực đại thỏa mãn $|\Delta \phi| < \pi$, hoàn toàn phù hợp với bài toán phát hiện tổn thương giai đoạn sớm có kích thước nhỏ và độ tương phản vận tốc âm thấp so với môi trường nền.

2. Dữ liệu thực nghiệm kiểm chứng trong Chương 2 được thu thập với cấu hình như thế nào?
Trả lời: Dữ liệu thực nghiệm do nhóm nghiên cứu tại Đại học Illinois (UIUC, Hoa Kỳ) cung cấp, đo tín hiệu áp suất tán xạ từ một quả bóng cao su đường kính $0,3\text{ mm} \pm 0,05\text{ mm}$ đặt trong môi trường nước ở $18,1^\circ\text{C}$ ($c_0 = 1484\text{ m/s}$). Hệ đo có khoảng cách từ tâm tới máy phát là $104\text{ mm}$, khoảng cách tới máy thu là $91\text{ mm}$, hoạt động ở tần số $0,64\text{ MHz}$, đầu thu thực hiện quét với bước dịch $10^\circ$ trong phạm vi góc $120^\circ$.

3. Ý nghĩa của việc xác định số vòng lặp chuyển giao $N_{f1}$ trong thuật toán DF-DBIM là gì?
Trả lời: Việc xác định chính xác số vòng lặp $N_{f1}$ là điểm mấu chốt để thuật toán đạt hiệu quả cao nhất. Tần số thấp $f_1$ giúp giải thuật hội tụ nhanh về vùng nghiệm thực tế nhưng cho độ phân giải thấp; tần số cao $f_2$ giúp tăng độ phân giải không gian. Nếu chuyển giao quá sớm khi tần số thấp chưa ổn định, giải thuật sẽ mất hội tụ; nếu chuyển giao quá muộn, số vòng lặp còn lại không đủ để tần số cao làm sắc nét ảnh. Thực nghiệm mô phỏng chỉ ra giá trị tối ưu là $N_{f1} = 3$ hoặc $4$ vòng lặp trên tổng số 8 vòng lặp.

4. Kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS) giải quyết vấn đề gì của lấy mẫu nén ngẫu nhiên (CS) truyền thống?
Trả lời: Kỹ thuật CS ngẫu nhiên truyền thống đòi hỏi các máy thu phải bố trí ngẫu nhiên trên hệ đo, điều này yêu cầu hệ thống phải tích hợp các bộ phát số ngẫu nhiên phần cứng (HRNG) hoặc giả ngẫu nhiên (PRNG) rất phức tạp khi chế tạo. Kỹ thuật DCS sử dụng chuỗi ánh xạ logistic có tính chất tất định, cho phép xác định cố định trước vị trí hình học của các đầu dò trên mảng đo, giúp việc thiết kế mạch phần cứng trở nên khả thi và đơn giản hơn mà vẫn giữ được hiệu quả giảm số lượng phép đo.

5. Cơ sở nào để luận án lựa chọn kích thước lưới $N = 32$ trong các mô phỏng?
Trả lời: Trong Chương 2, tác giả đã tiến hành đối sánh dữ liệu tán xạ mô phỏng ở ba mức lưới $N = 16, 32, 64$ với dữ liệu đo thực nghiệm từ UIUC. Kết quả cho thấy mức lưới $N = 16$ có sai số lớn, trong khi mức lưới $N = 32$ và $N = 64$ đều cho kết quả tương đồng cao với thực nghiệm sau khi lọc trung bình. Do chi phí tính toán của phương pháp nghịch đảo tỷ lệ với $N^2$ ($O(N_t N_r N^2)$), tác giả lựa chọn $N = 32$ nhằm đảm bảo độ chính xác của ảnh tái tạo đồng thời giảm thiểu tối đa khối lượng tính toán.

Kết luận

Bản luận án tiến sĩ của nghiên cứu sinh Trần Quang Huy đã xây dựng và kiểm chứng thành công các giải pháp xử lý tín hiệu nâng cao chất lượng tạo ảnh siêu âm cắt lớp dựa trên phương pháp lặp vi phân Born. Bằng việc đề xuất thuật toán kết hợp hai tần số DF-DBIM với thời điểm chuyển giao tối ưu và thuật toán lấy mẫu nén giả ngẫu nhiên DCS-DBIM sử dụng chuỗi logistic kết hợp chuẩn hóa $l_1$-LSP, công trình đã giải quyết đồng thời hai nút thắt lớn về tốc độ tính toán và độ phức tạp phần cứng. Các kết quả mô phỏng số và đối sánh dữ liệu thực nghiệm đã xác nhận tính đúng đắn của mô hình, mở ra khả năng ứng dụng thực tế cho các hệ thống chẩn đoán hình ảnh y sinh thế hệ mới.