Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Dữ liệu ảnh viễn thám quang học đóng vai trò quan trọng trong nhiều lĩnh vực kinh tế - xã hội như đo đạc bản đồ, quản lý đất đai, giám sát tài nguyên và môi trường. Người sử dụng thông thường tiếp cận dữ liệu ảnh ở các mức xử lý thương mại (mức 1A, 2A, 3A, 3B) và đánh giá chất lượng dựa trên các tiêu chí độ phân giải danh định (không gian, phổ, bức xạ, thời gian). Tuy nhiên, các thông số thiết kế ban đầu trong phòng thí nghiệm thường bị suy giảm khi vệ tinh hoạt động trong môi trường không gian khắc nghiệt do rung lắc khi phóng, biến dạng nhiệt, suy thoái linh kiện quang điện tử và vận động động lực học của vệ tinh.

Việc đánh giá tình trạng hoạt động thực tế của thiết bị chụp ảnh trên quỹ đạo thông qua dữ liệu ảnh thô (mức 0 và mức 1A) là bài toán kỹ thuật phức tạp, thường do nhà sản xuất hoặc cơ quan quản lý vệ tinh thực hiện. Tại Việt Nam, sau khi tiếp nhận và vận hành hệ thống vệ tinh quan sát Trái Đất đầu tiên (VNREDSat-1), công tác kiểm soát chất lượng ảnh chủ yếu dựa trên tài liệu hướng dẫn chuyển giao kỹ thuật của nhà sản xuất (Astrium, Pháp) với các thông số hạn chế, hoặc tập trung vào xử lý nâng cao ở các mức sản phẩm cuối. Việt Nam còn thiếu hệ thống phòng thí nghiệm chuyên dụng và các mô hình mô phỏng cơ khí - quang học phức tạp. Từ năm 2017, khi bãi kiểm định nhân tạo tại thành phố Buôn Ma Thuột (tỉnh Đắk Lắk) được đưa vào khai thác, điều kiện thực tế đã cho phép thiết lập cơ sở khoa học và quy trình định lượng độc lập để đánh giá chất lượng ảnh viễn thám quang học cho các thế hệ vệ tinh của Việt Nam.

Mục tiêu nghiên cứu

Mục tiêu nghiên cứu của luận án được xác định gồm hai nội dung:

  1. Đề xuất thông số để đánh giá chất lượng ảnh viễn thám quang học phù hợp với điều kiện của Việt Nam.
  2. Lựa chọn phương pháp và đề xuất quy trình đánh giá chất lượng ảnh viễn thám quang học đối với điều kiện của Việt Nam.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: 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 của dữ liệu ảnh viễn thám quang học.
  • Phạm vi không gian:
    • Bãi kiểm định nhân tạo: Bãi kiểm định tại Salon de Provence (Pháp) và bãi kiểm định nhân tạo tại thành phố Buôn Ma Thuột (tỉnh Đắk Lắk, Việt Nam).
    • Bãi kiểm định tự nhiên: Các khu vực sa mạc tại Châu Phi (Libya, Algeria) và vùng biển sâu tại Đại Tây Dương.
  • Phạm vi dữ liệu và thời gian: Dữ liệu ảnh vệ tinh VNREDSat-1 của Việt Nam ở mức 0 và mức 1A; số liệu đo quang phổ phản xạ bề mặt tại bãi kiểm định Buôn Ma Thuột thu thập trong hai đợt năm 2017 và năm 2018.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Công tác đánh giá chất lượng hệ thống chụp ảnh vệ tinh trên thế giới được chia thành ba giai đoạn: trước khi phóng (trong phòng thí nghiệm), khi mới vào quỹ đạo (hiệu chỉnh ban đầu) và hoạt động định kỳ trên quỹ đạo. Các phương pháp đánh giá gián tiếp qua dữ liệu ảnh tập trung vào hai nhóm thông số không gian và bức xạ:

  • Đánh giá chất lượng bức xạ và thông số SNR:

    • Nghiên cứu đánh giá chéo bức xạ giữa Landsat 8 (OLI/TIR) và Landsat 7 trên khu vực sa mạc Libya cho thấy độ sai khác bức xạ khoảng 3% [72].
    • Nghiên cứu hiệu chỉnh bức xạ và khử nhiễu xạ giữa hai kênh phổ VNIR và SWIR của bộ cảm ASTER [9].
    • Hiệu chỉnh độ tương phản phi tuyến tính trên 8 kênh ảnh WorldView-2 tại 7 điểm kiểm tra thực địa [96].
    • Đánh giá SNR trên hệ thống vệ tinh SPOT dựa trên phân tách mô hình nhiễu hàng và nhiễu cột qua ma trận $50 \times 50$ điểm ảnh trên vùng sa mạc hoặc băng tuyết [62].
    • Cơ quan Hàng không và Vũ trụ Hoa Kỳ (NASA) ứng dụng phương pháp địa thống kê (geostatistics), sử dụng phép lọc notch loại trừ nhiễu chu kỳ và biểu đồ bán phương sai (semivariogram) để tách nhiễu cảm biến cho dữ liệu phổ siêu cao AVIRIS [29].
    • Hệ thống Sentinel-2 (Cơ quan Vũ trụ Châu Âu - ESA) đánh giá SNR hàng tháng qua ảnh có khuếch tán mặt trời và ảnh dòng tối, đánh giá định kỳ 6 tháng qua các bề mặt đồng nhất tự nhiên (sa mạc, Greenland, Nam Cực) theo mô hình tán xạ Lambert [67].
  • Đánh giá chất lượng không gian và thông số MTF:

    • Boreman [20] và Ronald cùng cộng sự [85] chuẩn hóa lý thuyết MTF là giá trị tuyệt đối của biến đổi Fourier từ hàm lan truyền điểm (Point Spread Function - PSF) trong hệ thống dịch chuyển tuyến tính bất biến (Linear Shift Invariant - LSI).
    • Paul W. Scott [90] (DigitalGlobe) đánh giá MTF của QuickBird qua 16 cảnh ảnh (2001–2003) bằng thuật toán cạnh ("gờ") với vật chuẩn sơn trên mặt đất, so sánh sự khác biệt giữa phép nội suy song tuyến và nội suy 8-Point Sinc, chỉ ra giá trị MTF tại tần số Nyquist theo hai hướng dọc/ngang quỹ đạo đạt yêu cầu thiết kế ($> 0$).
    • Đánh giá MTF cho dữ liệu OrbView-3 (11-bit) tại các bãi kiểm định kẻ ô đen trắng và mép các tòa nhà cao tầng [56].
    • Vệ tinh Hyperion [76] và Landsat-7 ETM+ (James C. Storey, Trung tâm EROS/USGS [98]) sử dụng thuật toán cạnh trên cấu trúc nhân tạo dạng xung/cạnh là cầu qua vịnh Destin (Florida, rộng 13,02 m) và cầu hồ Pontchartrain Causeway (Louisiana) để tính toán hàm lan truyền đường (Line Spread Function - LSF) và trích xuất MTF tại tần số Nyquist (đạt 0,39 – 0,42).
    • Tiêu chuẩn quốc tế ISO 12233:2000 và ISO 12233:2017 xác lập phương pháp cạnh nghiêng (slanted-edge) tính MTF cho cảm biến CMOS/CCD [39]. Hang Li, Changxiang Yan và Jianbing Chao [64] cải tiến phương pháp cạnh nghiêng cho hệ thống hồng ngoại kết hợp kỹ thuật chiết tách góc cạnh và mô hình hàm lan truyền cạnh (Edge Spread Function - ESF) thích ứng.
    • Wu Zhaocong và cộng sự [109] áp dụng phương pháp đánh giá toàn diện mờ (fuzzy comprehensive evaluation) kết hợp 2 thuộc tính MTF và 7 chỉ số bức xạ (Nyquist, MTF 0.5, entropy, peak SNR, sai khác trung bình, cường độ cạnh, gradient trung bình, độ tương phản, khoảng cách lấy mẫu mặt đất - GSD) để so sánh định lượng ảnh GF2 và SPOT7.
    • Jin Li và cộng sự [65] đề xuất quy trình 5 bước ước tính MTF không phụ thuộc vật chuẩn trên ảnh SPOT5 và QuickBird. Latry, Despringre và Valorge [57] phát triển phương pháp bình phương nhỏ nhất dựa trên cặp ảnh độ phân giải kép (tỉ lệ phân giải ~3).

Các nghiên cứu tại Việt Nam và khoảng trống khoa học

Tại Việt Nam, các nghiên cứu trước đây chủ yếu tập trung vào dữ liệu đã qua xử lý mức cao:

  • Lương Chính Kế [3] (2011) nghiên cứu chuẩn hóa bức xạ phổ đa thời gian và kỹ thuật trộn ảnh.
  • Chu Hải Tùng [5] nghiên cứu mô hình hóa hình học, nắn chỉnh và cân bằng dàn màu phổ cho ảnh VNREDSat-1.
  • Nguyễn Văn Hùng [2] và Đặng Trường Giang [1] khảo sát chất lượng phổ, độ sắc nét và quản lý chất lượng mảnh ghép bình đồ ảnh vệ tinh (SPOT 6/7, Pleiades, VNREDSat-1, Sentinel-2, Landsat-8) ở mức sản phẩm 3A, 3B.
  • Các văn bản quy phạm kỹ thuật của Bộ Tài nguyên và Môi trường (Thông tư số 10/2015/TT-BTNMT [6] và Thông tư số 52/2016/TT-BTNMT [7]) quy định kiểm tra độ phủ mây, độ đồng đều phổ, độ chính xác tọa độ hình học và tỉ lệ lỗi bit (Bit Error Rate - BER).
  • Nghiêm Văn Tuấn và cộng sự [4] bước đầu đề xuất đánh giá các thông số hình học (MTF, GSD) và bức xạ (SNR, tín hiệu tối DS, hồi đáp không đồng đều PRNU, dải động bức xạ) cho VNREDSat-1 so với thiết kế xuất xưởng.

Khoảng trống nghiên cứu: Chưa có công trình nào tại Việt Nam xây dựng quy trình định lượng hoàn chỉnh để đánh giá độc lập chất lượng ảnh thô (mức 0 và mức 1A) thông qua các thông số vật lý cốt lõi (SNR và MTF) có tính đến điều kiện cơ sở hạ tầng bãi kiểm định trong nước, đồng thời liên kết kết quả đo đạc kỹ thuật với phân cấp mức chất lượng phục vụ từng mục đích sử dụng bản đồ cụ thể.


Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và các khái niệm kỹ thuật

Luận án dựa trên lý thuyết truyền xạ điện từ trong viễn thám quang học, lý thuyết xử lý tín hiệu và quang học Fourier:

  • Hệ thống tuyến tính bất biến (LSI): Giả thiết hệ thống thu nhận ảnh vệ tinh thỏa mãn tính chất xếp chồng tuyến tính và bất biến dịch chuyển: $$L{a f(x) + b g(x)} = a L{f(x)} + b L{g(x)} \quad \text{và} \quad L{f(x - c)} = g(x - c)$$
  • Hàm truyền điều biến (MTF): Được xác định từ độ biến thiên tương phản $M = \frac{A_{\max} - A_{\min}}{A_{\max} + A_{\min}}$. Về mặt toán học, MTF là biên độ của biến đổi Fourier 2 chiều của hàm lan truyền điểm $PSF(x, y)$: $$MTF(f_x, f_y) = |\mathcal{F}{PSF(x, y)}|$$ MTF tổng thể của hệ thống là tích của các thành phần suy giảm độc lập: $$MTF_{\text{hệ thống}} = MTF_{\text{quang học}} \times MTF_{\text{cảm biến}} \times MTF_{\text{điện tử}} \times MTF_{\text{chuyển động}} \times MTF_{\text{khí quyển}}$$
  • Tỉ lệ tín hiệu trên nhiễu (SNR): Biểu diễn tỉ số giữa mức năng lượng tín hiệu trung bình thu nhận từ vùng đồng nhất ($m$) và độ lệch chuẩn của nhiễu ($\sigma$): $$SNR = \frac{m}{\sigma}$$ Nhiễu bức xạ trên bộ cảm mảng CCD được tách biệt thành nhiễu thiết bị (nhiễu cột, phân bố Poisson) và nhiễu chuẩn hóa (nhiễu hàng).
  • Tiêu chuẩn Rayleigh và độ phân giải không gian: Quy định giới hạn phân giải khi điểm cực đại của mô hình nhiễu xạ điểm thứ nhất trùng với cực tiểu đầu tiên của mô hình nhiễu xạ điểm thứ hai.

Phương pháp nghiên cứu cụ thể

Luận án kết hợp các phương pháp nghiên cứu định tính và định lượng:

  • Phương pháp phân tích, tổng hợp lý thuyết: Phân tích tài liệu kỹ thuật, kế thừa các mô hình toán học và kết quả vận hành hệ thống vệ tinh VNREDSat-1.
  • Phương pháp thu thập và đo đạc thực địa: Đo phản xạ phổ bề mặt bãi kiểm định Buôn Ma Thuột (Đắk Lắk) bằng thiết bị đo quang phổ chuyên dụng trong hai năm 2017 và 2018; thu thập ảnh vệ tinh mức 0 và mức 1A tại các bãi kiểm định nhân tạo và tự nhiên.
  • Phương pháp xử lý ảnh viễn thám định lượng:
    • Đánh giá SNR: Áp dụng phương pháp độ lệch chuẩn cục bộ (Local Standard Deviation - LSD) chia ảnh thành các ô $4 \times 4$ pixel trên khu vực đồng nhất, phân loại 150 khoảng giá trị độ lệch chuẩn giữa giá trị nhỏ nhất và 1,2 lần giá trị trung bình để lọc giá trị nhiễu thực tế. Đồng thời phân tích so sánh với phương pháp tương quan không gian - phổ và phương pháp địa thống kê.
    • Đánh giá MTF: Sử dụng bãi kiểm định dạng cạnh Heaviside (yêu cầu khoảng chuyển tiếp $L_T > 3 \div 5 \times GSD$). Tích hợp thuật toán toán tử tách biên Canny để tự động dò tìm góc nghiêng và chiết tách đường biên cạnh chính xác trên ảnh phân giải cao, thay thế thuật toán xấp xỉ tuyến tính cổ điển; từ đó tính toán ESF, lấy đạo hàm thu được LSF và biến đổi Fourier ra đồ thị MTF.
  • Phương pháp thực nghiệm và lập trình: Xây dựng module phần mềm chuyên dụng để tự động hóa quá trình tính toán các tham số SNR, MTF, DS, PRNU từ tập dữ liệu ảnh thô.
  • Phương pháp chuyên gia: Tham vấn ý kiến chuyên gia trong lĩnh vực công nghệ vũ trụ và kỹ thuật viễn thám để xác lập ngưỡng chất lượng ảnh theo nhu cầu sử dụng thực tế.
Nhóm thông số Thông số kỹ thuật Khả năng đánh giá tại Việt Nam Lý do lựa chọn / Loại trừ trong luận án
Không gian MTF Thực hiện được Được chọn: Phản ánh độ sắc nét, tương phản; tính được qua bãi kiểm định Buôn Ma Thuột và Canny.
GSD Thực hiện được Không chọn: Hiệu chỉnh hình học trên mặt đất đã đảm bảo theo thiết kế.
Độ méo ảnh Không khả thi Loại trừ: Đòi hỏi mô hình mô phỏng thiết kế quang học chi tiết.
Khoảng cách lấy mẫu góc Không khả thi Loại trừ: Đòi hỏi mô hình mô phỏng động lực học vệ tinh.
Chất lượng định hướng Không khả thi Loại trừ: Chưa có khung cơ học mẫu thử mô phỏng trên mặt đất.
Độ rộng dải chụp Thực hiện được Loại trừ: Tốn kém kinh phí và không ảnh hưởng trực tiếp đến chất lượng từng điểm ảnh.
Bức xạ SNR Thực hiện được Được chọn: Đặc trưng cho mức độ nhiễu bức xạ; tính được trên vùng đồng nhất và bãi kiểm định.
Tín hiệu tối (DS) Thực hiện được Được chọn làm bước tiền xử lý: Đánh giá trên ảnh đại dương/sa mạc ở mức 0 để chuẩn hóa trước khi tính SNR.
Hồi đáp không đồng đều (PRNU) Thực hiện được Được chọn làm bước tiền xử lý: Đánh giá và cân bằng dải cảm biến ở mức 0 trước khi tính SNR/MTF.
Dải động bức xạ Không khả thi Loại trừ: Chưa có hệ thống mục tiêu chuẩn đạt độ trắng tuyệt đối tại thực địa.

Nội dung chính theo từng chương

Chương 1: Tổng quan về vấn đề nghiên cứu

Chương 1 hệ thống hóa khái niệm về chất lượng ảnh viễn thám quang học, bao gồm các chỉ tiêu kỹ thuật không gian và bức xạ, phân biệt chất lượng ảnh tuyệt đối và tương đối. Tác giả phân tích tính tất yếu của việc kiểm soát chất lượng ảnh qua 3 giai đoạn hoạt động của vệ tinh, đặc biệt nhấn mạnh giai đoạn vận hành trên quỹ đạo cần đánh giá gián tiếp định kỳ qua dữ liệu ảnh.

Chương này tổng thuật chi tiết các công trình quốc tế về đánh giá SNR (trên ảnh đơn, ảnh tổng hợp mây dày trên đại dương, phương pháp LSD, địa thống kê, tương quan phổ) và MTF (bãi kiểm định dạng cạnh, dạng xung, dạng chu kỳ, phương pháp cạnh nghiêng ISO 12233, thuật toán mờ, lọc hình thái học, phân giải kép). Tại Việt Nam, tác giả điểm lại các nghiên cứu xử lý ảnh của Lương Chính Kế, Chu Hải Tùng, Nguyễn Văn Hùng, Đặng Trường Giang, Nghiêm Văn Tuấn và hệ thống văn bản quy phạm kỹ thuật hiện hành. Trên cơ sở phân tích điều kiện trang thiết bị và sự sẵn có của bãi kiểm định Buôn Ma Thuột (hoạt động từ năm 2017), tác giả luận giải tính khoa học và thực tiễn của việc lựa chọn hai thông số trọng tâm là SNR và MTF (kết hợp hiệu chỉnh DS và PRNU) để đánh giá chất lượng ảnh viễn thám quang học tại Việt Nam.

Chương 2: Cơ sở khoa học và phương pháp đánh giá chất lượng ảnh viễn thám quang học

Chương 2 thiết lập cơ sở toán học và cơ chế vật lý hình thành chất lượng ảnh trên hệ thống thu nhận quang điện tử. Tác giả phân tích chi tiết ảnh hưởng của hiện tượng quang sai (quang sai màu, quang sai đơn sắc) và hiện tượng nhiễu xạ ánh sáng (đĩa sáng Airy, tiêu chuẩn tách Rayleigh).

Về mặt bức xạ, chương phân tích độ phân giải bức xạ qua các đại lượng $NE\Delta\rho$, $NE\Delta T$, $NE\Delta L$ và ảnh hưởng của nhiễu đối với dữ liệu mã hóa $\ge 11\text{-bit}$. Mô hình toán học của SNR được thiết lập chi tiết, phân tách hai thành phần: nhiễu cột (nhiễu thiết bị do dao động Poisson) và nhiễu hàng (nhiễu chuẩn hóa). Tác giả làm rõ vai trò của việc loại trừ sai số hệ thống thông qua hiệu chỉnh tín hiệu tối (DS) và mức độ hồi đáp không đồng đều của điểm ảnh (PRNU) từ dữ liệu ảnh mức 0 chụp trên các bãi kiểm định giả bất biến (sa mạc Libya, Algeria, biển sâu Đại Tây Dương).

Về mặt không gian, chương trình bày chi tiết lý thuyết hàm truyền điều biến MTF, hàm lan truyền điểm PSF, hàm lan truyền đường LSF và hàm lan truyền cạnh ESF. Tác giả phân tích đồ thị đường cong MTF suy giảm theo tần số không gian và sự phụ thuộc của MTF vào chuyển động vệ tinh và sol khí quyển. Phân tích so sánh các dạng bãi kiểm định cho thấy bãi kiểm định dạng cạnh Heaviside nhân tạo là giải pháp tối ưu cho vệ tinh quang học phân giải cao nhằm khắc phục hiệu ứng răng cưa (aliasing). Tác giả đề xuất ứng dụng toán tử Canny để chiết tách biên cạnh chính xác trên ảnh thực tế thay cho các thuật toán tuyến tính đơn giản.

Chương 3: Quy trình đánh giá chất lượng ảnh viễn thám quang học

Chương 3 xây dựng quy trình công nghệ gồm các bước đánh giá chất lượng ảnh viễn thám quang học phù hợp với điều kiện Việt Nam:

  1. Thu nhận và tiền xử lý dữ liệu ảnh thô mức 0 trên các bãi kiểm định tự nhiên (sa mạc, đại dương) để tính toán, hiệu chỉnh giá trị DS và PRNU.
  2. Xử lý dữ liệu ảnh mức 1A chụp bãi kiểm định Buôn Ma Thuột và Salon de Provence: áp dụng phương pháp độ lệch chuẩn cục bộ LSD để tính toán SNR trên các vùng đồng nhất với nhiều mức bức xạ khác nhau.
  3. Ứng dụng thuật toán Canny trích xuất cạnh bãi kiểm định, xây dựng đường cong ESF, đạo hàm lấy LSF và chuyển đổi Fourier để xác định giá trị MTF tại tần số Nyquist theo hai hướng dọc và ngang quỹ đạo bay.
  4. Tổng hợp các giá trị định lượng thu được, đối chiếu với chỉ tiêu thiết kế kỹ thuật ban đầu của nhà sản xuất.
  5. Xây dựng ma trận phân cấp chất lượng ảnh dựa trên giá trị MTF và SNR kết hợp với yêu cầu độ chính xác của người sử dụng đối với từng tỉ lệ thành lập bản đồ địa hình và ứng dụng chuyên ngành.

Chương 4: Thực nghiệm: đánh giá chất lượng ảnh viễn thám quang học VNREDSat-1 của Việt Nam

Chương 4 triển khai thực nghiệm toàn bộ quy trình trên hệ thống dữ liệu thực tế của vệ tinh VNREDSat-1:

  • Tính toán tín hiệu tối DS và mức độ hồi đáp không đồng đều PRNU từ các cảnh ảnh mức 0 chụp khu vực biển Đại Tây Dương và sa mạc Bắc Phi (Libya, Algeria).
  • Thực nghiệm đo đạc phổ phản xạ thực địa tại bãi kiểm định Buôn Ma Thuột (Đắk Lắk) vào các năm 2017 và 2018 để đối soát tham số phản xạ phổ.
  • Ứng dụng module phần mềm xây dựng trong luận án để trích xuất cạnh bãi kiểm định bằng thuật toán Canny trên ảnh toàn sắc (Panchromatic - Pan, độ phân giải 2,5 m) và ảnh đa phổ (Multispectral - MS, độ phân giải 10 m) của VNREDSat-1 mức 1A chụp bãi kiểm định Buôn Ma Thuột và bãi kiểm định Salon de Provence (Pháp).
  • Tính toán giá trị MTF tại tần số Nyquist và giá trị SNR trên các kênh phổ.
  • Kết quả thực nghiệm khẳng định chất lượng hình ảnh và trạng thái hoạt động của hệ thống quang điện tử trên vệ tinh VNREDSat-1 vẫn duy trì ổn định, đáp ứng các tiêu chuẩn kỹ thuật thiết kế ban đầu sau nhiều năm vận hành trên quỹ đạo. Đồng thời, chương đưa ra bảng khuyến nghị phân mức chất lượng ảnh VNREDSat-1 cho các bài toán ứng dụng thực tế tại Việt Nam.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Đề xuất thành công cơ sở khoa học và lựa chọn bộ thông số đại diện cốt lõi gồm tỉ lệ tín hiệu trên nhiễu (SNR - đại diện chất lượng bức xạ) và hàm truyền điều biến (MTF - đại diện chất lượng không gian) phù hợp với điều kiện cơ sở hạ tầng kỹ thuật và thực tiễn dữ liệu viễn thám tại Việt Nam.
  2. Xây dựng hoàn chỉnh quy trình công nghệ đánh giá chất lượng ảnh viễn thám quang học từ mức dữ liệu thô (mức 0 và mức 1A), kết hợp tiền xử lý chuẩn hóa DS, PRNU với đánh giá SNR và MTF.
  3. Ứng dụng thành công phương pháp chiết tách cạnh Canny thay thế cho các thuật toán tuyến tính truyền thống trong quy trình tính toán MTF từ ảnh chụp bãi kiểm định, nâng cao độ chính xác và tính thuận tiện khi đánh giá chất lượng ảnh vệ tinh độ phân giải cao.

Đóng góp mới về mặt thực tiễn và giải pháp đề xuất

  1. Chứng minh bằng số liệu thực nghiệm định lượng rằng chất lượng ảnh của hệ thống vệ tinh quan sát Trái Đất đầu tiên của Việt Nam (VNREDSat-1) vẫn đảm bảo các chỉ tiêu kỹ thuật thiết kế xuất xưởng cho đến thời điểm đánh giá, cung cấp cơ sở tin cậy cho công tác vận hành và khai thác dữ liệu.
  2. Quy trình và công cụ phần mềm do luận án đề xuất có tính mở, có thể chuyển giao áp dụng trực tiếp cho các trạm thu và xử lý ảnh viễn thám hiện có của Việt Nam (đang thu nhận ảnh VNREDSat-1, SPOT 6/7) cũng như sẵn sàng áp dụng để kiểm định chất lượng dữ liệu của các hệ thống vệ tinh tương lai (như KOMPSAT-3 và các vệ tinh quang học siêu cao, siêu phổ tiếp theo).
  3. Đề xuất khung phân cấp mức chất lượng ảnh theo ngưỡng SNR và MTF tương ứng với từng nhu cầu sử dụng cụ thể trên thị trường và trong các ngành quản lý tài nguyên, đo đạc bản đồ.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu

  • Do điều kiện kỹ thuật và hạ tầng thiết bị trong nước, luận án chưa thể đánh giá toàn bộ các thông số liên quan đến mô hình cơ học - quang học phức tạp như độ méo ảnh, khoảng cách lấy mẫu góc, chất lượng định hướng tuyệt đối của vệ tinh và dải động bức xạ tuyệt đối.
  • Nghiên cứu chỉ tập trung thực nghiệm trên tập dữ liệu ảnh quang học mức 0 và mức 1A của vệ tinh VNREDSat-1, chưa mở rộng thực nghiệm đồng thời trên nhiều hệ cảm biến siêu phổ hoặc cảm biến radar khẩu độ tổng hợp (SAR).

Hướng nghiên cứu tiếp theo

  • Mở rộng quy trình đánh giá và chuẩn hóa thuật toán cho các nguồn dữ liệu vệ tinh độ phân giải siêu cao và ảnh siêu phổ mà Việt Nam chuẩn bị tiếp nhận và khai thác (như KOMPSAT-3).
  • Nghiên cứu hoàn thiện các giải pháp đánh giá MTF và SNR không phụ thuộc hoàn toàn vào bãi kiểm định nhân tạo, kết hợp thuật toán đánh giá tự động trên các đối tượng tự nhiên đặc thù trên lãnh thổ Việt Nam.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu đối với:

  • Cơ quan quản lý và vận hành vệ tinh viễn thám (như Viện Công nghệ vũ trụ, Trung tâm Điều khiển và Khai thác vệ tinh nhỏ, Cục Viễn thám Quốc gia): Cung cấp quy trình định lượng và phần mềm thực nghiệm để kiểm định định kỳ chất lượng dữ liệu ảnh mức 0 và 1A trước khi sản xuất các sản phẩm thương mại.
  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Kỹ thuật Trắc địa - Bản đồ, Viễn thám và Xử lý ảnh: Cung cấp cơ sở lý thuyết toán học - vật lý chặt chẽ về MTF, SNR, LSI, PSF, thuật toán Canny và phương pháp đánh giá chất lượng ảnh quang điện tử.
  • Cơ quan xây dựng tiêu chuẩn kỹ thuật (Bộ Tài nguyên và Môi trường, Bộ Khoa học và Công nghệ): Cung cấp luận cứ khoa học để hoàn thiện các thông tư, quy chuẩn kỹ thuật quốc gia về kiểm định chất lượng ảnh viễn thám quang học tại Việt Nam.

Câu hỏi thường gặp

1. Tại sao luận án lại lựa chọn hai thông số MTF và SNR thay vì các thông số chất lượng không gian và bức xạ khác?

Hệ thống vệ tinh có nhiều thông số chất lượng, tuy nhiên các thông số như độ méo ảnh, khoảng cách lấy mẫu góc, chất lượng định hướng đòi hỏi các khung mô phỏng quang học - cơ học phức tạp mà Việt Nam chưa có; thông số GSD đã được đảm bảo qua quy trình hiệu chỉnh hình học; dải động bức xạ thiếu mục tiêu trắng tuyệt đối. Trong khi đó, SNR (đại diện bức xạ) và MTF (đại diện không gian) là hai chỉ tiêu cốt lõi phản ánh rõ nhất độ sắc nét và mức độ nhiễu của ảnh, hoàn toàn có thể tính toán chính xác và chủ động tại Việt Nam thông qua bãi kiểm định nhân tạo tại Buôn Ma Thuột và các bãi kiểm định tự nhiên.

2. Thuật toán Canny đóng vai trò gì trong việc tính toán hàm truyền điều biến MTF?

Trong phương pháp bãi kiểm định dạng cạnh Heaviside, việc xác định chính xác vị trí và độ nghiêng của đường biên giữa vùng sáng và vùng tối là yếu tố quyết định để xây dựng hàm lan truyền cạnh (ESF). Thay vì sử dụng các thuật toán xấp xỉ tuyến tính cổ điển dễ bị sai số do nhiễu và răng cưa, luận án sử dụng toán tử tách biên Canny giúp tự động hóa quá trình nhận dạng biên cạnh chính xác trên ảnh vệ tinh độ phân giải cao, từ đó nâng cao độ chính xác khi tính đạo hàm ra LSF và biến đổi Fourier xác định MTF.

3. Tín hiệu tối (DS) và mức độ hồi đáp không đồng đều (PRNU) được xử lý như thế nào trong quy trình của luận án?

Bản thân thiết bị chụp ảnh khi hoạt động có phát sinh nhiệt, tạo ra dòng tín hiệu ngay cả khi không có ánh sáng vào cảm biến (tín hiệu tối DS); đồng thời các phần tử cảm biến trên mảng CCD có độ nhạy khác nhau (PRNU). Luận án sử dụng dữ liệu ảnh mức 0 chụp các vùng đồng nhất tự nhiên (biển sâu Đại Tây Dương cho tín hiệu tối và sa mạc Libya, Algeria cho PRNU) để tính toán các hệ số hiệu chỉnh, loại trừ sai số hệ thống thô này trước khi đưa dữ liệu mức 1A vào tính toán SNR và MTF.

4. Dữ liệu thực nghiệm được sử dụng để kiểm chứng trong luận án gồm những nguồn nào?

Dữ liệu thực nghiệm bao gồm: ảnh vệ tinh VNREDSat-1 của Việt Nam ở mức 0 và mức 1A; số liệu đo quang phổ phản xạ bề mặt thực địa tại bãi kiểm định Buôn Ma Thuột (Đắk Lắk) thu thập trong hai năm 2017 và 2018; dữ liệu ảnh chụp bãi kiểm định nhân tạo Salon de Provence (Pháp); cùng dữ liệu ảnh chụp các bãi kiểm định tự nhiên tại sa mạc Châu Phi và Đại Tây Dương.


Kết luận

Luận án "Nghiên cứu đánh giá chất lượng ảnh viễn thám quang học của Việt Nam" đã giải quyết thành công bài toán xây dựng cơ sở khoa học và quy trình công nghệ định lượng độc lập để đánh giá chất lượng dữ liệu ảnh viễn thám mức thô (mức 0 và mức 1A) thông qua hai thông số then chốt là SNR và MTF. Việc ứng dụng toán tử tách biên Canny kết hợp khai thác bãi kiểm định Buôn Ma Thuột đã tối ưu hóa độ chính xác của phép đo MTF trên ảnh vệ tinh độ phân giải cao trong điều kiện thực tế của Việt Nam. Kết quả thực nghiệm trên vệ tinh VNREDSat-1 không chỉ khẳng định hệ thống thu nhận ảnh vẫn vận hành ổn định theo tiêu chuẩn thiết kế mà còn tạo tiền đề công nghệ để chủ động kiểm định chất lượng cho các thế hệ vệ tinh viễn thám tiếp theo của đất nước.