Tổng quan luận án

Luận án tiến sĩ chuyên ngành Hệ thống thông tin với đề tài "Nghiên cứu phân loại và xây dựng bộ dữ liệu lớp phủ đô thị tại Việt Nam" do nghiên cứu sinh Phạm Tuấn Dũng thực hiện tại Trường Đại học Công nghệ, Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của PGS. Doãn Minh Chung và TS. Bùi Quang Hưng (bảo vệ năm 2023).

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

Quá trình đô thị hóa tại Việt Nam diễn ra nhanh chóng cùng với sự tăng trưởng kinh tế và biến động sử dụng đất (tỷ lệ dân số đô thị chiếm khoảng 36,63% dân số cả nước, tốc độ tăng trưởng đô thị hóa đạt khoảng 2,9%/năm). Việc theo dõi biến động lớp phủ đô thị trên diện rộng và theo chuỗi thời gian dài là yêu cầu cấp thiết phục vụ quản lý quy hoạch, phát triển bền vững và đánh giá tác động môi trường.

Tuy nhiên, các nghiên cứu viễn thám về lớp phủ đô thị tại Việt Nam phần lớn mới triển khai ở quy mô cục bộ (từng thành phố hoặc vùng kinh tế), chỉ có duy nhất một công trình công bố trên phạm vi toàn quốc vào năm 2006 (Saksena và cộng sự) dựa trên dữ liệu thống kê dân số, nông nghiệp kết hợp viễn thám thô. Việc áp dụng trực tiếp các bộ dữ liệu lớp phủ toàn cầu vào Việt Nam bộc lộ nhiều hạn chế: độ chính xác suy giảm do thiếu dữ liệu kiểm chứng mặt đất đặc thù; độ phân giải không gian chưa tối ưu; tín hiệu quang phổ bị ảnh hưởng bởi mây bao phủ và tính biến thiên theo mùa của khí hậu nhiệt đới gió mùa; đồng thời các nguồn dữ liệu đa nguồn chưa được xử lý đồng bộ về độ phân giải và ngưỡng phân tách.

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

Luận án xác định 3 mục tiêu nghiên cứu cụ thể:

  1. Nghiên cứu, đánh giá các phương pháp biến đổi độ phân giải đối với dữ liệu đa nguồn trong bài toán phân loại lớp phủ đô thị tại Việt Nam.
  2. Nghiên cứu và cải tiến phương pháp phân loại lớp phủ đô thị của GLCNMO cho khu vực Việt Nam trên cơ sở lựa chọn dữ liệu và tính toán các ngưỡng phù hợp (bao gồm quy trình tính toán chỉ số theo chu kỳ 23 kỳ/năm của dữ liệu NDVI và thuật toán tách ngưỡng tự động).
  3. Xây dựng bộ dữ liệu lớp phủ đô thị Việt Nam qua các mốc thời gian và ứng dụng bộ dữ liệu này trong đánh giá sự ảnh hưởng của quá trình phát triển đô thị tới vấn đề ô nhiễm bụi mịn PM2.5 trong không khí tại Việt Nam giai đoạn 2004–2015.

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

  • Đối tượng nghiên cứu: Lớp phủ đô thị, các phương pháp biến đổi độ phân giải dữ liệu raster đa nguồn, thuật toán phân loại lớp phủ thích nghi theo ngưỡng và mô hình tương quan giữa mở rộng không gian đô thị với nồng độ bụi mịn PM2.5.
  • Phạm vi không gian: Toàn bộ lãnh thổ đất liền Việt Nam; đi sâu phân tích chi tiết tại các đô thị lớn (Hà Nội, Thành phố Hồ Chí Minh).
  • Phạm vi thời gian và chuỗi dữ liệu: Dữ liệu viễn thám và không gian đa nguồn trong giai đoạn 2004–2015, tập trung vào 4 mốc thời gian chính: năm 2004, 2008, 2012 và 2015.
  • Nguồn dữ liệu sử dụng: Dữ liệu lưới điểm ảnh mật độ dân số Worldpop (100m), dữ liệu bề mặt xây dựng không thấm nước EstISA (1km), dữ liệu ánh sáng ban đêm DMSP-OLS Version 4 (1km), dữ liệu chỉ số thực vật MOD13Q1 (250m, 16 ngày/kỳ), dữ liệu bề mặt nước MOD44W (250m) và dữ liệu ranh giới hành chính dạng vector GADM.

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

Tác giả hệ thống hóa các hướng nghiên cứu phân loại lớp phủ mặt đất và lớp phủ đô thị trên thế giới qua hai nhánh tiếp cận chính:

  1. Nhánh xây dựng các bộ dữ liệu lớp phủ toàn cầu: Các công trình kinh điển từ dữ liệu AVHRR như GLCC của USGS (1993, 17 lớp, độ chính xác 66,9%), UMd của Đại học Maryland (Hansen và Reed, 2000, 14 lớp, độ chính xác 69%); GLC2000 của EC (Bartholomé và Belward, 2005, 23 lớp, độ chính xác 68,6%); MODIS LC của Đại học Boston (Friedl và cộng sự, 2002, độ phân giải 500m/1km, 17 lớp, độ chính xác 75%); GlobCover của ESA (Defourny và cộng sự, 2009, độ phân giải 300m, 22 lớp, độ chính xác 67,1%); GLCNMO các phiên bản 2003, 2008, 2013 của Dự án Bản đồ Toàn cầu (Tateishi và cộng sự, 2011, 2014; Kobayashi và cộng sự, 2017, độ phân giải 500m, 20 lớp, độ chính xác 77,9%); MCD12Q1 của NASA (Sulla-Menashe và cộng sự, 2019, 500m); FROM-GLC của Đại học Thanh Hoa (Gong và cộng sự, 2013, độ phân giải 30m, 10 lớp, độ chính xác 64,9%); GlobLand30 (Chen và cộng sự, 2014, 30m, 10 lớp, độ chính xác 80%); CGLS-LC100 của Copernicus (Buchhorn và cộng sự, 2020, 100m, 20 lớp, độ chính xác 80,2%); iMap World 1.0 (Liu và cộng sự, 2021, 30m, 33 lớp, độ chính xác 73%); và các bộ dữ liệu độ phân giải 10m từ Sentinel năm 2020 gồm Dynamic World (Google, 2022), WorldCover 2020 (ESA, Zanaga và cộng sự, 2021, độ chính xác 74,4%) và Land Cover 2020 (Esri, Karra và cộng sự, 2021, độ chính xác 63%).

  2. Nhánh trích xuất đô thị bằng các chỉ số quang phổ và học máy: Các chỉ số chuyên biệt cho bề mặt xây dựng như NDBI (Zha và cộng sự, 2003; He và cộng sự, 2010), IBI (Xu, 2008), BUI (Kavzoglu và Tonbul, 2018), NBI (Jieli Chen và cộng sự, 2010), BAEI (Bouzekri và cộng sự, 2015), EBBI (As-syakur và cộng sự, 2012), NDCCI (Naili và cộng sự, 2015). Các nghiên cứu cấp quốc gia tiêu biểu gồm nghiên cứu quá trình đô thị hóa tại Hoa Kỳ giai đoạn 1990–2010 (Balk và cộng sự, 2018); phát triển đô thị tại Trung Quốc (Zhang và cộng sự, 2018; Gong và cộng sự, 2019 với độ chính xác bề mặt không thấm nước trên 90%); nghiên cứu tại Nhật Bản (Bagan và Yamagata, 2012) và Ấn Độ (Gorsevski và cộng sự, 2020 với SVM, CART, RF đạt độ chính xác 87%).

Các nghiên cứu tại Việt Nam

  • Quy mô toàn quốc: Hoàng Thanh Tùng và cộng sự (2020) phân loại 12 lớp phủ bằng ảnh radar và quang học (JAXA) với phương pháp KDE đạt độ chính xác 85,6%; Meyfroidt và Lambin (2008) phân tích biến động rừng thập niên 1990; Saksena và cộng sự (2014, công bố nghiên cứu thực hiện năm 2006) phân loại 4 vùng chuyển đổi đô thị (nông thôn, ngoại ô, đô thị, lõi đô thị) cấp xã/phường đạt độ chính xác 90%.
  • Khu vực Hà Nội: Zhao và Liu (2010, SVM, 86,25%); Mẫn Đức Chức và cộng sự (2018, Landsat 8 ensemble learning, 84%); Nguyễn Đình Dương (2016, phân tích phổ tự động); Lasko và cộng sự (2018, Sentinel-1A RF phân loại lúa >90%); Bùi Bá Nam và cộng sự (2019, Sentinel 1/2 XGBoost trên GEE, 86%); Bùi Quang Thành và cộng sự (2019, SPOT-7 object-based, 89%; 2016, SPOT-7 kết hợp GOA và ANN, 87,19%); Phạm Minh Hải và Yamaguchi (2011, Landsat/ASTER 1975–2003, MLC); Ha và Weng (2016, Landsat bề mặt không thấm nước DT, 74%–82%); Trần Đức Thiên và cộng sự (2019, ngưỡng NDVI); Nguyễn Thị Mai và cộng sự (2019, MLC, 90,75% và phân tích UHI).
  • Khu vực TP. Hồ Chí Minh và các tỉnh thành khác: Schaefer và Nguyễn Xuân Thịnh (2019, Landsat/SPOT MLC, 87,01%); Nguyễn Vân Anh và cộng sự (2021, Sentinel-2 vs VNREDSat-1, 83% và 86%); Goldblatt và cộng sự (2018, Landsat/Sentinel RF trên GEE, 80%); Vũ Tường Thụy và cộng sự (2018, SVM, 82,96%); Trần Thị Vân và cộng sự (2011, UHI từ NDVI/ASTER); Nguyễn Thanh Sơn và cộng sự (2017, Landsat NDVI/NDCI, >86,1%); Phan Thành Nội và Kappas (2018, Đồng bằng sông Hồng, Sentinel-2 RF/kNN/SVM, 90%–95%); Nguyễn Thị Thanh Hương và cộng sự (2018, Đắk Lắk Landsat 8 RF, 90,32%); Hà Văn Tuyên và cộng sự (2018, Thái Nguyên Landsat/DEM RF, 93,4%); Phạm Thị Mai Thy và cộng sự (2018, Cần Thơ Landsat/ASTER SVM 71,38%, DT 65,86%; 2021, hạn hán miền Trung); Bùi Quang Ngọc và cộng sự (2017, LiDAR Quảng Xương, 97%); Phan Văn Thọ và Matos (2022, Quy Nhơn Landsat SVM, 77%).

Khoảng trống nghiên cứu được lựa chọn

Luận án tập trung giải quyết khoảng trống: thiếu một phương pháp phân loại lớp phủ đô thị chuyên biệt cho toàn bộ lãnh thổ Việt Nam có khả năng tích hợp hiệu quả dữ liệu viễn thám đa nguồn có độ phân giải khác nhau, tự động hóa việc xác định ngưỡng phân loại thích nghi với đặc thù thổ nhưỡng và khí hậu Việt Nam, xây dựng bộ dữ liệu lịch sử nhiều mốc thời gian để làm cơ sở định lượng tác động của đô thị hóa tới ô nhiễm không khí (bụi mịn PM2.5).

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

Hệ thống phân loại và định nghĩa sử dụng

Luận án kế thừa hệ thống phân loại lớp phủ quốc tế LCCS (FAO) và CORINE (Ủy ban Châu Âu), đồng thời phân tích định nghĩa đô thị của Liên Hợp Quốc (mật độ >= 1.500 người/km2 cho thành phố; >= 300 người/km2 cho đô thị vệ tinh) và tiêu chuẩn phân loại đô thị tại Việt Nam (Nghị định số 42/2009/NĐ-CP và Nghị quyết số 1210/2016/UBTVQH13).

Trên cơ sở đó, luận án đưa ra định nghĩa tường minh về lớp phủ đô thị:

Định nghĩa 1.2: "Lớp phủ đô thị là các khu vực trên mặt đất có mật độ dân cư lớn (mật độ dân số tính trên diện tích đất xây dựng đô thị đạt từ 1.500 người/km2 trở lên), các công trình xây dựng chiếm tỷ lệ chủ yếu (trên 50% tổng diện tích khu vực), không chứa các diện tích cây xanh, bề mặt nước, hoặc ít ánh sáng ban đêm và có diện tích tối thiểu là 1 km2."

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

Luận án kết hợp các nhóm phương pháp:

  • Phương pháp kế thừa và cải tiến thuật toán: Kế thừa khung phân loại lớp phủ đô thị toàn cầu của GLCNMO (Dự án Bản đồ Toàn cầu) và cải tiến thuật toán trích xuất ngưỡng cho phù hợp với dữ liệu tại Việt Nam.
  • Phương pháp xử lý ảnh số và biến đổi độ phân giải (Resampling):
    • Nhóm phương pháp hạ độ phân giải (Downsampling/Aggregation): Luật đa số (Majority), lấy giá trị ngẫu nhiên (Random), lấy giá trị điểm trung tâm (Center), giá trị trung bình (Mean), trung bình dựa trên trọng số (Weighted Mean), trung vị (Median), giá trị lớn nhất (Maximum), giá trị nhỏ nhất (Minimum), tính tổng giá trị (Sum).
    • Nhóm phương pháp tăng độ phân giải (Upsampling/Interpolation): Nội suy láng giềng gần nhất (Nearest Neighbor), nội suy song tuyến tính (Bilinear Interpolation), nội suy xoắn bậc ba (Cubic Convolution).
  • Phương pháp đo lường chất lượng dữ liệu biến đổi: Sử dụng ba chỉ số toán học:
    • Sai số bình phương trung bình ($MSE$)
    • Tỷ số tín hiệu cực đại trên nhiễu ($PSNR$)
    • Chỉ số tương đồng cấu trúc ($SSIM$)
  • Phương pháp thống kê và tách ngưỡng tự động: Phân tích biểu đồ tần suất (Histogram) từ tập mẫu huấn luyện để xác định ngưỡng phân tách tối ưu giữa lớp đô thị và các lớp không phải đô thị.
  • Phương pháp chuyên gia và kiểm chứng thực địa: Kết hợp kiểm tra chéo các điểm mẫu với ảnh độ phân giải cao trên Google Earth và ý kiến chuyên gia viễn thám.
  • Phương pháp tương quan không gian: Mô hình hóa mối tương quan giữa diện tích mở rộng đô thị và phân bố nồng độ bụi mịn PM2.5 (ước lượng từ độ dày quang học sol khí MODIS AOD/AOT kết hợp dữ liệu quan trắc trạm mặt đất).

Bảng tổng hợp các nguồn dữ liệu đầu vào

Tên dữ liệu Mô tả chi tiết Độ phân giải gốc Mốc thời gian Nguồn cung cấp
Worldpop Dữ liệu phân bố và mật độ dân số Việt Nam (dựa trên Tổng điều tra dân số 2009) 100 m 2004, 2008, 2012, 2015 WorldPop Hub
DMSP-OLS Dữ liệu ánh sáng ban đêm Version 4 (ảnh tổng hợp trung bình năm đã khử mây) 1 km 2004, 2008, 2012, 2013 NOAA / NGDC
MOD13Q1 Chuỗi dữ liệu chỉ số thực vật NDVI chu kỳ 16 ngày (23 chu kỳ/năm, các cảnh h27v06, h27v07, h28v06, h28v07, h28v08) 250 m 2004, 2008, 2012, 2015 NASA LP DAAC
EstISA Dữ liệu ước tính mật độ bề mặt xây dựng không thấm nước toàn cầu 1 km 2001, 2010 NOAA / NGDC
MOD44W Bản đồ mặt nước tĩnh MODIS 250 m 2004, 2012, 2015 NASA LP DAAC
GADM Dữ liệu ranh giới địa giới hành chính Việt Nam dạng vector Vector 2015 GADM

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

CHƯƠNG 1: TỔNG QUAN VỀ LỚP PHỦ VÀ CÁC PHƯƠNG PHÁP PHÂN LOẠI LỚP PHỦ

Chương 1 thiết lập nền tảng lý luận toàn diện về viễn thám và phân loại lớp phủ:

  • Phân tích chi tiết các hệ thống phân loại chuẩn quốc tế: LCCS (FAO) với 3 mức chi tiết phân biệt bề mặt thực vật và phi thực vật; hệ thống CORINE với 3 cấp độ phân loại từ chuỗi đô thị đến các mảng không gian xanh và hạ tầng.
  • Làm rõ bản chất vật lý của dữ liệu viễn thám: nguyên lý thu nhận năng lượng phản xạ/bức xạ điện từ; phân loại cảm biến chủ động (SAR, LiDAR) và cảm biến thụ động (quang học đa phổ, nhiệt); 4 tham số độ phân giải cơ bản (độ phân giải không gian, độ phân giải quang phổ, độ phân giải bức xạ, độ phân giải thời gian).
  • Khảo sát các kỹ thuật biến đổi độ phân giải ảnh raster (Downsampling và Upsampling) và chỉ ra rằng việc chuyển đổi tỷ lệ không gian giữa các bộ dữ liệu khác nguồn luôn gây ra sai lệch thông tin và mất mát chi tiết hình học.
  • Đặt ra bài toán phân loại dữ liệu đa nguồn: sự kết hợp giữa các lớp dữ liệu vật lý (quang phổ NDVI, mặt nước), dữ liệu hoạt động kinh tế - xã hội (ánh sáng ban đêm DMSP-OLS, bề mặt không thấm nước EstISA) và dữ liệu xã hội học (mật độ dân số Worldpop).

CHƯƠNG 2: NGHIÊN CỨU PHƯƠNG PHÁP PHÂN LOẠI LỚP PHỦ ĐÔ THỊ TẠI VIỆT NAM

Chương 2 là trọng tâm phương pháp luận của luận án, tập trung cải tiến mô hình phân loại GLCNMO cho phù hợp với Việt Nam:

[Dữ liệu đa nguồn: Worldpop (100m), DMSP-OLS (1km), MOD13Q1 (250m), EstISA (1km), MOD44W (250m)]
             [Tiền xử lý & Biến đổi về cùng độ phân giải chuẩn (500m)]
                  (Đánh giá kiểm soát qua chỉ số MSE, PSNR, SSIM)
               [Lấy mẫu ngẫu nhiên & Tách tập Học / Kiểm tra]
              (Kiểm chứng chuyên gia & Ảnh phân giải cao Google Earth)
               [Tính toán ngưỡng tự động dựa trên Histogram]
             (Worldpop_th, DMSP_th, EstISA_th, NDVI_MAX_th)
                 [Cây quyết định phân loại lớp phủ đô thị]
             [Khử nhiễu, tích hợp MOD44W và Đánh giá độ chính xác]

  1. Đánh giá phương pháp phân loại lớp phủ đô thị của GLCNMO: Phương pháp gốc của GLCNMO chia thế giới thành các nhóm quốc gia theo thu nhập (Bảng 2.1) và gán các ngưỡng tĩnh cho từng khu vực (Phụ lục 4). Tại khu vực Châu Á (Việt Nam thuộc nhóm thu nhập thấp/trung bình thấp trong các phiên bản cũ), các ngưỡng áp đặt cố định không phản ánh đúng thực tế tập trung dân cư dày đặc và hình thái đô thị xen cài nông nghiệp tại Việt Nam.

  2. Đánh giá biến đổi độ phân giải dữ liệu đầu vào: Thực nghiệm so sánh các tổ hợp biến đổi độ phân giải trên dữ liệu Worldpop, DMSP-OLS, MOD13Q1. Kết quả kiểm chứng qua các chỉ số $MSE$, $PSNR$ và $SSIM$ chỉ ra sự suy giảm chất lượng dữ liệu đối với từng thuật toán nội suy và tổng hợp, từ đó lựa chọn được tổ hợp phương pháp tối ưu giúp bảo toàn cấu trúc không gian và độ phản xạ phổ của các đối tượng mặt đất.

  3. Thuật toán cải tiến phân loại đô thị tại Việt Nam:

    • Xử lý chu kỳ NDVI: Dữ liệu NDVI được tổng hợp từ 23 chu kỳ trong năm của MOD13Q1 để trích xuất giá trị chỉ số thực vật cực đại ($NDVI_MAX$), giúp loại bỏ ảnh hưởng của mùa vụ nông nghiệp và mây mù che phủ.
    • Xây dựng tập mẫu huấn luyện và kiểm tra: Đề xuất phương pháp lấy mẫu ngẫu nhiên kết hợp đối soát đa tầng với ảnh vệ tinh độ phân giải cao trên Google Earth và bản đồ hiện trạng, được thẩm định bởi chuyên gia viễn thám.
    • Xác định ngưỡng phân loại tự động: Xây dựng thuật toán phân tích histogram của tập mẫu học để tự động xác định các giá trị ngưỡng phân tách cho từng biến: ngưỡng dân số ($Worldpop$), ngưỡng ánh sáng đêm ($DMSP$), ngưỡng bề mặt không thấm nước ($EstISA$) và ngưỡng thực vật ($NDVI_MAX$).
    • Quy trình phân loại: Ứng dụng cây quyết định dựa trên các ngưỡng động đã tối ưu để gán nhãn điểm ảnh đô thị và phi đô thị, sau đó kết hợp lớp mặt nước MOD44W và ranh giới hành chính GADM để loại trừ các điểm ảnh nhiễu.

  4. Đánh giá độ chính xác: Ma trận kiểm định độ chính xác trên tập mẫu kiểm tra độc lập (Bảng 2.5 và Bảng 2.6) cho thấy phương pháp cải tiến đạt độ chính xác tổng thể (Overall Accuracy) và hệ số Kappa cao hơn đáng kể so với việc áp dụng trực tiếp bản đồ gốc của GLCNMO tại lãnh thổ Việt Nam.

CHƯƠNG 3: XÂY DỰNG BỘ DỮ LIỆU LỚP PHỦ ĐÔ THỊ TẠI VIỆT NAM VÀ ỨNG DỤNG

Chương 3 trình bày việc triển khai hệ thống phần mềm, xây dựng cơ sở dữ liệu không gian và ứng dụng thực tiễn trong phân tích môi trường:

  1. Thiết kế và cài đặt hệ thống xử lý dữ liệu lớn:

    • Xây dựng kiến trúc hệ thống thu thập, xử lý, lưu trữ và hiển thị dữ liệu địa không gian dựa trên nền tảng SEAP (Big Spatial Data Exploration and Analysis Platform).
    • Thiết kế mô hình lưu trữ phân tán sử dụng HDFS và hệ quản trị cơ sở dữ liệu không gian, hỗ trợ các ca sử dụng tải dữ liệu raster/vector tự động, xử lý chuỗi thời gian song song và trực quan hóa bản đồ tương tác.

  2. Xây dựng bộ dữ liệu lớp phủ đô thị 2004–2015:

    • Thành lập thành công 4 bản đồ lớp phủ đô thị toàn quốc tại các năm 2004, 2008, 2012 và 2015.
    • Định lượng diện tích đô thị Việt Nam qua từng thời kỳ: diện tích đô thị tăng trưởng liên tục qua các giai đoạn 2004–2008, 2008–2012 và 2012–2015, tập trung mạnh nhất tại hai cực tăng trưởng là vùng Đồng bằng sông Hồng (trọng tâm Hà Nội) và vùng Đông Nam Bộ (trọng tâm TP. Hồ Chí Minh).

  3. Ứng dụng đánh giá tác động của mở rộng đô thị tới ô nhiễm bụi mịn PM2.5:

    • Phương pháp ước lượng PM2.5: Sử dụng dữ liệu độ dày quang học của sol khí (AOD/AOT) từ vệ tinh kết hợp dữ liệu khí tượng và dữ liệu đo đạc thực tế tại 6 trạm quan trắc mặt đất để xây dựng mô hình ước lượng nồng độ bụi mịn PM2.5 trung bình năm.
    • Xây dựng bản đồ phân bố PM2.5: Thành lập bản đồ phân bố nồng độ PM2.5 trung bình toàn quốc và khu vực Hà Nội cho các năm 2004, 2008, 2012, 2015.
    • Phân tích tương quan: Đánh giá mối quan hệ không gian giữa tốc độ mở rộng diện tích lớp phủ đô thị và sự gia tăng nồng độ ô nhiễm PM2.5 tại Việt Nam nói chung, thành phố Hà Nội và TP. Hồ Chí Minh nói riêng trong giai đoạn 2004–2015. Kết quả chứng minh sự gia tăng diện tích xây dựng và mật độ hoạt động đô thị có mối tương quan thuận rõ rệt với mức độ gia tăng ô nhiễm bụi mịn trong không khí.

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

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

  1. Định lượng tác động của quá trình biến đổi độ phân giải: Đánh giá có hệ thống mức độ suy giảm chất lượng ảnh viễn thám đa nguồn thông qua bộ ba chỉ số $MSE$, $PSNR$, $SSIM$, đồng thời chứng minh mức độ ảnh hưởng trực tiếp của sai số biến đổi độ phân giải đối với độ chính xác phân loại lớp phủ mặt đất tại Việt Nam.
  2. Cải tiến phương pháp phân loại lớp phủ toàn cầu của GLCNMO:
    • Đề xuất quy trình lấy mẫu ngẫu nhiên phân tầng kết hợp ảnh độ phân giải cao và tri thức chuyên gia để xây dựng tập mẫu học và tập mẫu kiểm tra tin cậy.
    • Đề xuất thuật toán tự động trích xuất ngưỡng phân tách từ biểu đồ tần suất (histogram) của các tập mẫu, thay thế hoàn toàn việc áp đặt ngưỡng tĩnh theo phân nhóm thu nhập quốc gia của GLCNMO.
    • Xây dựng quy trình tích hợp chuỗi thời gian NDVI 23 kỳ/năm để khử triệt để ảnh hưởng của mây và tính thời vụ của thảm thực vật nhiệt đới.

Đóng góp về mặt thực tiễn và ứng dụng

  1. Bộ dữ liệu lớp phủ đô thị Việt Nam (2004–2015): Xây dựng thành công bộ dữ liệu không gian số dạng raster về lớp phủ đô thị bao phủ toàn quốc tại 4 mốc thời gian (2004, 2008, 2012, 2015), cung cấp tư liệu lịch sử quan trọng cho các nghiên cứu đa ngành.
  2. Hệ thống phần mềm xử lý dữ liệu không gian lớn: Thiết kế hoàn chỉnh kiến trúc hệ thống thu thập, lưu trữ, xử lý song song và hiển thị trực quan bản đồ đô thị trên nền tảng công nghệ phân tán.
  3. Định lượng mối liên hệ môi trường đô thị: Cung cấp bằng chứng định lượng về sự tương quan giữa biến động diện tích không gian đô thị và tình trạng ô nhiễm bụi mịn PM2.5 tại Việt Nam trong chuỗi thời gian 11 năm.

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

Hạn chế của luận án

  • Dữ liệu đầu vào sử dụng trong mô hình phân loại toàn quốc có độ phân giải không gian ở mức trung bình (250m đến 1km), do đó chưa phản ánh được các đối tượng đô thị quy mô nhỏ hoặc cấu trúc vi mô bên trong các ngõ hẻm đô thị.
  • Chuỗi thời gian nghiên cứu của luận án dừng lại ở mốc năm 2015 do tính sẵn có và sự đồng bộ của một số chuỗi dữ liệu vệ tinh lịch sử (đặc biệt là dữ liệu ánh sáng ban đêm DMSP-OLS và dữ liệu bề mặt không thấm nước EstISA toàn cầu).

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

  • Tích hợp và thử nghiệm phương pháp phân loại trên các nguồn dữ liệu thế hệ mới có độ phân giải không gian và bức xạ cao hơn như dữ liệu ảnh vệ tinh Sentinel-1/2, dữ liệu ánh sáng ban đêm Suomi-NPP VIIRS.
  • Mở rộng chuỗi dữ liệu lớp phủ đô thị Việt Nam cập nhật liên tục từ năm 2016 đến thời điểm hiện tại.
  • Ứng dụng các kiến trúc mạng nơ-ron học sâu (Deep Learning) và điện toán đám mây địa không gian để nâng cao mức độ tự động hóa và độ chính xác phân loại lớp phủ mặt đất chi tiết.

Giá trị tham khảo

Luận án là tài liệu tham khảo khoa học có giá trị chuyên môn cao đối với:

  • Nghiên cứu sinh, học viên và giảng viên: Chuyên ngành Hệ thống thông tin, Công nghệ thông tin, Địa tin học (GIS), Viễn thám và Khoa học môi trường; đặc biệt hữu ích trong phương pháp luận xử lý dữ liệu địa không gian đa nguồn, kỹ thuật biến đổi độ phân giải raster và quy trình tách ngưỡng tối ưu.
  • Các cơ quan quản lý nhà nước và viện nghiên cứu quy hoạch: Bộ Xây dựng, Bộ Tài nguyên và Môi trường, Viện Quy hoạch Đô thị và Nông thôn Quốc gia trong việc sử dụng bộ dữ liệu lịch sử đô thị 2004–2015 để giám sát tốc độ đô thị hóa, đánh giá thực trạng sử dụng đất và xây dựng mô hình cảnh báo ô nhiễm không khí đô thị.

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

1. Luận án định nghĩa lớp phủ đô thị dựa trên những tiêu chí định lượng nào?

Theo Định nghĩa 1.2 trong luận án, lớp phủ đô thị được xác định dựa trên 4 tiêu chí định lượng kết hợp:

  • Mật độ dân số tính trên diện tích đất xây dựng đô thị đạt từ 1.500 người/km2 trở lên (dựa theo tiêu chuẩn phân loại đô thị tại Việt Nam và quy định của Liên Hợp Quốc).
  • Tỷ lệ diện tích các công trình xây dựng chiếm trên 50% tổng diện tích khu vực.
  • Không chứa các diện tích thảm thực vật cây xanh lớn, diện tích mặt nước hoặc khu vực không có/ít ánh sáng ban đêm.
  • Diện tích khu vực tập trung đạt tối thiểu từ 1 km2 trở lên.

2. Phương pháp của luận án có điểm gì khác biệt và cải tiến so với phương pháp gốc của GLCNMO?

Phương pháp gốc của GLCNMO gán các ngưỡng cố định dựa trên phân nhóm thu nhập bình quân đầu người của các quốc gia năm 2008 (Việt Nam thuộc nhóm nước thu nhập thấp/trung bình thấp với ngưỡng mật độ dân số, ánh sáng đêm và bề mặt không thấm nước rất thấp). Luận án đã cải tiến 3 điểm cốt lõi:

  • Đề xuất phương pháp lấy mẫu ngẫu nhiên kết hợp ảnh vệ tinh độ phân giải cao và kiểm chứng chuyên gia để tạo tập mẫu học và tập kiểm tra riêng cho Việt Nam.
  • Đề xuất thuật toán tự động tính toán ngưỡng tối ưu dựa trên phân tích biểu đồ tần suất (histogram) của các tập mẫu thực tế tại Việt Nam thay vì áp đặt ngưỡng tĩnh.
  • Xử lý chuỗi dữ liệu $NDVI$ theo 23 chu kỳ trong năm từ MOD13Q1 để lấy giá trị cực đại ($NDVI_MAX$), loại trừ ảnh hưởng của thời tiết nhiều mây và chu kỳ mùa vụ cây trồng.

3. Những chỉ số nào được sử dụng để đánh giá chất lượng các phương pháp biến đổi độ phân giải ảnh viễn thám?

Luận án sử dụng 3 chỉ số định lượng toán học để đánh giá mức độ biến dạng và sai lệch thông tin khi thực hiện chuyển đổi độ phân giải (tăng/giảm độ phân giải) trên các ảnh raster:

  • $MSE$ (Mean Squared Error): Sai số bình phương trung bình giữa ảnh gốc và ảnh sau biến đổi.
  • $PSNR$ (Peak Signal to Noise Ratio): Tỷ số giữa tín hiệu cực đại và nhiễu (tính bằng decibel - dB).
  • $SSIM$ (Structural Similarity Index): Chỉ số so sánh sự tương đồng cấu trúc không gian giữa hai ảnh.

4. Luận án đã sử dụng những nguồn dữ liệu viễn thám và không gian đa nguồn nào để phân loại đô thị?

Luận án tích hợp đồng thời 6 nguồn dữ liệu khác nhau:

  • Dữ liệu mật độ dân số Worldpop độ phân giải 100m.
  • Dữ liệu ảnh ánh sáng ban đêm DMSP-OLS Version 4 độ phân giải 1km.
  • Dữ liệu bề mặt xây dựng không thấm nước EstISA độ phân giải 1km.
  • Dữ liệu chỉ số khác biệt thực vật MOD13Q1 NDVI độ phân giải 250m (chu kỳ 16 ngày/kỳ, 23 kỳ/năm).
  • Dữ liệu bản đồ bề mặt nước MOD44W độ phân giải 250m.
  • Dữ liệu ranh giới địa giới hành chính GADM dạng vector.

5. Kết quả ứng dụng bộ dữ liệu lớp phủ đô thị vào bài toán ô nhiễm bụi mịn PM2.5 giai đoạn 2004–2015 chỉ ra điều gì?

Luận án đã ước lượng nồng độ PM2.5 từ dữ liệu độ dày quang học sol khí (AOD/AOT) của MODIS kết hợp dữ liệu hiệu chỉnh từ 6 trạm quan trắc mặt đất, xây dựng bản đồ phân bố PM2.5 của Việt Nam và Hà Nội qua các năm 2004, 2008, 2012, 2015. Kết quả phân tích không gian chứng minh sự mở rộng diện tích lớp phủ đô thị có mối tương quan thuận rõ nét với mức độ gia tăng nồng độ ô nhiễm bụi mịn PM2.5, đặc biệt nghiêm trọng tại hai đô thị trung tâm là Hà Nội và TP. Hồ Chí Minh trong giai đoạn 2004–2015.

Kết luận

Luận án của tác giả Phạm Tuấn Dũng đã giải quyết hoàn chỉnh bài toán phân loại và xây dựng cơ sở dữ liệu lớp phủ đô thị trên toàn lãnh thổ Việt Nam giai đoạn 2004–2015 từ nguồn dữ liệu viễn thám đa nguồn. Bằng việc cải tiến mô hình GLCNMO thông qua quy trình lấy mẫu tin cậy, thuật toán tách ngưỡng histogram tự động và xử lý chuỗi thời gian NDVI cực đại, luận án đã nâng cao độ chính xác phân loại lớp phủ đô thị so với các bộ dữ liệu toàn cầu hiện hành. Bên cạnh đóng góp về thuật toán và hệ thống phần mềm lưu trữ phân tán SEAP, công trình đã chứng minh tính ứng dụng thực tiễn cao khi lượng hóa thành công tác động của quá trình mở rộng không gian đô thị đến tình trạng gia tăng ô nhiễm bụi mịn PM2.5 tại Việt Nam.