Tổng quan nghiên cứu

Quá trình đô thị hóa tại Việt Nam trong hai thập kỷ qua diễn ra với tốc độ nhanh chóng, kéo theo sự chuyển dịch mạnh mẽ từ đất nông nghiệp sang đất công trình với tỷ lệ biến động ước tính khoảng 5% đến 8% mỗi năm tại các vùng kinh tế trọng điểm. Sự gia tăng mật độ xây dựng và mở rộng không gian đô thị đặt ra yêu cầu cấp thiết về việc theo dõi, giám sát và phân loại lớp phủ bề mặt đất đai. Tuy nhiên, các nguồn dữ liệu viễn thám vệ tinh thô như MODIS (độ phân giải 250m, 500m, 1000m), DMSP-OLS (1000m) và WorldPop (100m) khi thu thập ban đầu đều có sự sai lệch lớn về tỷ lệ không gian, chứa nhiều nhiễu khí quyển và sai số hình học do góc quét cảm biến.

Mục tiêu cốt lõi của nghiên cứu là xây dựng, đánh giá toàn diện các phương pháp tiền xử lý và tổng hợp dữ liệu không gian (Spatial Data Aggregation), nhằm chuẩn hóa toàn bộ các tập dữ liệu viễn thám đa nguồn về cùng một độ phân giải lưới 500m thống nhất. Nghiên cứu tập trung giải quyết bài toán phân loại lớp phủ đô thị trên quy mô toàn lãnh thổ Việt Nam, đồng thời tiến hành kiểm chứng chi tiết tại khu vực đô thị trọng điểm Thành phố Hồ Chí Minh dựa trên dữ liệu năm 2015 và khung phân loại theo Nghị định số 42/2009/NĐ-CP.

Ý nghĩa khoa học của luận văn thể hiện ở việc thiết lập quy trình chuẩn hóa dữ liệu viễn thám quang học và số liệu không gian với độ tin cậy cao, hạn chế tối đa suy hao thông tin phổ. Về mặt thực tiễn, mô hình mang lại giải pháp công nghệ giúp nâng chỉ số F1-Score phân loại lớp phủ đạt mức trên 84%, hỗ trợ đắc lực cho các cơ quan quản lý nhà nước trong công tác quy hoạch đô thị bền vững.

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 lý thuyết cảm thám quang học và nguyên lý tương tác bức xạ điện từ. Cảm biến MODIS (trên hai vệ tinh Terra và Aqua) ghi nhận năng lượng bức xạ phản xạ từ bề mặt Trái Đất qua 36 kênh phổ trong dải bước sóng từ 0,405 µm đến 14,385 µm, với chu kỳ chụp lặp lại từ 1 đến 2 ngày.

Mô hình phân loại cốt lõi áp dụng trong luận văn là thuật toán GLCNMO (Global Land Cover by National Mapping Organizations) phiên bản 2 mở rộng, được điều chỉnh linh hoạt theo các tiêu chuẩn đô thị tại Việt Nam. Ba khái niệm nền tảng bao gồm:

  • Chỉ số thực vật chuẩn hóa (NDVI) và chỉ số thực vật cải tiến (EVI): Phản ánh mật độ sinh khối và độ che phủ cây xanh dựa trên tương quan phản xạ giữa kênh đỏ và cận hồng ngoại.
  • Bề mặt không thấm nước (Impervious Surface Area - EstISA): Đo lường mật độ các công trình nhân tạo như bê tông, nhựa đường, mái nhà đặc trưng cho không gian xây dựng.
  • Tổng hợp dữ liệu không gian (Spatial Aggregation): Phân chia thành tổng hợp dữ liệu phân lớp (Categorical Aggregation) áp dụng cho dữ liệu định tính và tổng hợp dữ liệu số (Numerical Aggregation) cho các dữ liệu trường liên tục.

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

Nghiên cứu sử dụng nguồn dữ liệu đa nguồn phong phú thu thập trong giai đoạn 2010 - 2015, bao gồm:

  • Ảnh chỉ số thực vật MODIS MOD13Q1 phiên bản 6 (độ phân giải 250m, lưới 4800 x 4800 điểm ảnh, chu kỳ tổng hợp 16 ngày).
  • Bản đồ mật độ dân số WorldPop năm 2015 (độ phân giải 100m x 100m, hệ tọa độ WGS84) được lựa chọn thay thế cho bộ dữ liệu LandScan 1km nhờ độ chi tiết không gian vượt trội và tính mở của dữ liệu.
  • Ảnh vệ tinh ánh sáng ban đêm DMSP-OLS phiên bản 4 (độ phân giải 1km, năm 2013) đã được khử nhiễu ánh sáng đi lạc.
  • Ảnh bề mặt không thấm nước toàn cầu EstISA năm 2010 (độ phân giải 1km) và lớp mặt nạ nước MOD44W năm 2015 (độ phân giải 250m).

Quy trình phân tích dữ liệu thực hiện qua 6 bước khép kín:

  • Bước 1: Tiền xử lý hiệu chỉnh bức xạ, hiệu chỉnh hình học và đưa về hệ quy chiếu đồng nhất Sinusoidal.
  • Bước 2: Tái mẫu dữ liệu (Resampling) đưa tất cả 5 lớp dữ liệu về lưới đồng nhất 500m. Áp dụng phương pháp nội suy song tuyến (Bilinear Interpolation) cho dữ liệu ánh sáng DMSP-OLS và EstISA; phương pháp tính tổng (Sum) cho dữ liệu dân số WorldPop; thuật toán luật đa số (Majority Rule) cho lớp mặt nạ nước MOD44W; và thử nghiệm 4 thuật toán tổng hợp riêng biệt cho lớp NDVI gồm: luật đa số (MRB), ngẫu nhiên (RRB), cửa sổ dịch chuyển trọng số khoảng cách (PDW) và lấy giá trị điểm ảnh trung tâm.
  • Bước 3: Đánh giá chất lượng tổng hợp ảnh thông qua hai thước đo toán học định lượng: Chỉ số tương đồng cấu trúc (SSIM - Structural Similarity Index) và Tỷ số tín hiệu cực đại trên nhiễu (PSNR - Peak Signal-to-Noise Ratio).
  • Bước 4 và 5: Xây dựng tập dữ liệu huấn luyện, tập kiểm tra độc lập và tiến hành phân loại lớp phủ đô thị bằng thuật toán GLCNMO mở rộng với ngưỡng mật độ dân số tối thiểu 2000 người/km² và đơn vị diện tích đô thị tối thiểu 1 km².
  • Bước 6: Kiểm định độ chính xác bằng ma trận nhầm lẫn (Confusion Matrix) và chỉ số F1-Score.

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

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

Phân tích định lượng trên các tập dữ liệu tổng hợp ảnh MODIS MOD13Q1 từ độ phân giải 250m về 500m mang lại các kết quả nổi bật:

  • Về bảo toàn thông tin phổ và cấu trúc không gian: Phương pháp cửa sổ dịch chuyển trọng số khoảng cách (PDW) và phương pháp lấy điểm ảnh trung tâm (Central Pixel) đạt chỉ số tương đồng cấu trúc SSIM trung bình lên tới 0,94, cao hơn 21,5% so với phương pháp lấy mẫu ngẫu nhiên (RRB chỉ đạt SSIM khoảng 0,77).
  • Về kiểm soát độ nhiễu: Phương pháp PDW đạt tỷ số PSNR vượt mức 36,8 dB và sai số bình phương trung bình (MSE) giảm xuống dưới 12,4 đơn vị số học, cho thấy sự ổn định tuyệt đối trong việc duy trì dải phổ thực vật liên tục.
  • Về độ chính xác phân loại lớp phủ đô thị: Bản đồ phân loại đô thị Việt Nam năm 2015 khi sử dụng dữ liệu NDVI được tổng hợp bằng phương pháp PDW kết hợp thuật toán GLCNMO mở rộng đạt độ chính xác tổng thể 88,6% và chỉ số F1-Score đạt 0,84, vượt trội hơn 6,2% so với việc sử dụng phương pháp tổng hợp theo luật đa số thông thường (F1-Score đạt 0,79).

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phương pháp PDW và phương pháp lấy điểm ảnh trung tâm đạt hiệu suất cao là nhờ khả năng bảo tồn độ dốc phổ tại các khu vực ranh giới chuyển tiếp giữa đô thị và nông thôn. Tại các vùng ven đô ở Việt Nam, cấu trúc sử dụng đất rất phân mảnh với sự xen kẽ giữa nhà ở, vườn cây và kênh rạch nhỏ. Khi hạ độ phân giải từ 250m xuống 500m, phương pháp ngẫu nhiên (RRB) làm đứt gãy tính liên tục của cấu trúc địa vật, dẫn đến việc mất mát các mảng xanh đô thị.

So sánh với các nghiên cứu trước đây tại khu vực Thành phố Hồ Chí Minh chỉ dựa trên ảnh Landsat 30m đơn lẻ hoặc phân loại nhiệt độ bề mặt cục bộ, việc kết hợp đa nguồn giữa WorldPop, DMSP-OLS và MODIS mang lại khả năng bao quát toàn quốc với chu kỳ cập nhật liên tục. Tiêu chí mật độ dân số tối thiểu 2000 người/km² theo Nghị định 42/2009/NĐ-CP kết hợp với ngưỡng ánh sáng ban đêm giúp loại bỏ hoàn toàn hiện tượng phân loại nhầm trên 32% diện tích các cồn cát ven biển miền Trung và các vùng đất trống đồi trọc vốn có giá trị phản xạ phổ rất giống với bê tông đô thị.

Toàn bộ dữ liệu phân tích có thể được mô hình hóa trực quan thông qua:

  • Biểu đồ phân phối tần suất Histogram của chỉ số NDVI trước và sau khi tổng hợp: Cho thấy đường cong phân bố của phương pháp PDW gần như trùng khít hoàn hảo với ảnh gốc 250m ban đầu.
  • Bản đồ chuyên đề phân loại lớp phủ đô thị Thành phố Hồ Chí Minh năm 2015: Minh họa sự phân tách ranh giới rõ rệt giữa vùng lõi trung tâm có mật độ xây dựng cao với các vùng đệm nông thôn lân cận.

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

Dựa trên kết quả thực nghiệm, bốn nhóm giải pháp kỹ thuật và chính sách được đề xuất cụ thể như sau:

  • Tự động hóa quy trình tiền xử lý và tích hợp dữ liệu viễn thám đa nguồn: Xây dựng các module tiền xử lý ảnh vệ tinh tự động tích hợp thuật toán PDW và phép nội suy song tuyến, nhằm rút ngắn thời gian chuẩn hóa dữ liệu viễn thám từ 10 ngày xuống dưới 24 giờ cho mỗi chu kỳ quan sát. Chủ thể thực hiện: Trung tâm Vũ trụ Việt Nam và các viện nghiên cứu chuyên ngành; thời gian triển khai: giai đoạn 2026 - 2027.
  • Thiết lập cơ sở dữ liệu không gian lớp phủ đô thị quốc gia chuẩn hóa: Chuẩn hóa hệ thống dữ liệu viễn thám giám sát đô thị trên lưới không gian 100m và 500m thống nhất, đặt mục tiêu bao phủ 100% diện tích các đô thị loại đặc biệt đến loại III trên cả nước. Chủ thể thực hiện: Bộ Tài nguyên và Môi trường phối hợp cùng Bộ Xây dựng; thời gian hoàn thành: quý IV năm 2027.
  • Nâng cấp thuật toán phân loại tích hợp trí tuệ nhân tạo (AI/Deep Learning): Nghiên cứu tích hợp mạng nơ-ron tích chập (CNN) kết hợp cùng thuật toán GLCNMO mở rộng, nhằm đẩy chỉ số F1-Score phân loại lớp phủ tự động đạt mốc trên 92% trong các bài toán nhận diện hạ tầng ngầm và công trình cao tầng. Chủ thể thực hiện: Các nhóm nghiên cứu Hệ thống thông tin tại các trường đại học; thời gian nghiên cứu: giai đoạn 2027 - 2028.
  • Hoàn thiện khung kỹ thuật giám sát biến động đô thị theo luật định: Ban hành bộ chỉ dẫn kỹ thuật áp dụng chỉ số bề mặt không thấm nước EstISA và mật độ dân cư vệ tinh trong việc thẩm định quy hoạch đô thị theo Nghị định số 42/2009/NĐ-CP. Chủ thể thực hiện: Vụ Quy hoạch Kiến trúc - Bộ Xây dựng và Ủy ban nhân dân các tỉnh, thành phố; thời gian áp dụng: từ năm 2026.

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

Luận văn là tài liệu tham khảo giá trị cao cho bốn nhóm đối tượng chuyên môn:

  • Chuyên gia quy hoạch và phát triển đô thị: Khai thác phương pháp phân tích lớp phủ viễn thám để đánh giá tốc độ mở rộng không gian xây dựng, kiểm soát tỷ lệ đất không thấm nước và tối ưu hóa quy hoạch hạ tầng giao thông tại các đô thị loại I và loại II.
  • Học viên, nghiên cứu sinh ngành Hệ thống thông tin và Địa tin học (GIS/Remote Sensing): Nắm bắt trọn vẹn quy trình tiền xử lý ảnh thô, các phương pháp tái lấy mẫu (Resampling) và phương pháp luận đo lường độ suy giảm chất lượng dữ liệu bằng chỉ số SSIM, PSNR.
  • Kỹ sư phát triển phần mềm và hệ thống dữ liệu không gian: Ứng dụng các thuật toán tổng hợp dữ liệu số và phân lớp (PDW, MRB, Bilinear) vào việc phát triển các thư viện xử lý ảnh raster lớn trên các nền tảng điện toán đám mây.
  • Cơ quan quản lý tài nguyên, môi trường và khí tượng thủy văn: Sử dụng quy trình chiết xuất thông tin đa nguồn (MODIS, DMSP-OLS, WorldPop) để phục vụ công tác kiểm kê đất đai định kỳ 5 năm một lần và giám sát hiện tượng đảo nhiệt đô thị.

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

  • Tại sao cần phải chuẩn hóa các nguồn dữ liệu viễn thám về cùng độ phân giải 500m? Mỗi bộ cảm biến vệ tinh hoạt động ở một độ cao quỹ đạo và góc nhìn tức thời (IFOV) khác nhau, dẫn đến kích thước điểm ảnh chênh lệch lớn từ 100m (WorldPop), 250m (MOD13Q1) đến 1000m (DMSP-OLS, EstISA). Việc đồng nhất dữ liệu về độ phân giải 500m là bước bắt buộc để thực hiện các phép chồng lớp ma trận đại số không gian trong thuật toán GLCNMO mà không gây sai lệch vị trí điểm ảnh.

  • Chỉ số SSIM và PSNR mang ý nghĩa gì trong đánh giá chất lượng ảnh tổng hợp? SSIM đo lường mức độ bảo toàn độ sáng, độ tương phản và cấu trúc địa vật của ảnh sau khi tổng hợp so với ảnh gốc (giá trị càng tiệm cận 1,0 càng hoàn hảo). Trong khi đó, PSNR (tính bằng decibel - dB) thể hiện tỷ lệ tín hiệu hữu ích trên lượng nhiễu phát sinh; giá trị PSNR vượt trên 35 dB khẳng định thuật toán không làm biến dạng giá trị phản xạ phổ nguyên bản.

  • Lý do luận văn lựa chọn bộ dữ liệu WorldPop 2015 thay vì LandScan? Bộ dữ liệu LandScan có độ phân giải không gian 1km và thuộc quyền sở hữu có giới hạn truy cập từ cơ quan quản lý năng lượng quốc tế. Ngược lại, WorldPop năm 2015 cung cấp độ phân giải cao tới 100m x 100m hoàn toàn mở, sử dụng thuật toán phân bổ dân số chi tiết dựa trên dữ liệu điều tra thực tế tại Việt Nam, mang lại độ chính xác cao gấp 10 lần tại các khu vực dân cư nông thôn chuyển dịch.

  • Tiêu chuẩn xác định một "khu đô thị" trong bài toán phân loại này là gì? Theo thuật toán GLCNMO mở rộng áp dụng trong luận văn, một khu vực được định nghĩa là đô thị khi đáp ứng đồng thời 3 tiêu chí: mật độ dân số tối thiểu đạt 2000 người/km², tỷ lệ bề mặt không thấm nước EstISA và cường độ ánh sáng ban đêm DMSP-OLS vượt qua các ngưỡng quy định, đồng thời đơn vị diện tích bản đồ tối thiểu phải đạt từ 1 km² trở lên.

  • Dữ liệu ánh sáng ban đêm DMSP-OLS đóng vai trò gì trong việc phân loại? Dữ liệu DMSP-OLS năm 2013 thu nhận bức xạ ánh sáng nhân tạo ban đêm, giúp phân biệt chính xác giữa các bề mặt đất trống không có người ở với các khu dân cư xây dựng có hoạt động kinh tế. Việc kết hợp dữ liệu này giúp mô hình triệt tiêu hoàn toàn sai số phân loại nhầm giữa đất cát, đất hoang và bê tông đô thị.

Kết luận

  • Chuẩn hóa thành công quy trình tiền xử lý và tích hợp 5 nguồn dữ liệu viễn thám không gian (MODIS 250m, WorldPop 100m, DMSP-OLS 1km, EstISA 1km, MOD44W 250m) về lưới dữ liệu đồng nhất 500m trên phạm vi toàn quốc.
  • Chứng minh phương pháp cửa sổ dịch chuyển trọng số khoảng cách (PDW) đạt hiệu quả vượt trội trong bảo toàn cấu trúc phổ với SSIM đạt 0,94 và PSNR vượt 36,8 dB.
  • Xây dựng hoàn chỉnh bản đồ lớp phủ đô thị Việt Nam năm 2015 bằng thuật toán GLCNMO mở rộng với độ chính xác tổng thể 88,6% và F1-Score đạt 0,84.
  • Cung cấp luận cứ khoa học thực tiễn phục vụ theo dõi biến động diện tích xây dựng và phân loại đô thị theo Nghị định 42/2009/NĐ-CP.
  • Định hình lộ trình nâng cấp thuật toán ứng dụng học sâu (Deep Learning) và dữ liệu vệ tinh phân giải cao Sentinel/Landsat trong giai đoạn 2026 - 2030.

Nghiên cứu mở ra hướng tiếp cận chuẩn mực cho công tác xử lý ảnh viễn thám đa nguồn tại Việt Nam. Độc giả, nhà nghiên cứu và cơ quan quản lý quan tâm có thể tiếp tục khai thác khung phương pháp này để ứng dụng trong các bài toán quan trắc môi trường và phát triển đô thị thông minh trong kỷ nguyên chuyển đổi số.