Tổng quan nghiên cứu

Quá trình đô thị hóa tại thành phố Hà Nội diễn ra với tốc độ nhanh chóng, kéo theo sự chuyển đổi mạnh mẽ của khoảng 30% đến 40% diện tích đất nông nghiệp và tự nhiên sang đất xây dựng và các bề mặt không thấm. Sự gia tăng các công trình bê tông kiên cố cùng với sự xuất hiện của các khu vực đất trống xen kẹt, đất dự án chậm triển khai đã đặt ra bài toán phức tạp cho công tác quản lý tài nguyên và quy hoạch không gian. Mật độ dân cư cao và sự phân bố dày đặc của các tòa nhà cao tầng gây hạn chế lớn cho phương pháp đo vẽ trắc địa trực tiếp trên thực địa.

Mục tiêu trọng tâm của nghiên cứu là đánh giá định lượng độ chính xác của kỹ thuật chiết tách đất xây dựng và đất trống khu vực đô thị bằng các chỉ số quang phổ viễn thám, thực nghiệm trên địa bàn thành phố Hà Nội bằng dữ liệu vệ tinh Landsat 8 năm 2018. Nghiên cứu tập trung so sánh hiệu năng của bốn chỉ số chuyên biệt gồm NDBI, EBBI, IBI và UI nhằm tìm ra giải pháp tối ưu cho việc phân loại lớp phủ.

Về mặt ý nghĩa thực tiễn, việc ứng dụng các chỉ số ảnh viễn thám giúp giảm thiểu hơn 50% chi phí khảo sát thực địa và rút ngắn thời gian thành lập bản đồ hiện trạng từ nhiều tuần xuống chỉ còn vài giờ xử lý trên máy tính. Kết quả nghiên cứu bảo đảm độ tin cậy phân loại tổng thể đạt trên 85%, cung cấp luận cứ khoa học tin cậy hỗ trợ các cơ quan chức năng giám sát chuyển đổi mục đích sử dụng đất, phục vụ công tác quy hoạch phát triển bền vững và đánh giá tác động của hiện tượng đảo nhiệt đô 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 trên nền tảng lý thuyết tương tác giữa bức xạ điện từ và các đối tượng tự nhiên, nhân tạo trên bề mặt Trái Đất. Năng lượng phản xạ phổ của vật thể biến thiên theo các dải bước sóng khác nhau: dải ánh sáng nhìn thấy (0,4 đến 0,76 μm), dải cận hồng ngoại NIR (0,77 đến 1,34 μm), dải hồng ngoại trung SWIR (1,55 đến 2,4 μm) và dải hồng ngoại nhiệt TIR (3 đến 22 μm). Các đối tượng lớp phủ được phân định qua ba khái niệm trọng tâm:

  • Đất xây dựng: Toàn bộ các bề mặt không thấm bao gồm nhà ở, hạ tầng giao thông, công trình công cộng và cơ sở sản xuất công nghiệp.
  • Đất trống: Các bề mặt không có thực vật che phủ như đất nông nghiệp sau thu hoạch, đất bỏ hoang hoặc mặt bằng dự án đang san lấp.
  • Điểm ảnh hỗn hợp (pixel): Các ô lưới không gian chứa nhiều loại vật thể khác nhau, gây ra hiện tượng pha trộn phản xạ phổ.

Hệ thống mô hình quang phổ được áp dụng bao gồm bốn chỉ số chính: Chỉ số đất xây dựng chuẩn hóa (NDBI), Chỉ số đất xây dựng và đất trống tăng cường (EBBI), Chỉ số đất xây dựng dựa trên chỉ số thành phần (IBI), và Chỉ số đô thị (UI).

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

Nghiên cứu sử dụng tư liệu ảnh vệ tinh Landsat 8 thu nhận năm 2018, bao gồm bộ cảm quang học OLI (độ phân giải không gian 30 m ở các kênh đa phổ, 15 m ở kênh toàn sắc, độ phân giải bức xạ 16-bit) và bộ cảm nhiệt TIRS (độ phân giải 100 m).

Quy trình tiền xử lý ảnh được thực hiện nghiêm ngặt qua các bước: hiệu chỉnh bức xạ chuyển đổi giá trị số (DN) sang giá trị phản xạ phổ đỉnh khí quyển và bề mặt, hiệu chỉnh ảnh hưởng khí quyển, cùng với hiệu chỉnh hình học dựa trên mạng lưới các điểm khống chế mặt đất (GCPs) nhằm duy trì sai số vị trí không gian nhỏ hơn 1 pixel (dưới 30 m).

Về phương pháp chọn mẫu, nghiên cứu áp dụng kỹ thuật lấy mẫu ngẫu nhiên phân tầng với cỡ mẫu gồm 280 điểm kiểm chứng thực tế phân bố đồng đều trên các dạng địa hình đặc trưng của Hà Nội. Lý do lựa chọn phương pháp phân tích chỉ số ảnh trên phần mềm ENVI 4.6 là nhờ khả năng tính toán đại số ma trận nhanh chóng, tự động hóa quy trình phân tách lớp phủ và hạn chế tối đa tính chủ quan so với phương pháp giải đoán ảnh truyền thống bằng mắt. Độ chính xác được kiểm định chặt chẽ thông qua việc lập ma trận sai số và tính toán hệ số Kappa.

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

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

Quá trình thực nghiệm và phân tích dữ liệu trên địa bàn thành phố Hà Nội mang lại các phát hiện cốt lõi sau:

  • Đối với nhiệm vụ chiết tách đất xây dựng, chỉ số NDBI mang lại hiệu quả cao nhất với độ chính xác toàn cục đạt 88,4% và hệ số tương đồng Kappa đạt 0,79. Chỉ số này phân định rất rõ nét ranh giới giữa các khối công trình bê tông, khu dân cư tập trung và hệ thống đường giao thông chính.
  • Đối với nhiệm vụ chiết tách đất trống, hai chỉ số IBI và NDBI thể hiện tính năng phân loại vượt trội, đạt độ chính xác tương ứng khoảng 85,2% và 86,7%. Sự khác biệt phản xạ phổ tại các kênh SWIR và NIR giúp nhận diện chính xác các diện tích đất canh tác bỏ hoang và đất công trình dở dang.
  • Chỉ số EBBI chứng minh khả năng thích ứng tốt tại các khu vực ven đô có cấu trúc lớp phủ phức tạp, đạt độ chính xác phân tách tổng hợp khoảng 84,6%, nhờ tích hợp kênh phổ nhiệt để giảm thiểu sự nhầm lẫn giữa đất xây dựng mật độ thấp và đất trống khô cằn.
  • Chỉ số UI cho tốc độ phản ánh nhanh các vùng đô thị hóa quy mô lớn nhưng tỷ lệ phân loại nhầm giữa đất trống và vật liệu ngói cũ vẫn ở mức 13,5%, đòi hỏi phải có bước lọc phân ngưỡng bổ sung.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự khác biệt giữa các chỉ số xuất phát từ tính chất phản xạ phổ của các loại vật liệu xây dựng và trạng thái độ ẩm của đất. Bê tông, nhựa đường và đất cát khô đều có xu hướng phản xạ mạnh ở dải sóng hồng ngoại trung (SWIR) và hấp thụ ở dải sóng cận hồng ngoại (NIR), gây ra hiện tượng đồng phổ dị vật. Đặc thù đô thị cổ Hà Nội với các khu phố cũ lợp mái ngói rêu phong nằm xen kẽ hệ thống sông hồ và cây xanh tạo ra các điểm ảnh hỗn hợp có đặc trưng quang phổ rất phức tạp.

So với các nghiên cứu tiền nhiệm tại khu vực Đông Á và Đông Nam Á, việc ứng dụng dữ liệu Landsat 8 với dải phổ thu hẹp đã giúp tăng độ chính xác phân loại thêm từ 4% đến 6% so với các thế hệ vệ tinh Landsat 5 hay Landsat 7 trước đây. Kết quả phân tích được biểu diễn trực quan thông qua bảng ma trận sai số đa chiều và biểu đồ tần số phân bố ngưỡng phổ, cho phép xác lập chính xác các khoảng giá trị ngưỡng tối ưu để tự động hóa việc trích xuất thông tin chuyên đề cho từng khu vực cụ thể.

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

  • Xây dựng quy trình tự động hóa chiết tách đất xây dựng định kỳ bằng chỉ số NDBI trên nền tảng ảnh vệ tinh đa phổ, nhằm mục tiêu giảm 70% thời gian cập nhật bản đồ hiện trạng sử dụng đất cho Sở Tài nguyên và Môi trường thành phố Hà Nội trong giai đoạn 2024 - 2026.
  • Ứng dụng tích hợp bộ đôi chỉ số IBI và EBBI trong công tác giám sát các dự án chậm tiến độ và diện tích đất nông nghiệp bỏ hoang, hướng đến chỉ tiêu phát hiện sớm 90% các biến động đất trống trên địa bàn các quận, huyện ngoại thành theo chu kỳ quét 16 ngày của vệ tinh.
  • Chuẩn hóa cơ sở dữ liệu phổ phản xạ mặt đất đối với các loại vật liệu xây dựng đặc thù tại các đô thị nhiệt đới, giao cho các trường đại học và viện nghiên cứu chuyên ngành thực hiện trong vòng 12 tháng tới để nâng độ chính xác phân loại lên mức trên 92%.
  • Ban hành tài liệu hướng dẫn kỹ thuật liên ngành giữa ngành Tài nguyên Môi trường và ngành Xây dựng về việc khai thác tư liệu viễn thám kết hợp hệ thống thông tin địa lý (GIS), triển khai thí điểm tại 5 đô thị trực thuộc trung ương trong vòng 2 năm tới.

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

  • Cơ quan quản lý nhà nước về tài nguyên, đất đai và quy hoạch đô thị: Khai thác quy trình phân tích chỉ số ảnh để giám sát biến động lớp phủ mặt đất, thanh tra hiện trạng xây dựng và kiểm soát các dự án sử dụng đất chậm tiến độ một cách khách quan.
  • Giảng viên, nghiên cứu viên ngành Trắc địa - Bản đồ và GIS: Sử dụng các công thức tính toán chỉ số quang phổ (NDBI, EBBI, IBI, UI) và phương pháp đánh giá độ chính xác qua ma trận sai số làm học liệu giảng dạy và mở rộng các hướng nghiên cứu chuyên sâu.
  • Học viên cao học và sinh viên các ngành Khoa học Trái Đất: Nắm bắt phương pháp tiền xử lý dữ liệu viễn thám Landsat 8, kỹ thuật nâng cao chất lượng ảnh và kỹ năng thao tác trên phần mềm chuyên dụng ENVI.
  • Doanh nghiệp tư vấn khảo sát và thành lập bản đồ: Ứng dụng giải pháp chiết tách tự động nhằm tiết kiệm 40% đến 60% chi phí nhân lực khảo sát thực địa trong các dự án quy hoạch vùng và thành lập bản đồ chuyên đề quy mô lớn.

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

1. Tại sao chỉ số NDBI lại cho kết quả chiết tách đất xây dựng tối ưu nhất tại Hà Nội?

Chỉ số NDBI khai thác triệt để sự tương phản giữa kênh hồng ngoại trung (SWIR) và kênh cận hồng ngoại (NIR). Với độ chính xác 88,4% và hệ số Kappa 0,79, NDBI phản ánh rất nhạy các bề mặt không thấm có mật độ xây dựng cao tại đô thị nhiệt đới.

2. Chỉ số EBBI có điểm gì vượt trội so với các chỉ số viễn thám truyền thống?

Chỉ số EBBI kết hợp đồng thời cả dải sóng hồng ngoại trung và hồng ngoại nhiệt (TIR). Sự kết hợp này giúp phân biệt rõ rệt bức xạ nhiệt giữa đất trống khô và các khối bê tông nhân tạo, mang lại độ chính xác phân tách đạt khoảng 84,6%.

3. Dữ liệu vệ tinh Landsat 8 mang lại lợi thế kỹ thuật nào cho nghiên cứu?

Landsat 8 sở hữu bộ cảm biến OLI và TIRS với độ phân giải bức xạ 16-bit và các dải sóng được tinh chỉnh hẹp hơn, giúp giảm thiểu đáng kể hiện tượng nhiễu khí quyển và tăng độ nhạy phân biệt các lớp phủ bề mặt lên hơn 5% so với thế hệ ảnh cũ.

4. Yếu tố nào gây nhầm lẫn nhiều nhất khi phân biệt đất trống và đất xây dựng?

Sự tương đồng về đặc tính phản xạ phổ ở các kênh quang học giữa đất khô cằn và bê tông cũ, kết hợp với các điểm ảnh hỗn hợp tại khu vực ngoại thành có diện tích nhỏ hơn 30 m, là nguyên nhân chính tạo ra khoảng 12% đến 15% sai số phân loại.

5. Phương pháp chiết tách bằng ảnh chỉ số có thể thay thế hoàn toàn đo đạc thực địa không?

Ảnh chỉ số không thay thế hoàn toàn mà đóng vai trò hỗ trợ đắc lực. Phương pháp này giúp tự động hóa khâu phân loại ban đầu với độ tin cậy trên 85%, trong khi kiểm chứng thực địa với mẫu chọn lọc là bước bắt buộc để hiệu chỉnh và hoàn thiện bản đồ.

Kết luận

  • Luận văn đã hoàn thiện quy trình khoa học chuẩn mực về xử lý ảnh viễn thám Landsat 8 trên phần mềm ENVI 4.6 để chiết tách thông tin không gian đô thị.
  • Đánh giá và lượng hóa thành công hiệu năng của bốn chỉ số quang phổ NDBI, EBBI, IBI và UI trong điều kiện cảnh quan đô thị nhiệt đới đặc thù tại Hà Nội.
  • Khẳng định chỉ số NDBI đạt hiệu quả tối ưu cho đất xây dựng với độ chính xác trên 88%, trong khi IBI và NDBI là lựa chọn hàng đầu cho đất trống với độ chính xác trên 85%.
  • Xác lập tầm quan trọng của công tác tiền xử lý bức xạ và hình học với sai số kiểm soát chặt chẽ dưới 1 pixel nhằm đảm bảo tính chuẩn xác cho các nghiên cứu lớp phủ.
  • Cung cấp cơ sở lý luận và thực tiễn vững chắc phục vụ công tác quy hoạch, điều tra tài nguyên đất và giám sát hiện tượng đảo nhiệt đô thị.

Trong giai đoạn tiếp theo, các nghiên cứu mở rộng nên tích hợp dữ liệu vệ tinh đa thời gian có độ phân giải cao hơn như Sentinel-2 hoặc ảnh siêu phân giải kết hợp thuật toán máy học. Các nhà khoa học và nhà quản lý được khuyến khích áp dụng ngay quy trình chỉ số phổ này để tối ưu hóa công tác giám sát biến động không gian đô thị một cách nhanh chóng và chuẩn xác.