Tổng quan về luận án

Sự phát triển mạnh mẽ của công nghệ vũ trụ đã đưa các hệ thống vệ tinh nhỏ viễn thám (Small Earth Observation Satellites) trở thành công cụ quan trọng trong quản lý tài nguyên thiên nhiên, quan trắc môi trường và bảo đảm quốc phòng - an ninh. Tuy nhiên, trong môi trường không gian khắc nghiệt, các hệ thống quang học độ phân giải cao trên vệ tinh nhỏ phải đối mặt với sự suy giảm hiệu năng nghiêm trọng do bức xạ vũ trụ, chu kỳ nhiệt, hiện tượng thoát khí và sự lão hóa của cảm biến quang điện tử (CCD/CMOS). Khác với các vệ tinh quan sát Trái đất cỡ lớn vốn được trang bị hệ thống hiệu chuẩn chủ động trên bo mạch (on-board calibration systems) đắt tiền và cồng kềnh, các vệ tinh nhỏ bị giới hạn nghiêm ngặt về kích thước, khối lượng và công suất tiêu thụ (SWaP - Size, Weight, and Power).

Khoảng trống nghiên cứu (Research Gap) cốt lõi tồn tại trong y văn viễn thám hiện nay là sự thiếu vắng một quy trình khép kín kết hợp giữa phương pháp kiểm định gián tiếp trên quỹ đạo (in-flight vicarious calibration) không phụ thuộc phần cứng chuyên dụng với các thuật toán toán học phục hồi chất lượng dữ liệu ảnh (image restoration) về gần trạng thái ban đầu lúc mới phóng. Các công trình quốc tế như của Storey (2001) đối với Landsat-7 hay Léger và cộng sự (2002) đối với SPOT-5 phần lớn chỉ tập trung vào việc định lượng suy giảm chất lượng mà chưa tích hợp mô hình phục hồi bức xạ và không gian thích nghi theo chuỗi thời gian hoạt động.

Luận án tiến sĩ của Nghiên cứu sinh Thái Thị Bích Hồng, chuyên ngành Vật lí vô tuyến và điện tử (Mã số: 9440130.03) tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của TS. Bùi Trọng Tuyên và PGS. Vũ Anh Phi, đã giải quyết trọn vẹn bài toán khoa học này. Luận án đặt ra hai câu hỏi nghiên cứu then chốt:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để định lượng chính xác sự suy giảm chất lượng không gian và chất lượng bức xạ của thiết bị chụp ảnh quang học trên vệ tinh nhỏ viễn thám thông qua dữ liệu ảnh mặt đất mà không cần thiết bị hiệu chuẩn trên tàu?
  2. Câu hỏi nghiên cứu 2 (RQ2): Thuật toán toán học nào cho phép phục hồi đồng thời độ sắc nét không gian và bù trừ suy giảm phổ bức xạ của ảnh viễn thám về tiệm cận thời điểm kiểm tra ban đầu khi vệ tinh mới vào quỹ đạo (In-Orbit Test - IOT)?

Hai giả thuyết khoa học tương ứng được kiểm chứng:

  • Giả thuyết 1 (H1): "Chất lượng của thiết bị chụp ảnh quang học trên vệ tinh nhỏ viễn thám được kiểm định gián tiếp qua thông số MTF (chất lượng không gian) và tỷ lệ tín hiệu trên nhiễu SNR (chất lượng bức xạ) là phù hợp và đảm bảo độ chính xác."
  • Giả thuyết 2 (H2): "Chất lượng ảnh được phục hồi về gần thời điểm vệ tinh mới phóng lên quỹ đạo dựa trên kết quả kiểm định thiết bị chụp ảnh là khả thi đối với hệ thống vệ tinh nhỏ viễn thám."

Khung lý thuyết của nghiên cứu được xây dựng vững chắc trên nền tảng Lý thuyết hệ thống quang học dịch chuyển tuyến tính bất biến (Linear Shift-Invariant - LSI), Quang học Fourier và Lý thuyết xử lý tín hiệu ngẫu nhiên. Phạm vi thực nghiệm của luận án tập trung vào hệ thống vệ tinh nhỏ viễn thám đầu tiên của Việt Nam - VNREDSat-1 (khối lượng ~120 kg, kích thước $1000 \times 600 \times 600\text{ mm}$, hoạt động trên quỹ đạo đồng bộ mặt trời SSO ở độ cao 680 km, độ nghiêng $98,128^\circ$, độ phân giải mặt đất GSD 2,5 m ở kênh toàn sắc PAN và 10 m ở 4 kênh đa phổ MS gồm Blue, Green, Red, NIR, lượng tử hóa 10-bit). Tập dữ liệu thực nghiệm bao gồm chuỗi 66 cảnh ảnh đa thời gian từ năm 2013 đến 2022 chụp các bãi kiểm chuẩn quốc tế Salon de Provence (Pháp), bãi kiểm định Buôn Ma Thuột (Đắk Lắk, Việt Nam) và hàng chục khu vực địa hình đặc thù trên lãnh thổ Việt Nam.

Literature Review và Positioning

Nghiên cứu về đánh giá và phục hồi chất lượng ảnh viễn thám quang học trên thế giới được định hình bởi hai dòng tiếp cận chính: đánh giá chất lượng hệ thống trên quỹ đạo (In-flight Image Quality Assessment) và kỹ thuật phục hồi ảnh ngược (Inverse Image Restoration).

Trong trường phái đánh giá chất lượng không gian, hàm truyền điều biến (Modulation Transfer Function - MTF) tại tần số Nyquist ($f_{Nyquist} = 0.5\text{ cycles/pixel}$) là chỉ số chuẩn mực xác định khả năng phân tách chi tiết biên cạnh. Storey (2001) đã tiên phong sử dụng cây cầu Lake Pontchartrain Causeway (Louisiana, Hoa Kỳ) như một bãi kiểm định dạng xung để ước tính MTF cho cảm biến ETM+ trên vệ tinh Landsat-7, chỉ ra sự suy thoái hiệu suất quang học theo thời gian. Kohm (2004) áp dụng thuật toán cạnh gờ trên các mục tiêu nhân tạo đen trắng để đánh giá vệ tinh OrbView-3. Đối với dòng vệ tinh SPOT, Léger và cộng sự (2002) đã kiểm tra tái hội tụ và đánh giá MTF của camera HRG và HRS trên SPOT-5 bằng phương pháp bãi chuẩn cạnh. Tuy nhiên, tiêu chuẩn quốc tế ISO 12233 về phương pháp cạnh nghiêng (slanted-edge method) áp dụng cho cảm biến ảnh đòi hỏi góc nghiêng của cạnh mục tiêu phải nằm trong khoảng hẹp dưới $10^\circ$ (khuyến nghị $5^\circ$), gây khó khăn lớn cho việc ứng dụng trên các vệ tinh nhỏ có góc quay chụp linh hoạt (lên tới $35^\circ$ dọc và $30^\circ$ ngang hướng bay như VNREDSat-1).

Đối với chất lượng bức xạ, tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR) được xem là thước đo hiệu suất thu nhận photon. Curran và Dungan (1989), Roger và Arnold (1996) đã phát triển các phương pháp ước tính nhiễu và SNR trên dữ liệu phổ AVIRIS. Hệ thống SPOT (Meygret et al., 2000) tách biệt nhiễu thành thành phần hàng và cột trên các ảnh đồng nhất như sa mạc hoặc băng tuyết. Gần đây hơn, Bouvet và cộng sự (2019) thiết lập mạng lưới hiệu chuẩn tự động RadCalNet nhằm cung cấp thông số phản xạ bề mặt đo đạc tại chỗ phục vụ hiệu chuẩn bức xạ gián tiếp.

Trong lĩnh vực phục hồi ảnh, Aouinti và cộng sự (2016) đề xuất phương pháp giải chập Wiener can thiệp kết hợp giải thuật di truyền (Genetic Algorithm) nhằm khôi phục ảnh vệ tinh bị mờ, song phương pháp này đòi hỏi giả định khắt khe về mật độ phổ công suất (Power Spectral Density - PSD) của ảnh gốc. Shao và cộng sự (2012) áp dụng phương pháp đồng dư pha (phase congruency) và chuẩn biến phân toàn phần (Total Variation - TV) để phục hồi ảnh vệ tinh độ phân giải cao dựa trên MTF ước tính từ hàm lan truyền cạnh (Edge Spread Function - ESF).

Tại Việt Nam, các nghiên cứu trước đây chủ yếu tập trung vào bình diện xử lý tăng cường ảnh ở mức sản phẩm thương mại 3A, 3B nhằm phục vụ mục đích phân loại lớp phủ hoặc chiết xuất chuyên đề (Lương Chính Kế, 2011; Chu Hải Tùng, 2014; Đặng Trường Giang, 2017). Công trình của Nghiêm Văn Tuấn (2018) và Nguyễn Minh Ngọc (2023) đã bước đầu đặt nền móng đánh giá MTF và SNR cho vệ tinh VNREDSat-1. Tuy nhiên, luận án của tác giả Nguyễn Minh Ngọc (2023) mới chỉ dừng lại ở việc tính toán các chỉ số kỹ thuật và tăng cường độ tương phản hiển thị phục vụ khai thác thông thường, hoàn toàn chưa xây dựng mô hình toán học để giải bài toán nghịch đảo phục hồi chất lượng ảnh nguyên thủy về trạng thái quang phổ ban đầu.

Luận án của Thái Thị Bích Hồng định vị chính xác tại điểm giao thoa giữa vật lý vô tuyến, quang học điện tử và khoa học xử lý ảnh viễn thám. Bằng cách thiết lập cầu nối trực tiếp giữa giá trị MTF/SNR đo đạc thực nghiệm trên bãi chuẩn và các ma trận giải chập không gian cùng hàm hồi quy bù trừ bức xạ đa phổ, nghiên cứu này đại diện cho bước tiến tiên phong tại Việt Nam và tiệm cận các chuẩn mực xử lý dữ liệu vệ tinh tiên tiến của Cơ quan Vũ trụ Châu Âu (ESA) và NASA.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết tạo ảnh tuyến tính bất biến (LSI) và quang học Fourier trong điều kiện suy thoái vật lý của môi trường không gian quỹ đạo thấp (Low Earth Orbit - LEO). Mô hình toán học tổng quát của quá trình tạo ảnh quang học viễn thám bị suy giảm chất lượng được mô tả qua phương trình tích chập hai chiều:

$$g(x, y) = f(x, y) * h(x, y) + n(x, y)$$

Trong đó:

  • $g(x, y)$ là ma trận cường độ bức xạ của ảnh thu nhận bị suy thoái tại trạm mặt đất.
  • $f(x, y)$ là phân bố bức xạ thực tế của đối tượng cảnh trên bề mặt Trái đất.
  • $h(x, y)$ là hàm lan truyền điểm (Point Spread Function - PSF) đại diện cho toàn bộ suy giảm quang sai của thấu kính, gương phản xạ Cassegrain/Gregorian, độ rung động tư thế vệ tinh (jitter) và vết mờ do vận tốc bay.
  • $n(x, y)$ là thành phần nhiễu ngẫu nhiên cộng tính phát sinh từ nhiễu lượng tử hóa (quantization noise), nhiễu nhiệt (thermal noise) và tín hiệu tối (Dark Signal - DS) của cảm biến CCD.

Về mặt lý thuyết, biến đổi Fourier 2D của hàm $h(x, y)$ chính là hàm truyền quang học (Optical Transfer Function - OTF):

$$OTF(u, v) = \mathcal{F}{h(x, y)} = MTF(u, v) \cdot e^{i\Phi(u, v)}$$

Trong đó phần thực biên độ $MTF(u, v)$ phản ánh mức độ suy giảm độ tương phản của hệ thống ở từng tần số không gian $(u, v)$, và $\Phi(u, v)$ là hàm truyền pha (Phase Transfer Function - PTF).

Luận án đóng góp về mặt lý thuyết mô hình suy thoái bức xạ đa thời gian: Độ phản xạ và độ xám $DN_k(t)$ tại kênh phổ thứ $k$ không giữ nguyên giá trị theo thời gian $t$ mà suy thoái theo một hàm tuyến tính do hiện tượng giảm độ nhạy quang phổ và suy thoái lớp phủ quang học:

$$DN_k(t) = \alpha_k \cdot DN_k(t_0) + \beta_k(t) + \epsilon_k$$

Việc chứng minh bằng thực nghiệm rằng hệ số suy giảm bức xạ của vệ tinh nhỏ tuân theo quy luật suy giảm tuyến tính theo chuỗi ngày hoạt động là một đóng góp lý thuyết có giá trị, cho phép tham số hóa quá trình lão hóa cảm biến trong không gian.

Khung phân tích độc đáo

Khung phân tích của luận án được thiết kế đồng bộ theo chu trình khép kín gồm hai phân hệ: Phân hệ Kiểm định gián tiếp (Indirect Verification) và Phân hệ Phục hồi thích nghi (Adaptive Restoration).

+-------------------------------------------------------------------------------+
|                       DỮ LIỆU ẢNH THÔ MỨC 1A (VNREDSat-1)                     |
+---------------------------------------+---------------------------------------+
                                        |
     +----------------------------------+----------------------------------+
     |                                                                     |
     v                                                                     v
[PHÂN HỆ KIỂM ĐỊNH GIÁN TIẾP]                               [PHÂN HỆ PHỤC HỒI CHẤT LƯỢNG]
  1. Trích xuất ROI bãi kiểm chuẩn                            1. Phục hồi suy giảm không gian:
     (Salon de Provence / Buôn Ma Thuột)                         - Ước tính ma trận nghịch đảo PSF từ MTF
  2. Tính toán MTF (Chất lượng không gian):                      - Tích chập ngược / Wiener Deconvolution
     Otsu Threshold -> Cubic Spline Edge Detection            2. Phục hồi suy giảm bức xạ:
     -> ESF -> Savitzky-Golay Filter -> LSF -> FFT               - Trích xuất điểm sáng nhất trên ô mẫu chuẩn
     -> $MTF(f_{Nyquist})$                                       - Xây dựng mô hình hồi quy đa thời gian
  3. Tính toán SNR (Chất lượng bức xạ):                          - Tính toán sai lệch bù $\Delta_{Band}$
     Local Standard Deviation (LSD) trên ô mẫu                   - Bù giá trị bức xạ trực tiếp vào mảng Raster
     $SNR = \mu / \sigma$                                        (GDAL/NumPy 10-bit Engine)
     |                                                                     |
     +----------------------------------+----------------------------------+
                                        v
+-------------------------------------------------------------------------------+
|            SẢN PHẨM ẢNH PHỤC HỒI TIỆM CẬN THỜI ĐIỂM BAN ĐẦU (IOT)            |
+-------------------------------------------------------------------------------+

Khung phân tích này tích hợp ba nền tảng toán học:

  1. Thuật toán cạnh nghiêng cải tiến: Thay vì dùng đạo hàm vi phân trực tiếp (vốn khuếch đại nhiễu bậc cao), quy trình áp dụng phân ngưỡng Otsu tự động xác định đường biên cạnh, nội suy bậc ba (cubic spline) để tái tạo tọa độ dưới mức pixel (sub-pixel resolution), phân dải binning tạo hàm ESF và áp dụng bộ lọc Savitzky-Golay (cửa sổ 51 điểm, đa thức bậc 3) để làm mịn trước khi vi phân thành hàm LSF và biến đổi Fourier rời rạc (FFT 2048 điểm).
  2. Phương pháp độ lệch chuẩn cục bộ (LSD): Ước tính SNR thông qua việc tính toán giá trị trung bình $\mu$ và độ lệch chuẩn $\sigma$ trên các cửa sổ trượt đồng nhất của các ô mẫu kiểm chuẩn (ô màu trắng, xám nhạt, xám, đen).
  3. Đại số ma trận bù phổ đa kênh: Sử dụng công cụ lập trình Python kết hợp các thư viện xử lý dữ liệu không gian mã nguồn mở GDAL, NumPy, SciPy để can thiệp trực tiếp vào mảng giá trị số thực (DN array) của ảnh mức 1A, thực hiện phép dịch chuyển histogram bức xạ độc lập cho từng kênh phổ (Blue, Green, Red, NIR).

Điều kiện biên (Boundary Conditions) của khung phân tích được xác định rõ: Phương pháp áp dụng cho các hệ thống chụp ảnh quang học viễn thám có đáp ứng tuyến tính hoặc gần tuyến tính, làm việc trong dải phổ nhìn thấy và cận hồng ngoại (VNIR), trên các sản phẩm ảnh sơ cấp mức 1A chưa qua nắn chỉnh hình học làm méo mó giá trị bức xạ gốc.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng (Positivism) kết hợp phương pháp luận thực nghiệm toán - lý định lượng nghiêm ngặt (Quantitative Experimental Physics). Thiết kế nghiên cứu được tổ chức theo mô hình đánh giá đa tầng (Multi-level validation):

  • Tầng 1 - Mô phỏng và Lý thuyết: Thiết lập các mô hình toán học về quang sai thấu kính, lỗi đầu sóng (Wave-Front Error - WFE), hàm phân bố phản xạ hai chiều (BRDF), và mô hình suy giảm chất lượng do mờ và nhiễu.
  • Tầng 2 - Lập trình thuật toán: Xây dựng bộ công cụ phần mềm độc lập bằng ngôn ngữ Python với giao diện tương tác đồ họa (GUI) cho phép trích xuất vùng quan tâm (ROI), tự động hóa toàn bộ quy trình tính toán MTF, SNR và phục hồi ảnh.
  • Tầng 3 - Thực nghiệm kiểm chứng trên dữ liệu vệ tinh thực tế: Ứng dụng trên toàn bộ chuỗi dữ liệu ảnh vệ tinh VNREDSat-1 qua các mốc thời gian từ khi phóng (2013) đến năm thứ 9 hoạt động trên quỹ đạo (2022).

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm tuân thủ các tiêu chuẩn nghiêm ngặt về kiểm chuẩn viễn thám quốc tế:

  1. Chiết xuất vùng kiểm chuẩn: Sử dụng các bãi kiểm chuẩn nhân tạo có tọa độ địa lý chính xác tại Salon de Provence (Pháp) và Buôn Ma Thuột (Đắk Lắk). Cảnh ảnh được trích xuất ở mức 1A để bảo toàn nguyên vẹn giá trị đo bức xạ thô (Raw DN) mà không bị ảnh hưởng bởi quá trình tái lấy mẫu hình học (resampling).
  2. Tính toán hàm lan truyền cạnh (ESF): Đo đạc gradient độ xám qua biên cạnh phân cách giữa các ô có độ tương phản cao. Tọa độ cạnh được xác định chính xác ở mức dưới pixel bằng phép xấp xỉ hàm bậc ba:

$$f(y) = a y^3 + b y^2 + c y + d$$

  1. Làm mịn thích nghi và vi phân thành LSF: Ứng dụng bộ lọc Savitzky-Golay để loại bỏ nhiễu ngẫu nhiên trên ESF mà không làm suy giảm độ dốc thực tế của biên cạnh:

$$LSF(x) = \frac{d}{dx} \left[ ESF_{smooth}(x) \right]$$

  1. Biến đổi Fourier và chuẩn hóa MTF: Áp dụng thuật toán FFT 2048 điểm để chuyển đổi LSF sang miền tần số không gian và chuẩn hóa về biên độ đơn vị tại tần số 0:

$$MTF(f) = \frac{\left| \int_{-\infty}^{\infty} LSF(x) e^{-i 2\pi f x} dx \right|}{\int_{-\infty}^{\infty} LSF(x) dx}$$

Giá trị MTF tại tần số Nyquist ($f = f_{Nyquist}$) được trích xuất riêng biệt theo cả hai phương: dọc hướng bay (along-track) và ngang hướng bay (cross-track).

  1. Đo lường SNR bằng phương pháp LSD: Tính toán trên các ô mẫu đồng nhất có kích thước $N \times M$ pixel:

$$SNR = \frac{\mu_{ROI}}{\sigma_{ROI}} = \frac{\frac{1}{NM} \sum_{i=1}^N \sum_{j=1}^M DN(i,j)}{\sqrt{\frac{1}{NM-1} \sum_{i=1}^N \sum_{j=1}^M \left( DN(i,j) - \mu_{ROI} \right)^2}}$$

Data và phân tích

Tập dữ liệu nghiên cứu được chọn lọc bao gồm 66 cảnh ảnh viễn thám quang học VNREDSat-1 mức 1A:

  • 34 cảnh ảnh đa phổ và toàn sắc tại bãi kiểm định Salon de Provence (Pháp) chụp định kỳ từ 10/07/2015 đến 21/05/2022.
  • 2 cảnh ảnh toàn sắc tại bãi kiểm định Buôn Ma Thuột chụp ngày 14/11/2017 và 02/11/2018.
  • 30 cảnh ảnh thử nghiệm trên các vùng sinh thái đặc trưng của Việt Nam: Hà Nội, Đà Nẵng, Thủy điện Hòa Bình, Rừng quốc gia Chư Yang Sin, Sân bay Tân Sơn Nhất, Đảo Thổ Chu, Vùng biển Trường Sa (Scaborough).

Phần mềm và môi trường tính toán: Toàn bộ thuật toán được phát triển trên nền tảng Python 3.8, sử dụng các thư viện chuyên dụng: osgeo.gdal (xử lý tệp raster không gian đa kênh định dạng GeoTIFF), numpy (tính toán đại số ma trận đa chiều), scipy.interpolate và scipy.signal (bộ lọc Savitzky-Golay, nội suy spline bậc ba), cv2 (thuật toán phân ngưỡng Otsu) và matplotlib (trực quan hóa đồ thị phổ và raster histogram).

Độ tin cậy của mô hình được đối chuẩn (cross-validation) với kết quả công bố của tác giả Nguyễn Minh Ngọc (2023) trên cùng tập dữ liệu VNREDSat-1 chụp bãi kiểm định Buôn Ma Thuột các năm 2017 và 2018.

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án mang lại 4 phát hiện khoa học mang tính đột phá với các bằng chứng thực nghiệm định lượng rõ ràng:

  1. Sự ổn định vượt bậc của chất lượng không gian (MTF) trên hệ thống quang học VNREDSat-1 qua 9 năm hoạt động: Kết quả đo đạc MTF tại tần số Nyquist theo phương pháp cạnh nghiêng cải tiến cho thấy: Vào năm 2015, MTF đạt giá trị $0,121$ (dọc hướng bay) và $0,125$ (ngang hướng bay). Đến ngày 21/05/2022 (sau 9 năm hoạt động), giá trị MTF vẫn duy trì ở mức trên $0,082$. Dù có sự suy giảm tự nhiên so với thời điểm kiểm tra nghiệm thu trên quỹ đạo (IOT - In-Orbit Test: $MTF \approx 0,20$), toàn bộ các giá trị MTF đo đạc qua chuỗi thời gian đều vượt tiêu chuẩn thiết kế ban đầu của nhà sản xuất (yêu cầu kỹ thuật thiết kế: $MTF \ge 0,080$).

  2. Hiệu suất đo bức xạ (SNR) thực tế vượt trội so với yêu cầu thiết kế: Thực nghiệm tính toán SNR trên các ô mẫu phản xạ chuẩn tại bãi kiểm định Buôn Ma Thuột cho thấy các giá trị thực đo đạt mức rất cao:

    • Ô màu trắng (năm 2017): $SNR = 146$ (năm 2018: $SNR = 145$).
    • Ô màu xám nhạt (năm 2017): $SNR = 161$ (năm 2018: $SNR = 117$).
    • Ô màu xám (năm 2017): $SNR = 175$ (năm 2018: $SNR = 149$).
    • Ô màu đen (năm 2017): $SNR = 380$ (năm 2018: $SNR = 410$). Toàn bộ các giá trị này đều vượt xa ngưỡng thiết kế ban đầu ($SNR_{thiết\ kế} \ge 100$), khẳng định cảm biến của vệ tinh VNREDSat-1 duy trì chất lượng đo bức xạ tốt sau nhiều năm làm việc trong môi trường bức xạ vũ trụ.

  3. Quy luật suy giảm bức xạ không đồng đều giữa các kênh phổ: Phân tích chuỗi số liệu độ xám của điểm ảnh sáng nhất trên các ô mẫu kiểm chuẩn tại Salon de Provence qua các năm cho thấy tốc độ lão hóa quang học có sự phân hóa rõ rệt giữa các bước sóng: Kênh cận hồng ngoại (NIR) chịu mức suy giảm bức xạ mạnh nhất, tiếp theo là kênh Đỏ (Red), kênh Xanh lam (Blue) và ít suy giảm nhất là kênh Xanh lục (Green). Cụ thể, trên cảnh ảnh ngày 21/05/2022, độ lệch suy giảm bức xạ cần bù đắp cho từng kênh được xác định chính xác là:

    • Kênh Blue (Band 1): $\Delta_{Blue} = -46\text{ DN}$
    • Kênh Green (Band 2): $\Delta_{Green} = -31\text{ DN}$
    • Kênh Red (Band 3): $\Delta_{Red} = -48\text{ DN}$
    • Kênh NIR (Band 4): $\Delta_{NIR} = -61\text{ DN}$

  4. Khả năng phục hồi dữ liệu ảnh về trạng thái ban đầu bằng thuật toán toán học: Sau khi áp dụng mô hình bù bức xạ tuyến tính đa kênh và giải chập không gian dựa trên ma trận PSF thực tế, các biểu đồ phân bố độ xám (Raster Histogram) của 32 cảnh ảnh thử nghiệm trên toàn lãnh thổ Việt Nam đã được tái cấu trúc hoàn chỉnh. Các dải giá trị bị bão hòa hoặc nén dải động được mở rộng trở lại tương đương với phân bố histogram của các cảnh ảnh chụp năm 2013-2014, loại bỏ hiệu ứng mờ nhòe và tăng cường đáng kể độ tương phản chi tiết bề mặt.

Bảng đối sánh kết quả thực nghiệm với nghiên cứu tiền nhiệm

Dưới đây là bảng trích xuất dữ liệu đối sánh trực tiếp giữa kết quả của luận án với công trình của tác giả Nguyễn Minh Ngọc (2023) trên cùng tập dữ liệu ảnh VNREDSat-1 chụp bãi kiểm chuẩn Buôn Ma Thuột:

Thông số đo đạc Đối tượng / Ô mẫu Mốc thời gian Kết quả đề xuất của Luận án Kết quả của N. M. Ngọc (2023) Sai số khác biệt
MTF tại Nyquist Dọc hướng bay Năm 2017 0,118 0,121 2,48%
Ngang hướng bay Năm 2017 0,120 0,122 1,64%
Dọc hướng bay Năm 2018 0,115 0,119 3,36%
Ngang hướng bay Năm 2018 0,117 0,119 1,68%
SNR Bức xạ Ô màu trắng Năm 2017 146 148 1,35%
Ô màu trắng Năm 2018 145 147 1,36%
Ô màu xám nhạt Năm 2017 161 156 +3,21%
Ô màu xám nhạt Năm 2018 117 111 +5,41%
Ô màu xám Năm 2017 175 170 +2,94%
Ô màu xám Năm 2018 149 139 +7,19%
Ô màu đen Năm 2017 380 372 +2,15%
Ô màu đen Năm 2018 410 402 +1,99%

Sự tương đồng cao giữa hai phương pháp độc lập (sai số dưới 7,2%) đã chứng minh độ tin cậy khoa học tuyệt đối của thuật toán đề xuất, đồng thời phương pháp của luận án vượt trội hơn nhờ khả năng chuyển tiếp trực tiếp sang khâu phục hồi ảnh.

Implications đa chiều

  • Về mặt lý thuyết và phương pháp luận: Luận án khẳng định tính khả thi của việc sử dụng các bãi kiểm chuẩn mặt đất kết hợp xử lý thuật toán tinh vi để thay thế hoàn toàn các hệ thống đo lường phần cứng phức tạp trên vệ tinh nhỏ, mở ra hướng đi tiết kiệm chi phí cho các chương trình không gian của các nước đang phát triển.
  • Về mặt ứng dụng thực tiễn: Giải pháp phần mềm được phát triển cho phép tích hợp trực tiếp vào quy trình xử lý dữ liệu tại Trạm thu nhận ảnh viễn thám mặt đất của Việt Nam, giúp kéo dài thời gian khai thác hữu ích của các vệ tinh quan sát Trái đất đã vượt quá tuổi thọ thiết kế lý thuyết (VNREDSat-1 có tuổi thọ thiết kế 5 năm nhưng đã hoạt động hiệu quả hơn 9 năm).
  • Về mặt chính sách: Cung cấp cơ sở khoa học xác đáng cho Cục Viễn thám Quốc gia, Viện Hàn lâm Khoa học và Công nghệ Việt Nam (VAST) trong việc ban hành quy chuẩn kỹ thuật quốc gia về kiểm định định kỳ chất lượng dữ liệu ảnh vệ tinh quang học.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra 3 hạn chế chính:

  1. Sự phụ thuộc vào bãi kiểm định nhân tạo: Quy trình đo MTF và SNR yêu cầu phải có dữ liệu chụp các bãi kiểm chuẩn chuyên dụng (như Salon de Provence hoặc Buôn Ma Thuột) với các điều kiện phản xạ bề mặt được đo đạc chính xác. Điều này làm giới hạn tần suất kiểm định khi vệ tinh không thể chụp bãi chuẩn do mây che phủ hoặc ưu tiên các nhiệm vụ khẩn cấp.
  2. Giả định tuyến tính của hệ thống quang học: Mô hình phục hồi bức xạ giả định suy giảm tuyến tính đồng nhất trên toàn dải cảm biến, chưa mô hình hóa đầy đủ các phi tuyến cục bộ do hiệu ứng nhiệt độ cảm biến thay đổi đột ngột giữa vùng sáng và vùng tối của quỹ đạo.
  3. Chưa xử lý triệt để ảnh hưởng của khí quyển: Các hệ số bù bức xạ hiện tại phản ánh tổng hòa của suy giảm thiết bị và suy giảm do tán xạ khí quyển trung bình, chưa tích hợp mô hình đo đạc độ dày quang học sol khí (AOD) thời gian thực tại thời điểm chụp.

Chương trình nghiên cứu tương lai (Future Research Agenda) bao gồm 4 hướng đi cụ thể:

  • Hướng 1: Mở rộng thuật toán kiểm định gián tiếp trên các bãi kiểm định giả bất biến tự nhiên (Pseudo-Invariant Calibration Sites - PICS) như sa mạc Sahara, lòng chảo muối Uyuni để tăng tần suất đánh giá mà không phụ thuộc bãi nhân tạo.
  • Hướng 2: Ứng dụng mạng nơ-ron tích chập sâu (Deep Learning CNN/Transformer) và kỹ thuật giải chập thời gian thực tăng tốc bằng GPU (đã được tác giả bước đầu công bố quốc tế) để phục hồi tự động ảnh vệ tinh mà không cần ước tính trước PSF.
  • Hướng 3: Xây dựng mô hình mô phỏng quang sai và chất lượng ảnh cho các dự án chùm vệ tinh nhỏ viễn thám thế hệ mới của Việt Nam (như VNREDSat-2, chùm vệ tinh CubeSat) ngay từ giai đoạn thiết kế tiền khả thi.
  • Hướng 4: Tích hợp mô hình hiệu chỉnh bức xạ khí quyển 6S (Second Simulation of a Satellite Signal in the Solar Spectrum) kết hợp đồng thời với thuật toán bù suy thoái thiết bị.

Tác động và ảnh hưởng

Luận án tạo ra những tác động đa chiều và sâu rộng:

  • Tác động học thuật: Thiết lập chuẩn mực nghiên cứu mới trong lĩnh vực vật lý vô tuyến và viễn thám quang học tại Việt Nam; cung cấp tài liệu tham khảo nền tảng cho các chương trình đào tạo sau đại học về công nghệ vũ trụ và kỹ thuật địa tin học.
  • Chuyển đổi công nghiệp không gian: Giảm thiểu chi phí chế tạo vệ tinh nhỏ bằng cách chuyển dịch gánh nặng kiểm định từ phần cứng (on-board hardware) sang phần mềm xử lý mặt đất (ground software algorithms), thúc đẩy ngành công nghiệp vệ tinh thương mại giá rẻ (NewSpace).
  • Tác động an ninh - quốc phòng và giám sát quốc gia: Phục hồi độ sắc nét của các cảnh ảnh vệ tinh chụp khu vực biên giới, hải đảo (như Trường Sa, Thổ Chu), giúp nâng cao độ chính xác trong việc phát hiện, nhận dạng tàu thuyền và biến động công trình ven biển.
  • Lợi ích kinh tế - xã hội: Tiết kiệm ngân sách hàng triệu USD thông qua việc tối ưu hóa và kéo dài tuổi thọ sử dụng của nguồn dữ liệu vệ tinh quốc gia, hạn chế việc phải mua dữ liệu ảnh thương mại từ nước ngoài.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giảng viên ngành Vật lý Vô tuyến, Kỹ thuật Không gian, Trắc địa - Viễn thám: Tiếp cận một khung lý thuyết hoàn chỉnh và mã nguồn thuật toán Python minh bạch để phát triển các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư vận hành hệ thống vệ tinh (Satellite Operations Engineers): Ứng dụng trực tiếp quy trình kiểm định MTF/SNR định kỳ để giám sát "sức khỏe" của tải trọng quang học trên quỹ đạo.
  • Các Trung tâm R&D Vũ trụ (như VNSC, Viện Công nghệ Vũ trụ): Sử dụng các bài học kinh nghiệm và mô hình suy thoái bức xạ để tối ưu hóa thiết kế quang - cơ - điện cho các thế hệ vệ tinh tiếp theo.
  • Cơ quan Quản lý Nhà nước (Bộ Tài nguyên và Môi trường, Bộ Công an, Bộ Quốc phòng): Được cung cấp nguồn dữ liệu ảnh viễn thám đã phục hồi với độ chính xác và độ tin cậy cao, phục vụ hiệu quả công tác quản lý tài nguyên, quy hoạch lãnh thổ và an ninh quốc gia.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Luận án đã mở rộng Lý thuyết Hệ thống Quang học Dịch chuyển Tuyến tính Bất biến (LSI) và Quang học Fourier sang miền ứng dụng suy thoái vật lý dài hạn ngoài không gian. Bằng cách thiết lập mối quan hệ hàm số giữa sự suy giảm độ nhạy phổ đa kênh $\Delta_{Band}(t)$ với thời gian hoạt động trên quỹ đạo và phản ánh trực tiếp sự suy giảm này vào mô hình phân rã ma trận suy biến tích chập quang học, luận án đã biến các chỉ số kiểm định MTF/SNR từ vai trò đánh giá thụ động thành tham số điều khiển chủ động cho thuật toán phục hồi chất lượng ảnh ngược.

2. Đột phá phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế tiền nhiệm là gì?

So với phương pháp cạnh nghiêng tiêu chuẩn ISO 12233 (vốn bị giới hạn ở góc nghiêng $<10^\circ$ như trong nghiên cứu của Magna, 2010) và phương pháp bãi chuẩn dạng xung của Storey (2001 trên Landsat-7), luận án đã đột phá khi tích hợp thuật toán phân ngưỡng Otsu tự động kết hợp nội suy bậc ba (cubic spline) và bộ lọc làm mịn Savitzky-Golay bậc 3 trên miền gradient ESF. Quy trình này triệt tiêu hiệu ứng răng cưa (aliasing) do lấy mẫu thiếu và loại bỏ nhiễu vi sai bậc cao, cho phép đo MTF chính xác ngay cả trên các cảnh ảnh chụp nghiêng góc lớn ($>30^\circ$) của vệ tinh nhỏ.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất từ dữ liệu ảnh VNREDSat-1?

Phát hiện bất ngờ nhất là sự duy trì hiệu suất quang học vượt bậc của hệ thống sau 9 năm hoạt động (2013-2022). Mặc dù MTF suy giảm từ mức IOT ($0,20$) xuống khoảng $0,082 - 0,121$, nó vẫn nằm trên ngưỡng thiết kế ban đầu ($0,080$). Đồng thời, phát hiện về sự suy giảm bức xạ mạnh nhất ở kênh NIR ($\Delta_{NIR} = -61\text{ DN}$) so với kênh Green ($\Delta_{Green} = -31\text{ DN}$) đã cung cấp bằng chứng thực nghiệm quan trọng về sự lão hóa quang phổ không đồng nhất của lớp kính lọc và đi-ốt quang silicon trong môi trường bức xạ LEO.

4. Luận án có cung cấp quy trình và mã nguồn tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết mã nguồn lập trình trên ngôn ngữ Python trong các phụ lục chuyên sâu, sử dụng các thư viện chuẩn như GDAL, NumPy, SciPy, OpenCV. Quy trình trích xuất ROI, tính toán ESF, LSF, FFT MTF, tính toán SNR theo độ lệch chuẩn cục bộ (LSD) và thuật toán bù bức xạ ma trận raster đều được cấu trúc hóa theo dạng module rõ ràng, cho phép tái lập kết quả $100%$ trên bất kỳ tập dữ liệu ảnh viễn thám nào có định dạng GeoTIFF chuẩn.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tới được tác giả định hướng theo 3 trụ cột: (1) Tự động hóa kiểm định thời gian thực trên quỹ đạo (On-orbit Real-time Edge Computing) tích hợp vi xử lý AI nhúng; (2) Phát triển các mô hình học sâu bán giám sát (Semi-supervised Deep Learning) để giải bài toán siêu phân giải (Super-resolution) và khử mờ thích nghi không cần bãi chuẩn nhân tạo; (3) Mở rộng khung kiểm định cho các chùm vệ tinh siêu nhỏ (CubeSat Constellations) và vệ tinh siêu phổ (Hyperspectral Satellites).

Kết luận

Luận án tiến sĩ của tác giả Thái Thị Bích Hồng là một công trình khoa học công phu, mẫu mực, kết hợp nhuần nhuyễn giữa lý thuyết vật lý vô tuyến - quang điện tử hiện đại với kỹ thuật lập trình xử lý ảnh không gian tiên tiến. Sáu đóng góp cốt lõi mang tính bước ngoặt của luận án bao gồm:

  1. Hoàn thiện cơ sở khoa học và quy trình kiểm định gián tiếp trên quỹ đạo cho thiết bị chụp ảnh quang học vệ tinh nhỏ thông qua hai thông số trụ cột: MTF (chất lượng không gian) và SNR (chất lượng bức xạ).
  2. Phát triển thuật toán cạnh nghiêng cải tiến ứng dụng phân ngưỡng Otsu, nội suy spline bậc ba và bộ lọc Savitzky-Golay, giải quyết triệt để bài toán khử nhiễu vi sai khi tính toán ESF và LSF.
  3. Xây dựng thành công bộ công cụ phần mềm độc lập bằng Python/GDAL, cho phép tự động hóa quy trình đánh giá và xử lý dữ liệu ảnh vệ tinh mức 1A.
  4. Chứng minh định lượng sự ổn định chất lượng của vệ tinh VNREDSat-1 sau 9 năm hoạt động, khẳng định hệ thống luôn duy trì các chỉ số MTF và SNR vượt tiêu chuẩn thiết kế ban đầu ($MTF \ge 0,080$; $SNR \ge 100$).
  5. Thiết lập mô hình phục hồi bức xạ đa kênh thích nghi, bù trừ chính xác sự suy thoái độ xám cho từng kênh phổ riêng biệt ($\Delta_{Blue} = -46$, $\Delta_{Green} = -31$, $\Delta_{Red} = -48$, $\Delta_{NIR} = -61$), tái tạo thành công phân bố histogram tiệm cận thời điểm mới phóng cho 32 cảnh ảnh thực nghiệm.
  6. Đặt nền móng phương pháp luận vững chắc cho công tác vận hành, khai thác và kiểm định chất lượng các thế hệ vệ tinh viễn thám tiếp theo của Việt Nam, đóng góp trực tiếp vào sự nghiệp phát triển khoa học công nghệ vũ trụ và bảo đảm an ninh - quốc phòng của đất nước.