Hình ảnh radar (radar imaging) là kỹ thuật tạo lập ảnh bề mặt hoặc mục tiêu thông qua việc phát và thu nhận các chùm sóng điện từ trong dải tần số vi ba (microwave), đo lường độ phản xạ tán xạ ngược (backscatter) và độ trễ pha thời gian để tái tạo bản đồ không gian hai chiều hoặc ba chiều với độ phân giải cao độc lập với điều kiện chiếu sáng và thời tiết. Khác với cảm biến quang học thụ động phụ thuộc vào ánh sáng mặt trời và dễ bị che khuất bởi mây mù, radar là hệ thống cảm biến chủ động (active remote sensing) có khả năng chụp ảnh xuyên mây, hoạt động liên tục cả ngày lẫn đêm và nhạy cảm với cấu trúc hình thái học, độ nhám bề mặt và hàm lượng độ ẩm điện môi của vật thể. Kể từ công trình lịch sử của Cutrona và cộng sự (1966) trên Proceedings of the IEEE đặt nền móng xử lý đồng pha cho radar khẩu độ tổng hợp (Synthetic Aperture Radar - SAR), công nghệ tạo ảnh radar đã có những bước nhảy vọt. Tổng quan chuẩn mực của Moreira và cộng sự (2013) trên IEEE Geoscience and Remote Sensing Magazine và các công trình kinh điển về giao thoa radar (InSAR) của Rosen và cộng sự (2000) cùng Bamler & Hartl (1998) trên Inverse Problems đã hệ thống hóa toàn diện lý thuyết SAR, kỹ thuật phân cực Polarimetric SAR, và quan trắc biến dạng mặt đất độ chính xác milimet. Bài viết này trình bày toàn diện về nguyên lý vật lý vi ba, kỹ thuật tạo ảnh SAR hai chiều, nén xung tín hiệu, giao thoa InSAR, phân cực đa kênh, cùng các ứng dụng viễn thám thực tiễn tại Việt Nam.
Nguyên lý Vật lý Vi ba và Cơ chế Tán xạ Ngược
Hình ảnh radar được hình thành dựa trên sự tương tác giữa sóng điện từ vi ba và các phần tử vật chất trên bề mặt:
1. Các Băng tần Radar Vi ba Phổ biến
- Băng X (bước sóng ngắn xentimet): Độ phân giải không gian rất cao, tán xạ chủ yếu tại bề mặt tán lá trên cùng, bề mặt băng tuyết và các công trình nhân tạo. Ứng dụng phổ biến trên các vệ tinh TerraSAR-X, COSMO-SkyMed.
- Băng C (bước sóng xentimet trung bình): Băng tần tiêu chuẩn cho quan sát biển, giám sát lũ lụt, băng biển và thảm phủ thực vật vừa phải. Đại diện tiêu biểu là chùm vệ tinh Sentinel-1 của Cơ quan Vũ trụ Châu Âu (ESA) và RADARSAT.
- Băng L (bước sóng đeximet dài): Bước sóng dài có khả năng xuyên thấu qua tán rừng nhiệt đới dày đặc và lớp đất mặt khô cằn để tương tác với thân cây và nền đất, lý tưởng cho ước tính sinh khối rừng và quan trắc đứt gãy kiến tạo (ALOS-2 PALSAR).
- Băng P (bước sóng dài nhất): Khả năng xuyên thấu sâu nhất qua tán rừng rậm và lớp đất dưới sâu, phục vụ nghiên cứu sinh khối rừng toàn cầu (nhiệm vụ BIOMASS).
2. Cơ chế Tán xạ Ngược (Backscattering Mechanisms)
Cường độ tín hiệu phản xạ lại máy thu radar phụ thuộc vào cơ chế tán xạ đặc trưng của từng loại bề mặt:
- Tán xạ bề mặt (Surface Scattering / Rough Surface): Chi phối bởi độ nhám bề mặt và hằng số điện môi. Bề mặt nhẵn (mặt nước phẳng, đường băng) đóng vai trò gương phản xạ hướng sóng ra xa máy thu (cho ảnh tối đen); bề mặt nhám gồ ghề tạo tán xạ khuếch tán đa hướng (cho ảnh sáng vừa).
- Tán xạ thể (Volume Scattering): Xảy ra khi sóng vi ba thâm nhập vào môi trường không đồng nhất như tán lá rừng ngập mặn, tuyết xốp hoặc thảm thực vật dày, tán xạ qua lại nhiều lần giữa các cành lá trước khi phản xạ về máy thu.
- Tán xạ góc đôi (Double-bounce Scattering / Corner Reflector): Sóng phản xạ liên tiếp giữa hai bề mặt trực giao (mặt đất và tường nhà, hoặc mặt nước và thân cây đước ngập mặn), tạo ra tín hiệu tán xạ ngược cực mạnh (các điểm sáng rực trên ảnh radar đô thị).
Nguyên lý Radar Khẩu độ Tổng hợp (SAR) theo Moreira (2013) và Cutrona (1966)
Để đạt độ phân giải không gian cao mà không cần dùng ăng-ten vật lý khổng lồ dài hàng kilômét, kỹ thuật SAR đã tổng hợp khẩu độ ảo thông qua chuyển động của bệ mang vệ tinh hoặc máy bay:
1. Độ phân giải cự ly (Range Resolution) và Kỹ thuật Nén Xung
Được xác định bởi thời gian tồn tại của xung phát hoặc độ rộng băng tần của tín hiệu điều tần tuyến tính (Chirp pulse):
$$\delta_r = \frac{c}{2 B \sin \theta}$$
Trong đó $c$ là vận tốc ánh sáng, $B$ là độ rộng băng tần của xung chirp, và $\theta$ là góc tới (incidence angle). Kỹ thuật nén xung (matched filtering) theo Cutrona và cộng sự (1966) cho phép nén một xung dài có công suất phát lớn thành một xung cực ngắn, đạt độ phân giải cự ly cao mà vẫn đảm bảo tỷ số tín hiệu trên tạp âm (SNR).
2. Độ phân giải phương vị (Azimuth Resolution) và Xử lý Lọc Doppler
Theo công trình nền tảng của Cutrona và cộng sự (1966) và Moreira và cộng sự (2013), bằng cách lưu giữ lịch sử pha kết hợp của tín hiệu phản xạ Doppler khi bệ mang di chuyển dọc theo quỹ đạo, thuật toán xử lý SAR tổng hợp một ăng-ten ảo có chiều dài tỷ lệ thuận với cự ly và bước sóng. Độ phân giải phương vị lý thuyết đạt được là:
$$\delta_{az} = \frac{L}{2}$$
Độ phân giải phương vị hoàn toàn độc lập với khoảng cách cự ly và bước sóng, chỉ phụ thuộc vào chiều dài vật lý $L$ của ăng-ten thực. Ăng-ten vật lý càng ngắn thì độ phân giải hình ảnh SAR thu được càng mịn.
Bảng đối chiếu Các Chế độ Hoạt động của Hệ thống SAR Hiện đại
| Chế độ chụp ảnh | Độ phân giải không gian | Độ rộng dải quét (Swath) | Ứng dụng và Mục tiêu điển hình |
|---|---|---|---|
| Spotlight | Rất cao (cấp dưới mét) | Hẹp | Giám sát công trình hạ tầng, mục tiêu quân sự, tàu thuyền |
| Stripmap | Trung bình (vài mét) | Vừa phải | Lập bản đồ địa hình, giám sát đứt gãy địa chất |
| ScanSAR | Thấp (hàng chục mét) | Rất rộng | Giám sát băng biển toàn cầu, quan sát hải dương học |
| TOPSAR (Sentinel-1) | Trung bình | Rộng | Quan sát bề mặt Trái Đất định kỳ, can thiệp biến dạng InSAR |
Radar Khẩu độ Tổng hợp Đa Phân cực (PolSAR)
Kỹ thuật phân cực radar (Polarimetric SAR) khai thác đầy đủ thông tin vector điện trường của sóng vi ba:
- Ma trận Tán xạ Sinclair: Mỗi pixel ảnh radar đa phân cực được biểu diễn bằng một ma trận tán xạ phức $2 \times 2$ bao gồm 4 kênh phân cực cơ bản: $HH$ (phát ngang - thu ngang), $HV$ (phát ngang - thu dọc), $VH$ (phát dọc - thu ngang), và $VV$ (phát dọc - thu dọc).
- Phân rã Tán xạ Pauli và Freeman-Durden: Theo tổng quan của Moreira và cộng sự (2013), các giải thuật phân rã mục tiêu (target decomposition) phân tách ma trận hiệp phương sai thành ba thành phần tán xạ vật lý độc lập: tán xạ bề mặt (Odd-bounce), tán xạ góc đôi (Double-bounce), và tán xạ thể ngẫu nhiên (Volume scattering), cho phép phân loại lớp phủ bề mặt và rừng tự động với độ chính xác cao.
Kỹ thuật Giao thoa Radar (InSAR) và Đo Biến dạng theo Rosen (2000) & Bamler (1998)
Giao thoa radar InSAR (SAR Interferometry) khai thác sự sai biệt pha giữa hai ảnh SAR chụp cùng một khu vực:
- Giao thoa đồ (Interferogram): Được tạo ra bằng cách nhân phức ảnh thứ nhất với ảnh liên hợp phức của ảnh thứ hai. Pha giao thoa $\Delta \phi$ chứa thông tin độ cao địa hình và độ dịch chuyển biến dạng: $$\Delta \phi = \phi_{topo} + \phi_{def} + \phi_{atm} + \phi_{noise}$$
- Mở pha giao thoa (Phase Unwrapping): Giải quyết bài toán mơ hồ pha tuần hoàn trong khoảng $[-\pi, \pi]$ để tái tạo giá trị pha liên tục tuyệt đối theo mô hình toán học nghịch đảo của Bamler & Hartl (1998) thông qua các thuật toán nhánh cắt (Branch Cut) hoặc dòng chi phí tối thiểu (Minimum Cost Flow).
- Giao thoa vi sai DInSAR và PSInSAR: Loại bỏ thành phần địa hình từ mô hình số độ cao (DEM) chuẩn để đo lường biến dạng lún sụt do khai thác nước ngầm hoặc trồi sụt do động đất với độ chính xác cấp milimet (Rosen và cộng sự, 2000).
- Kỹ thuật Chuỗi Thời gian Điểm Tán xạ Ổn định (PSInSAR / SBAS): Tận dụng các điểm phản xạ mạnh và ổn định theo thời gian (mái nhà bê tông, kết cấu thép, mỏm đá lộ thiên) để loại trừ triệt để ảnh hưởng của tầng đối lưu khí quyển, cung cấp biểu đồ lún sụt liên tục nhiều năm.
Ranh giới, Giới hạn Kỹ thuật và Hiện tượng Nhiễu trong Ảnh Radar
Việc diễn giải và xử lý ảnh radar đòi hỏi phải nhận thức rõ các hạn chế và hiện tượng hình học đặc thù:
- Nhiễu đốm (Speckle Noise): Hiện tượng giao thoa kết hợp của nhiều tán xạ nhỏ trong một pixel gây ra biến thiên ngẫu nhiên cường độ sáng tối giống như hạt muối tiêu, làm giảm khả năng phân biệt chi tiết biên ảnh. Xử lý giảm nhiễu đốm đòi hỏi các bộ lọc không gian thích nghi như bộ lọc Lee, Frost hoặc Gamma-MAP.
- Hiện tượng méo hình học do góc chiếu bên: Bao gồm hiện tượng dồn bóng (Foreshortening), bóng đè đỉnh núi (Layover) và vùng khuất bóng radar (Radar Shadow) tại các khu vực địa hình đồi núi dốc đứng.
- Hiện tượng mất tương quan pha (Decorrelation): Mất tương quan thời gian do thảm thực vật thay đổi hoặc mất tương quan không gian do khoảng cách đường đáy (baseline) quá lớn, làm giảm độ tin cậy của phép đo InSAR.
Ứng dụng Thực tiễn của Viễn thám Radar tại Việt Nam
Tại Việt Nam, công nghệ hình ảnh radar đóng vai trò chiến lược trong giám sát tài nguyên thiên nhiên và phòng chống thiên tai:
- Giám sát Sụt lún Đô thị và Đồng bằng sông Cửu Long: Sử dụng chuỗi dữ liệu Sentinel-1 và kỹ thuật PSInSAR để định lượng tốc độ sụt lún nền đất tại Thành phố Hồ Chí Minh, Hà Nội và Cần Thơ do khai thác nước ngầm quá mức và áp lực tải trọng xây dựng.
- Cảnh báo và Bản đồ hóa Vùng ngập Lũ lụt: Tận dụng khả năng chụp ảnh xuyên mây của sóng vi ba trong mùa mưa bão để khoanh vùng diện tích ngập úng tại các tỉnh duyên hải miền Trung và hạ lưu sông Hồng, hỗ trợ cứu hộ cứu nạn kịp thời.
- Giám sát Nông nghiệp và Ước tính Năng suất Lúa: Theo dõi sự biến thiên tán xạ ngược theo chu kỳ sinh trưởng của cây lúa nước tại vựa lúa Đồng bằng sông Cửu Long.
- Giám sát An ninh Biển đảo và Tàu thuyền: Ứng dụng thuật toán phát hiện tàu thuyền (CFAR) trên ảnh SAR băng C và băng X phục vụ quản lý giao thông hàng hải và giám sát vùng đặc quyền kinh tế.