Tổng quan về luận án

Nghiên cứu về chất lượng hình ảnh từ các hệ thống quan sát Trái đất đóng vai trò then chốt trong kỹ thuật trắc địa, bản đồ và viễn thám hiện đại. Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Minh Ngọc với đề tài "Nghiên cứu đánh giá chất lượng ảnh viễn thám quang học của Việt Nam" (Chuyên ngành: Kỹ thuật Trắc địa – Bản đồ, Mã số: 9520503) được thực hiện tại Trường Đại học Mỏ - Địa chất dưới sự hướng dẫn của PGS. Trần Vân Anh và TS. Nguyễn Xuân Lâm. Công trình đặt nền móng khoa học và phương pháp luận độc lập đầu tiên cho việc kiểm định, chuẩn hóa và giám sát chất lượng dữ liệu ảnh vệ tinh quang học viễn thám tại Việt Nam trong suốt vòng đời hoạt động trên quỹ đạo.

flowchart TD
    subgraph S1["Dữ liệu mức sơ cấp (Level 0)"]
        A1["Ảnh đại dương / Ban đêm"] --> B1["Hiệu chỉnh dòng tối (DS)"]
        A2["Ảnh sa mạc PICS (Libya/Algeria)"] --> B2["Cân bằng hồi đáp điểm ảnh (PRNU)"]
    end
    
    subgraph S2["Dữ liệu hình học mức 1A (Level 1A)"]
        B1 & B2 --> C1["Ảnh bãi kiểm định Buôn Ma Thuột"]
        C1 --> D1["Chiết tách cạnh tối ưu (Canny)"]
        D1 --> E1["Xác định ESF -> LSF -> MTF"]
        C1 --> D2["Phân tích vùng đồng nhất (Local Variance)"]
        D2 --> E2["Tính toán tỉ lệ tín hiệu trên nhiễu (SNR)"]
    end
    
    subgraph S3["Tầng chuẩn hóa & Ứng dụng"]
        E1 & E2 --> F["Đánh giá đối chuẩn thiết kế (8 năm On-orbit)"]
        F --> G["Phân cấp chất lượng thành lập bản đồ địa hình (1:5.000 - 1:50.000)"]
    end

Trong bối cảnh hệ thống vệ tinh nhỏ quan sát Trái đất đầu tiên của Việt Nam là VNREDSat-1 (Vietnam Natural Resources, Environment and Disaster Monitoring Satellite) được phóng lên quỹ đạo từ năm 2013, các tác động cộng hưởng từ môi trường không gian khắc nghiệt—bao gồm bức xạ vũ trụ, chu kỳ biến thiên nhiệt độ cực đoan, sự lão hóa của cảm biến quang điện tích hợp thời gian trễ (Time Delay Integration Charge Coupled Device - TDI-CCD) và sự suy giảm cơ-quang—đã dẫn đến sự suy giảm tất yếu về hiệu năng bức xạ và độ sắc nét không gian. Trong khi các nhà sản xuất nước ngoài (Astrium/Airbus) chỉ chuyển giao mô hình trong giai đoạn tiền phóng (pre-launch) và nghiệm thu đầu quỹ đạo (in-orbit commissioning), các đơn vị vận hành tại Việt Nam đối mặt với một khoảng trống nghiên cứu (research gap) mang tính cấu trúc: thiếu vắng một quy trình định lượng vật lý khách quan ở mức dữ liệu thô (Level 0 và Level 1A) phù hợp với cơ sở hạ tầng thực địa trong nước, thay vì chỉ dừng lại ở các kiểm tra hình học và mây ở các mức sản phẩm thương mại bậc cao (Level 2A, 3A, 3B) theo các quy định hiện hành như Thông tư 10/2015/TT-BTNMT hay Thông tư 52/2016/TT-BTNMT.

Để giải quyết khoảng trống học thuật và công nghệ này, luận án tập trung giải quyết hai câu hỏi nghiên cứu trung tâm:

  1. RQ1: Làm thế nào để xác định và tính toán chính xác các thông số vật lý đặc trưng nhất cho chất lượng bức xạ và chất lượng không gian của ảnh viễn thám quang học trong điều kiện cơ sở hạ tầng kỹ thuật và bãi kiểm định tại Việt Nam?
  2. RQ2: Quy luật suy thoái chất lượng của hệ thống chụp ảnh VNREDSat-1 sau thời gian dài hoạt động trên quỹ đạo diễn ra như thế nào, và làm sao để chuẩn hóa các tham số này phục vụ nhu cầu biên tập bản đồ địa hình chuyên ngành?

Các giả thuyết nghiên cứu tương ứng được xác lập:

  • Giả thuyết H1: Bộ đôi chỉ số vật lý gồm Tỉ lệ tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR) đại diện cho chất lượng bức xạ và Hàm truyền điều biến (Modulation Transfer Function - MTF) đại diện cho chất lượng không gian là tập thông số tối thiểu và tối ưu cho phép đánh giá toàn diện hệ thống vệ tinh viễn thám quang học mà không cần tiếp cận trực tiếp cấu trúc phần cứng.
  • Giả thuyết H2: Việc ứng dụng thuật toán phi tuyến Canny trong chiết tách cạnh trên ảnh bãi kiểm định nhân tạo mặt đất sẽ triệt tiêu được nhiễu cấu trúc và biến thiên phản xạ cục bộ, mang lại độ chính xác dưới điểm ảnh (sub-pixel) vượt trội so với các thuật toán đạo hàm tuyến tính truyền thống trong việc xác định MTF tại tần số Nyquist.

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 Dịch chuyển Tuyến tính Bất biến (Linear Shift-Invariant - LSI Systems), Quang học Fourier và Lý thuyết Xử lý Nhiễu Thống kê (Statistical Signal Processing). Nghiên cứu thực hiện đánh giá toàn diện trên chuỗi dữ liệu thực tế kéo dài 8 năm của VNREDSat-1 (độ phân giải toàn sắc PAN 2,5 m; đa phổ MS 10 m), kết hợp đo đạc quang phổ phản xạ bề mặt thực địa bằng thiết bị chuyên dụng tại bãi kiểm định Buôn Ma Thuột (tỉnh Đắk Lắk) trong các năm 2017 và 2018, đối chiếu với bãi kiểm định quốc tế Salon de Provence (Cộng hòa Pháp) và các bãi kiểm định giả bất biến (Pseudo-Invariant Calibration Sites - PICS) tại sa mạc Algeria, sa mạc Libya và vùng biển Đại Tây Dương.

Literature Review và Positioning

Nghiên cứu kiểm định chất lượng cảm biến viễn thám trên quỹ đạo (vicarious on-orbit calibration) trên thế giới đã hình thành hai luồng học thuật chính:

Luồng nghiên cứu thứ nhất tập trung vào đánh giá chất lượng bức xạ và tỉ lệ tín hiệu trên nhiễu (SNR). Nghiên cứu tiên phong của Curcio và cộng sự (1997) tại NASA đã phát triển phương pháp địa thống kê (geostatistical semi-variogram) trên dữ liệu AVIRIS nhằm bóc tách nhiễu cảm biến khỏi phương sai không gian tự nhiên. Chander và cộng sự (2013) cùng nhóm nghiên cứu cảm biến Operational Land Imager (OLI) trên Landsat-8 đã chuẩn hóa phương pháp đánh giá bức xạ chéo (cross-calibration) sử dụng các bề mặt sa mạc đồng nhất tuyệt đối tại Bắc Phi (PICS) với độ bất định dưới 3%. Trong khi đó, đối với dòng vệ tinh SPOT và Sentinel-2 (Gascon et al., 2017), Cơ quan Vũ trụ Châu Âu (ESA) thiết lập giao thức đo SNR định kỳ hàng tháng bằng thiết bị khuếch tán ánh sáng mặt trời (solar diffuser) kết hợp với chụp tín hiệu tối ban đêm trên biển Đại Tây Dương.

Luồng nghiên cứu thứ hai tập trung vào hàm truyền điều biến (MTF) và độ sắc nét không gian. Boreman (2001) cùng Campbell và Wynne (2011) đã hệ thống hóa nền tảng quang học Fourier, định nghĩa MTF là đáp ứng biên độ không gian của hệ thống đối với các kích thích hình sin ở các tần số khác nhau. Để ước tính MTF trên quỹ đạo, các phương pháp vật chuẩn cạnh (edge method hay knife-edge) được chuẩn hóa trong tiêu chuẩn quốc tế ISO 12233:2000 và bản cập nhật ISO 12233:2017. Scott (2003) tại DigitalGlobe đã chứng minh tính hiệu quả của phương pháp cạnh nghiêng trên 16 cảnh ảnh vệ tinh QuickBird, chỉ ra rằng việc tái chia mẫu nội suy (re-sampling) ảnh hưởng trực tiếp đến giá trị MTF tại tần số Nyquist ($f_N$). Viallefont-Robinet và Léger (2010) áp dụng thuật toán cạnh gờ trên bãi thử Salon de Provence để giám sát độ suy giảm thấu kính HRG trên SPOT-5. Storey (2001) và Barsi và cộng sự (2014) khai thác các vật thể dạng xung (pulse targets) như cây cầu vượt hồ Pontchartrain để kiểm tra đáp ứng đường (Line Spread Function - LSF) cho Landsat-7 ETM+.

graph LR
    subgraph "Trường phái Lý thuyết & Tiền phóng"
        A["Mô hình hóa tham số quang sai<br/>(Boreman, 2001)"] --> C["Giới hạn: Không phản ánh<br/>suy thoái thực tế trên quỹ đạo"]
    end
    
    subgraph "Trường phái Thực nghiệm Quốc tế"
        B1["Hiệu chuẩn chéo PICS<br/>(Chander et al., 2013)"]
        B2["Phương pháp cạnh Heaviside<br/>(Scott, 2003 - QuickBird)"]
        B3["Tiêu chuẩn ISO 12233:2017<br/>(Góc nghiêng hạn chế < 10°)"]
    end
    
    subgraph "Định vị Luận án (Nguyễn Minh Ngọc, 2022)"
        D["Tích hợp Hiệu chỉnh DS/PRNU thô<br/>+ Thuật toán Canny phi tuyến<br/>+ Bãi thử thực địa Buôn Ma Thuột"]
    end

    C -.-> D
    B1 & B2 & B3 --> D

Tuy nhiên, các tranh luận học thuật sâu sắc vẫn tồn tại xung quanh phương pháp luận:

  1. Tranh cãi về mô hình tham số trực tiếp đối lập với phương pháp gián tiếp vicarious: Các nhà sản xuất vệ tinh thường ưu tiên mô hình quang học nội tại kết hợp cảm biến nhiệt độ trên bo mạch. Tuy nhiên, các nhà nghiên cứu độc lập (Helder et al., 2004; Wu et al., 2019) chứng minh rằng các mô hình lý thuyết không thể phản ánh chính xác hiện tượng biến dạng dẻo nhiệt cấu trúc vệ tinh và sự rung lắc vi mô (attitude jitter) trong môi trường bay thực tế, do đó phương pháp đánh giá gián tiếp từ dữ liệu ảnh thực địa là thước đo chân thực nhất.
  2. Tranh luận về giải thuật chiết tách cạnh: Thuật toán cạnh gờ tuyến tính cổ điển nhạy cảm cao với nhiễu đốm và sự gồ ghề của vật liệu sơn bãi chuẩn, dẫn đến sai số trong việc khớp hàm phân bố bức xạ Heaviside. Các phương pháp biến đổi mờ hay độ phân giải kép (Latry et al., 2004) đòi hỏi các cặp ảnh chụp đồng thời với tỉ số phân giải bằng 3, vốn bất khả thi đối với các hệ thống vệ tinh đơn lẻ.

Luận án của NCS. Nguyễn Minh Ngọc đã định vị nghiên cứu chính xác tại giao điểm này: kế thừa khung lý thuyết vicarious calibration quốc tế, nhưng tạo bước đột phá phương pháp luận bằng cách kết hợp thuật toán lọc biên phi tuyến Canny với chuỗi dữ liệu đa thời gian tại bãi thử Buôn Ma Thuột, giải quyết trọn vẹn bài toán tự chủ công nghệ đánh giá chất lượng vệ tinh quang học cho Việt Nam.

Đó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 Hệ thống Dịch chuyển Tuyến tính Bất biến (LSI) trong quang học Fourier ứng dụng cho hệ thống cảm biến đẩy (pushbroom sensor) trên quỹ đạo thấp (Low Earth Orbit - LEO). Cốt lõi của lý thuyết này xem toàn bộ chuỗi thu nhận ảnh—từ khí quyển, hệ quang học viễn vọng Korsch, bộ chuyển đổi quang điện TDI-CCD, mạch điện tử tiền khuếch đại đến quá trình nén và lượng tử hóa—như một bộ lọc thông thấp tuyến tính nhiều tầng.

Hàm lan truyền điểm tổng thể của hệ thống $PSF_{sys}(x,y)$ là tích chập của các PSF thành phần: $$PSF_{sys}(x,y) = PSF_{optics}(x,y) * PSF_{detector}(x,y) * PSF_{motion}(x,y) * PSF_{atm}(x,y)$$

Tương ứng trong miền tần số không gian $(f_x, f_y)$, hàm truyền quang học $OTF(f_x, f_y)$ được biểu diễn thông qua biến đổi Fourier hai chiều: $$OTF(f_x, f_y) = \mathcal{F}{PSF_{sys}(x,y)} = MTF(f_x, f_y) \cdot e^{i\Phi(f_x, f_y)}$$

Trong đó, Hàm truyền điều biến $MTF(f)$ là môđun độ lớn thể hiện sự suy giảm độ tương phản của mục tiêu dạng sóng hình sin khi truyền qua hệ thống: $$MTF(f) = |OTF(f)| = \frac{M_{image}(f)}{M_{object}(f)}$$ với $M = \frac{I_{max} - I_{min}}{I_{max} + I_{min}}$ là độ tương phản (modulation).

flowchart LR
    A["Vật chuẩn Cạnh Mặt Đất<br/>f(x)"] -->|"Lấy mẫu qua cảm biến"| B["Hàm lan truyền cạnh<br/>ESF(x)"]
    B -->|"Đạo hàm bậc nhất d/dx"| C["Hàm lan truyền đường<br/>LSF(x)"]
    C -->|"Biến đổi Fourier rời rạc"| D["Hàm truyền quang học<br/>OTF(f)"]
    D -->|"Lấy môđun biên độ"| E["Hàm truyền điều biến<br/>MTF(f)"]
    E -->|"Chuẩn hóa"| F["MTF(Nyquist) >= 0.10"]

Luận án thiết lập hai mệnh đề lý thuyết nền tảng:

  • Mệnh đề 1 (Proposition 1): Đối với hệ cảm biến TDI-CCD quét đẩy, sự suy giảm MTF theo phương dọc hướng bay (Along-track) phản ánh tổng hòa của quang sai thấu kính, kích thước phần tử cảm biến và độ trôi vận tốc góc đồng bộ (TDI synchronization error), trong khi MTF theo phương ngang hướng bay (Across-track) phản ánh độc lập đặc tính quang sai và khẩu độ cảm biến.
  • Mệnh đề 2 (Proposition 2): Nhiễu bức xạ trong ảnh thô là sự kết hợp trực giao giữa nhiễu ngẫu nhiên Poisson của dòng photon/dòng tối và nhiễu sọc hệ thống phát sinh từ hệ số hồi đáp không đồng đều của các phần tử cảm biến (PRNU), có thể mô hình hóa độc lập thông qua ma trận phương sai cục bộ hàng và cột.

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

Khung phân tích của luận án tích hợp liên hoàn ba miền kiến trúc lý thuyết:

  1. Lý thuyết Hiệu chuẩn Bức xạ Tương đối & Tuyệt đối Gián tiếp: Khử bỏ hoàn toàn độ dịch dòng tối (Dark Signal - DS) và chuẩn hóa phản ứng điểm ảnh (Pixel Response Non-Uniformity - PRNU) trên dữ liệu thô Level 0 bằng các cảnh ảnh bề mặt biển đêm và sa mạc đồng nhất Lambert.
  2. Kỹ thuật Tối ưu hóa Gradient Biên Phi tuyến (Canny-based Sub-pixel Edge Profiling): Khắc phục triệt để nhược điểm của toán tử sai phân bậc nhất cổ điển bằng việc tích hợp bộ lọc Gauss khử nhiễu, thuật toán triệt tiêu không cực đại (Non-Maximum Suppression) và ngưỡng trễ kép (Hysteresis Thresholding) để trích xuất chính xác vị trí và góc nghiêng của cạnh mục tiêu $\theta$ trên bãi kiểm định.
  3. Mô hình Khớp Đa thức Bậc cao dựng hàm ESF và LSF: Chiếu các điểm ảnh hai chiều lên trục trực giao với đường biên cạnh để xây dựng Hàm lan truyền cạnh liên tục (Edge Spread Function - ESF), thực hiện vi phân để thu được Hàm lan truyền đường (Line Spread Function - LSF) và tính biến đổi Fourier rời rạc chuẩn hóa: $$MTF(f) = \frac{|\mathcal{F}{LSF(x)}|}{|\mathcal{F}{LSF(x)}|_{f=0}}$$

Điều kiện biên xác định: Khung phân tích áp dụng nghiêm ngặt cho cảm biến quang học vùng phổ nhìn thấy và cận hồng ngoại (VNIR: 0,45 – 0,89 µm), góc chụp nghiêng không quá 10°, độ che phủ mây 0% trên khu vực ô mẫu và bãi kiểm định có độ tương phản phản xạ bề mặt $\Delta\rho \ge 30%$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng luận (Positivism), dựa trên các phép đo lường vật lý khách quan, khả năng tái lập quy trình toán học và kiểm định thống kê chuẩn mực. Thiết kế nghiên cứu được tổ chức theo mô hình đa tầng (multi-level experimental design):

  • Tầng vật lý cảm biến: Phân tích dữ liệu mức thô (Level 0) để xác định các đặc tính dòng tối và sai số hồi đáp bức xạ.
  • Tầng tham số quang - hình học: Xử lý dữ liệu mức 1A để đo đạc và tính toán chính xác MTF và SNR tại các tần số không gian quy chuẩn.
  • Tầng ứng dụng kỹ thuật: Đối chiếu các thông số vật lý với tiêu chuẩn sai số bản đồ địa hình quốc gia nhằm thiết lập ngưỡng chất lượng ứng dụng.
graph TD
    subgraph "Tầng 1: Vật lý cảm biến (Level 0)"
        M0["Ảnh thô VNREDSat-1"]
        DS["Đo dòng tối DS ban đêm<br/>(Đại Tây Dương)"]
        PRNU["Tính hệ số PRNU<br/>(Sa mạc Algeria/Libya)"]
        M0 --> DS & PRNU
    end

    subgraph "Tầng 2: Tham số hình học & Bức xạ (Level 1A)"
        M1["Ảnh hiệu chỉnh Level 1A"]
        DS & PRNU --> M1
        CANNY["Chiết tách biên Canny<br/>(Sub-pixel Edge Fitting)"]
        NOISE["Phân tích ma trận 50x50<br/>(Độ lệch chuẩn sigma)"]
        M1 --> CANNY & NOISE
        CANNY --> RES_MTF["MTF tại f_Nyquist<br/>(Dọc & Ngang hướng bay)"]
        NOISE --> RES_SNR["Tỉ lệ tín hiệu/nhiễu SNR<br/>(Kênh PAN & Kênh MS)"]
    end

    subgraph "Tầng 3: Ứng dụng Địa lý & Bản đồ"
        RES_MTF & RES_SNR --> CLASS["Mô hình phân hạng chất lượng ảnh"]
        CLASS --> MAP1["Mức 1: Bản đồ 1:5.000 / 1:10.000"]
        CLASS --> MAP2["Mức 2: Bản đồ 1:25.000"]
        CLASS --> MAP3["Mức 3: Bản đồ 1:50.000"]
        CLASS --> MAP4["Mức 4: Giám sát chuyên đề định tính"]
    end

Quy trình nghiên cứu rigorous

Quy trình đánh giá chất lượng ảnh được xây dựng qua 4 giai đoạn chuẩn hóa:

[Dữ liệu thô Mức 0] 
       │
       ▼
[GIAI ĐOẠN 1: HIỆU CHỈNH BỨC XẠ]
 ├─ Khử Tín hiệu tối (Dark Signal - DS) ──> DN_DS_corr(p) = DN(p) - DS(p)
 └─ Cân bằng phản ứng điểm ảnh (PRNU)  ──> DN_PRNU_corr(p) = g(p) * DN_DS_corr(p)
       │
       ▼
[Dữ liệu chuẩn Mức 1A]
       │
       ├─────────────────────────────────────────┐
       ▼                                         ▼
[GIAI ĐOẠN 2: TÍNH TOÁN SNR]             [GIAI ĐOẠN 3: TÍNH TOÁN MTF]
 ├─ Chọn ô mẫu đồng nhất bãi thử         ├─ Trích xuất góc nghiêng biên (Toán tử Canny)
 ├─ Tính giá trị trung bình m            ├─ Chiếu mẫu 1D tái tạo đường cong ESF
 ├─ Tách phương sai hàng & cột           ├─ Vi phân số xác định đáp ứng xung LSF
 └─ SNR = m / sqrt(σ_col² + σ_row²)      └─ Biến đổi Fourier chuẩn hóa tìm MTF(f_Nyquist)
       │                                         │
       └────────────────────┬────────────────────┘
                            │
                            ▼
           [GIAI ĐOẠN 4: PHÂN HẠNG ỨNG DỤNG BẢN ĐỒ]
  1. Hiệu chỉnh Tín hiệu tối (Dark Signal - DS): Trích xuất giá trị dòng trôi nhiệt của từng điểm ảnh $p$ trên toàn bộ 7.000 phần tử mảng cảm biến bằng việc chụp khu vực đại dương hoàn toàn không có ánh sáng vào ban đêm. $$DN_{corr}(p) = DN(p) - \overline{DS}(p)$$
  2. Hiệu chỉnh mức độ hồi đáp không đồng đều (PRNU): Thu nhận các cảnh ảnh trên bề mặt sa mạc đồng nhất tuyệt đối tại Algeria và Libya trong điều kiện góc chiếu cao. Hệ số khuếch đại tương đối $g(p)$ của phần tử cảm biến $p$ được chuẩn hóa qua giá trị trung bình toàn mảng: $$g(p) = \frac{\frac{1}{N}\sum_{i=1}^N \overline{DN}{homo}(i)}{\overline{DN}{homo}(p)}$$
  3. Quy trình trích xuất MTF qua thuật toán Canny:
    • Lọc phẳng nhiễu bằng nhân Gauss kích thước $5 \times 5$ với độ lệch chuẩn $\sigma_g = 1.4$.
    • Tính toán ma trận độ lớn Gradient và hướng góc cạnh $\theta = \arctan(G_y / G_x)$.
    • Triệt tiêu các điểm không phải cực đại cục bộ dọc theo hướng Gradient để định vị chính xác tâm biên độ ở mức dưới điểm ảnh.
    • Tái chia mẫu các điểm lân cận với khoảng cách $\Delta x = 0.1 \text{ pixel}$ vuông góc với đường hồi quy biên cạnh để xây dựng đường cong ESF.
    • Tính đạo hàm bậc nhất bằng sai phân trung tâm để thu được LSF, áp dụng cửa sổ Hanning giảm rò rỉ phổ trước khi tính biến đổi Fourier nhanh (FFT) để xác định giá trị MTF tại tần số Nyquist ($f_N = 0.5 \text{ cycle/pixel}$).
  4. Quy trình tính toán SNR: Chọn các ô mẫu đồng nhất kích thước $50 \times 50$ điểm ảnh tại bãi kiểm định Buôn Ma Thuột. Tính toán giá trị tín hiệu trung bình $\overline{L}$ và độ lệch chuẩn nhiễu tổng hợp $\sigma = \sqrt{\sigma_{col}^2 + \sigma_{row}^2}$ để loại trừ tác động riêng rẽ của nhiễu thiết bị và nhiễu chuẩn hóa: $$SNR = \frac{\overline{L}}{\sigma}$$

Data và phân tích

Dữ liệu thực nghiệm của luận án bao gồm:

  • Chuỗi dữ liệu ảnh viễn thám VNREDSat-1: Gồm 18 cảnh ảnh thô Level 0 chụp sa mạc Algeria, Libya, biển Đại Tây Dương; 24 cảnh ảnh Level 1A chụp bãi kiểm định Buôn Ma Thuột và Salon de Provence; cùng các cảnh ảnh chụp các khu vực địa hình đặc trưng tại Việt Nam (Đức Thọ - Hà Tĩnh, Nha Trang - Khánh Hòa, Mường Lát - Thanh Hóa, Bình Đại - Bến Tre) từ năm 2013 đến năm 2021.
  • Dữ liệu thực địa quang phổ kế: Đo đạc phản xạ phổ mặt đất tại 4 ô mẫu bê tông đen, xám, trắng tại Buôn Ma Thuột bằng máy đo quang phổ phân giải cao trường nhìn hẹp trong hai đợt đo tháng 11/2017 và tháng 11/2018.
  • Công cụ phần mềm: Nghiên cứu sinh tự phát triển bộ công cụ chuyên dụng trên môi trường lập trình bậc cao tích hợp xử lý ảnh ma trận, cho phép tự động đọc tệp thông số hiệu chỉnh Calibration Parameter File (CPF), trích xuất vùng quan tâm (ROI), khớp toán học biên Canny và xuất đồ thị MTF/SNR.

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

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

Sau quá trình thực nghiệm nghiêm ngặt và xử lý thống kê, luận án đã mang lại 4 phát hiện mang tính bước ngoặt:

pie title Phân bố suy giảm MTF theo thành phần hệ thống VNREDSat-1
    "Quang sai thấu kính Korsch" : 35
    "Kích thước phần tử TDI-CCD" : 30
    "Nhiễu trôi vận tốc góc & Rung lắc" : 20
    "Tán xạ khí quyển & Khác" : 15
  1. Hiệu năng MTF của hệ quang học VNREDSat-1 duy trì xuất sắc sau 8 năm bay: Kết quả đo đạc thực nghiệm trên bãi thử Buôn Ma Thuột và Salon de Provence cho thấy giá trị MTF tại tần số Nyquist ($f_N$) của kênh toàn sắc PAN đạt giá trị trung bình $MTF_{Nyquist} = 0.142 \pm 0.015$ (theo phương dọc hướng bay - Along-track) và $0.168 \pm 0.012$ (theo phương ngang hướng bay - Across-track). Các giá trị này đều vượt ngưỡng tiêu chuẩn thiết kế ban đầu của nhà sản xuất Astrium ($MTF \ge 0.10$). Đối với 4 kênh đa phổ (B1-Blue, B2-Green, B3-Red, B4-NIR), giá trị MTF Nyquist duy trì ổn định trong khoảng $0.22 - 0.28$. Điều này khẳng định hệ thống thấu kính viễn vọng Korsch không bị suy thoái quang học nghiêm trọng do môi trường vũ trụ.
  2. Tỉ lệ SNR thực tế vượt ngưỡng thiết kế ở mọi kênh phổ: Dữ liệu thực nghiệm tại các ô mẫu chuẩn Buôn Ma Thuột ghi nhận tỉ lệ SNR của kênh PAN đạt $112.4$ (ở mức chiếu sáng tiêu chuẩn), cao hơn đáng kể so với mức yêu cầu thiết kế danh định ($SNR_{req} \ge 90$). Trên các kênh đa phổ, SNR đạt từ $165.8$ đến $214.5$ (vượt xa ngưỡng yêu cầu $\ge 150$).
  3. Phát hiện quy luật dị hướng không gian của MTF do chuyển động quét: Luận án chỉ ra sự chênh lệch có hệ thống giữa MTF dọc hướng bay và ngang hướng bay ($\Delta MTF \approx 0.026$). Phân tích vật lý chứng minh nguyên nhân xuất phát từ sự tích hợp thời gian trễ (TDI) của cảm biến mảng kết hợp với độ dịch chuyển góc cực nhỏ (pitch/roll micro-vibrations) của thân vệ tinh trong quá trình phơi sáng khi di chuyển với vận tốc quỹ đạo xấp xỉ $7.5 \text{ km/s}$.
  4. Chứng minh tính tương đương chuẩn mực của bãi kiểm định Buôn Ma Thuột: Kết quả so sánh chéo giá trị phản xạ phổ và tính toán MTF giữa bãi thử Buôn Ma Thuột (Việt Nam) và bãi thử chuẩn quốc tế Salon de Provence (Pháp) cho thấy độ lệch sai số tương đối $< 3.8%$, khẳng định bãi thử Buôn Ma Thuột hoàn toàn đáp ứng các tiêu chí kỹ thuật khắt khe của một bãi kiểm định mặt đất có trang thiết bị (Land Equipped Site - LES) đạt chuẩn quốc tế.
Kênh phổ VNREDSat-1 Bước sóng ($\mu$m) Độ phân giải (GSD) SNR Thiết kế SNR Thực nghiệm (Buôn Ma Thuột) MTF Nyquist Thiết kế MTF Nyquist Thực nghiệm (Trung bình)
PAN (Toàn sắc) 0.45 – 0.75 2.5 m $\ge 90$ 112.4 $\ge 0.10$ 0.155
B1 (Blue) 0.45 – 0.52 10.0 m $\ge 150$ 165.8 $\ge 0.10$ 0.245
B2 (Green) 0.53 – 0.59 10.0 m $\ge 150$ 188.2 $\ge 0.10$ 0.262
B3 (Red) 0.62 – 0.69 10.0 m $\ge 150$ 174.6 $\ge 0.10$ 0.258
B4 (NIR) 0.77 – 0.89 10.0 m $\ge 150$ 214.5 $\ge 0.10$ 0.278

Implications đa chiều

  • Ý nghĩa lý thuyết: Nghiên cứu đã làm phong phú lý thuyết đánh giá cảm biến viễn thám trong điều kiện khí hậu nhiệt đới gió mùa ẩm đặc trưng của Việt Nam—nơi có độ ẩm không khí và sol khí (aerosol) cao thường xuyên làm suy giảm độ tương phản quang học.
  • Đổi mới phương pháp luận: Việc thay thế thuật toán sai phân tuyến tính bằng thuật toán Canny tích hợp nội suy dưới điểm ảnh đã cung cấp một công cụ toán học mạnh mẽ, có thể chuyển giao áp dụng cho bất kỳ hệ thống vệ tinh quang học nào khác (như SPOT-6/7, Pleiades, WorldView, KOMPSAT-3).
  • Ứng dụng thực tiễn và chính sách: Luận án đã cụ thể hóa việc chuyển đổi các giá trị MTF và SNR thành bộ tiêu chí phân hạng chất lượng ảnh phục vụ các tỷ lệ bản đồ địa hình quốc gia:
    • Chất lượng Mức 1 ($MTF \ge 0.12; SNR \ge 90$): Đủ điều kiện thành lập và cập nhật cơ sở dữ liệu nền địa lý, bản đồ địa hình tỷ lệ 1:5.000 và 1:10.000.
    • Chất lượng Mức 2 ($0.08 \le MTF < 0.12; 60 \le SNR < 90$): Phục vụ thành lập bản đồ địa hình tỷ lệ 1:25.000 và 1:50.000.
    • Chất lượng Mức 3 ($MTF < 0.08; SNR < 60$): Chỉ sử dụng cho mục đích nhận biết chuyên đề định tính, giám sát biến động lớp phủ bề mặt không đòi hỏi độ chính xác biên giới thửa đất.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn khoa học mang tính khách quan:

  1. Hạn chế về mô hình phản xạ góc BRDF: Nghiên cứu chưa đo đạc liên tục Hàm phân bố phản xạ hai chiều (Bidirectional Reflectance Distribution Function - BRDF) đa góc của bề mặt bê tông bãi thử Buôn Ma Thuột theo mọi góc tà mặt trời trong năm, dẫn đến sai số nhỏ khi vệ tinh chụp với góc nghiêng lớn.
  2. Hạn chế về mô hình hiệu chỉnh khí quyển tức thời: Các phép tính MTF chủ yếu thực hiện trong các ngày thời tiết tối ưu (quang mây, tầm nhìn xa $> 10 \text{ km}$), chưa tích hợp hoàn toàn mô hình truyền xạ bức xạ khí quyển thời gian thực (như 6S hoặc MODTRAN) gắn kèm trạm đo thông số sol khí tự động (AERONET).
  3. Phạm vi loại cảm biến: Luận án mới tập trung vào cảm biến viễn thám quang học dải sóng nhìn thấy và cận hồng ngoại (VNIR), chưa mở rộng sang cảm biến hồng ngoại nhiệt (TIR) hay viễn thám siêu phổ (Hyperspectral) và Ra-đa khẩu độ tổng hợp (SAR).
  4. Kích thước vật chuẩn nhân tạo: Kích thước các ô mẫu tại bãi thử Buôn Ma Thuột hiện hữu còn hạn chế đối với các vệ tinh có độ phân giải không gian trung bình và thấp ($GSD > 10 \text{ m}$).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển thuật toán tự động hóa trích xuất MTF từ các đối tượng nhân tạo tự nhiên (cạnh mái tòa nhà, đường băng sân bay, ranh giới đê kè) kết hợp mạng thần kinh tích chập (Deep CNN) để giám sát liên tục chất lượng ảnh mà không phụ thuộc vào bãi kiểm định chuyên dụng.
  • Nghiên cứu thiết kế mạng lưới các mục tiêu hiệu chuẩn di động (sơn phản xạ đặc biệt dạng nan quạt, thang độ xám phản xạ chuẩn) phục vụ kiểm định các chùm vệ tinh siêu nhỏ (CubeSat/MicroSat) trong tương lai.
  • Xây dựng mô hình số đánh giá suy giảm MTF theo thời gian thực (Real-time MTF tracking model) phục vụ tự động bù trừ độ sắc nét trong phần mềm tiền xử lý trạm mặt đất.

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

  • Tác động học thuật (Academic Impact): Luận án là tài liệu tham khảo chuẩn mực, cung cấp cơ sở phương pháp luận cho các nghiên cứu tiếp theo về kỹ thuật viễn thám, đo ảnh và xử lý tín hiệu không gian tại các viện nghiên cứu và trường đại học. Dự báo công trình sẽ nhận được sự trích dẫn rộng rãi trong các nghiên cứu về hiệu chuẩn vệ tinh tại khu vực Đông Nam Á.
  • Chuyển đổi công nghiệp và vận hành vệ tinh: Kết quả nghiên cứu đã được chuyển giao trực tiếp cho Trung tâm Điều khiển và Khai thác Vệ tinh nhỏ (Viện Công nghệ Vũ trụ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam), giúp đội ngũ kỹ sư làm chủ hoàn toàn công tác giám sát sức khỏe cảm biến VNREDSat-1 mà không cần phụ thuộc vào chuyên gia nước ngoài, tiết kiệm hàng trăm nghìn USD chi phí bảo trì chuyên gia hàng năm.
  • Tác động chính sách quốc gia: Cung cấp cơ sở khoa học định lượng chuẩn xác giúp Bộ Tài nguyên và Môi trường cập nhật, hoàn thiện các Thông tư quy định kỹ thuật và định mức kinh tế kỹ thuật trong sản xuất, xử lý và nghiệm thu dữ liệu ảnh viễn thám quang học độ phân giải cao và siêu cao.

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

mindmap
  root((Hệ sinh thái hưởng lợi))
    Nghiên cứu sinh & Giới học thuật
      Kế thừa phương pháp luận vicarious calibration
      Khai thác mã nguồn thuật toán Canny MTF
      Mở rộng nghiên cứu cho vệ tinh thế hệ mới
    Cơ quan Vận hành Vệ tinh
      Viện Công nghệ Vũ trụ VAST
      Trung tâm Điều khiển VNREDSat-1
      Giám sát suy thoái cảm biến định kỳ
    Doanh nghiệp & Đơn vị Đo đạc Bản đồ
      Cục Đo đạc, Bản đồ và Thông tin địa lý
      Các công ty khảo sát trắc địa
      Phân loại chính xác ảnh cho từng tỷ lệ bản đồ
    Cơ quan Quản lý & Hoạch định Chính sách
      Bộ Tài nguyên và Môi trường
      Cục Viễn thám Quốc gia
      Xây dựng Tiêu chuẩn Quốc gia TCVN
  1. Nghiên cứu sinh và Giới học thuật Viễn thám: Tiếp cận một khung phương pháp luận toàn diện từ tiền xử lý mức 0 đến đánh giá tham số vật lý mức 1A, mở ra các hướng nghiên cứu xử lý ảnh vệ tinh chuyên sâu.
  2. Cơ quan quản lý và khai thác hệ thống vệ tinh (VAST, Cục Viễn thám Quốc gia): Nắm giữ quy trình kiểm chuẩn độc lập, chủ động kiểm tra chất lượng các nguồn ảnh vệ tinh mua từ nước ngoài (như SPOT-6/7, Pleiades, WorldView, KOMPSAT-3) trước khi nhập kho dữ liệu quốc gia.
  3. Đơn vị sản xuất bản đồ và cơ sở dữ liệu không gian: Đảm bảo độ tin cậy tuyệt đối về chất lượng hình học và độ sắc nét quang phổ khi sử dụng ảnh vệ tinh để thành lập bản đồ địa chính, bản đồ địa hình và bản đồ hiện trạng rừng.
  4. Nhà hoạch định chính sách: Có căn cứ khoa học vững chắc để xây dựng các tiêu chuẩn quốc gia (TCVN) và quy chuẩn kỹ thuật quốc gia (QCVN) trong lĩnh vực công nghệ vũ trụ và viễn thám.

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ì?

Luận án đã hình thành mô hình định lượng hóa chuỗi suy thoái cảm biến quang điện TDI-CCD trên quỹ đạo thấp trong điều kiện khí quyển nhiệt đới bằng cách liên kết hệ thống các lý thuyết quang học Fourier (mô hình LSI), lý thuyết phân tích thống kê nhiễu trực giao (phương sai hàng - cột) và lý thuyết xử lý biên phi tuyến. Đóng góp độc đáo nhất là chứng minh sự liên kết giữa chỉ số suy giảm MTF vật lý tại tần số Nyquist với khả năng phân giải hình học thực tế phục vụ các cấp tỷ lệ bản đồ địa hình.

2. Đổi mới phương pháp luận nổi bật khi so sánh với các nghiên cứu quốc tế?

So với phương pháp cạnh nghiêng truyền thống của ISO 12233:2000 (vốn bị giới hạn nghiêm ngặt ở góc nghiêng $< 10^\circ$ và nhạy cảm với nhiễu) hay phương pháp lọc tuyến tính trên SPOT-5 (Viallefont-Robinet, 2010), luận án đã tích hợp thuật toán lọc biên tối ưu Canny với cơ chế làm mịn Gauss và ngưỡng trễ kép. Giải pháp này cho phép định vị tim cạnh mục tiêu với độ chính xác dưới điểm ảnh ($\pm 0.1 \text{ pixel}$), triệt tiêu hoàn toàn hiện tượng răng cưa (aliasing) do lấy mẫu dưới ngưỡng và không bị ảnh hưởng bởi độ gồ ghề của bề mặt bãi thử bê tông.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sau 8 năm hoạt động trên quỹ đạo (vượt quá tuổi thọ thiết kế danh định 5 năm), hệ thống quang học VNREDSat-1 vẫn duy trì chất lượng vượt trội: MTF kênh PAN đạt $0.155$ (cao hơn $55%$ so với ngưỡng thiết kế tối thiểu $\ge 0.10$), và SNR đạt $112.4$ (vượt ngưỡng yêu cầu $\ge 90$). Đồng thời, nghiên cứu phát hiện sự bất đối xứng MTF có quy luật giữa phương dọc và ngang hướng bay ($\Delta \approx 0.026$) do hiệu ứng tích hợp thời gian trễ TDI.

4. Quy trình nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không?

Luận án cung cấp một giao thức thực nghiệm hoàn chỉnh, chi tiết từng bước: từ việc lựa chọn tọa độ cảnh ảnh kiểm định PICS trên thế giới, thông số hiệu chỉnh dòng tối DS và hệ số PRNU trên dữ liệu mức 0, tọa độ góc các ô mẫu tại bãi thử Buôn Ma Thuột, đến các phương trình toán học trích xuất cạnh và mã nguồn công cụ tính toán chi tiết trong phần phụ lục. Quy trình này có thể tái lập hoàn toàn trên bất kỳ trạm thu ảnh viễn thám nào.

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

Chương trình nghiên cứu 10 năm tới tập trung vào: (1) Tự động hóa kiểm định MTF/SNR trên quỹ đạo sử dụng trí tuệ nhân tạo và các mục tiêu tự nhiên ngẫu nhiên; (2) Thiết lập mạng lưới bãi chuẩn đa phổ và viễn thám siêu phổ quốc gia; (3) Nghiên cứu cơ chế suy thoái và giải thuật hiệu chuẩn cho chùm vệ tinh quan sát Trái đất thế hệ mới của Việt Nam (VNREDSat-2 và chùm vệ tinh nhỏ).

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Minh Ngọc là một công trình nghiên cứu khoa học công phu, nghiêm túc, giải quyết trọn vẹn bài toán thực tiễn cấp bách và đóng góp những giá trị học thuật xuất sắc cho ngành Trắc địa - Bản đồ và Viễn thám Việt Nam.

Tóm lược 5 đóng góp đột phá cốt lõi của công trình:

  1. Xác lập cơ sở khoa học và bộ thông số chuẩn hóa: Chứng minh luận điểm MTF và SNR là hai thông số đại diện tối ưu, đầy đủ nhất để đánh giá chất lượng bức xạ và không gian của ảnh viễn thám quang học tại Việt Nam.
  2. Đổi mới giải thuật xử lý ảnh kiểm định: Đề xuất thành công việc ứng dụng thuật toán chiết tách cạnh Canny trong tính toán MTF từ bãi kiểm định, nâng cao độ chính xác nội suy dưới điểm ảnh và triệt tiêu nhiễu bề mặt.
  3. Xây dựng quy trình công nghệ hoàn chỉnh: Thiết lập quy trình đánh giá khép kín từ khâu hiệu chỉnh dòng tối DS, cân bằng phản ứng điểm ảnh PRNU trên ảnh mức 0 đến tính toán MTF/SNR trên ảnh mức 1A.
  4. Kiểm định thực nghiệm toàn diện hệ thống vệ tinh VNREDSat-1: Cung cấp bộ dữ liệu lịch sử vô giá chứng minh hệ thống chụp ảnh VNREDSat-1 vẫn hoạt động hoàn hảo, vượt các chỉ tiêu thiết kế sau 8 năm bay trên quỹ đạo.
  5. Gắn kết kiểm định vật lý với tiêu chuẩn bản đồ quốc gia: Đề xuất bảng phân hạng chất lượng ảnh viễn thám thành 4 mức phục vụ trực tiếp công tác thành lập bản đồ địa hình từ tỷ lệ 1:5.000 đến 1:50.000, mở ra bước tiến tự chủ công nghệ đánh giá vệ tinh viễn thám cho Việt Nam.