Phương pháp dựa trên cường độ khác phương pháp dựa trên đặc trưng như thế nào?

Phương pháp dựa trên cường độ so khớp trực tiếp toàn bộ giá trị mức xám của các voxel mà không cần phân vùng trước, trong khi phương pháp dựa trên đặc trưng yêu cầu trích xuất các điểm mốc, đường biên hoặc bề mặt hình học tường minh.

Tại sao thông tin tương hỗ được ưu tiên cho đăng ký ảnh đa phương thức?

Thông tin tương hỗ đo lường mức độ phụ thuộc thống kê chung giữa các phân bố xác suất mức xám, hoạt động hiệu quả ngay cả khi mối quan hệ cường độ giữa các phương thức thu nhận như CT và MRI mang tính phi tuyến phức tạp.

Chiến lược kim tự tháp đa độ phân giải giải quyết vấn đề gì?

Chiến lược đa độ phân giải làm mịn và giảm kích thước ảnh ở nhiều tầng để tối ưu hóa trước trên ảnh độ phân giải thấp, giúp thuật toán tránh mắc kẹt vào các điểm cực trị cục bộ và tăng tốc độ hội tụ.

Hiện tượng trường biến thiên cường độ (bias field) ảnh hưởng thế nào đến độ chính xác?

Trường biến thiên cường độ tạo ra dải thay đổi mức xám chậm không đồng nhất trên ảnh chụp cộng hưởng từ, làm sai lệch biểu đồ tần suất mức xám chung và làm suy giảm độ chính xác của các hàm đo tương đồng nếu không được khử nhiễu tiền xử lý.

Phương pháp dựa trên cường độ là gì? Nguyên lý và ứng dụng

Phương pháp dựa trên cường độ (intensity-based method) là nhóm thuật toán tối ưu hóa trong xử lý ảnh số và thị giác máy tính, thực hiện đăng ký hoặc phân tích ảnh bằng cách so khớp trực tiếp toàn bộ giá trị mức xám của các điểm ảnh hoặc phần tử thể tích (voxel) giữa các tập dữ liệu mà không cần trích xuất trước các đặc trưng hình học tường minh. Kỹ thuật này đóng vai trò nền tảng trong việc đồng chỉnh không gian giữa các ảnh chụp khác thời điểm hoặc từ các thiết bị thu nhận khác nhau trong chẩn đoán y tế và viễn thám. Bài viết này trình bày bản chất toán học, các hàm đo độ tương đồng cường độ, thuật toán tối ưu hóa, phạm vi ứng dụng lâm sàng và các giới hạn kỹ thuật của phương pháp.

Bản chất và nguyên lý của phương pháp dựa trên cường độ

Trong bài toán đăng ký ảnh (image registration), mục tiêu là tìm kiếm một phép biến đổi hình học không gian phù hợp nhằm ánh xạ tọa độ của một ảnh nguồn (thường gọi là ảnh di động) sang hệ tọa độ của một ảnh đích (ảnh tham chiếu). Khác với các phương pháp dựa trên đặc trưng hình học vốn đòi hỏi phải phát hiện ranh giới, đường bao giải phẫu hoặc các điểm mốc đánh dấu nhân tạo, phương pháp dựa trên cường độ tiếp cận dữ liệu ảnh dưới góc độ các phân bố giá trị đo lường quang phổ hoặc độ cản xạ vật lý.

Phương pháp này tính toán một chỉ số định lượng độ tương đồng hoặc độ sai khác trên toàn bộ miền giao nhau của hai ảnh. Hàm mục tiêu được định dạng như một bài toán tối ưu hóa toán học liên tục, trong đó các tham số của phép biến đổi hình học được tinh chỉnh lặp qua các bước nhằm cực đại hóa độ tương đồng hoặc cực tiểu hóa sai số cường độ. Do tận dụng toàn bộ thông tin mức xám sẵn có của các voxel, phương pháp dựa trên cường độ hạn chế tối đa sai số tích lũy phát sinh từ giai đoạn phân vùng hoặc trích xuất đặc trưng hình học ban đầu.

Các độ đo tương đồng cường độ phổ biến

Việc lựa chọn hàm đo độ tương đồng phụ thuộc chặt chẽ vào việc hai ảnh thuộc cùng một phương thức thu nhận (đơn phương thức - monomodal) hay từ các nguồn vật lý khác biệt (đa phương thức - multimodal).

Độ đo cho đăng ký ảnh đơn phương thức

Khi hai ảnh được ghi nhận từ cùng một loại thiết bị thu nhận (ví dụ hai lần chụp cộng hưởng từ ở các mốc thời gian khác nhau), mối quan hệ giữa cường độ của các điểm ảnh tương ứng thường được giả định là đồng nhất hoặc có quan hệ tuyến tính.

Tổng bình phương sai phân (Sum of Squared Differences - SSD): Giả định rằng mức xám tại các vị trí tương ứng chỉ khác nhau bởi một thành phần nhiễu Gauss độc lập. Hàm sai số có dạng:
$SSD = \sum_{x} (I_1(x) - I_2(T(x)))^2$
Trong đó $I_1$ là cường độ ảnh tham chiếu, $I_2$ là cường độ ảnh di động, và $T(x)$ là phép biến đổi không gian tác động lên tọa độ $x$ . Phương pháp tìm kiếm giá trị tham số sao cho $SSD$ đạt cực tiểu.
Hệ số tương quan chéo chuẩn hóa (Normalized Cross-Correlation - NCC): Được áp dụng khi tồn tại sự suy giảm hoặc khuếch đại tuyến tính về cường độ giữa hai lần chụp:
$NCC = \frac{\sum_{x} (I_1(x) - \mu_1)(I_2(T(x)) - \mu_2)}{\sqrt{\sum_{x} (I_1(x) - \mu_1)^2 \sum_{x} (I_2(T(x)) - \mu_2)^2}}$
Trong đó $\mu_1$ và $\mu_2$ lần lượt biểu diễn cường độ trung bình của các ảnh trên miền không gian chồng lấn.

Độ đo cho đăng ký ảnh đa phương thức

Trong các ứng dụng y tế phối hợp hình ảnh chụp cắt lớp vi tính (CT) và cộng hưởng từ (MRI), giá trị cường độ biểu diễn các tính chất vật lý hoàn toàn khác biệt: CT đo mật độ cản tia X của mô xương, trong khi MRI phản ánh mật độ proton và thời gian hồi phục từ tính của mô mềm. Mối quan hệ mức xám giữa hai phương thức này mang tính phi tuyến phức tạp, khiến các độ đo tuyến tính như tương quan chéo hoàn toàn không còn khả năng hoạt động.

Năm 1995, Maes và cộng sự cùng Wells và cộng sự đã độc lập đề xuất ứng dụng độ đo thông tin tương hỗ (Mutual Information - MI) dựa trên lý thuyết thông tin Shannon trong xử lý ảnh y tế. Độ đo này xem xét mức xám của hai ảnh như các biến ngẫu nhiên và đo lường mức độ phụ thuộc thống kê giữa chúng.

Hàm thông tin tương hỗ giữa hai ảnh được tính toán thông qua entropy biên và entropy đồng thời:

I(A, B) = H(A) + H(B) - H(A, B)

Trong đó $H(A)$ và $H(B)$ là entropy biên biểu diễn mức độ bất định của phân bố mức xám trong từng ảnh, còn $H(A, B)$ là entropy đồng thời được xây dựng từ biểu đồ tần suất mức xám chung (joint histogram). Trong điều kiện miền không gian giao nhau cố định, khi hai ảnh được căn chỉnh chính xác, phân bố mức xám chung trở nên tập trung nhất, dẫn tới entropy đồng thời $H(A, B)$ đạt giá trị cực tiểu. Tuy nhiên, do diện tích vùng giao nhau thực tế luôn biến thiên trong quá trình biến đổi hình học, việc tối đa hóa thông tin tương hỗ $I(A, B)$ giúp bù trừ sự thay đổi của các entropy biên thay vì chỉ tối thiểu hóa entropy đồng thời đơn thuần.

Năm 1999, Studholme và cộng sự đề xuất biến thể thông tin tương hỗ chuẩn hóa (Normalized Mutual Information - NMI) nhằm khắc phục triệt để sự nhạy cảm của hàm đo trước sự thay đổi diện tích vùng giao nhau giữa hai ảnh:

NMI(A, B) = \frac{H(A) + H(B)}{H(A, B)}

Chỉ số $NMI$ này cung cấp tính ổn định cao trước sự thay đổi của miền chồng lấn giữa các khối dữ liệu thể tích.

Độ đo tương đồng	Ký hiệu	Giả định về mối quan hệ mức xám	Phạm vi áp dụng phù hợp
Tổng bình phương sai phân	SSD	Cường độ giống nhau cộng nhiễu Gauss ngẫu nhiên	Ảnh đơn phương thức (cùng loại thiết bị)
Tương quan chéo chuẩn hóa	NCC	Mối quan hệ biến thiên tuyến tính giữa các mức xám	Ảnh đơn phương thức có biến thiên độ sáng tuyến tính
Thông tin tương hỗ	MI	Phụ thuộc thống kê chung giữa các phân bố xác suất	Ảnh đa phương thức (CT - MRI, MRI - PET)
Thông tin tương hỗ chuẩn hóa	NMI	Phụ thuộc thống kê không phụ thuộc tỷ lệ vùng giao nhau	Ảnh đa phương thức có diện tích chồng lấn biến thiên

Không gian biến đổi và kỹ thuật tối ưu hóa

Quá trình đăng ký ảnh dựa trên cường độ liên tục cập nhật các tham số của mô hình hình học để tìm vị trí cực trị của hàm mục tiêu.

Mô hình biến đổi hình học

Biến đổi cứng (Rigid transformation): Bao gồm phép tịnh tiến và phép quay trong không gian ba chiều, bảo toàn khoảng cách giữa các điểm. Thích hợp cho đăng ký ảnh sọ não hoặc các cấu trúc xương cứng không bị biến dạng cơ học.
Biến đổi afin (Affine transformation): Mở rộng biến đổi cứng bằng cách bổ sung hệ số co giãn và biến dạng trượt, bảo toàn tính song song của các đường thẳng.
Biến đổi phi tuyến (Non-rigid / Deformable transformation): Sử dụng các mô hình trường dịch chuyển tự do dựa trên hàm spline (B-spline Free-Form Deformation) hoặc mô hình vật lý đàn hồi để mô phỏng sự chuyển động và biến dạng của các tạng mềm như gan, phổi hoặc tim trong chu kỳ thở.

Chiến lược tối ưu hóa đa độ phân giải

Do hàm tương đồng cường độ thường chứa nhiều điểm cực trị cục bộ (local extrema), các thuật toán tối ưu hóa chuẩn như dốc giảm dần (Gradient Descent), thuật toán Gauss-Newton hoặc phương pháp Powell rất dễ bị mắc kẹt tại các vị trí sai lệch nếu khởi tạo ban đầu cách xa điểm hội tụ toàn cục.

Để giải quyết thách thức này, phương pháp dựa trên cường độ thường kết hợp chiến lược kim tự tháp đa độ phân giải (multiresolution pyramid). Hình ảnh ban đầu được làm mịn và giảm kích thước ở nhiều mức độ khác nhau. Quá trình tối ưu hóa được thực thi trước tiên trên các phiên bản ảnh có độ phân giải thấp nhằm xác định gần đúng vị trí tổng thể, sau đó kết quả được chuyển tiếp làm điểm khởi tạo cho các tầng có độ phân giải cao hơn. Kỹ thuật này vừa giúp thuật toán thoát khỏi các bẫy cực trị cục bộ vừa tiết kiệm đáng kể thời gian tính toán.

Bước phát triển hiện đại: Đăng ký ảnh dựa trên cường độ bằng học sâu

Theo tổng quan của Zou và cộng sự (2022), sự kết hợp giữa phương pháp dựa trên cường độ truyền thống và kiến trúc mạng nơ-ron tích chập (CNN) đã tạo ra bước đột phá trong tốc độ xử lý. Trong các mô hình học sâu không giám sát hiện đại, mạng nơ-ron nhận đầu vào là cặp ảnh cần đăng ký và trực tiếp dự đoán trường dịch chuyển không gian dày đặc.

Thay vì cần nhãn chuẩn giải phẫu tốn kém, hàm mất mát (loss function) huấn luyện mạng được thiết lập trực tiếp từ các hàm đo tương đồng cường độ kinh điển như hệ số tương quan chéo cục bộ (Local Cross-Correlation) hoặc thông tin tương hỗ, kết hợp với thành phần phạt độ trơn của trường biến dạng. Khi mô hình đã hoàn tất quá trình huấn luyện ngoại tuyến, thời gian thực thi đăng ký cho một cặp ảnh thể tích ba chiều giảm từ hàng chục phút xuống chỉ còn vài giây, mở ra tiềm năng ứng dụng trực tiếp trong dẫn đường phẫu thuật thời gian thực.

Ứng dụng thực tế trong khoa học và kỹ thuật

Phương pháp dựa trên cường độ giữ vị trí trung tâm trong nhiều quy trình công nghệ cao:

Định vị và dẫn đường trong phẫu thuật thần kinh: Căn chỉnh ảnh chụp cộng hưởng từ mô mềm trước phẫu thuật với ảnh chụp cắt lớp vi tính xương trong phòng mổ, giúp bác sĩ phẫu thuật xác định chính xác vị trí ranh giới khối u mà không cần gắn khung cố định cơ học xâm lấn lên đầu bệnh nhân.
Lập kế hoạch xạ trị định vị thân (SBRT): Tích hợp dữ liệu từ ảnh chụp CT mô phỏng với ảnh cắt lớp phát xạ positron (PET) để vạch rõ thể tích đích sinh học của khối u, giảm thiểu liều phóng xạ tác động lên các mô lành lân cận.
Theo dõi diễn tiến bệnh mạn tính: Đăng ký chuỗi ảnh cộng hưởng từ não thu nhận định kỳ để định lượng mức độ teo hồi hải mã hoặc sự tiến triển của các tổn thương mất myelin trong bệnh thoái hóa thần kinh.
Xử lý ảnh viễn thám và trắc địa: Ghép nối và đồng chỉnh các khung ảnh vệ tinh đa phổ chụp ở các dải tần quang học và radar khẩu độ tổng hợp khác nhau nhằm theo dõi biến động lớp phủ mặt đất và xói lở bờ biển.

Ưu điểm và các hạn chế kỹ thuật

Ưu điểm nổi bật

Lợi thế lớn nhất của phương pháp dựa trên cường độ là tính tự động hóa hoàn toàn. Phương pháp loại bỏ sự can thiệp thủ công vốn tiềm ẩn nhiều sai số chủ quan của con người trong khâu chấm điểm mốc giải phẫu. Nghiên cứu thực nghiệm kinh điển của Wells và cộng sự (1996) trên ảnh thể tích sọ não đa phương thức đã chứng minh phương pháp cho phép đạt độ chính xác căn chỉnh ở mức dưới một voxel (subvoxel), tạo độ tin cậy cao cho các can thiệp ngoại khoa tinh vi.

Hạn chế và điều kiện áp dụng

Mặc dù có nhiều ưu điểm, phương pháp dựa trên cường độ vẫn tồn tại các giới hạn kỹ thuật cần được cân nhắc kỹ lưỡng:

Chi phí tính toán lớn: Do phải duyệt qua hàng triệu điểm ảnh hoặc voxel trong mỗi bước lặp của hàm mục tiêu, khối lượng phép tính là rất lớn, đặc biệt đối với các mô hình biến đổi phi tuyến nhiều bậc tự do.
Độ nhạy với trường biến thiên cường độ (Bias field): Trong chụp cộng hưởng từ, sự không đồng nhất của từ trường thường tạo ra một dải biến thiên mức xám chậm trên toàn bộ thể tích ảnh, làm sai lệch phân bố xác suất của biểu đồ tần suất chung và giảm độ chính xác của hàm thông tin tương hỗ nếu không qua bước chuẩn hóa tiền xử lý.
Yêu cầu về vùng chồng lấn tối thiểu: Khi hai ảnh ban đầu có khoảng cách lệch quá xa hoặc vùng giao nhau quá nhỏ, các hàm đo thống kê không có đủ lượng mẫu đại diện để ước lượng chính xác mật độ xác suất, dẫn tới thất bại hoàn toàn trong việc tìm hướng hội tụ.

Phân biệt với các khái niệm cùng tên trong các ngành khác

Cần lưu ý rằng cụm từ "phương pháp dựa trên cường độ" còn xuất hiện trong các ngữ cảnh học thuật khác với bản chất khoa học hoàn toàn tách biệt:

Trong quản lý môi trường và kinh tế carbon: Phương pháp dựa trên cường độ phát thải là cơ chế phân bổ hạn ngạch hoặc đặt mục tiêu giảm phát thải khí nhà kính dựa trên tỷ lệ phát thải trên một đơn vị sản lượng sản xuất (ví dụ tấn khí nhà kính trên một tấn sản phẩm), trái ngược với phương pháp đặt hạn mức trần phát thải tuyệt đối.
Trong kỹ thuật cảm biến quang học: Phương pháp đo dựa trên cường độ là kỹ thuật đo lường tín hiệu thông qua sự thay đổi biên độ công suất quang tán xạ ngược hoặc truyền qua, phân biệt với phương pháp đo giao thoa pha hoặc biến đổi bước sóng.

Mục từ này chỉ tập trung phân tích phương pháp dựa trên cường độ trong chuyên ngành xử lý ảnh và thị giác máy tính thuộc nhóm ngành kỹ thuật và công nghệ.