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

Quá trình đô thị hóa và phát triển kinh tế - xã hội diễn ra nhanh chóng tại các đô thị lớn kéo theo sự gia tăng mật độ dân số, mở rộng các khu công nghiệp, công viên tư nhân và hệ thống hạ tầng giao thông. Thành phố Đài Bắc (Taipei City) là trung tâm thương mại và công nghiệp lớn nhất của Đài Loan. Đến năm 2010, dân số thành phố đã đạt 2.618.772 người, tăng 11.344 người so với năm 2009 (theo số liệu từ Taipei lookbook, 2010). Sự gia tăng dân số nhanh chóng cùng các hoạt động phát triển kinh tế - xã hội đã thúc đẩy quá trình mở rộng đô thị mạnh mẽ ra các vùng ngoại ô, gây áp lực lớn lên tài nguyên đất và dẫn đến sự biến đổi sâu sắc của lớp phủ mặt đất (Land Cover Change - LCC).

Sự thay đổi lớp phủ bề mặt, đặc biệt là quá trình chuyển đổi thảm thực vật tự nhiên thành các bề mặt nhân tạo (bê tông, mái nhà, đường nhựa), làm thay đổi đặc tính phản xạ và phát xạ nhiệt, tác động tiêu cực đến thời tiết và khí hậu địa phương, làm gia tăng nhiệt độ bề mặt và tạo nên hiện tượng đảo nhiệt đô thị. Việc ứng dụng công nghệ viễn thám (Remote Sensing - RS) kết hợp với Hệ thống Thông tin Địa lý (GIS) cho phép rút ngắn thời gian thành lập bản đồ so với các phương pháp truyền thống, đồng thời cung cấp công cụ giám sát chính xác hiện trạng và biến động môi trường đô thị.

Khóa luận cử nhân ngành Khoa học và Quản lý Môi trường của tác giả Nguyễn Hương Giang tập trung giải quyết các mục tiêu và nhiệm vụ nghiên cứu sau:

  1. Thành lập bản đồ lớp phủ mặt đất thành phố Đài Bắc cho hai mốc thời gian năm 2004 và năm 2014 từ dữ liệu ảnh vệ tinh Landsat.
  2. Phát hiện, lượng hóa và phân tích biến động lớp phủ mặt đất (LCC) của thành phố Đài Bắc trong giai đoạn 10 năm (2004 – 2014) bằng thuật toán so sánh sau phân loại (post-classification comparison).
  3. Thành lập bản đồ nhiệt độ bức xạ bề mặt (Land Surface Brightness Temperature - LSBT) năm 2004 và 2014 từ kênh nhiệt của ảnh Landsat nhằm đánh giá tác động của quá trình đô thị hóa lên môi trường nhiệt.
  4. Đánh giá mối quan hệ tương quan giữa sự thay đổi phân bố không gian các loại hình lớp phủ và sự biến đổi nhiệt độ bề mặt đất, từ đó đề xuất định hướng quản lý đất đai và bảo vệ môi trường đô thị.

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

  • Đối tượng nghiên cứu: Bản đồ biến động lớp phủ và nhiệt độ bức xạ bề mặt đất giai đoạn 2004 – 2014 tại thành phố Đài Bắc.
  • Phạm vi không gian: Toàn bộ ranh giới hành chính thành phố Đài Bắc, Đài Loan (nằm trọn trong khung cảnh vệ tinh Landsat Path 117 / Row 43).
  • Phạm vi thời gian thực hiện đề tài: 03 tháng (từ ngày 20/03/2015 đến ngày 29/06/2015). Giai đoạn dữ liệu quan trắc phân tích kéo dài 10 năm (2004 – 2014).

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

Cơ sở lý thuyết và khung phân loại

Khóa luận dựa trên các khái niệm nền tảng về lớp phủ mặt đất (Land Cover), biến động lớp phủ mặt đất (Land-Cover Change - LCC) và Hệ thống Phân loại Lớp phủ Mặt đất (Land Cover Classification System - LCCS). Lớp phủ mặt đất được xác định là trạng thái vật lý và sinh học bao phủ bề mặt trái đất, bao gồm nước, thực vật, đất trống và các cấu trúc nhân tạo. Quá trình biến động lớp phủ do con người thúc đẩy là một trong những nguyên nhân hàng đầu dẫn đến suy thoái đa dạng sinh học và biến đổi khí hậu cục bộ.

Về phân loại lớp phủ, tác giả áp dụng khung phân loại 3 cấp dựa trên việc điều chỉnh hệ thống phân loại của Anderson (1976). Hệ thống phân loại trong nghiên cứu bao gồm 3 đối tượng chính:

  • Đất xây dựng (Building): Khu dân cư, thương mại, công nghiệp, công trình giao thông và khu vui chơi giải trí (thể hiện bằng màu hồng trên bản đồ).
  • Thực vật (Vegetation): Rừng rụng lá, rừng thường xanh, rừng hỗn giao, đất trồng trọt, đồng cỏ và đất trống có cây bụi (thể hiện bằng màu xanh lá đậm).
  • Mặt nước (Water Bodies): Sông ngòi, hồ chứa, ao đầm và các vùng nước đọng (thể hiện bằng màu xanh dương).

Về nguyên lý nhiệt độ bức xạ bề mặt, nghiên cứu tiếp cận khái niệm nhiệt độ bức xạ (Brightness Temperature) - thông số đo lường năng lượng bức xạ vi sóng/hồng ngoại nhiệt truyền từ đỉnh khí quyển đến cảm biến vệ tinh, biểu thị dưới dạng nhiệt độ của một vật đen tương đương. Sự gia tăng diện tích mặt đệm nhân tạo với vật liệu hấp thụ nhiệt như mái tôn, nhựa đường làm gia tăng phát xạ nhiệt, trực tiếp nâng cao nhiệt độ bề mặt đất.

Nguồn dữ liệu sử dụng

Nghiên cứu sử dụng hai cảnh ảnh vệ tinh Landsat đa thời gian được tải về từ cổng dữ liệu Earth Explorer của Cục Khảo sát Địa chất Hoa Kỳ (USGS):

Thuộc tính dữ liệu Cảnh ảnh năm 2004 Cảnh ảnh năm 2014
Vệ tinh / Cảm biến Landsat 5 TM Landsat 8 OLI/TIRS Pre-WRS2
Mã định danh cảnh ảnh (Scene ID) LT51170432004194BKT02 LC81170432014237LGN00
Ngày chụp 12/07/2004 25/08/2014
Tọa độ cảnh (Path / Row) Path 117 / Row 43 Path 117 / Row 43
Tỷ lệ mây che phủ (Cloud Cover) 1,56% 7,xx%
Chất lượng ảnh (Image Quality) 7/10 9/10
Các kênh tổ hợp phân loại (Spatial Res: 30m) Kênh 1, 2, 3, 4, 5, 7 Kênh 2, 3, 4, 5, 6, 7
Kênh hồng ngoại nhiệt (Thermal Band) Kênh 6 (Độ phân giải nhiệt 120m) Kênh 10 và 11 (Độ phân giải nhiệt 100m)

Phương pháp xử lý ảnh và thành lập bản đồ

Nghiên cứu ứng dụng hai phần mềm chuyên dụng là ENVI (xây dựng trên nền tảng ngôn ngữ lập trình IDL) và ArcGIS (bao gồm các mô-đun ArcMap, ArcCatalog, ArcToolbox) để thực hiện quy trình xử lý qua các bước:

[Dữ liệu ảnh Landsat (2004, 2014)]
[Tiền xử lý ảnh (Pre-processing)]
[Phân loại ảnh (Image Classification)]
[Hậu phân loại & Đánh giá (Post-Classification)]
[Phát hiện biến động & Chồng xếp bản đồ (Change Detection & Overlay)]

Đối với quy trình tính toán nhiệt độ bức xạ bề mặt (LSBT):

  1. Chuyển đổi giá trị số số học (Digital Number - DN) sang phổ bức xạ thực tế ($L_\lambda$): $$L_\lambda = M_L \cdot Q_{cal} + A_L$$ Trong đó: $L_\lambda$ là giá trị bức xạ phổ tại khẩu độ cảm biến ($\text{W}/(\text{m}^2 \cdot \text{sr} \cdot \mu\text{m})$); $M_L$ là hệ số nhân tỷ lệ kênh ($Rad_mult_band$); $A_L$ là hệ số cộng tỷ lệ kênh ($Rad_add_band$); $Q_{cal}$ là giá trị pixel số hóa (DN).
  2. Chuyển đổi phổ bức xạ sang nhiệt độ bề mặt tại vệ tinh (Kelvin): $$T = \frac{K_2}{\ln\left(\frac{K_1}{L_\lambda} + 1\right)}$$ Trong đó: $T$ là nhiệt độ hữu hiệu tại vệ tinh tính bằng Kelvin (K); $K_1$ và $K_2$ là các hằng số hiệu chuẩn nhiệt cụ thể của từng cảm biến (Landsat 5 kênh 6 có $K_1 = 607.76$, $K_2 = 1260.56$; Landsat 8 kênh 10 và 11 sử dụng bảng hằng số tương ứng từ metadata).
  3. Chuyển đổi đơn vị nhiệt độ sang độ Celsius ($^\circ\text{C}$): $$T_{^\circ\text{C}} = T - 273$$

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

Tổng quan cơ sở lý thuyết và các công trình nghiên cứu liên quan

Nội dung chương tổng hợp hệ thống lý luận về công nghệ viễn thám và GIS, vai trò của các chỉ số thực vật (NDVI), chỉ số tổng bức xạ phản xạ (TRRI), cũng như phân tích các công trình nghiên cứu ứng dụng viễn thám trên thế giới và tại Việt Nam:

  • Nghiên cứu quốc tế: Điển hình là nghiên cứu biến động lớp phủ tại Pulivendula Taluk (Ấn Độ) sử dụng ảnh IRS-P6; nghiên cứu tích hợp dữ liệu dân số và ảnh Landsat TM tại Đông Bắc Latvia giai đoạn 1992 - 2007 (Fonji & Taff); nghiên cứu phân loại Maximum Likelihood tại Rize (Thổ Nhĩ Kỳ) của Selçuk Reis; hệ thống phân loại lớp phủ của Roach & Anderson (1976); và quy trình xử lý bù trừ ảnh hưởng địa hình, mây mù tại vùng núi hồ Kivu (Trung Phi) của Basnet & Vodacek (2011).
  • Nghiên cứu tại Việt Nam: Phân tích biến động đất đô thị thành phố Vinh giai đoạn 1992 - 2005 kết hợp Landsat và SPOT (Nguyễn Ngọc Phi, 2009); đánh giá biến động lớp phủ huyện Tiên Yên (Quảng Ninh) bằng ảnh SPOT và phần mềm eCognition (Nguyễn Thị Thu Hiền, 2013); xác định hệ số xói mòn C tại Tam Nông (Phú Thọ) của Trần Quốc Vinh (2009); xây dựng bản đồ lớp phủ Việt Nam chuỗi 2001 - 2003 bằng dữ liệu MODIS 32 ngày của Nguyễn Đình Dương (2006); và giám sát biến động rừng tại huyện Tánh Linh (Bình Thuận) kết hợp chỉ số NDVI và TRRI trên dữ liệu đa thời gian Landsat TM (Nguyễn Đình Dương, Kim Thoa, Nguyễn Thanh Hoàn, 2005).

Dữ liệu và phương pháp nghiên cứu (Materials and Methods)

Chương này trình bày chi tiết quy trình thực nghiệm tại Trung tâm Nghiên cứu Không gian và Viễn thám (CSRSR) thuộc Đại học Quốc lập Trung ương Đài Loan (NCU). Trọng tâm bao gồm việc thiết lập không gian địa lý, các bước xử lý bức xạ, hiệu chỉnh hình học đa cấp độ (Registration, Rectification, Geocoding, Orthorectification), tạo vùng quan tâm (ROI), và thiết lập các tham số giải đoán thuật toán phân loại có kiểm định Maximum Likelihood dựa trên hàm phân phối xác suất Gaussian. Đồng thời, chương làm rõ các tham số hằng số nhiệt ($K_1, K_2$) và các bước đại số raster để tính toán nhiệt độ bề mặt đất trên kênh 6 (Landsat 5) và kênh 10, 11 (Landsat 8).

Kết quả phân loại lớp phủ và đánh giá độ chính xác (Classification and Accuracy Assessment)

Tác giả trình bày kết quả giải đoán ảnh lớp phủ năm 2004 và 2014 sau khi áp dụng công cụ lọc trung vị $3 \times 3$ để loại bỏ hiện tượng nhiễu hạt "muối tiêu". Độ chính xác của kết quả phân loại được kiểm định nghiêm ngặt thông qua ma trận nhầm lẫn (Error Matrix):

  • Độ chính xác của nhà sản xuất (Producer's Accuracy) đạt trên 89,4%.
  • Độ chính xác của người sử dụng (User's Accuracy) đạt từ 85,x%.
  • Độ chính xác tổng thể (Overall Accuracy) cho cả hai thời điểm phân loại năm 2004 và 2014 đạt 94,4%, với hệ số Kappa đạt mức tương đồng rất cao (gần như hoàn hảo theo thang phân cấp của Landis & Koch, 1977).

Biến động lớp phủ và nhiệt độ bức xạ bề mặt giai đoạn 2004 - 2014 (Land Cover Change & Brightness Temperature Analysis)

Nội dung chương tập trung lượng hóa các chuyển dịch không gian giữa các lớp phủ và đối chiếu với biến động nhiệt độ:

  • Phân tích thống kê cho thấy sự chuyển đổi mạnh mẽ từ thảm thực vật sang đất xây dựng nhân tạo.
  • Bản đồ biến động nhiệt độ bề mặt (BTM) phản ánh trực quan sự dịch chuyển các vùng có nhiệt độ cao theo các trục mở rộng công trình đô thị và mạng lưới giao thông.

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

Kết quả nghiên cứu định lượng chính

Quá trình phân tích dữ liệu đa thời gian giai đoạn 2004 – 2014 tại thành phố Đài Bắc đem lại các kết quả định lượng cụ thể:

Chỉ số biến động (2004 – 2014) Giá trị định lượng Tỷ lệ so với tổng diện tích
Diện tích đất xây dựng (Building) gia tăng +4.995 ha +18,39%
Diện tích thảm thực vật (Vegetation) suy giảm -5.075 ha -18,xx%
Diện tích thực vật chuyển đổi trực tiếp sang đất xây dựng 2.883 ha 10,61%
Mức tăng nhiệt độ bức xạ bề mặt trung bình (LSBT) +3,69 $^\circ\text{C}$ Tăng rõ rệt trong 10 năm
Độ chính xác tổng thể phân loại (Overall Accuracy) 94,4% Hệ số Kappa mức tin cậy cao
Biến động diện tích các đối tượng lớp phủ chính (2004 - 2014):
Đất xây dựng (Building):  [████████████████████] +4.995 ha (+18,39%)
Thảm thực vật (Vegetation): [░░░░░░░░░░░░░░░░░░░░] -5.075 ha (-18,xx%)
Chuyển đổi Thực vật -> Xây dựng: 2.883 ha (10,61% tổng diện tích)

Kết quả nghiên cứu chứng minh mối quan hệ tỷ lệ thuận giữa mật độ xây dựng và nhiệt độ bức xạ bề mặt: các khu vực có tỷ lệ công trình nhân tạo cao luôn đi kèm với chỉ số bức xạ nhiệt lớn. Mức tăng nhiệt độ trung bình 3,69 $^\circ\text{C}$ trong vòng một thập kỷ là minh chứng rõ ràng cho tác động của quá trình mở rộng đô thị đối với vi khí hậu địa phương tại Đài Bắc.

Kiến nghị và giải pháp đề xuất

Từ các phát hiện nghiên cứu, tác giả đề xuất các định hướng nhằm hỗ trợ các nhà quản lý đô thị Đài Bắc:

  1. Quy hoạch không gian đô thị hợp lý, kiểm soát tốc độ bê tông hóa tại các vùng ven và khu vực chuyển tiếp sinh thái.
  2. Thiết lập các vùng đệm xanh, duy trì và phục hồi diện tích thảm thực vật công cộng nhằm giảm thiểu hiệu ứng đảo nhiệt đô thị.
  3. Ứng dụng tích hợp công nghệ viễn thám và GIS như một công cụ thường trực trong việc giám sát, dự báo biến động môi trường và hỗ trợ ra quyết định quản lý tài nguyên đất đai.

Đóng góp của đề tài

  • Thiết lập chuỗi cơ sở dữ liệu không gian định lượng về biến động lớp phủ và nhiệt độ bề mặt đất tại thành phố Đài Bắc trong giai đoạn đô thị hóa mạnh mẽ (2004 – 2014).
  • Chuẩn hóa quy trình xử lý dữ liệu viễn thám kết hợp giữa ảnh Landsat 5 TM và Landsat 8 OLI/TIRS, minh chứng tính hiệu quả của phương pháp Maximum Likelihood kết hợp bộ lọc sau phân loại và giải đoán nhiệt bức xạ trong nghiên cứu môi trường đô thị.

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

  • Hạn chế: Nghiên cứu chỉ tập trung phân tích tại 2 mốc thời gian đại diện (tháng 7/2004 và tháng 8/2014) vào mùa hè; chưa phân tích chuỗi biến thiên nhiệt độ liên tục theo mùa trong năm. Độ phân giải không gian của kênh nhiệt Landsat 5 (120m) và Landsat 8 (100m) còn ở mức trung bình, chịu ảnh hưởng nhất định bởi độ che phủ mây tại khu vực nghiên cứu.
  • Hướng nghiên cứu tiếp: Mở rộng nghiên cứu bằng việc sử dụng chuỗi dữ liệu vệ tinh đa thời gian có độ phân giải không gian và thời gian cao hơn; tích hợp các thông số khí tượng quan trắc mặt đất thực tế; mở rộng phân tích các chỉ số vi khí hậu khác như chỉ số lớp phủ thực vật điều chỉnh theo đất (SAVI) và chỉ số chênh lệch xây dựng (NDBI).

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên môn hữu ích cho sinh viên, học viên và nghiên cứu viên thuộc các chuyên ngành Khoa học Môi trường, Quản lý Đất đai, Địa lý học, Hệ thống Thông tin Địa lý (GIS) và Viễn thám.

Đặc biệt, tài liệu cung cấp giá trị thực hành cao ở các nội dung:

  • Quy trình chuẩn các bước tiền xử lý và cắt lọc dữ liệu ảnh Landsat trên môi trường ENVI và ArcGIS.
  • Công thức toán học và bảng hằng số hiệu chuẩn chi tiết ($K_1, K_2, M_L, A_L$) phục vụ trích xuất nhiệt độ bề mặt từ kênh hồng ngoại nhiệt của hai thế hệ vệ tinh Landsat 5 và Landsat 8.
  • Phương pháp luận đánh giá độ chính xác phân loại ảnh viễn thám bằng ma trận sai số và hệ số Kappa.

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

1. Khóa luận sử dụng nguồn dữ liệu viễn thám nào để phân tích biến động lớp phủ và nhiệt độ bề mặt tại Đài Bắc?
Nghiên cứu sử dụng hai cảnh ảnh vệ tinh Landsat tải từ USGS Earth Explorer thuộc khung Path 117 / Row 43: ảnh Landsat 5 TM chụp ngày 12/07/2004 (độ phủ mây 1,56%) và ảnh Landsat 8 OLI/TIRS chụp ngày 25/08/2014 (độ phủ mây 7,xx%).

2. Thuật toán nào được áp dụng để phân loại lớp phủ mặt đất và các lớp đối tượng bao gồm những gì?
Nghiên cứu sử dụng thuật toán phân loại có kiểm định cực đại xác suất (Maximum Likelihood Algorithm) trên phần mềm ENVI, kết hợp bộ lọc trung vị $3 \times 3$ sau phân loại. Hệ thống lớp phủ gồm 3 nhóm chính: Đất xây dựng (Building), Thảm thực vật (Vegetation) và Mặt nước (Water bodies).

3. Độ chính xác giải đoán của bản đồ lớp phủ trong nghiên cứu đạt mức bao nhiêu?
Độ chính xác của nhà sản xuất (Producer's accuracy) đạt từ 89,4%, độ chính xác của người dùng (User's accuracy) đạt từ 85,x%, và độ chính xác tổng thể (Overall accuracy) cho cả hai năm 2004 và 2014 đạt 94,4% với hệ số Kappa ở mức tương đồng rất cao.

4. Biến động diện tích đất xây dựng, thảm thực vật và nhiệt độ bề mặt đất tại Đài Bắc giai đoạn 2004 - 2014 diễn ra như thế nào?
Trong giai đoạn 2004 – 2014, diện tích đất xây dựng tăng thêm 4.995 ha (tương đương 18,39% tổng diện tích), diện tích thảm thực vật giảm 5.075 ha, trong đó có 2.883 ha thực vật chuyển trực tiếp thành đất xây dựng. Nhiệt độ bức xạ bề mặt đất (LSBT) trung bình tăng thêm 3,69 $^\circ\text{C}$.

5. Công thức nào được dùng để tính toán nhiệt độ bức xạ bề mặt từ kênh nhiệt Landsat?
Nhiệt độ bề mặt được xác định qua bước tính phổ bức xạ $L_\lambda = M_L \cdot Q_{cal} + A_L$, chuyển sang nhiệt độ Kelvin qua công thức $T = \frac{K_2}{\ln\left(\frac{K_1}{L_\lambda} + 1\right)}$, và đổi sang độ Celsius bằng biểu thức $T_{^\circ\text{C}} = T - 273$.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Hương Giang đã ứng dụng thành công công nghệ viễn thám và GIS để xây dựng bản đồ biến động lớp phủ và nhiệt độ bức xạ bề mặt đất tại thành phố Đài Bắc giai đoạn 2004 – 2014 với độ chính xác tổng thể đạt 94,4%. Kết quả nghiên cứu phản ánh rõ nét sự suy giảm của thảm thực vật (giảm 5.075 ha) và sự mở rộng của đất xây dựng (tăng 4.995 ha), đồng thời xác nhận mức gia tăng nhiệt độ bề mặt trung bình 3,69 $^\circ\text{C}$ do tác động của đô thị hóa. Công trình là tài liệu khoa học có giá trị thực tiễn cao, đóng góp cơ sở dữ liệu định lượng phục vụ công tác quy hoạch đô thị và quản lý môi trường bền vững.