Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Bệnh ung thư vú là một trong những nguyên nhân gây tử vong hàng đầu ở phụ nữ. Thống kê trong văn bản cho thấy trên thế giới mỗi năm có khoảng 1,4 triệu ca mắc mới và 458.000 ca tử vong. Tại Việt Nam, mỗi năm ghi nhận khoảng 7.000 ca mắc mới với tỷ lệ tử vong khoảng 35%, đồng thời bệnh có xu hướng trẻ hóa. Việc phát hiện sớm các khối u có kích thước nhỏ dưới 5 mm đóng vai trò quyết định, giúp tỷ lệ điều trị thành công đạt mức 90%.

Phương pháp chụp X-quang tuyến vú (mammography) hiện tồn tại hạn chế lớn khi áp dụng cho phụ nữ dưới 50 tuổi do mật độ mô tuyến vú dày đặc, gây suy giảm độ tương phản hình ảnh. Phương pháp siêu âm B-mode truyền thống chỉ cung cấp thông tin định tính về trở kháng âm và khó nhận diện tổn thương kích thước nhỏ. Ngược lại, kỹ thuật siêu âm cắt lớp (Ultrasound Tomography - UST) dựa trên lý thuyết tán xạ sóng âm cho phép tái tạo định lượng các thông số âm học như tốc độ truyền âm và độ suy hao, hỗ trợ nhận diện các u nhỏ ở giai đoạn sớm.

Mặc dù có nhiều triển vọng, kỹ thuật siêu âm cắt lớp sử dụng phương pháp lặp vi phân Born (Distorted Born Iterative Method - DBIM) gặp phải các rào cản kỹ thuật lớn:

  • Độ phức tạp tính toán cao $O(N_t N_r N^2)$ và yêu cầu dung lượng bộ nhớ lớn, khiến thời gian tái tạo ảnh kéo dài.
  • Tốc độ hội tụ và độ phân giải không gian phụ thuộc chặt chẽ vào tần số sóng siêu âm; việc sử dụng đơn tần số hoặc đa tần số chưa tối ưu về điểm chuyển giao tần số dẫn đến chất lượng ảnh kém hoặc tăng gánh nặng tính toán.
  • Ứng dụng kỹ thuật lấy mẫu nén (Compressed Sensing - CS) ngẫu nhiên gặp khó khăn khi hiện thực hóa trên phần cứng do đòi hỏi các bộ phát số ngẫu nhiên vật lý (HRNG/PRNG) phức tạp để điều khiển vị trí đầu dò.

Mục tiêu nghiên cứu

Luận án xác định các mục tiêu cụ thể 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 quá trình lấy mẫu và nâng cao chất lượng tạo ảnh siêu âm cắt lớp.
  2. Đề xuất và khảo sát kỹ thuật kết hợp hai tần số (DF-DBIM) nhằm cải thiện độ phân giải không gian và xác định thời điểm chuyển giao tần số tối ưu giữa tần số thấp và tần số cao.
  3. Đề xuất và phát triển giải thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM) dựa trên chuỗi ánh xạ logistic nhằm giảm số lượng phép đo, giảm độ phức tạp tính toán và tạo khả năng hiện thực hóa trên phần cứng hệ thống đo.

Đố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 ($\Delta c$) trong môi trường không đồng nhất khi xuất hiện khối u; kỹ thuật tạo ảnh dựa trên bài toán tán xạ ngược sóng siêu âm; phương pháp lặp vi phân Born (DBIM); các kỹ thuật xử lý tín hiệu kết hợp (kỹ thuật hai tần số DF-DBIM và kỹ thuật lấy mẫu nén giả ngẫu nhiên DCS-DBIM).
  • Phạm vi không gian: Tập trung vào bài toán tạo ảnh lát cắt 2D trong hệ tọa độ trụ cho cấu hình đo tròn bao quanh đối tượng (mô phỏng mô tuyến vú), làm cơ sở để phát triển tạo ảnh 3D.
  • Phạm vi tham số: Đối tượng khảo sát là hình trụ tròn đặt trong môi trường nền đồng nhất (nước); kích thước đối tượng nhỏ hơn hoặc tương đương bước sóng ($\le \lambda$); độ chênh lệch tốc độ truyền âm nhỏ ($\Delta c$ khoảng vài phần trăm, bỏ qua hiện tượng phản xạ và khúc xạ mạnh); tần số siêu âm khảo sát trong dải chẩn đoán y sinh (từ 0,64 MHz đến các tần số khảo sát trong mô phỏng).

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

Các hướng nghiên cứu trước đây

Luận án hệ thống hóa các hướng tiếp cận chính trong lịch sử phát triển kỹ thuật siêu âm cắt lớp:

  1. Siêu âm cắt lớp dựa trên lý thuyết chiếu tia (Ray-based tomography):

    • Các hệ thống thực nghiệm từ cuối thập niên 1970 và đầu thập niên 1980 của Greenleaf và cộng sự (1981), Carson và cộng sự (1981), Glover (1977).
    • Hệ thống CURE do Viện Ung thư Karmanos phát triển (Nebeker và cộng sự, 2012; Duric và cộng sự, 2005, 2007) sử dụng phép chiếu tia thẳng và hiệu chỉnh khúc xạ để tái tạo tốc độ âm, độ suy hao và độ phản xạ.
    • Hệ thống HUTT tại Đại học Nam California (Huang và cộng sự, 2005, 2008) sử dụng kỹ thuật chiếu tia đa tần số để tái tạo hệ số suy hao.
    • Hệ thống siêu âm 3D tại Viện Công nghệ Karlsruhe (Ruiter và cộng sự, 2006; Gemmeke và cộng sự, 2007) sắp xếp đầu dò trong khẩu độ hình trụ 2D.
    • Các nghiên cứu của Zhang và cộng sự (2009, 2010) khôi phục tham số phi tuyến từ sóng hài bậc 2.
    • Hạn chế: Lý thuyết chiếu tia chỉ áp dụng chính xác cho cấu trúc có kích thước lớn hơn bước sóng (từ 2 đến 5 bước sóng), mất độ chính xác đối với cấu trúc tán xạ nhỏ (Mast, 1999; Simonetti và cộng sự, 2007).

  2. Siêu âm cắt lớp dựa trên lý thuyết tán xạ (Wave-based / Inverse scattering):

    • Phương pháp xấp xỉ Rytov bậc 1 (Devaney, 1981; Kaveh và cộng sự, 1982) không hạn chế kích thước đối tượng nhưng gặp hiện tượng gập pha (phase wrapping), gây sai số đột ngột khi giải bài toán nghịch đảo (Kak và Slaney, 1988).
    • Phương pháp xấp xỉ Born bậc 1 và phương pháp lặp Born (BIM - Wang và Chew, 1989), phương pháp lặp vi phân Born (DBIM - Chew và Wang, 1990). Điều kiện áp dụng yêu cầu độ dịch pha $|\Delta \phi| < \pi$.
    • Thiết bị thương mại SoftVue (Glide và cộng sự, 2012) áp dụng siêu âm cắt lớp với 2048 biến tử, 512 kênh phát, 512 kênh thu, tần số 3 MHz.

  3. Các nghiên cứu cải tiến thuật toán DBIM:

    • Kết hợp BIM và DBIM để cân bằng giữa khả năng kháng nhiễu và tốc độ hội tụ (Tran và cộng sự, 2013).
    • Tái tạo đồng thời thông tin mật độ và tốc độ truyền âm (Lavarello và cộng sự, 2006; Hesford và cộng sự, 2010).
    • Cấu hình đo phân bố ngẫu nhiên đầu dò và giải bài toán tối ưu $l_1$ (LSP) thay thế chuẩn tắc Tikhonov (Nguyen và cộng sự, 2013, 2014).
    • Kỹ thuật nội suy Nearest neighbor trong DBIM để giảm độ phức tạp (Tran và cộng sự, 2012, 2013).
    • Kỹ thuật đa tần số (Haddadin và Ebbini, 1998; Pratt và cộng sự, 2007; Lin và cộng sự, 2011; Roy và cộng sự, 2013) và kỹ thuật lấy mẫu nén CS (Candes và Tao, 2006; Donoho, 2006; Provost và cộng sự, 2009; Zhao và cộng sự, 2012).

Khoảng trống nghiên cứu được lựa chọn giải quyết

Luận án xác định hai khoảng trống kỹ thuật cốt lõi:

  • Trong kỹ thuật kết hợp tần số cho DBIM, các nghiên cứu trước chưa đưa ra tiêu chí xác định thời điểm chuyển giao ($N_{f1}$) từ tần số thấp sang tần số cao, dẫn đến nguy cơ thuật toán không hội tụ hoặc chất lượng ảnh kém hơn khi dùng đơn tần số.
  • Trong kỹ thuật lấy mẫu nén áp dụng cho DBIM, các nghiên cứu trước đây dùng phân bố ngẫu nhiên thuần túy (CS), gây cản trở lớn khi chế tạo phần cứng đầu dò. Luận án giải quyết bằng cách đề xuất cấu hình lấy mẫu nén giả ngẫu nhiên (DCS) có tính tiền định toán học.

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

Khung lý thuyết và mô hình toán học

Luận án xây dựng trên nền tảng lý thuyết truyền sóng âm trong môi trường không đồng nhất 2D:

  1. Phương trình sóng và biểu diễn trường áp suất âm: Tổng áp suất âm $p(\vec{r})$ tại vị trí $\vec{r}$ là tổng của áp suất sóng tới $p^{(i)}(\vec{r})$ và áp suất sóng tán xạ $p^{(sc)}(\vec{r})$: $$p(\vec{r}) = p^{(i)}(\vec{r}) + p^{(sc)}(\vec{r})$$ Phương trình vi phân sóng: $$\nabla^2 p(\vec{r}) + k_0^2 p(\vec{r}) = -O(\vec{r}) p(\vec{r})$$ Trong đó, $k_0 = \frac{2\pi f}{c_0}$ là số sóng trong môi trường nền đồng nhất, và hàm mục tiêu lý tưởng $O(\vec{r})$ biểu diễn sự thay đổi tốc độ truyền âm: $$O(\vec{r}) = k_0^2 \left( \frac{c_0^2}{c^2(\vec{r})} - 1 \right)$$ Với $c_0$ là tốc độ âm môi trường nền và $c(\vec{r})$ là tốc độ âm tại vị trí $\vec{r}$.

  2. Nghiệm tích phân và hàm Green: Biểu diễn trường áp suất tán xạ qua tích phân diện tích với hàm Green $G_0$: $$p^{(sc)}(\vec{r}) = \iint O(\vec{r}') p(\vec{r}') G_0(k_0 |\vec{r} - \vec{r}'|) d\vec{r}'$$ Trong môi trường nền 2D đồng nhất, hàm Green là hàm Hankel loại 1 bậc 0: $$G_0(k_0 |\vec{r} - \vec{r}'|) = \frac{i}{4} H_0^{(1)}(k_0 |\vec{r} - \vec{r}'|)$$

  3. Rời rạc hóa bằng phương pháp Moment (MoM): Vùng khảo sát được chia thành lưới gồm $N \times N$ điểm ảnh ($N^2$ ẩn số). Hệ phương trình ma trận được thiết lập:

    • Áp suất tổng tại các điểm bên trong vùng chia lưới: $$\bar{E} = (\bar{I} - \bar{C} \text{diag}(\bar{O}))^{-1} \bar{E}^{(i)}$$ Trong đó, $\bar{C}$ là ma trận hàm Green biểu thị tương tác giữa các điểm ảnh, $\bar{I}$ là ma trận đơn vị.
    • Áp suất tán xạ đo được tại các đầu thu ($N_t N_r \times 1$ phần tử): $$\bar{p}^{(sc)} = \bar{B} \text{diag}(\bar{E}) \bar{O}$$ Trong đó, $\bar{B}$ là ma trận hàm Green biểu thị tương tác từ các điểm ảnh tới các máy thu.

  4. Giải bài toán nghịch đảo trong DBIM: Sử dụng xấp xỉ Born để tuyến tính hóa mối liên hệ giữa sai lệch áp suất tán xạ $\Delta \bar{p}^{(sc)}$ và biến thiên hàm mục tiêu $\Delta \bar{O}$: $$\Delta \bar{p}^{(sc)} = \bar{M} \Delta \bar{O}$$ Hàm mục tiêu được cập nhật sau mỗi vòng lặp: $\bar{O}^{(n+1)} = \bar{O}^{(n)} + \Delta \bar{O}^{(n)}$. Nghiệm $\Delta \bar{O}$ được tìm qua phương pháp chuẩn tắc Tikhonov: $$\Delta \bar{O} = \arg\min_{\Delta \bar{O}} \left( |\Delta \bar{p}^{(sc)} - \bar{M} \Delta \bar{O}|_2^2 + \gamma |\Delta \bar{O}|_2^2 \right)$$ Với $\gamma$ là tham số điều chuẩn Tikhonov.

Bảng tham chiếu tốc độ sóng âm trong cơ thể

Văn bản cung cấp bảng đối sánh tốc độ truyền sóng âm trong các môi trường mô sinh học:

Cơ quan / Mô Tốc độ siêu âm ($m/s$)
Phổi 630
Mô mềm 1540
Não bộ 1550
Gan 1565
Máu 1570
Thận 1575
Cơ bắp 1600
Thần kinh thị giác 1615
Thấu kính mắt 1630
Xương 4090

Phương pháp nghiên cứu thực tế

  • Mô hình hóa và mô phỏng số: Thiết lập cấu hình hệ đo tròn bao quanh đối tượng 2D. Tín hiệu phát được mô phỏng bằng chùm Bessel bậc 0 trong hệ tọa độ trụ: $$p^{(i)}(r, \theta) = \sum_{m=-\infty}^{\infty} J_m(k_0 r) e^{im\theta}$$
  • Kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS): Sử dụng chuỗi ánh xạ phi tuyến (chuỗi logistic) $x_{k+1} = \mu x_k (1 - x_k)$ với tham số điều khiển $\mu$ để xác định tiền định vị trí các máy thu trên vòng tròn đo, sau đó giải bài toán tối ưu chuẩn $l_1$ ($l_1$-LSP): $$\min |\bar{x}|_1 \quad \text{thỏa mãn} \quad |\bar{y} - \bar{\Phi} \bar{\Psi} \bar{x}|_2 \le \epsilon$$
  • Nguồn dữ liệu thực nghiệm kiểm chứng: Dữ liệu đo áp suất tán xạ thực tế thu từ mô hình thực nghiệm:
    • Đối tượng mẫu: Quả bóng đặt trong môi trường nước, đường kính 0,3 mm (sai số 0,05 mm).
    • Tần số âm: 0,64 MHz.
    • Khoảng cách từ tâm đến máy phát: 104 mm; khoảng cách từ tâm đến máy thu: 91 mm.
    • Nhiệt độ môi trường nước: 18,1 °C.
    • Đầu thu quét từng bước dịch $1^\circ$ trên cung quét $120^\circ$.

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

Chương 1: Tổng quan

Chương 1 phân tích bức tranh tổng thể về dịch tễ ung thư vú và vai trò của các công nghệ chẩn đoán hình ảnh y sinh. Nội dung chương đối sánh hai nhóm kỹ thuật siêu âm cắt lớp chính: nhóm dựa trên lý thuyết chiếu tia và nhóm dựa trên lý thuyết tán xạ ngược.

Tác giả chỉ ra rằng phương pháp chiếu tia tuy đơn giản nhưng chỉ nhận diện được đối tượng lớn hơn 2 đến 5 lần bước sóng. Nhóm phương pháp tán xạ ngược (đặc biệt là DBIM với xấp xỉ Born bậc 1) giải quyết được bài toán nhận diện u kích thước nhỏ hơn bước sóng, song gặp hạn chế về chi phí tính toán lớn ($O(N_t N_r N^2)$). Chương 1 khảo sát các hướng cải tiến DBIM đã có như tích hợp thuật toán nội suy, kết hợp đa tần số và ứng dụng lấy mẫu nén, qua đó xác lập định hướng nghiên cứu chính của luận án là phát triển thuật toán hai tần số DF-DBIM và lấy mẫu nén giả ngẫu nhiên DCS-DBIM.

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

Chương 2 trình bày cơ sở toán học chi tiết của bài toán thuận và bài toán nghịch trong siêu âm cắt lớp sử dụng phương pháp lặp vi phân Born (DBIM).

Nội dung trọng tâm bao gồm:

  • Thiết lập hệ phương trình sóng tán xạ 2D sử dụng hàm Green trong môi trường nền đồng nhất dạng hàm Hankel loại 1 bậc 0.
  • Rời rạc hóa miền tính toán bằng phương pháp Moment (MoM) với hệ hàm cơ sở sinc và hàm delta để xây dựng các ma trận tương tác $\bar{C}$ và $\bar{B}$.
  • Áp dụng kỹ thuật điều chuẩn Tikhonov để giải hệ phương trình tuyến tính hóa tìm gia số hàm mục tiêu $\Delta \bar{O}$.
  • Xây dựng thuật toán tổng quát DBIM và thực hiện mô phỏng kiểm chứng với các tham số: tần số $f = 0,64\text{ MHz}$, đường kính vùng tán xạ 10 mm, độ tương phản tốc độ âm $\Delta c = 2%$, kích thước lưới $N = 16, 32, 64$.
  • Kiểm chứng mô hình mô phỏng thông qua tập dữ liệu thực nghiệm đo áp suất tán xạ từ quả bóng đường kính 0,3 mm trong môi trường nước ở nhiệt độ 18,1 °C với góc quét $120^\circ$. Kết quả so sánh cho thấy sự tương đồng giữa dạng sóng mô phỏng và thực nghiệm.

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 tập trung giải quyết bài toán tối ưu hóa quá trình tái tạo ảnh bằng thuật toán hai tần số DF-DBIM (Dual Frequency DBIM).

Các nội dung và kết quả chính:

  • Phương pháp luận: Quá trình tái tạo được chia làm hai giai đoạn trong tổng số 8 vòng lặp ($N_{iter} = 8$):
    • Giai đoạn đầu thực hiện $N_{f1}$ vòng lặp tại tần số thấp $f_1$ nhằm mở rộng miền hội tụ và ước tính biên dạng đối tượng thô.
    • Giai đoạn sau thực hiện $N_{f2} = 8 - N_{f1}$ vòng lặp tại tần số cao $f_2$ nhằm cải thiện độ phân giải không gian và độ sắc nét đường biên.
  • Khảo sát các kịch bản mô phỏng: Tác giả thiết lập 4 kịch bản khảo sát với số lượng máy phát $N_t$, máy thu $N_r$ và kích thước lưới $N$ khác nhau:
    • Kịch bản 1: $N_t = 16, N_r = 16, N = 22$.
    • Kịch bản 2: $N_t = 20, N_r = 20, N = 22$.
    • Kịch bản 3: $N_t = 24, N_r = 24, N = 22$.
    • Kịch bản 4: $N_t = 32, N_r = 32, N = 22$.
  • Xác định thời điểm chuyển giao tối ưu $N_{f1}$: Luận án tính toán lỗi chuẩn hóa (Normalized Error) qua các giá trị $N_{f1} \in [1, 7]$. Kết quả phân tích số liệu khẳng định:
    • Giá trị $N_{f1} = 2$ cho lỗi chuẩn hóa nhỏ nhất trên tất cả các kịch bản.
    • So với phương pháp DBIM đơn tần số truyền thống, thuật toán DF-DBIM với $N_{f1} = 2$ đạt tốc độ hội tụ nhanh hơn rõ rệt và giảm thiểu hiện tượng nhiễu đốm xung quanh vùng biên khối u.

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 nghiên cứu giải pháp giảm số lượng đầu dò và giảm độ phức tạp tính toán bằng cách kết hợp phương pháp DBIM với kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM).

Các đóng góp và số liệu quan trọng của chương:

  • Nguyên lý DCS: Thay vì sử dụng ma trận đo ngẫu nhiên phân bố chuẩn hoặc Bernoulli (CS thuần túy), thuật toán sử dụng chuỗi hỗn loạn sinh bởi ánh xạ logistic để xác định vị trí đặt máy thu một cách tiền định. Điều này cho phép cấu hình phần cứng cố định mà vẫn bảo đảm tính chất đẳng cự giới hạn (RIP) của ma trận nén.
  • Đánh giá hiệu suất DCS-DBIM so với DBIM truyền thống:
    • Cấu hình truyền thống (đo đều): $N_t = N_r = 20$ (tổng 400 phép đo) hoặc $N_t = N_r = 22$ (tổng 484 phép đo) hoặc $N_t = N_r = 30$ (tổng 900 phép đo).
    • Cấu hình đề xuất (DCS): Giảm xuống $N_t = N_r = 16$ (tổng 256 phép đo) hoặc $N_t = N_r = 20$ (tổng 400 phép đo).
    • Kết quả so sánh: Cấu hình đề xuất DCS-DBIM với 256 phép đo đạt lỗi chuẩn hóa tương đương hoặc thấp hơn cấu hình DBIM truyền thống 484 phép đo; cấu hình DCS-DBIM 400 phép đo cho chất lượng tái tạo vượt trội so với DBIM truyền thống 900 phép đo.
  • Thời gian tính toán và tỷ số nén: Thuật toán DCS-DBIM giải bài toán $l_1$-LSP giúp giảm đáng kể thời gian tạo ảnh sau 8 vòng lặp và duy trì độ ổn định khôi phục ngay cả khi giảm sâu số lượng đầu thu.
  • Đánh giá so sánh giữa DCS-DBIM và CS-DBIM: Ở các tỷ số nén khác nhau, DCS-DBIM cho sai số thặng dư tương đối (RRE) và lỗi chuẩn hóa tương đương với CS-DBIM ngẫu nhiên, nhưng vượt trội ở tính khả thi khi triển khai trên phần cứng thực tế.

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

Luận án đã đạt được các kết quả và đóng góp cụ thể sau:

Đóng góp về mặt lý luận và học thuật

  1. Phát triển thuật toán hai tần số DF-DBIM: Xác lập quy trình tạo ảnh kết hợp hai tần số cho siêu âm cắt lớp và giải quyết vấn đề xác định thời điểm chuyển giao tần số tối ưu ($N_{f1} = 2$ trong chu kỳ 8 vòng lặp), giúp thuật toán đạt tốc độ hội tụ cao và cải thiện độ phân giải không gian so với phương pháp đơn tần số.
  2. Đề xuất giải thuật lấy mẫu nén giả ngẫu nhiên DCS-DBIM: Ứng dụng thành công chuỗi ánh xạ logistic tiền định vào việc lấy mẫu không gian trong siêu âm cắt lớp, chứng minh tính tương đương về mặt toán học và chất lượng tái tạo ảnh giữa lấy mẫu giả ngẫu nhiên và lấy mẫu ngẫu nhiên thuần túy trong bài toán nghịch đảo sóng âm.

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

  1. Giảm thiểu độ phức tạp tính toán và số lượng đầu dò: Phương pháp DCS-DBIM cho phép giảm số lượng phép đo từ 484 xuống 256 phép đo (hoặc từ 900 xuống 400 phép đo) mà vẫn giữ nguyên độ chính xác tái tạo cấu trúc khối u, từ đó rút ngắn thời gian xử lý và giảm dung lượng bộ nhớ yêu cầu.
  2. Khả năng hiện thực hóa trên phần cứng: Việc thay thế phân bố ngẫu nhiên bằng chuỗi giả ngẫu nhiên tiền định loại bỏ yêu cầu sử dụng các mạch tạo số ngẫu nhiên phần cứng (HRNG/PRNG) phức tạp, tạo điều kiện thuận lợi cho việc thiết kế và chế tạo các hệ thống máy chụp siêu âm cắt lớp thế hệ mới.
  3. Kiểm chứng mô hình bằng dữ liệu thực nghiệm: Xây dựng mô hình mô phỏng tương thích và được kiểm chứng thông qua dữ liệu đo áp suất tán xạ thực tế với mẫu vật thể hình cầu kích thước 0,3 mm.

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

Hạn chế của luận án

  • Giới hạn không gian 2D: Nghiên cứu mới dừng lại ở việc mô phỏng và tái tạo lát cắt 2D trong hệ tọa độ trụ, chưa mở rộng sang mô hình tái tạo không gian 3D đầy đủ.
  • Giả thiết môi trường: Mô hình khảo sát lý tưởng hóa đối tượng là hình trụ tròn đặt trong môi trường nền đồng nhất (nước), bỏ qua sự phân bố phức tạp của các lớp mô tuyến vú xung quanh và các hiệu ứng phản xạ, khúc xạ mạnh.
  • Tham số tái tạo: Luận án mới chỉ tập trung khai thác tham số sự thay đổi tốc độ truyền âm ($\Delta c$), chưa tích hợp đồng thời các tham số âm học khác như hệ số suy hao âm và sự biến thiên mật độ mô vào giải thuật đề xuất.

Hướng nghiên cứu tiếp theo

  • Mở rộng thuật toán DF-DBIM và DCS-DBIM cho bài toán tái tạo ảnh 3D từ tập hợp các lát cắt 2D hoặc sử dụng khẩu độ hình trụ 3D.
  • Phát triển giải thuật tái tạo đa tham số đồng thời (tốc độ âm, độ suy hao và mật độ mô) để tăng cường thông tin chẩn đoán mô bệnh học.
  • Triển khai thuật toán trên các nền tảng phần cứng chuyên dụng (như FPGA hoặc GPU) và tích hợp vào hệ thống thu nhận dữ liệu siêu âm cắt lớp thực tế.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong các lĩnh vực:

  • Đối với nghiên cứu sinh và học viên cao học ngành Kỹ thuật điện tử, Kỹ thuật y sinh, Xử lý tín hiệu: Cung cấp khung toán học chi tiết về bài toán tán xạ ngược, phương pháp Moment (MoM), phương pháp điều chuẩn Tikhonov và kỹ thuật lấy mẫu nén ứng dụng trong bài toán ảnh y học.
  • Đối với các kỹ sư và nhà phát triển thiết bị y tế: Cung cấp giải pháp cấu hình phần cứng đầu dò theo chuỗi giả ngẫu nhiên (DCS) giúp tối ưu hóa số lượng kênh thu phát và giảm chi phí chế tạo phần cứng máy chụp siêu âm cắt lớp.
  • Phần nội dung đáng tham khảo nhất: Quy trình giải thuật DF-DBIM (Chương 3), phương pháp xây dựng ma trận lấy mẫu DCS dựa trên chuỗi logistic và kỹ thuật giải bài toán tối ưu $l_1$-LSP trong DCS-DBIM (Chương 4).

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

1. Tại sao kỹ thuật siêu âm cắt lớp (UST) dựa trên lý thuyết tán xạ lại có ưu thế hơn chụp X-quang tuyến vú và siêu âm B-mode trong phát hiện ung thư vú sớm?

Kỹ thuật chụp X-quang tuyến vú bị suy giảm hiệu quả ở phụ nữ dưới 50 tuổi do mô vú dày đặc không tạo đủ độ tương phản. Siêu âm B-mode chỉ cung cấp thông tin định tính về trở kháng âm và phụ thuộc vào kinh nghiệm chuyên gia. Ngược lại, kỹ thuật siêu âm cắt lớp dựa trên lý thuyết tán xạ ngược cho phép tái tạo định lượng các tham số cơ học (tốc độ truyền âm, độ suy hao), có khả năng phân biệt sự thay đổi vận tốc âm nhỏ và nhận diện các tổn thương, khối u có kích thước nhỏ dưới 5 mm (kích thước nhỏ hơn hoặc bằng bước sóng).

2. Phương pháp lặp vi phân Born (DBIM) có ưu điểm và nhược điểm gì so với phương pháp lặp Born (BIM)?

Trong phương pháp DBIM, hàm Green được cập nhật liên tục sau mỗi vòng lặp, giúp thuật toán có tốc độ hội tụ nhanh hơn so với BIM. Tuy nhiên, nhược điểm của DBIM là nhạy cảm hơn với nhiễu ở những vòng lặp đầu tiên và có độ phức tạp tính toán lớn ($O(N_t N_r N^2)$). Trong khi đó, BIM giữ cố định hàm Green nên có khả năng kháng nhiễu tốt hơn nhưng tốc độ hội tụ chậm và tốn nhiều vòng lặp hơn.

3. Trong giải thuật DF-DBIM, việc chuyển giao tần số được thực hiện như thế nào và giá trị vòng lặp tối ưu là bao nhiêu?

Giải thuật DF-DBIM chia tổng số 8 vòng lặp thành hai giai đoạn: giai đoạn đầu chạy $N_{f1}$ vòng lặp ở tần số thấp $f_1$ để định vị và đảm bảo hội tụ; giai đoạn sau chạy $N_{f2} = 8 - N_{f1}$ vòng lặp ở tần số cao $f_2$ để tăng độ phân giải biên ảnh. Qua khảo sát thực nghiệm số trên 4 kịch bản đo khác nhau, luận án xác định giá trị chuyển giao tối ưu là $N_{f1} = 2$. Tại điểm này, sai số chuẩn hóa của ảnh đạt giá trị nhỏ nhất và độ sắc nét cao nhất.

4. Kỹ thuật lấy mẫu nén giả ngẫu nhiên (DCS-DBIM) giải quyết bài toán phần cứng như thế nào so với kỹ thuật lấy mẫu nén ngẫu nhiên (CS-DBIM)?

Kỹ thuật CS truyền thống đòi hỏi phân bố vị trí máy thu ngẫu nhiên hoàn toàn, yêu cầu phải tích hợp các bộ sinh số ngẫu nhiên phần cứng (HRNG) hoặc bộ phát số giả ngẫu nhiên phức tạp, khó chế tạo hệ thống đầu dò cố định. Kỹ thuật DCS-DBIM sử dụng chuỗi ánh xạ logistic để tạo ra vị trí máy thu có tính chất giả ngẫu nhiên nhưng hoàn toàn tiền định về mặt toán học. Điều này cho phép cố định vị trí cơ khí của các đầu dò trên vành đo mà vẫn đáp ứng đầy đủ điều kiện phục hồi tín hiệu thưa với số lượng đầu thu giảm đáng kể.

5. Số lượng phép đo và thời gian tính toán được cải thiện như thế nào khi áp dụng DCS-DBIM?

Kết quả thực nghiệm số trong luận án chứng minh rằng phương pháp DCS-DBIM với 256 phép đo ($N_t = N_r = 16$) đạt độ chính xác khôi phục tương đương hoặc tốt hơn cấu hình DBIM truyền thống cần tới 484 phép đo ($N_t = N_r = 22$). Khi so sánh ở quy mô lớn hơn, DCS-DBIM với 400 phép đo cho chất lượng vượt trội so với DBIM truyền thống 900 phép đo. Việc cắt giảm hơn một nửa số phép đo giúp giảm mạnh kích thước ma trận hệ thống, tiết kiệm bộ nhớ và rút ngắn thời gian tái tạo ảnh sau 8 vòng lặp.

Kết luận

Luận án tiến sĩ của tác giả Trần Quang Huy đã nghiên cứu và phát triển thành công hai giải thuật xử lý tín hiệu cho kỹ thuật tạo ảnh siêu âm cắt lớp 2D: giải thuật hai tần số DF-DBIM và giải thuật lấy mẫu nén giả ngẫu nhiên DCS-DBIM. Công trình đã giải quyết được bài toán xác định thời điểm chuyển giao tần số tối ưu ($N_{f1} = 2$) để nâng cao tốc độ hội tụ và độ phân giải ảnh, đồng thời ứng dụng chuỗi logistic tiền định giúp giảm mạnh số lượng đầu dò thu phát và khối lượng tính toán. Các kết quả mô phỏng số và dữ liệu thực nghiệm kiểm chứng với vật thể mẫu 0,3 mm khẳng định tính đúng đắn của mô hình toán học, đóng góp cơ sở khoa học quan trọng cho việc phát triển các thế hệ thiết bị siêu âm cắt lớp trong chẩn đoán y sinh.