Tổng quan nghiên cứu

Năm 2013, sự kiện phóng thành công vệ tinh viễn thám quang học đầu tiên của Việt Nam là VNREDSat-1 với khối lượng 120 kg và độ phân giải không gian 10 mét đã mở ra bước ngoặt mới cho công nghệ quan sát trái đất nước nhà. Tuy nhiên, việc khai thác hiệu quả nguồn dữ liệu ảnh viễn thám phục vụ giám sát tài nguyên thiên nhiên, đặc biệt là tài nguyên nước tại các đô thị phức tạp, vẫn gặp nhiều rào cản do phụ thuộc vào các phần mềm thương mại ngoại nhập đắt đỏ hoặc phương pháp giải đoán thủ công tốn nhiều nhân lực. Nhằm giải quyết triệt để vấn đề này, đề tài nghiên cứu đã xây dựng thuật toán tự động trích xuất thông tin nước chuyên biệt cho ảnh đa phổ VNREDSat-1 khu vực đô thị.

Nghiên cứu tập trung vào không gian thử nghiệm tại Thủ đô Hà Nội với phạm vi tọa độ địa lý từ 20°53' đến 21°23' vĩ độ Bắc và 105°44' đến 106°02' kinh độ Đông, tổng diện tích nghiên cứu đạt 349 km² cùng chu vi đường bao 74,8 km. Khu vực này sở hữu hệ thống thủy văn đa dạng gồm 163 km chiều dài dòng sông Hồng chảy qua địa bàn thành phố, cùng mạng lưới hồ đầm dày đặc với điển hình là Hồ Tây có diện tích mặt nước lên tới 526 ha. Mục tiêu trọng tâm của luận văn là tự động hóa quy trình bóc tách các vùng nước mặt, phân biệt chính xác nước sông giàu phù sa, nước hồ lục địa và loại trừ triệt để bóng đổ của các khối nhà cao tầng. Về mặt giá trị ứng dụng, giải pháp giúp rút ngắn 75% thời gian xử lý dữ liệu ảnh viễn thám, nâng cao độ tin cậy của bản đồ hiện trạng nước mặt tỷ lệ 1:10.000, đóng góp thiết thực cho công tác quản lý tài nguyên và quy hoạch thoát nước đô thị.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự tích hợp của ba nền tảng lý thuyết cốt lõi trong lĩnh vực viễn thám và thị giác máy tính: Lý thuyết phản xạ phổ điện từ (Spectral Reflectance Theory), Mô hình phân lớp cây quyết định (Decision Tree Classification) và Lý thuyết phân tích hình thái không gian (Spatial-morphological Analysis). Hệ thống khái niệm trọng tâm bao gồm:

  • Độ phân giải phổ: Khả năng thu nhận tín hiệu sóng điện từ của bộ cảm quang học NAOMI trên 4 kênh phổ độc lập, bao gồm Kênh 1 Blue (0,45 - 0,52 µm), Kênh 2 Green (0,53 - 0,60 µm), Kênh 3 Red (0,62 - 0,69 µm) và Kênh 4 Cận hồng ngoại NIR (0,76 - 0,89 µm).
  • Độ phân giải không gian: Kích thước chi tiết của mỗi điểm ảnh trên mặt đất, đạt 10 mét đối với ảnh đa phổ và 2,5 mét đối với ảnh toàn sắc Panchromatic.
  • Chỉ số hình dạng hình học M: Tỷ số tương quan giữa diện tích và bình phương chu vi của cụm điểm ảnh nhằm nhận diện đặc trưng hình học của đối tượng.
  • Mô hình chuyển đổi Raster sang Vector: Kỹ thuật chuyển đổi cấu trúc ma trận điểm ảnh sang định dạng không gian vector gắn hệ tọa độ địa lý chuẩn.

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

Nguồn dữ liệu thực nghiệm sử dụng cảnh ảnh đa phổ VNREDSat-1 chụp ngày 22 tháng 5 năm 2013 tại khu vực Hà Nội. Cảnh ảnh gốc có kích thước 1750 x 1750 pixel với dung lượng 23,57 MB được nắn chỉnh hình học lên kích thước 2174 x 2145 pixel với dung lượng 36,42 MB tích hợp mô hình số hóa độ cao DEM. Nghiên cứu thực hiện chọn mẫu có kiểm soát trên 8 nhóm đối tượng bề mặt chính: nước Sông Hồng, nước lục địa, nhà ở, đường giao thông, thực vật, đất trống, mây và bóng đổ tòa nhà. Cỡ mẫu cho mỗi đối tượng được thu thập độc lập từ 30 đến 100 mẫu cục bộ không chồng lấn nhằm đảm bảo tính đại diện thống kê cao.

Phương pháp phân tích dựa trên mô hình vật lý chuyên dụng Imogen chạy trên phần mềm ERDAS IMAGINE 2014 để hiệu chỉnh hình học mức 2A về hệ tọa độ chuẩn Quốc gia VN-2000. Quy trình bố trí 13 điểm kiểm soát mặt đất GCP đạt sai số trung phương bình sai RMSE chỉ 1,01175 mét và sai số kiểm tra đạt 1,75914 mét, đáp ứng trọn vẹn tiêu chuẩn bản đồ tỷ lệ 1:10.000. Lý do lựa chọn cây quyết định kết hợp giải thuật tìm kiếm theo chiều rộng BFS 4 hướng kế cận là khả năng phân nhánh dữ liệu trực tiếp theo ngưỡng phổ, tính toán hình thái học nhanh chóng và dễ dàng tích hợp thành modul độc lập. Toàn bộ quy trình nghiên cứu, xây dựng mã nguồn Python và đối soát thực địa được hoàn thành liên tục trong giai đoạn 2013 - 2015.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, phân tích phổ phản xạ khẳng định kênh cận hồng ngoại NIR (Band 4: 0,76 - 0,89 µm) là dải phổ then chốt để cô lập nước. Giá trị bức xạ trung bình của nước lục địa tại kênh 4 chỉ đạt 63,8 và nước Sông Hồng là 81,4, trong khi đối tượng liền kề thấp nhất là đường giao thông đạt 201,4. Ngưỡng phân tách d1 = 141 được thiết lập giúp bóc tách hơn 98% đối tượng nước và bóng đổ ra khỏi nền bề mặt đô thị.

Thứ hai, nghiên cứu phát hiện kênh Blue (Band 1: 0,45 - 0,52 µm) cho phép phân tách tuyệt đối giữa nước Sông Hồng (giá trị trung bình 115,6) và nước hồ lục địa (giá trị trung bình 78,4) nhờ ngưỡng phổ d2 = 97, giải quyết triệt để bài toán nhận dạng nước có độ đục và hàm lượng phù sa cao.

Thứ ba, việc áp dụng ngưỡng diện tích giới hạn d3 = 22.000 m² kết hợp chỉ số hình M nằm trong dải từ 0,13 đến 0,23 và ngưỡng phổ phụ NIR > 108 đã loại bỏ chính xác trên 92% bóng đổ của các tòa nhà cao tầng có dạng hình chữ nhật hoặc dải hẹp song song.

Thứ tư, thuật toán được mở rộng thành công để trích xuất tự động thảm thực vật với ngưỡng kết hợp NIR > 292, Red < 144 và mây với ngưỡng Red > 426. Dữ liệu diện tích mặt nước sau khi bóc tách được biểu diễn qua bảng ma trận sai số và bản đồ vector đối sánh với diện tích Hồ Tây 526 ha và Hồ Gươm cho thấy độ tương đồng đạt trên 95% so với số liệu địa chính thực tế.

Thảo luận kết quả

Nguyên nhân của sự phân tách vượt trội này xuất phát từ đặc tính hấp thụ quang phổ cực mạnh của phân tử nước ở dải sóng cận hồng ngoại và sự tán xạ mạnh của các hạt vật chất lơ lửng trong nước sông phù sa ở dải sóng xanh lam. Mặt khác, bóng đổ công trình tuy có giá trị phổ hấp thụ thấp gần giống nước ở các kênh nhìn thấy nhưng lại có cấu trúc biên thẳng song song và kích thước diện tích nhỏ hơn các hồ thủy vực tự nhiên.

Khi so sánh với mô hình của các công trình quốc tế áp dụng trên vệ tinh SPOT5, thuật toán trên ảnh VNREDSat-1 thể hiện ưu thế rõ nét: độ phân giải kênh cận hồng ngoại của VNREDSat-1 đạt 10 mét đồng nhất (vượt trội hơn kênh SWIR 20 mét của SPOT5), đồng thời khoảng cách chênh lệch giá trị số giữa nước và nền đô thị rộng hơn tới 35%, giúp hạn chế hiện tượng nhầm lẫn điểm ảnh biên hỗn hợp. Kết quả phân tích có thể được tổng hợp trực quan qua các biểu đồ đường phổ đa kênh và ma trận phân loại điểm ảnh, mang lại độ tin cậy vượt bậc cho công nghệ xử lý ảnh viễn thám nội địa.

Đề xuất và khuyến nghị

  • Tích hợp module tự động vào trạm thu viễn thám quốc gia: Cục Viễn thám Quốc gia trực thuộc Bộ Tài nguyên và Môi trường cần đưa thuật toán cây quyết định vào quy trình xử lý tự động cấp độ 2A. Mục tiêu nhằm tự động hóa 95% công tác trích xuất mặt nước trong thời gian dưới 15 phút cho mỗi cảnh ảnh, triển khai trong lộ trình 12 tháng.
  • Xây dựng ngân hàng mẫu phổ chuẩn hóa đa thời gian: Viện Hàn lâm Khoa học và Công nghệ Việt Nam phối hợp cùng các trường đại học thu thập trên 5.000 mẫu phổ chuẩn cho hệ thống sông hồ tại 63 tỉnh thành phố. Chương trình cần thực hiện trong vòng 24 tháng nhằm thích ứng với sự thay đổi của thời tiết và các mùa mưa lũ.
  • Nâng cấp thuật toán kết hợp học máy và mô hình DEM độ phân giải cao: Nhóm nghiên cứu phát triển phần mềm cần bổ sung các mạng học sâu (Deep Learning) kết hợp dữ liệu độ cao địa hình chi tiết 2,5 mét. Mục tiêu giảm sai số nhận dạng tại các khu đô thị có mật độ cao tầng dày đặc xuống dưới 2% trong vòng 18 tháng tới.
  • Chuẩn hóa quy trình liên thông dữ liệu GIS phục vụ quản lý thoát nước: Chi cục Thủy lợi và các Công ty Thoát nước đô thị định kỳ hàng quý trích xuất bản đồ vector VN-2000 tỷ lệ 1:10.000 từ ảnh VNREDSat-1 để theo dõi dung tích trữ nước hồ đô thị và phát hiện kịp thời các hành vi lấn chiếm hành lang thoát lũ.

Đối tượng nên tham khảo luận văn

  • Kỹ sư và Chuyên gia Viễn thám - GIS: Nắm vững quy trình nắn chỉnh hình học ảnh VNREDSat-1 mức 2A bằng mô hình vật lý Imogen, kỹ thuật xử lý 13 điểm GCP đạt độ chính xác dưới 2 mét và quy chuẩn xuất bản đồ vector VN-2000.
  • Nhà nghiên cứu và Giảng viên Công nghệ Thông tin: Khai thác giải thuật cây quyết định tối ưu hóa theo đặc trưng phổ, kết hợp thuật toán duyệt đồ thị BFS 4 hướng để phân cụm không gian trên ma trận ảnh lớn 2174 x 2145 điểm ảnh.
  • Cán bộ Quản lý Tài nguyên Nước và Môi trường: Tiếp cận công cụ giám sát biến động diện tích mặt nước tự động cho các hồ lớn như Hồ Tây (526 ha) và 163 km sông Hồng, hỗ trợ điều hành cảnh báo lũ lụt và quy hoạch đô thị bền vững.
  • Học viên Cao học và Sinh viên ngành Kỹ thuật Phần mềm: Tài liệu học thuật chuẩn mực về quy trình lấy mẫu thống kê 30 đến 100 mẫu, kỹ thuật xử lý ảnh đa phổ 16-bit và phương pháp đối soát thực địa khoa học.

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

Thuật toán phân biệt nước lục địa và bóng đổ tòa nhà cao tầng bằng cách nào?
Thuật toán kết hợp đồng thời hai cơ chế: lọc theo ngưỡng diện tích d3 = 22.000 m² và phân tích chỉ số hình M = A / P². Các khối bóng đổ nhà thường có cạnh song song với chỉ số M từ 0,13 đến 0,23 kèm giá trị phổ NIR trung bình lớn hơn 108, giúp loại bỏ chính xác trên 92% bóng đổ công trình.

Tại sao kênh cận hồng ngoại NIR (Band 4) lại được chọn làm ngưỡng lọc đầu tiên?
Dải sóng cận hồng ngoại từ 0,76 đến 0,89 µm có đặc tính vật lý bị nước hấp thụ mạnh với giá trị phản xạ chỉ từ 63,8 đến 81,4. Mức chênh lệch lớn so với các bề mặt nhân tạo như đường nhựa (201,4) hay thảm thực vật (398,2) tạo ra khoảng cách phân tách an toàn d1 = 141, giúp loại bỏ 98% đối tượng không phải nước ngay từ bước đầu.

Độ chính xác nắn chỉnh hình học ảnh VNREDSat-1 trong nghiên cứu đạt mức nào?
Sử dụng mô hình vật lý Imogen trên nền tảng ERDAS IMAGINE 2014 cùng 13 điểm kiểm soát GCP chuẩn, kết quả bình sai khối ảnh đạt sai số trung phương RMSE mặt đất là 1,01175 mét và sai số kiểm tra đạt 1,75914 mét. Độ chính xác vị trí này đáp ứng hoàn hảo yêu cầu thành lập bản đồ địa hình tỷ lệ 1:10.000.

Thuật toán có phân biệt được nước sông phù sa đục với nước hồ đô thị không?
Có, thuật toán phân tách hiệu quả nhờ kiểm tra kênh Blue (Band 1: 0,45 - 0,52 µm). Phù sa trong nước Sông Hồng phản xạ mạnh khiến giá trị phổ đạt 115,6, cao hơn nhiều so với mức 78,4 của nước hồ lục địa. Việc áp dụng ngưỡng d2 = 97 giúp phân lập riêng biệt dòng chảy 163 km sông Hồng với độ chính xác cao.

Việc chuyển đổi dữ liệu kết quả từ dạng Raster sang Vector mang lại lợi ích gì?
Dữ liệu Raster dạng lưới pixel được vector hóa và gán trực tiếp hệ tọa độ Quốc gia VN-2000. Điều này giúp các cơ quan quản lý dễ dàng đo đạc trực tiếp chu vi, diện tích từng hồ chứa, tích hợp vào các hệ thống bản đồ số GIS chuyên ngành và đối chiếu tức thời với số liệu địa chính thực địa.

Kết luận

  • Xây dựng thành công thuật toán cây quyết định tự động trích xuất thông tin nước chuyên biệt cho ảnh vệ tinh đa phổ VNREDSat-1 đầu tiên của Việt Nam.
  • Thiết lập hệ thống 5 ngưỡng phổ và hình thái học tối ưu gồm d1 = 141, d2 = 97, d3 = 22.000 m², d4 = 0,13 và d5 = 108 giúp loại trừ triệt để bóng đổ đô thị.
  • Ứng dụng quy trình nắn chỉnh hình học Imogen mức 2A đạt sai số bình sai RMSE 1,01 mét, đảm bảo độ chính xác cho bản đồ tỷ lệ 1:10.000.
  • Ứng dụng giải thuật duyệt đồ thị BFS 4 hướng tính toán tự động diện tích, chu vi và vector hóa dữ liệu chính xác về hệ quy chiếu VN-2000.
  • Mở rộng năng lực trích xuất đa đối tượng gồm 163 km dòng chảy sông Hồng, thảm thực vật và độ phủ mây với độ tin cậy thực nghiệm trên 95%.

Đóng góp lớn nhất của luận văn là làm chủ giải pháp công nghệ tự động hóa xử lý ảnh viễn thám độ phân giải cao 10 mét, khẳng định năng lực tự chủ công nghệ không gian của Việt Nam. Trong lộ trình 2026 - 2028, mô hình cần tiếp tục được mở rộng sang các hệ thống học sâu và dữ liệu viễn thám siêu phổ. Hãy tham khảo ngay toàn văn công trình nghiên cứu để ứng dụng quy trình xử lý viễn thám tiên tiến này vào các dự án GIS và quản lý tài nguyên môi trường thực tế!