Nhiễu Rician khác với nhiễu Gauss thông thường ở điểm nào?

Nhiễu Gauss là nhiễu cộng tuyến tính, có phân bố đối xứng quanh kỳ vọng và hoàn toàn độc lập với cường độ tín hiệu gốc. Ngược lại, nhiễu Rician là nhiễu phi cộng xuất hiện khi trích xuất biên độ từ tín hiệu phức, có phân bố lệch phải và phụ thuộc phi tuyến vào giá trị tín hiệu thực.

Tại sao ảnh chụp cộng hưởng từ MRI lại xuất hiện nhiễu Rician?

Trong máy quét cộng hưởng từ, dữ liệu thô thu nhận trên miền tần số không gian k là tín hiệu phức với phần thực và phần ảo độc lập, mỗi phần đều chứa nhiễu nhiệt Gauss. Khi hiển thị ảnh lâm sàng, hệ thống thực hiện phép tính căn bậc hai của tổng bình phương phần thực và phần ảo, biến đổi nhiễu Gauss thành nhiễu Rician.

Khi nào nhiễu Rician suy biến thành phân bố Rayleigh?

Nhiễu Rician suy biến hoàn toàn thành phân bố Rayleigh khi cường độ tín hiệu thực tế bằng 0, điển hình là tại vùng nền không khí xung quanh cơ thể bệnh nhân trên ảnh cộng hưởng từ. Tại vùng này, giá trị đo được thuần túy là nhiễu và kỳ vọng của cường độ điểm ảnh bị nâng lên mức dương bằng khoảng 1,2533 lần độ lệch chuẩn nhiễu.

Làm thế nào để khử nhiễu Rician mà không làm mờ ranh giới cấu trúc giải phẫu?

Để khử nhiễu Rician hiệu quả, các nhà nghiên cứu thường áp dụng bộ lọc khối phi cục bộ Non-Local Means 3D thích ứng hoặc bộ lọc trong miền sóng con wavelet kết hợp hiệu chỉnh moment bậc hai. Các phương pháp này tận dụng tính tự tương đồng không gian và đặc tính đại số của moment bậc hai để triệt tiêu nhiễu mà vẫn giữ được độ sắc nét của các đường biên mô mềm.

Nhiễu Rician (Rician noise) là gì?

Nhiễu Rician (Rician noise) là mô hình phân bố xác suất mô tả sự biến thiên ngẫu nhiên của biên độ tín hiệu trong các hệ thống thu nhận dữ liệu phức chịu ảnh hưởng bởi nhiễu Gauss trắng cộng trên cả hai thành phần thực và ảo độc lập. Hiện tượng này đóng vai trò nền tảng trong kỹ thuật xử lý tín hiệu và tạo ảnh y tế, đặc biệt là ảnh chụp cộng hưởng từ khi tái tạo ảnh cường độ biên độ từ cuộn thu đơn kênh. Bài viết cung cấp phân tích toàn diện về nguồn gốc vật lý, biểu thức giải tích, các hành vi tiệm cận đặc thù, tác động gây sai lệch phép đo định lượng và các phương pháp lọc nhiễu thích ứng hiện đại.

Bản chất toán học và cơ chế hình thành

Trong các hệ thống đo lường vật lý và viễn thông hiện đại, tín hiệu thu nhận thường được biểu diễn dưới dạng số phức bao gồm hai thành phần vuông pha: thành phần đồng pha (phần thực) và thành phần vuông pha (phần ảo). Quá trình chụp ảnh cộng hưởng từ là ví dụ điển hình nhất cho cơ chế này. Dữ liệu thô thu nhận trên không gian k sau khi trải qua phép biến đổi Fourier ngược rời rạc sẽ tạo ra tín hiệu ảnh phức tại mỗi tọa độ không gian.

Trong quá trình thu nhận, các biến thiên nhiệt của cuộn thu và sự chuyển động nhiệt của các ion bên trong cơ thể người bệnh tạo ra nhiễu nhiệt ngẫu nhiên. Nhiễu nhiệt này được mô hình hoá chính xác bằng nhiễu Gauss trắng cộng (additive white Gaussian noise) với kỳ vọng bằng 0 và phương sai bằng $\sigma^2$ trên cả hai kênh đo lường độc lập. Giả sử tín hiệu phức thu được tại một điểm ảnh có dạng:

I = I_r + j I_i = (A \cos\theta + n_r) + j (A \sin\theta + n_i)

Trong biểu thức trên, đại lượng $A$ biểu diễn biên độ tín hiệu thực không nhiễu mang thông tin cấu trúc giải phẫu, đại lượng $\theta$ là pha của tín hiệu, trong khi $n_r$ và $n_i$ lần lượt là các thành phần nhiễu Gauss độc lập trên kênh thực và kênh ảo với cùng độ lệch chuẩn $\sigma$ . Do pha của tín hiệu thường chịu ảnh hưởng bởi sự không đồng nhất của từ trường tĩnh và các hiệu ứng sinh lý, các thuật toán hiển thị lâm sàng tiêu chuẩn thực hiện phép trích xuất độ lớn biên độ:

M = |I| = \sqrt{I_r^2 + I_i^2}

Phép lấy căn bậc hai của tổng bình phương hai biến ngẫu nhiên Gauss độc lập có kỳ vọng khác 0 là một phép biến đổi phi tuyến tính nghiêm ngặt. Phép biến đổi này làm biến đổi hoàn toàn bản chất thống kê của tín hiệu, biến thành phần nhiễu cộng tuyến tính ban đầu thành nhiễu phi cộng tuân theo phân bố Rician. Kết quả là giá trị đo lường $M$ không còn mang tính chất đối xứng và không còn độc lập với giá trị cường độ tín hiệu gốc $A$ .

Hàm mật độ xác suất và các tham số thống kê

Hàm mật độ xác suất (probability density function) của biến ngẫu nhiên biên độ $M$ tuân theo phân bố Rician được thiết lập dựa trên lý thuyết thống kê tín hiệu ngẫu nhiên do Stephen O. Rice xây dựng trong các công trình kinh điển về phân tích toán học nhiễu ngẫu nhiên. Đối với ảnh cộng hưởng từ, dạng giải tích tiêu chuẩn được xác định bởi công thức:

p(M | A, \sigma) = \frac{M}{\sigma^2} \exp\left(-\frac{M^2 + A^2}{2\sigma^2}\right) I_0\left(\frac{A M}{\sigma^2}\right), \quad M \ge 0

Trong phương trình xác suất này, các đại lượng toán học và vật lý được định nghĩa rõ ràng như sau:

Biến ngẫu nhiên $M$ là cường độ biên độ đo được tại điểm ảnh, luôn nhận giá trị không âm.
Tham số $A$ đại diện cho biên độ tín hiệu thực sự của mô sinh học khi không có nhiễu.
Tham số $\sigma$ là độ lệch chuẩn của nhiễu Gauss gốc trên mỗi kênh thu nhận thực và ảo.
Ký hiệu $I_0$ đại diện cho hàm Bessel biến đổi loại 1 bậc 0 (modified Bessel function of the first kind of order zero).

Hàm Bessel biến đổi loại 1 bậc 0 có thể được tính toán bằng chuỗi vô hạn hội tụ:

I_0(z) = \sum_{k=0}^{\infty} \frac{\left(\frac{1}{4} z^2\right)^k}{(k!)^2}

Sự hiện diện của hàm Bessel bậc 0 trong hàm mật độ xác suất khiến cho các phép tính kỳ vọng bậc 1 của $M$ trở nên rất phức tạp, đòi hỏi việc biểu diễn qua các hàm siêu hình học hợp lưu (confluent hypergeometric functions) hoặc đa thức Laguerre bậc nửa. Cụ thể, kỳ vọng toán học của biến ngẫu nhiên Rician được tính theo công thức:

\mathbb{E}[M] = \sigma \sqrt{\frac{\pi}{2}} L_{1/2}\left(-\frac{A^2}{2\sigma^2}\right)

Trong đó ký hiệu $L_{1/2}$ là đa thức Laguerre tổng quát. Trái ngược với kỳ vọng bậc 1 phức tạp, moment bậc 2 của phân bố Rician lại sở hữu một cấu trúc đại số tuyến tính vô cùng chuẩn mực và thuận tiện cho việc xử lý số học:

\mathbb{E}[M^2] = A^2 + 2\sigma^2

Hệ thức tuyến tính này cho thấy giá trị trung bình bình phương của biên độ quan sát được bằng tổng của bình phương tín hiệu thực cộng thêm 2 lần phương sai nhiễu kênh. Đây là một đặc tính mấu chốt được các nhà khoa học khai thác để thiết kế các thuật toán ước lượng phương sai nhiễu và hiệu chỉnh độ chệch mà không cần giải phương trình hàm Bessel phi tuyến phức tạp.

Hành vi tiệm cận và hiện tượng độ chệch Rician

Đặc tính thống kê của nhiễu Rician phụ thuộc sâu sắc vào tỷ số tín hiệu trên nhiễu cục bộ (SNR), được định nghĩa bằng tỷ số giữa biên độ tín hiệu thực và độ lệch chuẩn nhiễu $\text{SNR} = A/\sigma$ . Nghiên cứu thực nghiệm và lý thuyết của Gudbjartsson và Patz năm 1995 đã chỉ ra hai trạng thái tiệm cận cực hạn của phân bố này:

Trường hợp tỷ số tín hiệu trên nhiễu cao

Khi cường độ tín hiệu gốc vượt trội so với mức độ biến thiên của nhiễu nhiệt, cụ thể khi chỉ số SNR lớn hơn ngưỡng giá trị 3, hàm Bessel biến đổi loại 1 bậc 0 có thể xấp xỉ tiệm cận theo biểu thức:

I_0(z) \approx \frac{\exp(z)}{\sqrt{2\pi z}}

Khi thế biểu thức xấp xỉ này vào hàm mật độ xác suất ban đầu, phân bố Rician suy biến tiệm cận về phân bố chuẩn Gauss thông thường với kỳ vọng tiến gần về $A$ và phương sai xấp xỉ bằng $\sigma^2$ :

p(M | A, \sigma) \approx \frac{1}{\sqrt{2\pi\sigma^2}} \exp\left(-\frac{(M - A)^2}{2\sigma^2}\right)

Tại vùng mô có tín hiệu dồi dào trên ảnh cộng hưởng từ, nhiễu Rician có thể được xem xét tương đương với nhiễu Gauss cộng thông thường mà không gây ra sai lệch đáng kể cho các phép phân tích định tính của bác sĩ chẩn đoán hình ảnh.

Trường hợp nền không có tín hiệu

Tại các vùng không gian không có mô giải phẫu, chẳng hạn như vùng không khí bên ngoài cơ thể bệnh nhân, cường độ tín hiệu thực tế hoàn toàn triệt tiêu và đạt giá trị bằng 0. Khi tham số $A = 0$ , vì giá trị hàm Bessel $I_0(0) = 1$ , phân bố Rician suy biến hoàn toàn về phân bố Rayleigh:

p(M | 0, \sigma) = \frac{M}{\sigma^2} \exp\left(-\frac{M^2}{2\sigma^2}\right), \quad M \ge 0

Các đặc trưng thống kê của phân bố Rayleigh ở vùng nền được xác định tường minh bởi các đại lượng giải tích sau:

Giá trị kỳ vọng toán học tỷ lệ thuận với độ lệch chuẩn $\sigma$ theo hệ số căn bậc hai của pi chia 2 xấp xỉ bằng giá trị số 1.2533:
$\mathbb{E}[M] = \sigma \sqrt{\frac{\pi}{2}} \approx 1.2533 \sigma$
Phương sai của phân bố vùng nền bằng phương sai nhiễu kênh nhân với hiệu số giữa 2 và pi chia 2, tương ứng với hệ số số học xấp xỉ 0.4292:
$\text{Var}(M) = \sigma^2 \left(2 - \frac{\pi}{2}\right) \approx 0.4292 \sigma^2$

Hệ quả vật lý cực kỳ quan trọng của hiện tượng này là cường độ điểm ảnh tại vùng nền không bao giờ có kỳ vọng bằng 0, mà luôn bị đẩy lên một mức điện áp dương trung bình bằng 1.2533 lần độ lệch chuẩn $\sigma$ . Đây là nguyên nhân khiến nền đen của ảnh cộng hưởng từ luôn hiển thị lấm tấm các hạt xám sáng.

Vùng chuyển tiếp và độ chệch tín hiệu thấp

Trong khoảng chuyển tiếp quan trọng khi chỉ số SNR nhỏ hơn 2 hoặc nằm trong dải giá trị từ 1 đến 3, phân bố Rician thể hiện tính phi đối xứng mạnh mẽ. Vì biên độ là đại lượng không âm, các biến thiên âm của nhiễu bị triệt tiêu và phản xạ về phía dương qua phép lấy độ lớn. Điều này tạo nên độ chệch dương (positive bias) cố hữu, khiến giá trị trung bình đo được luôn lớn hơn cường độ thực $A$ . Hiện tượng độ chệch tín hiệu này làm giảm độ tương phản giữa các mô có cường độ tín hiệu yếu và làm sai lệch nghiêm trọng các phép đo động học định lượng.

So sánh nhiễu Rician với các mô hình phân bố thống kê lân cận

Để định vị chính xác vị trí của mô hình nhiễu Rician trong lý thuyết xác suất và xử lý ảnh, bảng tổng hợp dưới đây đối chiếu các đặc tính cốt lõi giữa phân bố Rician và các mô hình phân bố liên quan mật thiết:

Mô hình phân bố	Biểu thức toán học đặc trưng	Tính chất đối xứng	Bản chất phụ thuộc tín hiệu	Bối cảnh ứng dụng chủ yếu
Nhiễu Gauss (Gaussian)	Đối xứng hình chuông quanh kỳ vọng	Đối xứng hoàn toàn	Độc lập với cường độ tín hiệu gốc	Nhiễu nhiệt trên từng kênh phức, tín hiệu có SNR cao
Nhiễu Rician (Rice)	Phi đối xứng, chứa hàm Bessel bậc 0	Lệch phải ở SNR thấp, đối xứng ở SNR cao	Phụ thuộc phi tuyến vào biên độ tín hiệu gốc	Ảnh biên độ cộng hưởng từ cuộn đơn, kênh truyền sóng có tia trực tiếp
Nhiễu Rayleigh	Dạng suy biến của Rician khi tín hiệu bằng 0	Lệch dương mạnh, chặn dưới tại 0	Không có tín hiệu hữu ích, thuần túy là nhiễu	Vùng nền ảnh cộng hưởng từ, kênh fading đa đường che khuất hoàn toàn
Nhiễu Chi phi tâm (Non-central Chi)	Tổng bình phương đa kênh thu nhận	Phi đối xứng, phụ thuộc số bậc tự do	Phụ thuộc phi tuyến vào tín hiệu và số lượng cuộn thu	Hệ thống cộng hưởng từ đa cuộn thu ghép nối tổng bình phương

Tác động trong thực tiễn chẩn đoán hình ảnh và đo lường định lượng

Trong thực hành lâm sàng tại các bệnh viện và trung tâm nghiên cứu y sinh học, các máy quét chụp cộng hưởng từ phổ biến nhất hoạt động ở cường độ từ trường tĩnh tiêu chuẩn là 1.5 T và 3 T. Mặc dù từ trường cao giúp tăng cường đáng kể tỷ số tín hiệu trên nhiễu tổng thể, sự hiện diện của nhiễu Rician vẫn tạo ra nhiều rào cản kỹ thuật nghiêm trọng trong các ứng dụng chẩn đoán chuyên sâu:

Ảnh hưởng đến kỹ thuật chụp cộng hưởng từ khuếch tán (DWI và DTI): Trong kỹ thuật chụp khuếch tán sức căng, các gradient khuếch tán mạnh được áp dụng để khảo sát sự chuyển động nhiệt của phân tử nước trong mô não. Việc áp dụng gradient mạnh làm suy giảm mạnh mẽ biên độ tín hiệu, đẩy dữ liệu vào vùng SNR thấp dưới ngưỡng giá trị 2. Độ chệch dương Rician làm tăng nhân tạo cường độ tín hiệu bị suy giảm, dẫn đến việc ước lượng sai hệ số khuếch tán biểu kiến (apparent diffusion coefficient, ADC) và làm sai lệch hướng các bó sợi thần kinh chất trắng khi tái tạo đường dẫn truyền.
Sai lệch bản đồ thời gian hồi phục T1 và T2: Để định lượng chính xác các thông số vật lý mô học thông qua kỹ thuật lập bản đồ (mapping) định lượng, các mô hình toán học cần khớp đường cong suy giảm hàm mũ với dữ liệu đo được. Khi mức sàn của nhiễu bị nâng lên do độ chệch Rician ở các thời điểm trễ dài, đường cong suy giảm không còn tiến về 0 mà tiệm cận mức sàn dương, làm cho thời gian hồi phục ước lượng bị ngắn lại hoặc kéo dài bất thường.
Nhiễu hạt và suy giảm độ sắc nét đường biên: Tại các vùng chuyển tiếp ranh giới giữa mô tổn thương và mô lành có tỷ trọng tín hiệu thấp, nhiễu Rician tạo ra hiện tượng nhiễu hạt biến thiên phức tạp, che lấp các vi cấu trúc giải phẫu và làm giảm độ tin cậy của các thuật toán phân vùng tự động (image segmentation).

Các giải pháp lọc và khử nhiễu Rician

Do bản chất phi cộng và phụ thuộc tín hiệu của phân bố Rician, việc áp dụng trực tiếp các bộ lọc nhiễu Gauss truyền thống như bộ lọc trung bình tuyến tính hoặc bộ lọc Gauss không gian sẽ làm mờ nghiêm trọng các cạnh giải phẫu và bảo toàn nguyên vẹn độ chệch dương ở vùng nền. Để khắc phục triệt để thách thức này, cộng đồng nghiên cứu xử lý ảnh y tế quốc tế đã phát triển nhiều nhóm phương pháp lọc thích ứng chuyên biệt:

Phương pháp lọc trong miền Wavelet

Năm 1999, Robert D. Nowak đã công bố nghiên cứu bước ngoặt trên tạp chí quốc tế, phân tích cặn kẽ đặc tính thống kê của nhiễu Rician trong miền biến đổi sóng con (wavelet domain). Công trình của Nowak đã chứng minh rằng việc biến đổi dữ liệu bình phương biên độ cho phép tách biệt thành phần tín hiệu hữu ích khỏi thành phần nhiễu thông qua một hệ thức xấp xỉ tuyến tính hóa. Dựa trên cơ sở này, bộ lọc miền biến đổi wavelet thích ứng đã được thiết lập thành công, cho phép loại bỏ thành phần nhiễu hạt ngẫu nhiên đồng thời khôi phục lại độ sắc nét của các đường biên vi thể mà không làm biến dạng cấu trúc giải phẫu của não bộ.

Bộ lọc khối phi cục bộ Non-Local Means 3D

Khai thác tính chất tự tương đồng không gian trong thể tích hình học của cơ thể người, năm 2008, Pierrick Coupé và các cộng sự tại Pháp đã phát triển và tối ưu hóa giải thuật bộ lọc khối phi cục bộ (blockwise Non-Local Means filter) cho dữ liệu cộng hưởng từ ba chiều. Bằng cách so sánh độ tương đồng giữa các khối lân cận trong toàn bộ thể tích voxel thay vì chỉ xét từng điểm ảnh đơn lẻ, giải thuật của Coupé đã cải tiến trọng số lọc thích ứng với hàm mật độ xác suất Rician, cho phép triệt tiêu nhiễu vượt bậc ở cả các vùng mô có tỷ số SNR thấp mà vẫn giữ nguyên được độ sắc nét ranh giới mô mềm.

Mô hình hoá và ước lượng tham số phương sai

Năm 2016, hai nhà khoa học Santiago Aja-Fernández và Gonzalo Vegas-Sánchez-Ferrero đã xuất bản chuyên khảo học thuật toàn diện, hệ thống hóa toàn bộ các mô hình ước lượng tham số thống kê nhiễu trong ảnh cộng hưởng từ. Để giải quyết bài toán hiệu chỉnh độ chệch, các tác giả nhấn mạnh tầm quan trọng của việc ước lượng chính xác tham số phương sai $\sigma^2$ . Khi tham số độ lệch chuẩn $\sigma$ được trích xuất tin cậy từ vùng nền hoặc từ các phương pháp ước lượng moment cục bộ (LMMSE), giá trị tín hiệu thực $A$ có thể được tái cấu trúc chính xác thông qua phép hiệu chỉnh bình phương ngược:

\hat{A} = \sqrt{\max\left(M^2 - 2\sigma^2, 0\right)}

Điều kiện áp dụng và các hướng phát triển mở

Mặc dù mô hình phân bố Rician đã được kiểm chứng chuẩn xác cho các hệ thống thu nhận cuộn đơn truyền thống, sự phát triển vũ bão của công nghệ phần cứng cộng hưởng từ trong kỷ nguyên hiện đại đặt ra những giới hạn mới cho phạm vi áp dụng của mô hình này:

Hệ thống cuộn thu đa kênh song song: Trong các máy quét lâm sàng hiện đại, các kỹ sư sử dụng mảng cuộn thu gồm nhiều kênh độc lập để tăng tốc độ quét ảnh. Khi dữ liệu từ nhiều cuộn thu được kết hợp lại bằng kỹ thuật tổng bình phương (Sum of Squares, SoS), phân bố của ảnh biên độ tổng hợp không còn là phân bố Rician đơn thuần mà chuyển dịch thành phân bố Chi phi tâm (non-central Chi distribution) với số bậc tự do tỷ lệ thuận với số lượng kênh thu.
Tái tạo ảnh song song và nén cảm nhận: Khi sử dụng các kỹ thuật chụp ảnh tăng tốc cao cấp như SENSE, GRAPPA hoặc nén cảm nhận (compressed sensing), ma trận tái tạo kết hợp các phép nội suy phức tạp làm cho phương sai nhiễu không còn phân bố đồng đều mà biến thiên theo từng vị trí không gian (spatially variant noise). Trong bối cảnh này, việc áp dụng mô hình Rician đơn biến cố định trên toàn ảnh đòi hỏi phải được mở rộng thành các mô hình Rician thích ứng cục bộ theo bản đồ hình học g-factor.