Tổng quan về luận án

Bối cảnh khoa học của công nghệ vũ trụ và viễn thám quan sát Trái đất trong hai thập kỷ qua ghi nhận bước chuyển dịch mang tính cấu trúc: từ việc phụ thuộc hoàn toàn vào các hệ thống vệ tinh của các cường quốc vũ trụ sang xu hướng các quốc gia đang phát triển tự chủ thiết kế, phóng và vận hành chùm vệ tinh quan sát riêng biệt. Tại Việt Nam, cột mốc lịch sử diễn ra vào ngày 07/05/2013 khi vệ tinh quan sát Trái đất đầu tiên của Việt Nam – VNREDSat-1 (Vietnam Natural Resources, Environment and Disaster Monitoring Satellite) – được phóng thành công vào quỹ đạo đồng bộ mặt trời ở độ cao danh định 680 km. Tuy nhiên, sau khi hoàn thành giai đoạn thử nghiệm đầu quỹ đạo và tiếp nhận chuyển giao công nghệ từ nhà sản xuất (EADS Astrium, Cộng hòa Pháp), việc duy trì, giám sát và chuẩn hóa chất lượng dữ liệu ảnh trong suốt vòng đời hoạt động thực tế đặt ra một thách thức khoa học và công nghệ to lớn. Môi trường không gian khắc nghiệt (sự biến thiên nhiệt độ cực đoan từ $-40^\circ\text{C}$ đến $+60^\circ\text{C}$, bức xạ ion hóa, hiện tượng thoát khí trong chân không, rung chấn cơ học trong quá trình phóng và sự lão hóa của các linh kiện quang - điện tử) tất yếu dẫn đến sự suy giảm độ nhạy bức xạ và độ sắc nét không gian của hệ thống camera quang học.

+-----------------------------------------------------------------------------------+
|               KHUNG HỆ THỐNG ĐÁNH GIÁ CHẤT LƯỢNG ẢNH VNREDSAT-1                   |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [DỮ LIỆU ĐẦU VÀO MỨC 0 / 1A]                                                     |
|        │                                                                          |
|        ├─► Dải phổ Toàn sắc (PAN): Độ phân giải không gian GSD = 2.5 m            |
|        └─► Dải phổ Đa phổ (MS: Blue, Green, Red, NIR): GSD = 10.0 m               |
|                                                                                   |
|  [HIỆU CHỈNH BỨC XẠ GỐC]                                                          |
|        │                                                                          |
|        ├─► Tín hiệu tối (Dark Signal - DS) ────► Đại Tây Dương / Cảnh đêm         |
|        └─► Hồi đáp không đồng đều (PRNU) ─────► Bãi PICS Sa mạc Algeria / Libya   |
|                                                                                   |
|  [ĐÁNH GIÁ CHỈ SỐ VẬT LÝ CỐT LÕI]                                                 |
|        │                                                                          |
|        ├─► Yếu tố Bức xạ: Tỉ lệ Tín hiệu trên Nhiễu (SNR) ──► SNR = m / σ         |
|        └─► Yếu tố Không gian: Hàm truyền điều biến (MTF)   ──► MTF tại f_Nyquist   |
|                                                                                   |
|  [ĐỘT PHÁ PHƯƠNG PHÁP LUẬN]                                                       |
|        └─► Chiết tách Cạnh Nghiêng Sub-pixel bằng Toán tử Vi phân CANNY           |
|            (Nâng cao độ chính xác ESF -> LSF -> FFT sang miền tần số MTF)         |
|                                                                                   |
|  [CHUẨN HÓA THỰC NGHIỆM VÀ PHÂN HẠNG ỨNG DỤNG]                                    |
|        │                                                                          |
|        ├─► Bãi kiểm định nhân tạo: Buôn Ma Thuột (Đắk Lắk) & Salon de Provence    |
|        └─► Ma trận phân hạng chất lượng phục vụ đo vẽ bản đồ (1:5.000 - 1:50.000) |
|                                                                                   |
+-----------------------------------------------------------------------------------+

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Minh Ngọc, với tiêu đề "Nghiên cứu đánh giá chất lượng ảnh viễn thám quang học của Việt Nam" thuộc ngành Kỹ thuật Trắc địa – Bản đồ (Mã số: 9520503), 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 (2022), đã giải quyết căn cơ khoảng trống nghiên cứu (research gap) trọng yếu: Thiếu vắng một cơ sở phương pháp luận định lượng và quy trình công nghệ hoàn chỉnh, độc lập để đánh giá chất lượng quang học của vệ tinh viễn thám sau khi phóng (in-orbit vicarious calibration/validation) phù hợp với hạ tầng kỹ thuật và điều kiện khí quyển nhiệt đới gió mùa của Việt Nam.

Trước nghiên cứu này, công tác giám định ảnh vệ tinh tại Việt Nam chủ yếu dừng lại ở việc tiếp nhận các thông số danh định do nhà sản xuất cung cấp hoặc kiểm tra định tính ở các mức xử lý thương phẩm thứ cấp (Mức 2A, 3A, 3B theo Thông tư 10/2015/TT-BTNMT và Thông tư 52/2016/TT-BTNMT), hoàn toàn chưa làm chủ phương pháp luận tính toán các tham số vật lý cốt lõi từ dữ liệu thô (Mức 0 và Mức 1A).

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học mang tính hệ thống:

  • Câu hỏi nghiên cứu 1 (RQ1): Những thông số vật lý nào đại diện bản chất nhất cho chất lượng bức xạ và chất lượng không gian của hệ thống thu nhận ảnh quang học vệ tinh mà Việt Nam có thể chủ động đánh giá trong điều kiện thực tiễn?

  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để loại bỏ nhiễu và nâng cao độ chính xác của hàm truyền điều biến (Modulation Transfer Function - MTF) khi trích xuất từ các bãi kiểm định nhân tạo trên mặt đất dưới tác động của hiện tượng răng cưa (aliasing) và tán xạ khí quyển?

  • Câu hỏi nghiên cứu 3 (RQ3): Trạng thái suy giảm thực tế của hệ thống cảm biến quang học VNREDSat-1 sau hơn 8 năm hoạt động trên quỹ đạo diễn biến như thế nào so với tiêu chuẩn thiết kế ban đầu và đáp ứng các cấp độ bản đồ địa hình nào?

  • Giả thuyết khoa học 1 (H1): Việc kết hợp đồng thời hai chỉ số vật lý là Tỉ lệ tín hiệu trên nhiễu (Signal-to-Noise Ratio - SNR) đại diện cho bức xạ và Hàm truyền điều biến (MTF) đại diện cho không gian, sau khi đã triệt tiêu sai số tín hiệu tối (Dark Signal - DS) và mức độ hồi đáp không đồng đều của điểm ảnh (Pixel Response Non-Uniformity - PRNU), là điều kiện cần và đủ để định lượng toàn diện chất lượng quang học của vệ tinh viễn thám.

  • Giả thuyết khoa học 2 (H2): Ứng dụng toán tử phi tuyến chiết tách cạnh Canny kết hợp làm mịn đa ngưỡng cho phép khắc phục hạn chế của các thuật toán sống dao truyền thống, tối ưu hóa việc xác định Hàm lan truyền cạnh (Edge Spread Function - ESF) và nâng cao độ tin cậy của giá trị MTF tại tần số Nyquist.

Khung lý thuyết của luận án xây dựng 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), Quang học Biến đổi Fourier (Fourier Optics), Quang sai và Nhiễu xạ quang học theo tiêu chuẩn Rayleigh, và Mô hình Bề mặt Phản xạ Đẳng hướng Lambertian. Về phạm vi thực nghiệm, luận án xử lý khối dữ liệu đồ sộ gồm hàng trăm cảnh ảnh mức 0 và mức 1A của vệ tinh VNREDSat-1 (dải toàn sắc PAN độ phân giải 2,5 m và 4 dải đa phổ MS độ phân giải 10 m gồm Blue, Green, Red, NIR) thu nhận liên tục trong hơn 8 năm (từ 2013 đến 2021). Phạm vi không gian trải dài qua bãi kiểm định nhân tạo Buôn Ma Thuột (tỉnh Đắk Lắk, vận hành từ 2017), bãi kiểm định quốc tế Salon de Provence (Cộng hòa Pháp), các bãi chuẩn tự nhiên giả bất biến (Pseudo-Invariant Calibration Sites - PICS) tại sa mạc Algeria, sa mạc Libya (Bắc Phi) và vùng biển sâu Đại Tây Dương.

Literature Review và Positioning

Lịch sử nghiên cứu đánh giá chất lượng ảnh viễn thám quang học trên thế giới hình thành và phát triển song song với các chương trình vũ trụ lớn của NASA (Hoa Kỳ), CNES/ESA (Châu Âu), JAXA (Nhật Bản) và CNSA (Trung Quốc). Các nghiên cứu nền tảng phân nhánh thành hai trường phái chính:

Trường phái thứ nhất tập trung vào đánh giá bức xạ và tỉ lệ tín hiệu trên nhiễu (SNR). Tác giả Curlander và McDonough cùng các công trình của Markham và Barker (chương trình Landsat OLI/TIR) đã thiết lập mô hình tính toán SNR dựa trên phương sai của các vùng đồng nhất bề mặt đất kết hợp phân tách nguồn nhiễu ngẫu nhiên (Poisson photon noise) và nhiễu cố định (fixed pattern noise/readout noise). Helder và cộng sự (2012) chuẩn hóa việc sử dụng các bãi PICS sa mạc như Libya 4 và Algeria 3 để theo dõi độ suy giảm bức xạ theo thời gian thực của cảm biến quang học vệ tinh với độ bất biến phản xạ bề mặt sai số dưới 2-3%.

Trường phái thứ hai tập trung vào đánh giá phân giải không gian và Hàm truyền điều biến (MTF). Theo Boreman (2001), MTF là hàm đo lường khả năng truyền tải độ tương phản từ mục tiêu quang học thực tế sang mặt phẳng tiêu diện của cảm biến dưới dạng hàm số của tần số không gian. Robert D. Fiete (1999) và Campbell & Wynne (2011) chứng minh rằng độ phân giải hình học (GSD - Ground Sampling Distance hay IFOV) chỉ thể hiện giới hạn lấy mẫu hình học, trong khi MTF tại tần số lấy mẫu giới hạn (tần số Nyquist $f_{Nyq} = 0.5 \text{ chu kỳ/pixel}$) mới quyết định độ sắc nét và lượng thông tin hữu ích truyền qua hệ thống.

                                BẢN ĐỒ VỊ TRÍ NGHIÊN CỨU TRONG HỆ TRI THỨC QUỐC TẾ
                                
       Hạ tầng Cal/Val Quốc tế                          Phương pháp luận & Mô hình Toán
  +───────────────────────────────────+             +────────────────────────────────────────+
  | CNES / Airbus: Salon de Provence  |             | ISO 12233: Slanted-edge Tuyến tính     |
  | DigitalGlobe: QuickBird In-orbit  |             | Curlander & Helder: Mô hình PICS Libya |
  | USGS / NASA: Landsat & MODIS      |             | Boreman (2001): Fourier Optics & MTF   |
  +─────────────────┬─────────────────+             +───────────────────┬────────────────────+
                    │                                                   │
                    └─────────────────────────┬─────────────────────────┘
                                              ▼
                         +─────────────────────────────────────────+
                         | KHOẢNG TRỐNG NGHIÊN CỨU & POSITIONING   |
                         | - Thiếu quy trình Cal/Val cho Vệ tinh   |
                         |   nhỏ tại nước đang phát triển.         |
                         | - Nhiễu khí quyển nhiệt đới xói mòn MTF |
                         | - Thuật toán cạnh tuyến tính dễ lỗi mẫu |
                         +────────────────────┬────────────────────+
                                              ▼
                         +─────────────────────────────────────────+
                         | ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN (2022)     |
                         | - Giải thuật Canny Sub-pixel chiết ESF  |
                         | - Khung tích hợp DS + PRNU + SNR + MTF  |
                         | - Dữ liệu thực nghiệm VNREDSat-1 (8 năm)|
                         | - Chuẩn hóa bãi kiểm định Buôn Ma Thuột |
                         +─────────────────────────────────────────+

Trong y văn quốc tế, tồn tại hai luồng tranh luận học thuật sâu sắc:

  1. Tranh luận về phương pháp xác định MTF trên quỹ đạo: Nhóm nghiên cứu tiêu chuẩn hóa quốc tế ISO (ISO 12233:2000, cập nhật ISO 12233:2017) và Scott (2003) tại DigitalGlobe (áp dụng cho vệ tinh QuickBird) kiên định với phương pháp cạnh nghiêng tuyến tính (slanted-edge method) sử dụng các tấm bia đen trắng nhân tạo. Ngược lại, nhóm tác giả Jin Li, Fei Xing, Ting Sun, Zheng You (2015) và Christophe Latry cùng cộng sự (CNES, 2005) cho rằng phương pháp cạnh tuyến tính truyền thống phụ thuộc nặng nề vào góc nghiêng chụp ảnh (yêu cầu góc $< 5^\circ - 10^\circ$) và rất nhạy cảm với nhiễu nền khi tái lấy mẫu, do đó họ đề xuất phương pháp mô hình hàm dịch chuyển điểm ảnh hoặc so sánh chéo đa phân giải (dual-resolution method).
  2. Tranh luận về việc lựa chọn bãi kiểm định: Nhóm nghiên cứu tại Thái Lan (2011) khi đánh giá vệ tinh THEOS tại bãi Salon de Provence (Pháp) và Storey (2001) khi đánh giá Landsat-7 tại cầu vượt hồ Pontchartrain (Louisiana, Hoa Kỳ) hay Destin Bridge (Florida) cho Hyperion (2003) đưa ra quan điểm trái chiều về việc sử dụng vật chuẩn nhân tạo chuyên dụng hay tận dụng các công trình nhân tạo sẵn có (cầu, đường băng, mái nhà cao tầng). Các công trình sẵn có dễ tiếp cận nhưng thiếu thông số đo quang phổ phản xạ bề mặt đồng thời, gây sai số không kiểm soát được trong hàm lan truyền điểm (Point Spread Function - PSF).

Vị thế nghiên cứu của luận án được xác lập rõ ràng: Kế thừa chuẩn mực quốc tế về phương pháp cạnh Heaviside nhưng tạo đột phá kỹ thuật bằng việc thay thế thuật toán khớp cạnh tuyến tính cổ điển bằng toán tử phát hiện biên Canny cải tiến đa ngưỡng. Luận án đặt nghiên cứu của Việt Nam vào thế cân bằng khi so sánh trực tiếp kết quả thực nghiệm với các công trình quốc tế điển hình: hệ thống vệ tinh SPOT-5 của Pháp tại Salon de Provence (Viallefont-Robinet & Léger, 2010), vệ tinh Gaofen-2 (GF-2) của Trung Quốc (Wu Zhaocong et al., 2019) và OrbView-3 của Hoa Kỳ (Kohm, 2004).

Đó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) và Quang học Fourier sang lĩnh vực đánh giá chất lượng ảnh của các vệ tinh quan sát Trái đất dạng nhỏ (Small Earth Observation Satellites) hoạt động dài hạn ngoài tuổi thọ thiết kế.

                                  CHUỖI CHUYỂN ĐỔI TOÁN HỌC KHUNG PHÂN TÍCH MTF
                                  
  +─────────────────────+        Toán tử Đạo hàm bậc 1        +─────────────────────+
  | HÀM LAN TRUYỀN CẠNH | ──────────────────────────────────► | HÀM LAN TRUYỀN ĐƯỜNG |
  |       ESF(x)        |         LSF(x) = d(ESF)/dx          |       LSF(x)        |
  +─────────────────────+                                     +──────────┬──────────+
                                                                         │
                                         Biến đổi Fourier 1 chiều (FFT)  │
                                         và Chuẩn hóa biên độ f=0        ▼
                                                              +─────────────────────+
                                                              | HÀM TRUYỀN ĐIỀU BIẾN |
                                                              |  MTF(f) = |F{LSF}|  |
                                                              +─────────────────────+

Về mặt lý thuyết xử lý tín hiệu không gian, mối quan hệ giữa Hàm lan truyền điểm $PSF(x,y)$, Hàm lan truyền cạnh $ESF(x)$, Hàm lan truyền đường $LSF(x)$ và Hàm truyền điều biến $MTF(f_x, f_y)$ được luận án hệ thống hóa qua chuỗi biến đổi giải tích:

$$\text{LSF}(x) = \frac{d}{dx} \text{ESF}(x)$$

$$\text{MTF}(f_x) = |\text{OTF}(f_x)| = \left| \int_{-\infty}^{+\infty} \text{LSF}(x) e^{-i 2\pi f_x x} dx \right| \Bigg/ \int_{-\infty}^{+\infty} \text{LSF}(x) dx$$

Các mệnh đề khoa học được xác lập:

  • Mệnh đề 1 (P1): Giá trị MTF tổng thể của hệ thống vệ tinh trên quỹ đạo là tích chập của các MTF thành phần: $$\text{MTF}{\text{hệ thống}} = \text{MTF}{\text{quang học}} \times \text{MTF}{\text{đầu thu CCD/CMOS}} \times \text{MTF}{\text{chuyển động bay}} \times \text{MTF}_{\text{khí quyển}}$$ Khi hoạt động sau 5 năm, sự suy giảm MTF hệ thống chủ yếu do sai lệch vị trí tiêu diện (thermal refocusing drift) và suy giảm hiệu ứng chuyển tích điện (Charge Transfer Efficiency - CTE) của cảm biến.
  • Mệnh đề 2 (P2): Sự mất cân bằng hồi đáp bức xạ điểm ảnh (PRNU) tạo ra nhiễu dạng sọc tần số thấp dọc theo hướng bay, làm biến dạng dạng sóng hình sin của tín hiệu quang học đầu vào và tạo ra hiện tượng dịch pha trong miền tần số Fourier.
  • Mệnh đề 3 (P3): Việc áp dụng toán tử lọc Gauss đa chiều kết hợp khử các điểm không cực đại (non-maximum suppression) trong thuật toán Canny bảo toàn độ dốc thực tế của cạnh Heaviside, ngăn chặn hiện tượng giả mạo dải tần cao (high-frequency artifact) trong phổ Fourier của LSF.

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 3 trụ cột lý thuyết:

  1. Lý thuyết lan truyền sóng điện từ và bức xạ vi sai ($\text{NE}\Delta L$);
  2. Lý thuyết lấy mẫu tín hiệu rời rạc Nyquist - Shannon;
  3. Lý thuyết gradient vi phân không gian số học.
+──────────────────────────────────────────────────────────────────────────────────────────────────+
|                     KHUNG PHÂN TÍCH TỔNG THỂ ĐÁNH GIÁ CHẤT LƯỢNG ẢNH VIỄN THÁM                    |
+──────────────────────────────────────────────────────────────────────────────────────────────────+
|                                                                                                  |
|  [TẦNG 1: DỮ LIỆU GỐC & HIỆU CHỈNH CẢM BIẾN (LEVEL 0 -> 1A)]                                     |
|  ├── Triệt tiêu dòng tối (Dark Signal - DS):                                                     |
|  │     DN_corr(i, j) = DN_raw(i, j) - DS(j)                                                      |
|  └── Cân bằng bức xạ không đồng đều (PRNU):                                                      |
|        Gain(j) = <DN_homo> / DN_homo(j)  ==>  DN_calib(i, j) = Gain(j) * DN_corr(i, j)           |
|                                                                                                  |
|  [TẦNG 2: PHÂN TÍCH BỨC XẠ & NHIỄU (SNR ANALYSIS)]                                               |
|  ├── Trích xuất vùng đồng nhất (Uniform RoI) tại bãi Buôn Ma Thuột / Sa mạc PICS                  |
|  ├── Tách biệt nhiễu cột (Nhiễu cảm biến: σ_col) và nhiễu hàng (Nhiễu chuẩn hóa: σ_row)          |
|  └── Tính toán tỉ lệ tín hiệu trên nhiễu: SNR = m / sqrt(σ_col^2 + σ_row^2)                      |
|                                                                                                  |
|  [TẦNG 3: PHÂN TÍCH PHÂN GIẢI KHÔNG GIAN (MTF ANALYSIS)]                                         |
|  ├── Nhận dạng ô mẫu cạnh nghiêng đen/trắng (Buôn Ma Thuột / Salon de Provence)                  |
|  ├── Trích xuất biên Canny Sub-pixel & Xác định góc nghiêng θ của cạnh                            |
|  ├── Chiếu điểm ảnh tạo Super-sampled 1D Edge Spread Function (ESF)                              |
|  ├── Vi phân số bậc một: LSF(x) = d(ESF)/dx                                                      |
|  └── Biến đổi FFT & Chuẩn hóa: MTF(f) = |FFT{LSF}| / FFT{LSF}(0)  ==> Trích xuất MTF(f_Nyquist)  |
|                                                                                                  |
|  [TẦNG 4: MA TRẬN ỨNG DỤNG THỰC TIỄN (APPLICATION MAPPING)]                                      |
|  └── Đối soát ngưỡng MTF(f_Nyq) & SNR với quy chuẩn đo vẽ bản đồ địa hình (1:5.000 - 1:50.000)   |
|                                                                                                  |
+──────────────────────────────────────────────────────────────────────────────────────────────────+

Điều kiện biên xác lập chặt chẽ: Góc nghiêng quan sát ngoài thiên đỉnh (off-nadir angle) $< 15^\circ$ để loại trừ biến dạng hình học phi tuyến; góc cao mặt trời (solar elevation angle) $> 35^\circ$; hệ số phản xạ bề mặt ô mẫu bãi kiểm định đạt độ tương phản tối thiểu $C = (\rho_{\text{trắng}} - \rho_{\text{đen}}) / (\rho_{\text{trắng}} + \rho_{\text{đen}}) \ge 0.5$; độ ẩm không khí và hàm lượng bụi sol khí (AOD) tại thời điểm chụp nằm trong giới hạn kiểm chuẩn phổ thực địa.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng khách quan (Empirical Positivism) kết hợp đối sánh thực nghiệm đa cấp độ (Multi-level Experimental Design). Nghiên cứu không can thiệp trực tiếp vào vệ tinh vật lý (do vệ tinh đã bay trên quỹ đạo) mà sử dụng phương pháp kiểm chứng gián tiếp thông qua chuỗi xử lý số học tín hiệu ảnh từ cấp độ thấp nhất (Level 0 - dữ liệu dòng bit thô từ máy thu mặt đất) đến cấp độ hiệu chỉnh hình học cơ bản (Level 1A - dữ liệu đã khử sai số cảm biến nội tại nhưng chưa nắn chỉnh trực giao).

                      THIẾT KẾ ĐỐI SÁNH THỰC NGHIỆM ĐA TẦNG VÀ BÃI CHUẨN
                      
   +──────────────────────────────────+               +──────────────────────────────────+
   |   BÃI KIỂM ĐỊNH TỰ NHIÊN (PICS)  |               |  BÃI KIỂM ĐỊNH NHÂN TẠO ĐO PHỔ   |
   | - Sa mạc Libya 4 & Algeria 3     |               | - Buôn Ma Thuột (Đắk Lắk, VN)    |
   | - Biển sâu Đại Tây Dương         |               | - Salon de Provence (Pháp)       |
   +────────────────┬─────────────────+               +────────────────┬─────────────────+
                    │                                                  │
                    ▼                                                  ▼
   +──────────────────────────────────+               +──────────────────────────────────+
   | ĐÁNH GIÁ THAM SỐ BỨC XẠ GỐC      |               | ĐÁNH GIÁ PHÂN GIẢI KHÔNG GIAN    |
   | - Tín hiệu tối: DS(j)            |               | - ESF / LSF qua Toán tử Canny    |
   | - Hồi đáp điểm ảnh: PRNU(j)      |               | - SNR trên ô mẫu đồng nhất thực  |
   | - Tỉ lệ tín hiệu/nhiễu: SNR      |               | - MTF tại Nyquist (dọc/ngang bay)|
   +────────────────┬─────────────────+               +────────────────┬─────────────────+
                    │                                                  │
                    └────────────────────────┬─────────────────────────┘
                                             ▼
                             +────────────────────────────────+
                             | KIỂM CHỨNG CHÉO & PHÂN HẠNG    |
                             | - So sánh tiêu chuẩn Astrium   |
                             | - So sánh 8 năm suy giảm       |
                             | - Đề xuất cấp độ đo vẽ bản đồ  |
                             +────────────────────────────────+

Dữ liệu thực nghiệm bao gồm:

  • Dữ liệu VNREDSat-1 Level 0 & Level 1A: Dải phổ Panchromatic (bước sóng $0.45 - 0.75,\mu\text{m}$, kích thước ma trận $7000 \times 7000\text{ pixels}$, kích thước điểm ảnh $2.5\text{ m}$) và 4 dải phổ Multispectral gồm Blue ($0.45 - 0.52,\mu\text{m}$), Green ($0.53 - 0.59,\mu\text{m}$), Red ($0.62 - 0.69,\mu\text{m}$), NIR ($0.77 - 0.89,\mu\text{m}$) với kích thước ma trận $1750 \times 1750\text{ pixels}$, kích thước điểm ảnh $10.0\text{ m}$, mức lượng tử hóa 12-bit ghi nhận ở 10-bit truyền thông.
  • Dữ liệu đo quang phổ thực địa tại bãi Buôn Ma Thuột: Đo trực tiếp bằng máy đo phổ kế tán xạ (Field Spectroradiometer) đồng thời với thời điểm vệ tinh bay qua vào các ngày 14/11/2017 và 02/11/2018.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được chuẩn hóa qua 4 module thuật toán khép kín:

[Module 1: Hiệu chỉnh Dòng tối DS]
  │── Thu nhận cảnh ảnh biển sâu Đại Tây Dương / Cảnh đêm không có bức xạ
  │── Tính giá trị DN trung bình cho từng phần tử cảm biến j: DS(j) = (1/N) * Σ DN(i, j)
  └── Trừ nền dòng tối cho toàn bộ ảnh mức 0: DN_1(i, j) = DN_0(i, j) - DS(j)
        │
        ▼
[Module 2: Hiệu chỉnh Hồi đáp không đồng đều PRNU]
  │── Thu nhận cảnh ảnh đồng nhất sa mạc Algeria / Libya (PICS)
  │── Tính độ hồi đáp trung bình của toàn mảng cảm biến: <DN> = (1/M*N) * ΣΣ DN_1(i, j)
  │── Tính hệ số tăng ích tương đối: Gain(j) = <DN> / ((1/M) * Σ DN_1(i, j))
  └── Cân bằng mức xám: DN_2(i, j) = DN_1(i, j) * Gain(j)
        │
        ▼
[Module 3: Đánh giá Tỉ lệ Tín hiệu trên Nhiễu (SNR)]
  │── Trích xuất cửa sổ đồng nhất W (kích thước tối thiểu 50x50 pixels) trên bãi Buôn Ma Thuột
  │── Tách nhiễu cột (cảm biến) và nhiễu hàng (chuẩn hóa):
  │     σ_col = sqrt( (1/P) * Σ (DN_mean_col - DN_grand_mean)^2 )
  │     σ_row = sqrt( (1/Q) * Σ (DN_mean_row - DN_grand_mean)^2 )
  └── Tính toán chỉ số: SNR = m / sqrt(σ_col^2 + σ_row^2)
        │
        ▼
[Module 4: Đánh giá MTF bằng Thuật toán Canny Cải tiến]
  │── Cắt vùng quan tâm (ROI) chứa biên cạnh đen - trắng tại bãi Buôn Ma Thuột / Salon de Provence
  │── Lọc làm mịn ảnh bằng bộ lọc Gauss 2D: G(x, y, σ)
  │── Tính toán vector Gradient và hướng cạnh: ∇I = [Ix, Iy]^T, θ = arctan(Iy / Ix)
  │── Triệt tiêu các điểm không phải cực đại (Non-maximum suppression) và Phân ngưỡng Hysteresis
  │── Xác định đường bao cạnh chính xác ở mức dưới điểm ảnh (Sub-pixel Edge Fitting)
  │── Chiếu toàn bộ điểm ảnh lân cận vuông góc với đường cạnh để lập hàm ESF siêu mẫu (Super-sampled ESF)
  │── Làm mịn phân đoạn đa thức (Spline/Binning) và lấy vi phân để suy ra hàm LSF: LSF(x) = d(ESF)/dx
  └── Biến đổi Fourier nhanh (FFT) hàm LSF và lấy trị tuyệt đối: MTF(f) = |FFT{LSF}| / FFT{LSF}(0)

Data và phân tích

Toàn bộ quy trình thuật toán được lập trình phát triển thành công cụ phần mềm chuyên dụng (sử dụng ngôn ngữ C/C++ và Python, kết hợp thư viện OpenCV và tích hợp trong môi trường ENVI/IDL). Dữ liệu được kiểm tra tính vững (robustness checks) qua việc phân tích sai số giữa các chu kỳ đánh giá (2013, 2017, 2018, 2019, 2021), loại bỏ điểm dị biệt (outliers) bằng phương pháp 3-Sigma ($3\sigma$), và đối sánh chéo giữa hai hướng chuyển dịch của cạnh (chuyển từ đen sang trắng và từ trắng sang đen) cũng như giữa hai phương chuyển động (dọc hướng bay - Along-track và ngang hướng bay - Across-track).

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

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

Nghiên cứu của luận án cung cấp bức tranh định lượng chân thực, khách quan và chi tiết nhất từ trước đến nay về hiệu suất thực tế của hệ thống cảm biến vệ tinh VNREDSat-1:

+──────────────────────────────────────────────────────────────────────────────────────────────────+
|                    BẢNG TỔNG HỢP CÁC KẾT QUẢ ĐO ĐẠC VÀ ĐỐI SOÁT THỰC NGHIỆM                      |
+──────────────────────────┬──────────────────────┬──────────────────────┬─────────────────────────+
| THAM SỐ VẬT LÝ           | THIẾT KẾ DANH ĐỊNH   | THỰC NGHIỆM BÃI THỬ  | THỰC NGHIỆM BÃI THỬ     |
|                          | (ASTRIUM BAN ĐẦU)    | SALON DE PROVENCE    | BUÔN MA THUỘT (VIỆT NAM)|
+──────────────────────────┼──────────────────────┼──────────────────────┼─────────────────────────+
| SNR Kênh Toàn sắc (PAN)  | ≥ 90.0               | 112.4 – 118.6        | 108.5 – 115.2           |
| SNR Kênh Blue (B1)       | ≥ 150.0              | 185.3 – 192.0        | 178.4 – 186.7           |
| SNR Kênh Green (B2)      | ≥ 150.0              | 190.1 – 198.5        | 184.2 – 191.8           |
| SNR Kênh Red (B3)        | ≥ 150.0              | 182.4 – 189.7        | 176.9 – 183.5           |
| SNR Kênh NIR (B4)        | ≥ 150.0              | 175.6 – 184.2        | 169.8 – 177.3           |
| MTF Nyquist PAN (Dọc bay)| ≥ 0.150              | 0.198 – 0.205        | 0.189 – 0.196           |
| MTF Nyquist PAN (Ngang)  | ≥ 0.150              | 0.210 – 0.218        | 0.202 – 0.208           |
| MTF Nyquist MS (Dọc bay) | ≥ 0.200              | 0.235 – 0.248        | 0.228 – 0.239           |
| MTF Nyquist MS (Ngang)   | ≥ 0.200              | 0.252 – 0.264        | 0.241 – 0.255           |
+──────────────────────────┴──────────────────────┴──────────────────────┴─────────────────────────+

           DIỄN BIẾN HÀM TRUYỀN ĐIỀU BIẾN MTF TẠI TẦN SỐ NYQUIST CỦA VNREDSAT-1 (2013 - 2021)
  MTF
  0.30 ┼────────────────────────────────────────────────────────────────────────────────────
       │
  0.25 ┼     * Dải MS Ngang (0.264 -> 0.241)     * Dải MS Dọc (0.248 -> 0.228)
       │
  0.20 ┼     -------------------------------------------------------------------------------
       │     # Dải PAN Ngang (0.218 -> 0.202)    # Dải PAN Dọc (0.205 -> 0.189)
  0.15 ┼────────────────────────────────────────────────────────────────────────────────────
       │     ================== NGƯỠNG THIẾT KẾ TỐI THIỂU (MTF = 0.15) =====================
  0.10 ┼
       │
  0.05 ┼
       │
  0.00 ┴─────┬───────────────────────────┬───────────────────────────┬─────────────────────►
           Năm 2013                    Năm 2017                    Năm 2021 (Năm thứ 8)
  1. Hiệu suất hàm truyền điều biến (MTF) duy trì xuất sắc sau 8 năm: Tại tần số lấy mẫu giới hạn Nyquist ($f_{Nyq} = 0.5 \text{ cy/pixel}$), giá trị MTF thực nghiệm trên kênh PAN đo được tại bãi Buôn Ma Thuột đạt $0.189 - 0.196$ theo chiều dọc hướng bay và $0.202 - 0.208$ theo chiều ngang hướng bay. So sánh với ngưỡng thiết kế ban đầu ($\text{MTF} \ge 0.15$), chất lượng không gian của VNREDSat-1 sau 8 năm hoạt động liên tục trên không gian vẫn vượt chuẩn thiết kế hơn $26.0%$, chứng minh cơ cấu quang học phản xạ Korsch và hệ thống điều hòa nhiệt thụ động trên vệ tinh duy trì độ ổn định cơ-nhiệt tuyệt hảo.
  2. Tỉ lệ Tín hiệu trên Nhiễu (SNR) vượt mức thiết kế cơ sở: Giá trị SNR trên toàn bộ 4 dải đa phổ đều duy trì ở mức cao ($> 169.8$), vượt yêu cầu thiết kế danh định ($\text{SNR} \ge 150$). Kênh toàn sắc PAN đạt giá trị SNR $> 108.5$ (thiết kế $\ge 90$). Điều này khẳng định độ nhạy của linh kiện mảng cảm biến CCD và chuỗi mạch điện tử xử lý tín hiệu không bị suy thoái do bức xạ vũ trụ.
  3. Độ tương đồng cao giữa bãi kiểm định quốc nội và quốc tế: Sai số tính toán MTF và SNR giữa bãi Buôn Ma Thuột (Việt Nam) và bãi Salon de Provence (Pháp) nằm trong khoảng hẹp $3.2% - 4.5%$. Kết quả này khẳng định bãi kiểm định nhân tạo tại Buôn Ma Thuột hoàn toàn đáp ứng các tiêu chuẩn kỹ thuật khắt khe của quốc tế, đóng vai trò là cơ sở hiệu chuẩn chuẩn tắc cho các dòng vệ tinh quang học hiện tại và tương lai của Việt Nam.
  4. Hiệu quả vượt trội của thuật toán Canny trong việc khử nhiễu cạnh: So sánh với phương pháp đạo hàm bậc một trực tiếp của chuẩn ISO 12233, việc ứng dụng thuật toán Canny giúp giảm phương sai của hàm LSF tới $18.4%$, loại bỏ hoàn toàn các đỉnh nhiễu giả mạo do dị thường phản xạ vi mô trên bề mặt bê tông của ô mẫu, nâng độ chính xác của giá trị MTF tính toán lên mức $\pm 0.008$.

Implications đa chiều

+──────────────────────────────────────────────────────────────────────────────────────────────────+
|                                    MA TRẬN HỆ THỐNG IMPLICATIONS                                 |
+────────────────────────────────┬─────────────────────────────────────────────────────────────────+
| KHÍA CẠNH                      | NỘI DUNG VÀ TÁC ĐỘNG ĐỊNH LƯỢNG                                 |
+────────────────────────────────┼─────────────────────────────────────────────────────────────────+
| Lý thuyết & Học thuật          | - Mở rộng lý thuyết LSI Fourier trong quan trắc viễn thám thực địa.|
|                                | - Hoàn thiện mô hình suy giảm chất lượng cảm biến quang học quỹ |
|                                |   đạo thấp (LEO) trong điều kiện khí hậu nhiệt đới.             |
+────────────────────────────────┼─────────────────────────────────────────────────────────────────+
| Phương pháp luận               | - Thiết lập pipeline xử lý tự động Canny Sub-pixel ESF -> MTF.  |
|                                | - Khắc phục điểm yếu góc chụp của ISO 12233 cho vệ tinh linh hoạt.|
+────────────────────────────────┼─────────────────────────────────────────────────────────────────+
| Ứng dụng & Thực tiễn Kỹ thuật  | - Cung cấp bộ thông số hiệu chỉnh Calibration Parameter File (CPF)|
|                                |   chuẩn cho Trung tâm Khai thác Vệ tinh nhỏ (Viện Công nghệ vũ trụ).|
|                                | - Tối ưu hóa việc nâng cấp chất lượng ảnh thô trước khi thương mại.|
+────────────────────────────────┼─────────────────────────────────────────────────────────────────+
| Chính sách & Quy chuẩn Quốc gia| - Cơ sở khoa học sửa đổi Thông tư 10/2015/TT-BTNMT và Thông tư   |
|                                |   52/2016/TT-BTNMT của Bộ Tài nguyên và Môi trường.              |
|                                | - Thiết lập khung phân hạng ảnh theo tỉ lệ bản đồ địa hình.    |
+────────────────────────────────┴─────────────────────────────────────────────────────────────────+

Về mặt ứng dụng đo vẽ bản đồ, luận án đề xuất ma trận phân hạng chất lượng ảnh dựa trên MTF và SNR:

  • Mức chất lượng 1 (Rất tốt: $\text{MTF} \ge 0.18, \text{SNR} \ge 100$): Sử dụng trực tiếp để đo vẽ, thành lập và hiện chỉnh bản đồ địa hình tỷ lệ lớn $1:5.000$ (khi kết hợp mô hình số độ cao DSM chất lượng cao) và $1:10.000$.
  • Mức chất lượng 2 (Tốt: $0.12 \le \text{MTF} < 0.18, 80 \le \text{SNR} < 100$): Đáp ứng hoàn hảo yêu cầu thành lập bản đồ địa hình tỷ lệ $1:25.000$ và $1:50.000$, phân loại lớp phủ bề mặt và giám sát biến động tài nguyên rừng.
  • Mức chất lượng 3 (Trung bình: $0.08 \le \text{MTF} < 0.12, 50 \le \text{SNR} < 80$): Giới hạn sử dụng cho phân loại hiện trạng sử dụng đất tổng thể, cảnh báo thiên tai diện rộng, không dùng cho đo vẽ hình học chi tiết.
  • Mức chất lượng 4 (Không đạt: $\text{MTF} < 0.08$ hoặc $\text{SNR} < 50$): Dữ liệu bị loại bỏ hoặc phải chuyển sang quy trình xử lý tăng cường đặc biệt (Deconvolution/Restoration filtering).

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật mang tính khách quan:

  1. Rào cản thời tiết và mây mù đặc trưng nhiệt đới: Khí hậu Việt Nam với lượng mây bao phủ trung bình năm cao ($> 65%$) làm hạn chế tần suất thu nhận ảnh rõ nét tại bãi Buôn Ma Thuột, dẫn đến số lượng cảnh ảnh đạt chuẩn kiểm định trong năm không đồng đều giữa mùa mưa và mùa khô.
  2. Thiếu cơ cấu hiệu chỉnh nội tại chủ động trên vệ tinh: Vệ tinh VNREDSat-1 thuộc dòng vệ tinh nhỏ (khối lượng $\approx 115\text{ kg}$) nên không được trang bị hệ thống đèn chuẩn nội (on-board calibration lamps) hoặc tấm khuếch xạ mặt trời (solar diffuser) như các hệ thống vệ tinh cỡ lớn (Landsat, Sentinel-2), buộc công tác hiệu chuẩn bức xạ tuyệt đối phải phụ thuộc hoàn toàn vào phương pháp gián tiếp qua bãi PICS.
  3. Kích thước bãi kiểm định nhân tạo giới hạn: Kích thước các ô mẫu tại bãi Buôn Ma Thuột ($25\text{ m} \times 25\text{ m}$) được tối ưu cho dải PAN ($2.5\text{ m}$) nhưng đối với dải đa phổ MS ($10.0\text{ m}$), số lượng điểm ảnh thuần khiết (pure pixels) trên ô mẫu bị hạn chế, chịu ảnh hưởng nhất định của hiệu ứng biên (adjacency effect).

Chương trình nghiên cứu 10 năm tiếp theo (2022–2032) mở ra các hướng đột phá:

  • Hướng 1: Xây dựng mạng lưới bãi kiểm chuẩn tự động hóa hoàn toàn (Automated Vicarious Radiometric Calibration Network) tích hợp các trạm quan trắc bức xạ mặt trời tự động (AERONET) và camera toàn cảnh bầu trời tại nhiều vùng địa lý khác nhau của Việt Nam.
  • Hướng 2: Tích hợp mô hình chuyển xạ khí quyển nâng cao (Advanced Atmospheric Radiative Transfer Models: 6S code, MODTRAN-6) kết hợp dữ liệu sol khí thời gian thực từ vệ tinh Himawari để loại bỏ triệt để tán xạ khí quyển trong tính toán MTF/SNR.
  • Hướng 3: Mở rộng và chuẩn hóa quy trình đánh giá cho các hệ thống vệ tinh quan sát Trái đất thế hệ mới của Việt Nam (vệ tinh siêu phân giải không gian VNREDSat-2, vệ tinh siêu phổ Hyperspectral, vệ tinh radar khẩu độ tổng hợp SAR LOTUSat-1).

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

+──────────────────────────────────────────────────────────────────────────────────────────────────+
|                                    TÁC ĐỘNG VÀ LAN TỎA TOÀN DIỆN                                  |
+──────────────────────────┬───────────────────────────────────────────────────────────────────────+
| LĨNH VỰC TÁC ĐỘNG        | GIÁ TRỊ VÀ KẾT QUẢ ĐO LƯỜNG CỤ THỂ                                   |
+──────────────────────────┼───────────────────────────────────────────────────────────────────────+
| Học thuật & Xuất bản     | - Xuất bản các công trình trên tạp chí khoa học uy tín chuyên ngành    |
|                          |   Trắc địa - Bản đồ, Viễn thám và Công nghệ Vũ trụ.                  |
|                          | - Ước tính trích dẫn cao trong các nghiên cứu Cal/Val tại Đông Nam Á.  |
+──────────────────────────┼───────────────────────────────────────────────────────────────────────+
| Tự chủ Công nghệ Không gian| - Chấm dứt sự phụ thuộc vào phần mềm thương mại độc quyền nước ngoài. |
|                          | - Tiết kiệm hàng trăm nghìn USD chi phí dịch vụ hiệu chuẩn quốc tế.  |
+──────────────────────────┼───────────────────────────────────────────────────────────────────────+
| Hiệu quả Quản lý Nhà nước| - Chuẩn hóa nguồn dữ liệu đầu vào cho Hệ thống Dữ liệu Không gian      |
|                          |   Địa lý Quốc gia (NSDI).                                             |
|                          | - Nâng cao độ tin cậy trong giám sát sạt lở, lũ lụt và quy hoạch đô thị.|
+──────────────────────────┴───────────────────────────────────────────────────────────────────────+

Công trình đặt nền móng vững chắc cho việc quản trị chất lượng dữ liệu viễn thám cấp quốc gia, khẳng định vị thế của Việt Nam là một trong số ít các quốc gia trong khu vực ASEAN có khả năng tự chủ toàn diện quy trình kiểm chuẩn vệ tinh quan sát Trái đất sau phóng.

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

  • Nghiên cứu sinh và Giới nghiên cứu Viễn thám: Tiếp cận phương pháp luận mẫu mực về xử lý tín hiệu không gian Fourier, thuật toán Canny dưới điểm ảnh và mô hình tính toán SNR phân tách nguồn nhiễu.
  • Giảng viên và Viện nghiên cứu đầu ngành: Tài liệu tham khảo học thuật chuẩn mực để xây dựng giáo trình đào tạo sau đại học ngành Kỹ thuật Trắc địa - Bản đồ, Công nghệ Vũ trụ, Xử lý ảnh số tại các trường đại học kỹ thuật.
  • Đơn vị Vận hành và Khai thác Vệ tinh (Trung tâm Điều khiển và Khai thác Vệ tinh nhỏ - Viện Công nghệ Vũ trụ): Ứng dụng trực tiếp quy trình và mã nguồn phần mềm để cập nhật định kỳ tệp thông số CPF (Calibration Parameter File), giám sát sức khỏe cảm biến VNREDSat-1.
  • Cơ quan Quản lý Nhà nước (Cục Viễn thám Quốc gia, Cục Đo đạc, Bản đồ và Thông tin Địa lý Việt Nam - Bộ TN&MT): Cơ sở pháp lý và kỹ thuật để xây dựng định mức kinh tế kỹ thuật, chuẩn hóa quy định nghiệm thu và cấp phép dữ liệu ảnh viễn thám thương mại lưu hành tại Việt Nam.

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?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết Hệ thống Dịch chuyển Tuyến tính Bất biến (LSI) và Mô hình Phân giải Tần số Không gian Fourier cho cảm biến viễn thám hoạt động dài hạn ngoài không gian. Luận án đã chứng minh một cách tường minh rằng sự suy thoái chất lượng hình học không gian của ảnh vệ tinh không diễn ra đồng nhất trên toàn dải tần số mà tập trung làm sụt giảm biên độ điều biến ở các tần số cao ($f > 0.3 \text{ cy/pixel}$) do hiện tượng tán xạ vi mô trên bề mặt thấu kính/gương phản xạ và suy giảm hiệu quả tích nạp điện tích (TDI - Time Delay Integration). Nghiên cứu tích hợp thành công mô hình lan truyền bức xạ quang học với giải thuật gradient phi tuyến Canny để chuẩn hóa chuỗi chuyển đổi $PSF \rightarrow ESF \rightarrow LSF \rightarrow MTF$.

2. Sự đổi mới trong phương pháp nghiên cứu khi so sánh với ít nhất 2 công trình quốc tế?

  • So sánh với Chuẩn quốc tế ISO 12233 (slanted-edge method) vốn sử dụng phép hồi quy tuyến tính đơn giản để dò tìm cạnh, phương pháp của luận án vượt trội nhờ áp dụng toán tử lọc Gauss 2D kết hợp phân ngưỡng kép Hysteresis của Canny, giúp loại bỏ hiện tượng méo dạng hàm ESF khi chụp ở các góc nghiêng lớn ($> 10^\circ$) hoặc khi bề mặt bãi thử có độ nhám vi mô.
  • So sánh với công trình của Scott (2003) trên vệ tinh QuickBird (DigitalGlobe) vốn áp dụng phép nội suy 8-Point Sinc nhạy cảm với hiện tượng dao động Gibbs, thuật toán của luận án sử dụng kỹ thuật phân rã siêu mẫu (super-sampling binning) kết hợp làm mịn Spline bậc ba, giúp giảm thiểu phương sai ước tính MTF tại tần số Nyquist xuống dưới $2.1%$, mang lại độ ổn định tính toán cao hơn ngay cả trong điều kiện ảnh hưởng của tán xạ khí quyển nhiệt đới.

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ờ và có giá trị thực tiễn cao nhất là: Sau hơn 8 năm bay trên quỹ đạo (vượt 60% tuổi thọ thiết kế 5 năm), hệ thống quang học VNREDSat-1 hầu như không bị lệch tiêu cự cơ - quang học (optical defocusing). Giá trị MTF tại tần số Nyquist trên kênh toàn sắc vẫn duy trì ở mức $0.189 - 0.208$, cao hơn đáng kể so với ngưỡng thiết kế ban đầu ($\ge 0.15$). Hiện tượng suy giảm nhẹ ($< 6%$) chủ yếu diễn ra do sự tích lũy tín hiệu dòng tối (Dark Signal) và sự biến đổi nhỏ của độ hồi đáp không đồng đều (PRNU) trên các phần tử biên của mảng cảm biến CCD, hoàn toàn có thể triệt tiêu triệt để bằng việc cập nhật tệp thông số hiệu chuẩn CPF định kỳ 6 tháng một lần theo quy trình luận án đề xuất.

                    MÔ THỨC QUY TRÌNH TÁI LẬP HIỆU CHUẨN ĐỊNH KỲ (CPF PIPELINE)
                    
  +─────────────────────────────────────────────────────────────────────────────────────────+
  | KỲ QUAN TRẮC ĐỊNH KỲ (6 THÁNG / LẦN)                                                    |
  +─────────────────────────────────────────────────────────────────────────────────────────+
  |  1. Chụp Cảnh Đêm / Biển Sâu ────► Trích xuất DS Matrix [7000 x 1]                      |
  |  2. Chụp PICS Algeria/Libya   ────► Trích xuất Gain Normalization PRNU Vector [7000 x 1] |
  |  3. Chụp Bia Buôn Ma Thuột    ────► Đo phổ thực địa + Tính SNR & MTF Nyquist            |
  |  4. Cập nhật Tệp CPF Mới      ────► Nạp vào Trạm Xử lý Thu nhận Mặt đất (Ground Station)|
  +─────────────────────────────────────────────────────────────────────────────────────────+

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp bộ quy trình công nghệ chuẩn mực và khả thi 100% để tái lập thực nghiệm, bao gồm:

  • Toàn bộ sơ đồ khối và thuật giải chi tiết từng bước (Step-by-step algorithms) từ xử lý Level 0 sang Level 1A.
  • Mã nguồn công cụ phần mềm (C/C++ và Python) tại phần Phụ lục, cho phép tự động hóa hoàn toàn các thao tác: đọc tệp ảnh viễn thám dạng nhị phân thô, khử dòng tối DS, cân bằng PRNU, trích xuất ROI cạnh nghiêng, phát hiện biên Canny, khớp hàm ESF/LSF và thực hiện biến đổi Fourier FFT.
  • Bảng thông số tọa độ địa lý chính xác, quy cách kích thước hình học và hệ số phản xạ phổ quang học chuẩn của từng ô mẫu tại bãi kiểm định Buôn Ma Thuột và bãi Salon de Provence.

5. Khung lộ trình nghiên cứu 10 năm được phác thảo như thế nào?

Luận án xác lập lộ trình 3 giai đoạn:

  • Giai đoạn 1 (2022–2024): Hoàn thiện tích hợp phần mềm đánh giá MTF/SNR vào quy trình tiếp nhận dữ liệu tự động tại Đài Viễn thám Trung ương; công bố bộ tiêu chuẩn kỹ thuật quốc gia về đánh giá chất lượng ảnh vệ tinh quang học.
  • Giai đoạn 2 (2025–2028): Nâng cấp bãi kiểm định Buôn Ma Thuột thành Trạm kiểm chuẩn quang học tự động hóa (Automated Cal/Val Site); phát triển thuật toán đánh giá cho vệ tinh độ phân giải siêu cao VNREDSat-2 ($GSD \le 0.5\text{ m}$) và chùm vệ tinh Micro/CubeSat.
  • Giai đoạn 3 (2029–2032): Mở rộng mạng lưới kiểm chuẩn liên hoàn cho các dải phổ siêu cao tần (Hyperspectral) và radar khẩu độ tổng hợp (SAR), đưa hệ thống kiểm định của Việt Nam gia nhập Mạng lưới Hiệu chuẩn Quốc tế CEOS WGCV (Committee on Earth Observation Satellites - Working Group on Calibration and Validation).

Kết luận

Công trình nghiên cứu luận án tiến sĩ của tác giả Nguyễn Minh Ngọc đại diện cho một bước tiến dài của ngành Kỹ thuật Trắc địa – Bản đồ và Viễn thám Việt Nam. Luận án kết tinh 5 đóng góp cốt lõi:

  1. Đột phá phương pháp luận: Xác lập cơ sở khoa học và chứng minh tường minh rằng bộ đôi chỉ số Hàm truyền điều biến (MTF) và Tỉ lệ tín hiệu trên nhiễu (SNR) là tiêu chí đánh giá định lượng toàn diện nhất cho chất lượng ảnh quang học vệ tinh tại Việt Nam.
  2. Cải tiến kỹ thuật thuật toán: Khởi xướng và thực nghiệm thành công việc áp dụng toán tử vi phân Canny trong trích xuất cạnh nghiêng dưới điểm ảnh, khắc phục triệt để nhược điểm của các thuật toán tuyến tính cổ điển, nâng cao độ chính xác xác định hàm LSF và MTF trong điều kiện môi trường có nhiễu.
  3. Hiện thực hóa quy trình công nghệ độc lập: Xây dựng thành công quy trình đánh giá và hiệu chỉnh chất lượng ảnh viễn thám từ dữ liệu thô (Level 0 và Level 1A), giúp các đơn vị vận hành vệ tinh trong nước hoàn toàn làm chủ công nghệ mà không cần phụ thuộc vào nhà sản xuất nước ngoài.
  4. Kiểm định và khẳng định chất lượng VNREDSat-1: Cung cấp bằng chứng thực nghiệm vững chắc chứng minh vệ tinh VNREDSat-1 sau 8 năm bay trên quỹ đạo vẫn duy trì chất lượng bức xạ và hình học vượt chuẩn thiết kế ban đầu, dữ liệu hoàn toàn đáp ứng độ chính xác để thành lập bản đồ địa hình tỷ lệ lớn $1:5.000$ đến $1:10.000$.
  5. Tiêu chuẩn hóa bãi kiểm định quốc gia: Khẳng định giá trị khoa học và độ tin cậy tương đương quốc tế của bãi kiểm định nhân tạo Buôn Ma Thuột, mở ra kỷ nguyên tự chủ kiểm định và chuẩn hóa dữ liệu cho toàn bộ các thế hệ vệ tinh quan sát Trái đất tương lai của Việt Nam.