Bối cảnh và vấn đề nghiên cứu

Khóa luận tốt nghiệp của tác giả Kim Văn Phúc (chuyên ngành Quản lý Tài nguyên Rừng & Môi trường, Đại học Lâm nghiệp Việt Nam, dưới sự hướng dẫn của PGS. Trần Quang Bảo) tập trung vào việc ứng dụng công nghệ viễn thám và hệ thống thông tin địa lý (GIS) nhằm thành lập bản đồ lớp phủ mặt đất tại khu vực miền núi.

Xuất phát từ thực tiễn quản lý tài nguyên, bản đồ lớp phủ mặt đất đóng vai trò then chốt trong việc theo dõi biến động tài nguyên rừng, giám sát sử dụng đất và hỗ trợ quy hoạch kinh tế - xã hội tại địa phương. Tuy nhiên, các phương pháp điều tra, đo vẽ bản đồ truyền thống đòi hỏi nguồn kinh phí lớn, tốn nhiều nhân lực, thời gian thực hiện kéo dài và thông tin thu được thường có độ trễ cao. Tại huyện Ba Vì (thành phố Hà Nội), đặc biệt là xã Vân Hòa – khu vực nằm ở vùng đệm Vườn quốc gia Ba Vì với địa hình đồi núi phức tạp, việc áp dụng tư liệu ảnh vệ tinh thế hệ mới Landsat 8 kết hợp GIS để thành lập bản đồ lớp phủ chưa từng được triển khai đồng bộ.

Mục tiêu nghiên cứu của đề tài được xác định cụ thể như sau:

  • Mục tiêu chung: Góp phần tạo lập cơ sở khoa học để đề xuất quy trình thành lập bản đồ lớp phủ mặt đất, nâng cao hiệu quả ứng dụng ảnh vệ tinh Landsat 8 trong công tác quản lý và giám sát tài nguyên thiên nhiên.
  • Mục tiêu cụ thể:
    1. Đánh giá hiện trạng quản lý sử dụng đất tại xã Vân Hòa, huyện Ba Vì, thành phố Hà Nội.
    2. Xây dựng bản đồ lớp phủ mặt đất xã Vân Hòa bằng tư liệu viễn thám và GIS.
    3. Đề xuất các giải pháp kỹ thuật và quản lý nhằm nâng cao hiệu quả khai thác ảnh Landsat 8 trong xây dựng bản đồ lớp phủ và quản lý tài nguyên.

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

  • Đối tượng nghiên cứu: Toàn bộ các trạng thái lớp che phủ trên bề mặt đất thuộc địa bàn nghiên cứu (bao gồm thực vật rừng, đất nông nghiệp, khu dân cư, mặt nước, đất ngập nước và đất trống).
  • Phạm vi không gian: Toàn bộ diện tích tự nhiên 3.290,98 ha của xã Vân Hòa, huyện Ba Vì, thành phố Hà Nội.
  • Phạm vi thời gian: Sử dụng dữ liệu ảnh vệ tinh Landsat 8 thu nhận vào ngày 20 tháng 10 năm 2015, kết hợp số liệu điều tra thực địa và số liệu thống kê quản lý đất đai địa phương giai đoạn 2011–2016.

Cơ sở lý thuyết và phương pháp

Khóa luận thiết lập khung lý thuyết dựa trên các nguyên lý vật lý của viễn thám quang học, cấu trúc dữ liệu không gian của GIS và đặc tính quang phổ của các đối tượng tự nhiên:

  • Nguyên lý viễn thám và hệ thống Landsat 8: Viễn thám quang học ghi nhận năng lượng bức xạ mặt trời phản xạ từ bề mặt Trái Đất qua 5 mắt xích: nguồn năng lượng, tương tác khí quyển, phản xạ bề mặt, cảm biến thu nhận và xử lý ảnh số. Vệ tinh Landsat 8 (nhiệm vụ LDCM do NASA và USGS hợp tác phóng ngày 11/02/2013) bay ở độ cao 705 km, góc nghiêng quỹ đạo 98°, chu kỳ lặp 16 ngày, mang hai bộ cảm OLI (9 kênh phổ phản xạ, lượng tử hóa 12 bit) và TIRS (2 kênh nhiệt sóng dài).
  • Đặc trưng phản xạ phổ của các đối tượng tự nhiên:
    • Thực vật: Hấp thụ mạnh bức xạ ở dải sóng đỏ ($0,6 - 0,7,\mu\text{m}$) do sắc tố diệp lục quang hợp và phản xạ mạnh ở dải cận hồng ngoại NIR ($\lambda > 0,72,\mu\text{m}$) do cấu trúc tế bào mô giậu; tổng năng lượng phản xạ tuân theo biểu thức $E_R = E_1 + E_G + E_{IR}$.
    • Thổ nhưỡng: Phổ phản xạ tăng đều từ vùng nhìn thấy đến hồng ngoại, không có tính thấu quang; khả năng phản xạ biến thiên nghịch với độ ẩm đất, tỷ lệ chất hữu cơ ($0,5 - 5,0%$) và hàm lượng oxide sắt.
    • Mặt nước: Hấp thụ gần như toàn bộ bức xạ ở dải cận hồng ngoại và hồng ngoại sóng ngắn; phản xạ ánh sáng khả kiến chịu ảnh hưởng trực tiếp bởi độ sâu, độ đục và nồng độ chất lơ lửng.
  • Phân biệt Lớp phủ (Land cover) và Sử dụng đất (Land use): Lớp phủ mô tả trạng thái vật lý - sinh học quan sát được trên bề mặt; sử dụng đất thể hiện chức năng kinh tế - xã hội do con người thiết lập. Biến động lớp phủ gồm biến đổi loại hình (conversion) và biến đổi thuộc tính nội tại (modification).

Phương pháp nghiên cứu thực nghiệm bao gồm:

  1. Thu thập và kế thừa dữ liệu: Tải ảnh Landsat 8 cảnh ngày 20/10/2015 từ cổng thông tin EarthExplorer của USGS (hệ quy chiếu UTM múi 48 Bắc, kinh tuyến trục $105^\circ$, elipxoid WGS-84, kích thước pixel 30m). Bổ trợ tư liệu ảnh Google Earth độ phân giải cao có tỷ lệ che phủ mây dưới 10% phục vụ đối soát. Kế thừa các báo cáo kinh tế - xã hội và bản đồ kiểm kê đất đai xã Vân Hòa.
  2. Điều tra thực địa: Thiết lập mẫu phiếu điều tra ngoại nghiệp ghi nhận tọa độ GPS, ảnh chụp mặt đất, mô tả hiện trạng sử dụng. Khảo sát 30 điểm mẫu thực tế phân bố ngẫu nhiên đại diện cho 6 trạng thái lớp phủ bằng máy định vị GPS cầm tay.
  3. Xử lý và phân tích ảnh số trên phần mềm ArcGIS 10:
    • Tiền xử lý: Chuyển đổi giá trị số $DN$ sang hệ số phản xạ đỉnh khí quyển TOA (Top of Atmosphere) bằng công cụ Raster Calculator theo công thức chuẩn của cảm biến OLI. Tổ hợp các kênh phổ và tích hợp kênh toàn sắc Band 8 (15m) để gia tăng độ sắc nét; nắn chỉnh hình học và cắt ảnh theo đường ranh giới hành chính xã.
    • Tính chỉ số khác biệt thực vật (NDVI): $$\text{NDVI} = \frac{\text{Band 5 (NIR)} - \text{Band 4 (Red)}}{\text{Band 5 (NIR)} + \text{Band 4 (Red)}}$$
    • Phân loại ảnh không kiểm định (Unsupervised Classification): Sử dụng thuật toán IsoData để gom cụm các pixel có đặc tính phổ tương đồng.
    • Giải đoán bằng mắt và gộp lớp: Dựa trên các dấu hiệu giải đoán trực tiếp (màu sắc, độ đậm nhạt, hình dạng, cấu trúc, kích thước) và dấu hiệu chỉ định gián tiếp để lập bộ khóa giải đoán 6 lớp: Rừng, Đất nông nghiệp, Dân cư, Đất ngập nước, Mặt nước, Đất trống. Dùng lệnh Reclassify để chuẩn hóa bản đồ phân loại.
    • Đánh giá độ chính xác sau phân loại: Xây dựng ma trận sai số (Error Matrix), tính tỷ lệ phân loại đúng tổng thể và hệ số tin cậy Kappa ($\kappa$): $$\kappa = \frac{N\sum_{i=1}^r x_{ii} - \sum_{i=1}^r (x_{i+} \cdot x_{+i})}{N^2 - \sum_{i=1}^r (x_{i+} \cdot x_{+i})}$$

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

Phần I: Tổng quan về vấn đề nghiên cứu

Chương này hệ thống hóa lịch sử phát triển của công nghệ viễn thám từ các mốc chụp ảnh khinh khí cầu (1858), máy bay quân sự, đến chuỗi vệ tinh tài nguyên ERTS-1/Landsat 1 (1972) đến Landsat 8 (2013). Tác giả tổng kết các đặc trưng cảm biến OLI và TIRS của Landsat 8, làm rõ tính ưu việt của dải lượng tử hóa 12 bit và bổ sung kênh 1 (ven bờ/sol khí) cùng kênh 9 (phát hiện mây ti). Chương cũng đi sâu phân tích cơ sở vật lý quang phổ phản xạ của ba nhóm vật thể chính: thực vật, đất và nước, đồng thời tổng quan tình hình ứng dụng viễn thám trong giám sát tài nguyên tại Việt Nam và trên thế giới.

Phần II: Mục tiêu – Nội dung – Phương pháp nghiên cứu

Nội dung chương trình bày chi tiết mục tiêu nghiên cứu, phạm vi không gian, thời gian và các nhiệm vụ kỹ thuật. Tác giả xây dựng sơ đồ quy trình công nghệ 5 bước xử lý ảnh viễn thám: thu thập dữ liệu -> tiền xử lý ảnh -> phân loại không kiểm định kết hợp phân tích NDVI -> giải đoán bằng mắt và gộp lớp Reclassify -> đánh giá độ chính xác thực địa bằng ma trận sai số và chỉ số Kappa. Chương cũng mô tả chi tiết mẫu phiếu điều tra ngoại nghiệp và thông số kỹ thuật của cảnh ảnh Landsat 8 thu nhận ngày 20/10/2015.

Phần III: Điều kiện tự nhiên – Dân cư – Kinh tế – Xã hội

Chương cung cấp bức tranh thực tế về địa bàn xã Vân Hòa:

  • Tự nhiên: Diện tích tự nhiên 3.290,98 ha, địa hình dốc nghiêng từ Tây sang Đông, chuyển tiếp từ núi cao Tản Viên (1.200 m) xuống gò đồi trung du (15–40 m) và đồng bằng trũng (10 m). Khí hậu nhiệt đới gió mùa có mùa đông lạnh (nhiệt độ trung bình $17,9^\circ\text{C}$, cực đoan $6,5^\circ\text{C}$) và mùa hè nóng ẩm (trung bình $26,1^\circ\text{C}$, đỉnh điểm $38,2^\circ\text{C}$). Lượng mưa chân núi đạt trên $1.731,4,\text{mm/năm}$, tập trung trên 90% từ tháng 6 đến tháng 10. Hệ thống thủy lợi có trên 50 km kênh mương nội đồng nhưng mới kiên cố hóa được 12 km.
  • Dân sinh - Kinh tế: Xã có 14 thôn, 3.002 hộ với 12.943 nhân khẩu (người Mường chiếm 43% với 5.256 khẩu; theo số liệu tổng điều tra nông thôn 2016 là 10.999 người). Cơ cấu kinh tế chuyển biến mạnh với tỷ trọng chăn nuôi chiếm 70,81%, trồng trọt chiếm 16,22%, dịch vụ chiếm 12,97%. Nghề chăn nuôi bò sữa phát triển mạnh dẫn đến diện tích trồng cỏ đạt 180–230 ha.
  • Quản lý tài nguyên và môi trường: Toàn xã có 305 ha rừng trồng; năm 2016 xảy ra 5 vụ cháy rừng gây thiệt hại 2 ha tại thôn Đa Cuống và Hòa Trung. Công tác cấp giấy chứng nhận quyền sử dụng đất sau dồn điền đổi thửa mới đạt khoảng 50% (1.100/2.200 hồ sơ). Rác thải sinh hoạt do Công ty Minh Quân thu gom 2 lần/tuần (khoảng 180 tấn/năm), toàn xã tồn tại 3 bãi rác tự phát.

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

Chương tập trung công bố và phân tích các phát hiện thực nghiệm:

  • Hiện trạng quản lý đất đai năm 2011: Đất rừng chiếm 52%, đất nông nghiệp chiếm 33%, đất trồng cỏ phục vụ bò sữa chiếm 6%, đất chuyên dùng và phi nông nghiệp chiếm 9%.
  • Hệ thống văn bản quy phạm pháp luật làm căn cứ điều chỉnh: Luật Đất đai sửa đổi, các chính sách tài nguyên nước, bảo vệ môi trường, đo đạc và bản đồ.
  • Khảo sát 30 điểm ngoại nghiệp phân bố: 6 điểm rừng, 6 điểm đất nông nghiệp, 6 điểm mặt nước, 5 điểm dân cư, 4 điểm đất ngập nước, 3 điểm đất trống.
  • Bộ khóa giải đoán ảnh xây dựng dựa trên sự kết hợp giữa ảnh phản xạ Landsat 8, ảnh thực địa và tư liệu Google Earth độ phân giải cao.
  • Phân tích trường chỉ số thực vật NDVI: giá trị biến thiên từ 0,0178 đến 0,8125; phân định rõ ranh giới phản xạ của thảm thực vật rừng so với các đối tượng phi thực vật.
  • Thành lập bản đồ lớp phủ mặt đất tỷ lệ cấp xã và thực hiện kiểm chứng chéo với số liệu thống kê địa chính của UBND xã Vân Hòa.
  • Tính toán ma trận sai số trên 30 điểm mẫu khảo sát, xác định hệ số Kappa đạt 0,86.

Chương V: Kết luận, tồn tại, kiến nghị

Chương tổng kết toàn diện các kết quả kỹ thuật đạt được của khóa luận, chỉ ra những nguyên nhân gây nhầm lẫn quang phổ giữa các lớp phủ và đề xuất các giải pháp kỹ thuật, cơ chế phối hợp quản lý tài nguyên đất, bảo vệ rừng trên địa bàn xã Vân Hòa.

Kết quả và đóng góp

Nghiên cứu đã hoàn thành việc xây dựng bản đồ lớp phủ mặt đất xã Vân Hòa gồm 6 trạng thái và phân tích định lượng chỉ số thực vật NDVI:

Lớp đối tượng Khoảng giá trị NDVI Đặc điểm phản xạ phổ
Rừng 0,67 – 0,81 Độ che phủ cao, hấp thụ mạnh kênh Đỏ, phản xạ rất mạnh kênh NIR
Đất nông nghiệp 0,52 – 0,67 Thảm thực vật hoa màu, lúa nước, độ đồng nhất phụ thuộc mùa vụ
Dân cư 0,33 – 0,52 Hỗn hợp vật liệu nhân tạo, nhà cửa bê tông xen kẽ vườn cây
Đối tượng khác (Nước, đất trống...) 0,17 – 0,33 Đất trống phản xạ cao, mặt nước hấp thụ hầu hết phổ cận hồng ngoại

Kết quả đo vẽ diện tích 6 lớp phủ trên bản đồ giải đoán ảnh Landsat 8 khi đối chiếu với số liệu thống kê thực địa của UBND xã Vân Hòa cho thấy mức độ phù hợp cao:

STT Trạng thái lớp phủ Diện tích trên bản đồ (ha) Tỷ lệ bản đồ (%) Diện tích thực tế (ha) Tỷ lệ thực tế (%)
1 Rừng 1.385,28 42 1.440,47 44
2 Đất nông nghiệp 937,71 28 1.007,97 31
3 Dân cư 774,09 24 710,17 21
4 Đất ngập nước 74,88 2 54,68 1
5 Nước 62,64 2 42,67 1
6 Đất trống 52,74 2 35,25 1
Tổng 3.287,34 100 3.291,21 100

Độ chính xác phân loại được kiểm định độc lập thông qua ma trận sai số 30 điểm khảo sát thực địa:

Đối tượng giải đoán Rừng Dân cư Đất nông nghiệp Đất ngập nước Mặt nước Đất trống Tổng hàng (Bản đồ)
Rừng 6 0 0 0 0 0 6
Dân cư 0 4 1 0 0 0 5
Đất nông nghiệp 0 1 5 0 0 1 7
Đất ngập nước 0 0 0 3 0 0 3
Mặt nước 0 0 0 1 6 0 7
Đất trống 0 0 0 0 0 2 2
Tổng cột (Thực địa) 6 5 6 4 6 3 30

Hệ số Kappa tính toán đạt $\kappa = 0,86$, thể hiện độ chính xác phân loại ở mức cao và kết quả bản đồ hoàn toàn có thể chấp nhận được trong quản lý thực tiễn.

Đóng góp và giải pháp tác giả đề xuất:

  • Khẳng định tính khả thi của quy trình công nghệ tích hợp ảnh viễn thám Landsat 8 miễn phí với phần mềm ArcGIS 10 trong việc lập bản đồ lớp phủ cấp xã vùng bán sơn địa, giảm thiểu kinh phí đo đạc.
  • Đề xuất giải pháp kỹ thuật nâng cao chất lượng phân loại: tăng cường chất lượng ảnh bằng kênh toàn sắc Band 8 (đưa độ phân giải từ 30m về 15m), chuyển đổi giá trị số $DN$ sang độ phản xạ $TOA$ nhằm loại bỏ nhiễu khí quyển.
  • Kiến nghị chính quyền địa phương sử dụng tư liệu bản đồ viễn thám định kỳ để kiểm soát ranh giới rừng Vườn quốc gia Ba Vì, quy hoạch diện tích chuyển đổi đất trồng cỏ nuôi bò sữa, và phòng ngừa nguy cơ cháy rừng trong mùa khô hanh.

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

Khóa luận chỉ rõ những tồn tại kỹ thuật xuất phát từ đặc trưng tư liệu và điều kiện thực địa:

  • Hạn chế về độ phân giải không gian: Kích thước pixel 30m của ảnh Landsat 8 gây ra hiện tượng điểm ảnh hỗn hợp (mixed pixel) tại các khu vực ranh giới, đặc biệt là các điểm dân cư phân tán xen kẽ vườn cây ăn quả hoặc các dải đất ven suối.
  • Hiện tượng chồng lấn phổ: Phổ phản xạ giữa đất ngập nước và mặt nước nông có độ tương đồng lớn, dễ dẫn đến phân loại nhầm. Tương tự, một số vị trí đất dân cư và đất nông nghiệp/đất trống có tín hiệu phổ gần nhau.
  • Phụ thuộc tư liệu phụ trợ: Do ảnh hưởng của địa hình núi cao và mây mù bao phủ khu vực rừng Ba Vì, tác giả phải sử dụng ảnh Google Earth độ phân giải cao để hỗ trợ thiết lập khóa giải đoán.
  • Hướng nghiên cứu tiếp: Mở rộng ứng dụng các dòng ảnh vệ tinh có độ phân giải không gian và thời gian cao hơn (như Sentinel-2 hoặc ảnh máy bay không người lái UAV); thử nghiệm các thuật toán phân loại có kiểm định (Supervised Classification) như Maximum Likelihood, Random Forest để tối ưu hóa độ chính xác; thiết lập chuỗi quan trắc đa thời gian để đánh giá biến động lớp phủ qua các năm.

Giá trị tham khảo

Báo cáo khóa luận là tài liệu tham khảo hữu ích cho:

  • Sinh viên các ngành Quản lý Tài nguyên Rừng, Quản lý Đất đai, Khoa học Môi trường, Trắc địa – Bản đồ và Địa tin học (GIS/Remote Sensing) đang thực hiện các đồ án chuyên ngành hoặc khóa luận tốt nghiệp.
  • Cán bộ kỹ thuật thuộc Phòng Tài nguyên và Môi trường huyện Ba Vì, cán bộ địa chính xã Vân Hòa và lực lượng kiểm lâm Vườn Quốc gia Ba Vì trong công tác cập nhật hiện trạng sử dụng đất và kiểm kê rừng định kỳ.
  • Phần có giá trị tham khảo cao nhất là quy trình 5 bước xử lý ảnh số trên ArcGIS 10 (từ bước tính toán TOA, thuật toán tính NDVI, kỹ thuật phân nhóm IsoData đến phương pháp lập ma trận sai số và tính chỉ số Kappa).

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

1. Cảnh ảnh Landsat 8 sử dụng trong nghiên cứu có những thông số kỹ thuật nào?
Cảnh ảnh được tải từ cổng thông tin EarthExplorer của USGS, chụp ngày 20/10/2015, có kích thước pixel 30m, lượng tử hóa 12 bit, chiếu theo hệ quy chiếu UTM múi 48 Bắc, elipxoid WGS-84, định hướng theo hướng Bắc bản đồ.

2. Công thức tính chỉ số thực vật NDVI trong khóa luận sử dụng những kênh phổ cụ thể nào của Landsat 8?
Tác giả áp dụng công thức $\text{NDVI} = \frac{\text{Band 5} - \text{Band 4}}{\text{Band 5} + \text{Band 4}}$, trong đó Band 5 là kênh cận hồng ngoại (NIR, bước sóng $0,85 - 0,88,\mu\text{m}$) và Band 4 là kênh đỏ (Red, bước sóng $0,64 - 0,67,\mu\text{m}$).

3. Bản đồ thành lập phân chia lớp phủ mặt đất xã Vân Hòa thành mấy trạng thái và trạng thái nào chiếm diện tích lớn nhất?
Bản đồ phân chia thành 6 trạng thái: Rừng, Đất nông nghiệp, Dân cư, Đất ngập nước, Mặt nước và Đất trống. Trong đó, Rừng chiếm diện tích lớn nhất với 1.385,28 ha (tương ứng 42% diện tích giải đoán trên bản đồ).

4. Độ chính xác của bản đồ lớp phủ được định lượng bằng chỉ số nào và kết quả đạt được là bao nhiêu?
Độ chính xác được kiểm định qua ma trận sai số xây dựng từ 30 điểm mẫu khảo sát thực địa bằng thiết bị GPS cầm tay, cho hệ số tin cậy Kappa đạt $\kappa = 0,86$, phản ánh độ chính xác ở mức cao.

5. Tại sao tác giả phải sử dụng thêm ảnh Google Earth trong quá trình giải đoán?
Tác giả sử dụng ảnh Google Earth có độ phủ mây dưới 10% để thay thế và hỗ trợ nhận diện các đối tượng rõ nét hơn, khắc phục tình trạng che phủ mây tại khu vực đồi núi rừng Ba Vì, đảm bảo độ chính xác khi thành lập bộ khóa giải đoán.

Kết luận

Khóa luận của tác giả Kim Văn Phúc đã xây dựng thành công bản đồ lớp phủ mặt đất gồm 6 nhóm đối tượng tại xã Vân Hòa, huyện Ba Vì dựa trên tư liệu viễn thám Landsat 8 và công nghệ GIS. Kết quả giải đoán phản ánh sát thực tế quản lý với tổng diện tích 3.287,34 ha, đạt hệ số tin cậy Kappa 0,86 qua kiểm định 30 điểm ngoại nghiệp. Công trình cung cấp quy trình kỹ thuật tiền xử lý và phân loại ảnh số rõ ràng, đóng góp nguồn tư liệu khoa học hữu ích cho công tác giám sát tài nguyên rừng và quy hoạch đất đai tại địa phương.