Giới thiệu dự án

Sự biến động của đường bờ biển dưới tác động của biến đổi khí hậu, nước biển dâng và các hoạt động nhân sinh đang là thách thức nghiêm trọng đối với các vùng ven biển trên toàn cầu. Theo các nghiên cứu địa mạo quốc tế, hơn 70% các bãi biển cát trên thế giới đang bị xói lở với tốc độ trung bình từ 0.5 đến 1.0 m/năm. Tại Đài Loan – hòn đảo có tổng chiều dài bờ biển xấp xỉ 1.100 km – hơn 80% bờ biển cát đã ghi nhận hiện tượng xói lở nghiêm trọng trong 3 thập kỷ qua. Trong đó, huyện Nghi Lan (Yilan County) nằm ở phía Đông Bắc Đài Loan là khu vực chịu tổn thương nặng nề nhất do tần suất bão nhiệt đới dày đặc (hơn 350 cơn bão trong lịch sử quan trắc 100 năm) và sự xuất hiện của các công trình nhân tạo như đê chắn sóng và cảng biển.

Dự án khóa luận tốt nghiệp "Change Detection and Analysis of Shoreline in the Northeast of Taiwan Using Remotely Sensed Data" (Phát hiện và phân tích biến đổi đường bờ biển Đông Bắc Đài Loan bằng dữ liệu viễn thám) tập trung giải quyết bài toán giám sát không gian – thời gian động thái xói lở và bồi tụ bờ biển dọc dải cát dài 33 km từ Đầu Thành (Toucheng) đến Tô Áo (Su'ao) trong giai đoạn 15 năm (2000 – 2015).

+-----------------------------------------------------------------------------------+
|                              KHUNG BÀI TOÁN NGHIÊN CỨU                            |
+-----------------------------------------------------------------------------------+
|  [Dữ liệu Landsat 7/8 (2000-2015)] ---> [Hiệu chỉnh phổ & Hình học (20 GCPs)]    |
|                                                     |                             |
|                                                     v                             |
|  [Bản đồ Xói lở / Bồi tụ GIS]       <--- [Chỉ số nước MNDWI + Vector hóa tự động] |
+-----------------------------------------------------------------------------------+

Vấn đề thực tiễn và Pain Points

  • Hạn chế của phương pháp trắc địa truyền thống: Đo đạc thực địa bằng GPS hoặc toàn đạc điện tử trên bờ biển dài 103 km của Nghi Lan tốn kém chi phí, nguy hiểm trong mùa mưa bão và không thể tái hiện dữ liệu lịch sử liên tục.
  • Xung đột công trình nhân tạo: Cảng cá Wushi (Wushi Port) xây dựng đê chắn sóng nhô ra biển đã chặn đứng dòng bùn cát ven bờ (longshore sediment transport), gây bồi tụ cục bộ ở phía Bắc nhưng tạo ra hiện tượng "downdrift erosion" nghiêm trọng làm mất 2 km bãi biển Toucheng.
  • Biến động bồi tích sông Lanyang: Sông Lanyang dài 73 km bắt nguồn từ núi Nanhu (độ cao 3.740 m) mang lượng phù sa khổng lồ đổ ra biển Đông Bắc, tạo ra sự phân hóa phức tạp giữa bồi tụ cửa sông và xói lở vùng lân cận.

Mục tiêu dự án

  1. Trích xuất vị trí đường bờ biển tự động: Ứng dụng chỉ số khác biệt nước cải biên (Modified Normalized Difference Water Index - MNDWI) từ ảnh vệ tinh Landsat 7 ETM+ và Landsat 8 OLI để bóc tách ranh giới đất - nước với độ phân giải không gian 30 m.
  2. Định lượng tốc độ xói lở và bồi tụ: Xây dựng bản đồ phân vùng biến động đường bờ qua 15 năm (2000 – 2015) trên 7 bãi biển trọng điểm của huyện Nghi Lan.
  3. Phân tích nguyên nhân địa mạo học: Đánh giá tương quan giữa các yếu tố tự nhiên (bão, hoàn lưu gió mùa Đông Bắc, chế độ thủy triều, lưu lượng bồi tích sông Lanyang) và nhân sinh (xây dựng cảng Wushi, đê kè ven biển).
  4. Chuẩn hóa quy trình xử lý ảnh vệ tinh: Thiết lập chuỗi xử lý khép kín gồm hiệu chỉnh bức xạ tương đối (Relative Radiometric Rectification), trừ vật tối (Dark Object Subtraction - DOS) và nắn chỉnh hình học Image-to-Image với sai số RMS < 0.5 pixel.

Phương pháp tiếp cận và Kết quả kỳ vọng

  • Tiếp cận: Tích hợp công nghệ viễn thám đa thời gian (Multi-temporal Remote Sensing) và Hệ thông tin địa lý (GIS) trên nền tảng ENVI và ArcGIS.
  • Kết quả định lượng kỳ vọng: Bóc tách chính xác diện tích bồi tụ đạt ~1.04 km² và các vùng xói lở cục bộ; xác định vị trí biến động đường bờ với độ chính xác vị trí phụ thuộc pixel dưới 15 m sau nắn chỉnh.
  • Phạm vi nghiên cứu: 33 km bờ biển cát đồng bằng phù sa Lanyang (Lanyang alluvial plain); giới hạn thời gian từ ngày 10/08/2000 (Landsat 7 ETM+) đến ngày 16/11/2015 (Landsat 8 OLI).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí Khảo sát thực địa (GPS / RTK) Viễn thám quang học cổ điển (NDWI) Giải pháp đề xuất (MNDWI + Multi-temporal RS/GIS)
Chi phí triển khai Rất cao (nhân lực, thiết bị hiện trường) Thấp (ảnh miễn phí USGS) Rất thấp (tự động hóa xử lý dữ liệu mở)
Độ bao phủ không gian Rời rạc, giới hạn điểm đo Toàn diện (Swath width 185 km) Toàn diện (toàn bộ 103 km bờ biển Nghi Lan)
Khả năng triệt tiêu nhiễu Không áp dụng Kém (dễ nhầm lẫn giữa nước và đô thị/đất cát ẩm) Rất cao (kênh SWIR triệt tiêu triệt để nhiễu đô thị và đất trống)
Thời gian hồi quy (Revisit) Không liên tục (vài năm/lần) 16 – 18 ngày 16 ngày kết hợp đa cảm biến Landsat 7/8
Khả năng phân tích lịch sử Không khả thi trước thời điểm đo Khả thi từ 1984 Tối ưu với chuẩn hóa bức xạ đa thời gian

Yêu cầu người dùng theo mô hình MoSCoW

  • Must have (Bắt buộc): Bóc tách đường bờ biển năm 2000 và 2015; Tính toán chỉ số MNDWI trên ảnh Landsat; Bản đồ chuyên đề biến động xói lở/bồi tụ; Báo cáo định lượng diện tích thay đổi.
  • Should have (Cần có): Quy trình tiền xử lý ảnh tự động (DOS và Radiometric Calibration); Nắn chỉnh hình học ảnh đa thời gian với điểm khống chế mặt đất (GCPs).
  • Could have (Có thể có): Đối soát trực quan với ảnh vệ tinh độ phân giải cao Google Earth Historical Imagery; Phân tích tương quan theo từng phân khu hành chính (Toucheng, Zhuangwei, Wujie, Su'ao).
  • Won't have (Không làm trong phạm vi này): Mô hình hóa thủy động lực học sóng 3D thời gian thực; Hiệu chỉnh biến triều tự động theo trạm đo triều (Tidal gauge integration) do thiếu trạm quan trắc liên tục tại thời điểm vệ tinh chụp.
+------------------------------------------------------------------------------------+
|                      SƠ ĐỒ KIẾN TRÚC HỆ THỐNG XỬ LÝ DỮ LIỆU                        |
+------------------------------------------------------------------------------------+
|                                                                                    |
|   +---------------------------------+    +-------------------------------------+   |
|   | USGS EarthExplorer (Data Lake)  |    | ENVI 5.3 (Pre-processing Engine)    |   |
|   | - Landsat 7 ETM+ (2000-08-10)   |--->| - Radiometric Calibration (Gain/Off)|   |
|   | - Landsat 8 OLI/TIRS(2015-11-16)|    | - Dark Object Subtraction (DOS)     |   |
|   +---------------------------------+    | - Image-to-Image Registration (GCPs)|   |
|                                          +-------------------------------------+   |
|                                                             |                      |
|                                                             v                      |
|   +---------------------------------+    +-------------------------------------+   |
|   | ArcGIS 10.3 / Spatial Analysis  |    | Thuật toán MNDWI & Binarization     |   |
|   | - Vectorization & Topology Check|<---| - Green Band vs SWIR Band Ratio     |   |
|   | - Spatial Overlay & Union       |    | - Zero-Threshold Segmentation       |   |
|   | - Shoreline Change Hazard Map   |    +-------------------------------------+   |
|   +---------------------------------+                                              |
+------------------------------------------------------------------------------------+

Thiết kế hệ thống và Công nghệ

  • Technology Stack:
    • Xử lý viễn thám: ENVI 5.3, ERDAS Imagine 2015.
    • Phân tích không gian & Biên tập bản đồ: ESRI ArcGIS Desktop 10.3 (ArcMap, ArcCatalog, ArcToolbox).
    • Dữ liệu vệ tinh: Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI) do USGS cung cấp.
    • Hệ quy chiếu chuẩn: WGS_1984_UTM_Zone_51N (Dành cho khu vực Đài Loan).
+-----------------------------------------------------------------------------------+
|                CẤU TRÚC DỮ LIỆU ĐỊA KHÔNG GIAN (GEODATABASE SCHEMA)               |
+-----------------------------------------------------------------------------------+
|  Dataset: Shoreline_Dynamics_Yilan                                                |
|  |-- Raster: Landsat7_2000_MNDWI (Type: Float32, Resolution: 30m)                 |
|  |-- Raster: Landsat8_2015_MNDWI (Type: Float32, Resolution: 30m)                 |
|  |-- FeatureClass: Shoreline_2000_Vector (Type: Polyline, Fields: Length, Date)   |
|  |-- FeatureClass: Shoreline_2015_Vector (Type: Polyline, Fields: Length, Date)   |
|  |-- FeatureClass: Change_Polygons_2000_2015                                      |
|      |-- Attribute: Class (String: "Accretion" | "Erosion")                       |
|      |-- Attribute: Area_SqM (Double)                                             |
|      |-- Attribute: Area_Ha (Double)                                              |
|      |-- Attribute: Coastal_Sector (String: "Waiao", "Toucheng", "Zhuangwei", etc)|
+-----------------------------------------------------------------------------------+

Phương pháp luận nghiên cứu (Methodology)

Quy trình nghiên cứu được triển khai theo 4 giai đoạn logic chặt chẽ:

[Thu thập dữ liệu vệ tinh] ---> [Tiền xử lý ảnh phổ/hình học] ---> [Trích xuất đường bờ MNDWI] ---> [Phân tích biến động GIS]
  1. Giai đoạn 1 - Chuẩn bị và lọc dữ liệu: Thu thập ảnh Landsat có tỷ lệ che phủ mây dưới 10% (Cloud cover < 10%) vào cùng mùa trong năm để giảm thiểu sai số do chu kỳ thực vật và góc chiếu mặt trời.
  2. Giai đoạn 2 - Tiền xử lý hình ảnh:
    • Chuyển đổi giá trị số Digital Number (DN) sang bức xạ phổ (Top of Atmosphere Radiance).
    • Áp dụng thuật toán Dark Object Subtraction (DOS) để khử tán xạ khí quyển.
    • Nắn chỉnh hình học không gian (Geometric Registration) giữa ảnh năm 2000 (Base image) và ảnh năm 2015 (Warp image) thông qua 20 điểm GCPs phân bố đồng đều tại các điểm giao cắt giao thông, cửa sông và công trình kiên cố.
  3. Giai đoạn 3 - Trích xuất ranh giới mặt nước: Tính toán chỉ số MNDWI theo từng pixel, áp dụng ngưỡng phân ngưỡng nhị phân để bóc tách ranh giới đất/nước tự động.
  4. Giai đoạn 4 - Phân tích GIS đa thời gian: Vector hóa ranh giới nước thành polyline, đóng vùng không gian thành polygon và thực hiện toán tử Spatial Overlay (Union) để tính toán ma trận thay đổi diện tích bồi/xói.

Implementation và kết quả

Chi tiết thuật toán và Quy trình xử lý

Chỉ số khác biệt nước cải biên (Modified Normalized Difference Water Index - MNDWI) do Xu (2006) đề xuất được sử dụng thay thế cho chỉ số NDWI truyền thống của McFeeters (1996). MNDWI thay thế kênh cận hồng ngoại (NIR) bằng kênh hồng ngoại sóng ngắn (SWIR):

$$\text{MNDWI} = \frac{\text{Green} - \text{SWIR}}{\text{Green} + \text{SWIR}}$$

  • Đối với Landsat 7 ETM+ (2000):
    • $\text{Green} = \text{Band 2}\ (0.52 - 0.60\ \mu\text{m})$
    • $\text{SWIR 2} = \text{Band 7}\ (2.09 - 2.35\ \mu\text{m})$ $$\text{MNDWI}_{2000} = \frac{\text{Band 2} - \text{Band 7}}{\text{Band 2} + \text{Band 7}}$$
  • Đối với Landsat 8 OLI (2015):
    • $\text{Green} = \text{Band 3}\ (0.53 - 0.59\ \mu\text{m})$
    • $\text{SWIR 2} = \text{Band 7}\ (2.11 - 2.29\ \mu\text{m})$ $$\text{MNDWI}_{2015} = \frac{\text{Band 3} - \text{Band 7}}{\text{Band 3} + \text{Band 7}}$$

Script Python tự động hóa trích xuất và tính toán biến động đường bờ

Đoạn mã dưới đây minh họa quy trình tính toán MNDWI, phân ngưỡng nhị phân và bóc tách ranh giới đường bờ sử dụng rasterio, numpygeopandas:

import numpy as np
import rasterio
from rasterio.features import shapes
import geopandas as gpd
from shapely.geometry import shape

def compute_mndwi(green_path: str, swir_path: str, output_mndwi_path: str) -> np.ndarray:
    """
    Tính toán chỉ số MNDWI từ kênh Green và SWIR của Landsat.
    Triệt tiêu nhiễu đất và công trình đô thị hiệu quả.
    """
    with rasterio.open(green_path) as src_g, rasterio.open(swir_path) as src_s:
        green = src_g.read(1).astype(np.float32)
        swir = src_s.read(1).astype(np.float32)
        profile = src_g.profile

        # Tránh lỗi chia cho 0
        denominator = green + swir
        denominator[denominator == 0] = np.nan
        mndwi = (green - swir) / denominator

        # Ghi raster MNDWI
        profile.update(dtype=rasterio.float32, count=1, nodata=np.nan)
        with rasterio.open(output_mndwi_path, 'w', **profile) as dst:
            dst.write(mndwi, 1)

    return mndwi, profile

def extract_shoreline_vector(mndwi_path: str, threshold: float = 0.0, output_shp: str = None) -> gpd.GeoDataFrame:
    """
    Phân ngưỡng nhị phân mặt nước và vector hóa ranh giới bờ biển.
    """
    with rasterio.open(mndwi_path) as src:
        mndwi = src.read(1)
        transform = src.transform
        crs = src.crs

        # Phân ngưỡng: Mặt nước (MNDWI > threshold) = 1, Đất liền <= 0
        water_mask = np.where(mndwi > threshold, 1, 0).astype(np.int16)

        # Trích xuất polygon ranh giới
        results = (
            {'properties': {'raster_val': v}, 'geometry': s}
            for i, (s, v) in enumerate(shapes(water_mask, mask=(water_mask == 1), transform=transform))
        )

        geoms = list(results)
        gdf = gpd.GeoDataFrame.from_features(geoms, crs=crs)
        
        if output_shp:
            gdf.to_file(output_shp)
            
    return gdf

# Thực thi mẫu phân tích biến động
if __name__ == "__main__":
    # 1. Tính toán MNDWI cho Landsat 7 (2000) và Landsat 8 (2015)
    mndwi_2000, prof_2000 = compute_mndwi("L7_B2.tif", "L7_B7.tif", "mndwi_2000.tif")
    mndwi_2015, prof_2015 = compute_mndwi("L8_B3.tif", "L8_B7.tif", "mndwi_2015.tif")
    
    # 2. Bóc tách vector mặt nước
    water_2000_gdf = extract_shoreline_vector("mndwi_2000.tif", 0.0, "water_2000.shp")
    water_2015_gdf = extract_shoreline_vector("mndwi_2015.tif", 0.0, "water_2015.shp")
    
    # 3. Tính toán không gian giao thoa để xác định bồi tụ / xói lở
    # Bồi tụ (Accretion): Năm 2000 là nước, năm 2015 thành đất
    # Xói lở (Erosion): Năm 2000 là đất, năm 2015 thành nước
    accretion = gpd.overlay(water_2000_gdf, water_2015_gdf, how='difference')
    erosion = gpd.overlay(water_2015_gdf, water_2000_gdf, how='difference')
    
    print(f"Tổng diện tích bồi tụ: {accretion.geometry.area.sum() / 1e6:.2f} km²")
    print(f"Tổng diện tích xói lở: {erosion.geometry.area.sum() / 1e6:.2f} km²")

Kiểm định và Độ chính xác (Validation & Benchmarks)

  • Độ chính xác nắn chỉnh hình học (Geometric Accuracy): Sử dụng 20 điểm GCPs rải đều trên khu vực nghiên cứu, sai số toàn phương trung bình (Root Mean Square Error - RMSE) đạt giá trị ấn tượng 0.42 pixel (tương đương sai số mặt đất < 12.6 m trên độ phân giải 30 m), đáp ứng hoàn hảo tiêu chuẩn lập bản đồ biến động tỉ lệ trung bình.
  • Đối soát thực địa & Google Earth Imagery: Đường bờ biển trích xuất tự động qua ngưỡng MNDWI = 0.0 khớp hoàn toàn với đường viền mép nước trực quan trên ảnh màu thực RGB (Bands 7-4-2 và 3-2-1).
+------------------------------------------------------------------------------------+
|                KẾT QUẢ ĐỐI SOÁT VỊ TRÍ ĐIỂM KHỐNG CHẾ MẶT ĐẤT (GCPS)                |
+------------------------------------------------------------------------------------+
|  Số lượng điểm GCPs lựa chọn:               20 điểm                               |
|  Bậc mô hình chuyển đổi hình học:            Polynomial 1st Order (Affine)         |
|  Sai số RMS trục X:                          0.28 pixel (8.4 m)                   |
|  Sai số RMS trục Y:                          0.31 pixel (9.3 m)                   |
|  Sai số tổng hợp toàn hệ thống (Total RMS):  0.42 pixel (12.6 m)                  |
|  Tiêu chuẩn kỹ thuật chấp nhận:              RMS < 0.5 pixel  --> [ĐẠT YÊU CẦU]   |
+------------------------------------------------------------------------------------+

Kết quả đạt được theo từng khu vực trọng điểm

Sau 15 năm (2000 – 2015), đường bờ biển huyện Nghi Lan trải qua các biến động phức tạp với sự đan xen giữa bồi tụ mạnh và suy thoái xói lở cục bộ:

+---------------------------------------------------------------------------------------+
|                TỔNG HỢP BIẾN ĐỘNG ĐƯỜNG BỜ BIỂN NGHI LAN (2000 - 2015)                |
+------------------------------------+--------------------+-----------------------------+
| Phân khu ven biển (Sector)         | Xu hướng chính     | Hiện trạng & Nguyên nhân    |
+------------------------------------+--------------------+-----------------------------+
| Bãi biển Waiao & Cảng cá Wushi     | Bồi tụ phía Bắc /  | Đê chắn sóng chặn cát ven bờ|
|                                    | Xói lở phía Nam    | Cát tích tụ mạn đón dòng    |
| Phía Nam thị trấn Đầu Thành        | Xói lở nghiêm trọng| Bãi biển bị co hẹp dần      |
| Bờ biển Tráng Vi (Zhuangwei)       | Đan xen bồi / xói  | Tác động bão & sóng biển    |
| Cửa sông Lanyang & Ngũ Kết (Wujie) | Bồi tụ rất mạnh    | Nguồn phù sa bồi tích sông  |
| Phía Bắc Tô Áo (Su'ao)             | Bờ đá ổn định      | Cấu trúc địa chất vách đá   |
+------------------------------------+--------------------+-----------------------------+
  1. Khu vực Bãi biển Waiao và Cảng Wushi: Bồi tụ mở rộng ở phía Bắc đê chắn sóng cảng Wushi do cát bị chặn lại. Ngược lại, phía Nam cảng bị thiếu hụt nguồn bùn cát, dẫn đến hiện tượng xói lở sâu vào đất liền, đe dọa trực tiếp bãi tắm du lịch Toucheng vốn dài 2 km trước đây.
  2. Khu vực ven biển Zhuangwei và Wujie (Cửa sông Lanyang): Ghi nhận diện tích bồi tụ vượt trội đạt xấp xỉ 1.04 km². Lưu vực sông Lanyang đóng vai trò nguồn cung cấp trầm tích chính, bù đắp khối lượng cát bị bão biển cuốn trôi.
  3. Khu vực Su'ao: Do cấu trúc địa chất dạng vách đá sỏi kiên cố (cliff coast), đường bờ biển tại đây duy trì sự ổn định, biên độ dịch chuyển qua 15 năm là không đáng kể.

Đổi mới và đóng góp

Cải tiến kỹ thuật vượt trội

  1. Hiệu quả phân lập mép nước vượt trội của kênh SWIR: Việc khai thác kênh SWIR 2 (~2.2 $\mu\text{m}$) thay vì NIR giúp triệt tiêu hoàn toàn tín hiệu phản xạ giả từ các khu dân cư xây dựng tập trung ven biển Nghi Lan và lớp đất cát ẩm trên bãi triều, giảm thiểu tối đa hiện tượng phân loại nhầm (Misclassification rate giảm > 40% so với NDWI cổ điển).
  2. Quy trình nắn chỉnh tương đối không phụ thuộc trạm khí tượng: Kỹ thuật chuẩn hóa bức xạ tuyến tính (Relative Radiometric Rectification) kết hợp tập điểm đối chứng "Dark" và "Bright" giúp loại trừ độ lệch khí quyển giữa hai ảnh cách nhau 15 năm với sai số phản xạ tuyệt đối < 1% mà không cần dữ liệu độ đục khí quyển phức tạp tại trạm quan trắc.
  3. Phát hiện quy luật biến động do công trình nhân sinh: Nghiên cứu đã chứng minh bằng số liệu định lượng về tác động tiêu cực của đê kè cảng cá Wushi làm bẻ cong dòng chảy ven bờ một góc gần 90 độ ra biển khơi, tạo dòng rip current cuốn bùn cát ra xa bờ, cung cấp cơ sở khoa học xác thực cho các nhà hoạch định chính sách bảo vệ đới bờ.

Ứng dụng thực tế và triển khai

+------------------------------------------------------------------------------------+
|                         LỘ TRÌNH TRIỂN KHAI VÀ KHẢ NĂNG MỞ RỘNG                    |
+------------------------------------------------------------------------------------+
|  Phase 1: Chuẩn hóa dữ liệu viễn thám mở (Landsat / Sentinel)                     |
|     |                                                                              |
|  Phase 2: Tự động hóa Pipeline trích xuất ranh giới với Python / ArcPy / Cloud     |
|     |                                                                              |
|  Phase 3: Tích hợp hệ thống hỗ trợ ra quyết định quản lý tổng hợp đới bờ (ICZM)   |
|     |                                                                              |
|  Phase 4: Nhân rộng quy mô cho vùng Đồng bằng sông Cửu Long & Duyên hải Miền Trung|
+------------------------------------------------------------------------------------+
  • Ứng dụng quản lý tổng hợp đới bờ (Integrated Coastal Zone Management - ICZM): Cung cấp tài liệu địa không gian chính xác cho Cơ quan Xây dựng và Quy hoạch Đài Loan (CPAMI) trong việc cấp phép xây dựng công trình cảng biển và kè chắn sóng.
  • Quy hoạch du lịch bền vững: Cảnh báo nguy cơ biến mất bãi tắm tự nhiên tại Toucheng, hỗ trợ thiết kế các dự án nuôi bãi (Beach Nourishment) định kỳ sau mùa bão.
  • Khả năng chuyển giao và mở rộng: Quy trình xử lý viễn thám trong đồ án có khả năng áp dụng trực tiếp cho các vùng bờ biển nhạy cảm xói lở tại Việt Nam như dải duyên hải miền Trung, Cà Mau và khu vực Đồng bằng sông Cửu Long (Mekong River Delta).
  • Phân tích hiệu quả kinh tế (ROI): Tiết kiệm hơn 85% ngân sách so với các chiến dịch trắc địa định kỳ bằng tàu thuyền và thiết bị bay không người lái (UAV) trên quy mô vùng ven biển kéo dài hàng trăm km.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Độ phân giải không gian 30 m: Giới hạn của Landsat khiến các biến động đường bờ quy mô nhỏ hơn 15 – 30 m trong chu kỳ ngắn hạn chưa được thể hiện chi tiết.
  • Thiếu dữ liệu đồng bộ triều (Tidal Coordination): Do hạn chế về dữ liệu trạm quan trắc mực nước triều tức thời tại thời điểm vệ tinh bay qua, đường bờ bóc tách là đường viền mép nước tức thời (Instantaneous Waterline) thay vì đường triều cao trung bình nhiều năm (MHWL).
  • Độ bao phủ mây: Khí hậu Đông Bắc Á với gió mùa và mây mù dày đặc khiến số lượng cảnh ảnh quang học đạt tiêu chuẩn (Cloud < 10%) bị hạn chế.

Hướng nâng cấp đề xuất

  • Tích hợp ảnh viễn thám radar đa phân cực (SAR Sentinel-1): Khắc phục hoàn toàn trở ngại mây mù và hoạt động tốt trong điều kiện mưa bão với khả năng xuyên mây.
  • Nâng cấp độ phân giải với Sentinel-2 MSI (10 m) và PlanetScope (3 m): Giúp nâng độ chính xác vị trí đường bờ lên gấp 3 – 10 lần.
  • Tự động hóa với Machine Learning / Deep Learning: Ứng dụng các mạng nơ-ron tích chập (U-Net, SegNet) hoặc thuật toán Random Forest / Support Vector Machine (SVM) để tự động phân loại đa lớp địa mạo đới bờ.

Đối tượng hưởng lợi

+-------------------------------------------------------------------------------------+
|                              BẢNG PHÂN TÍCH NHÓM HƯỞNG LỢI                          |
+---------------------+---------------------------------------------------------------+
| Nhóm đối tượng      | Giá trị thực tiễn tiếp nhận                                   |
+---------------------+---------------------------------------------------------------+
| Sinh viên & Học viên| Nắm vững quy trình xử lý ảnh vệ tinh viễn thám chuẩn mực,     |
| ngành Môi trường/GIS| phương pháp tính chỉ số MNDWI và tích hợp ArcGIS/ENVI.        |
+---------------------+---------------------------------------------------------------+
| Kỹ sư Viễn thám /   | Sở hữu đoạn mã tự động hóa quy trình bóc tách bờ biển,        |
| GIS Developers      | tối ưu hóa pipeline xử lý raster địa không gian lớn.          |
+---------------------+---------------------------------------------------------------+
| Cơ quan Quản lý     | Tài liệu khoa học phục vụ quy hoạch đê kè ven biển, cảnh báo  |
| Đới bờ & Đô thị     | xói lở bãi biển du lịch và ứng phó biến đổi khí hậu.         |
+---------------------+---------------------------------------------------------------+
| Nhà nghiên cứu      | Khung phương pháp luận đối soát biến động không gian - thời   |
| Địa mạo bờ biển     | gian trên các dải cát chịu ảnh hưởng của bão và công trình.  |
+---------------------+---------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật phần cứng và phần mềm tối thiểu để triển khai giải pháp này là gì?

Hệ thống yêu cầu máy trạm cấu hình tối thiểu: CPU đa nhân (Intel Core i5 thế hệ 8 trở lên hoặc AMD tương đương), RAM tối thiểu 16 GB (khuyến nghị 32 GB để xử lý ảnh raster đa thời gian), ổ cứng SSD trống tối thiểu 50 GB. Về phần mềm: Hệ điều hành Windows/Linux, phần mềm ENVI 5.x / ArcGIS 10.x hoặc các công cụ nguồn mở thay thế như QGIS 3.x, Python 3.8+ cùng các thư viện chuyên dụng GDAL, Rasterio, GeoPandas, NumPy.

2. Làm thế nào để giải quyết bài toán sai lệch đường bờ do dao động thủy triều?

Để loại trừ sai số do triều, có hai giải pháp tối ưu: (1) Thu thập số liệu đo triều từng giờ từ các trạm hải văn lân cận để chọn các cảnh ảnh vệ tinh có cùng mực nước triều; (2) Áp dụng phương pháp mô hình hóa mặt triều 2D/3D (Tide-coordinated shoreline extraction) kết hợp thuật toán tam giác hóa Delaunay để ngoại suy đường triều chuẩn đồng nhất.

3. Tại sao chỉ số MNDWI lại vượt trội hơn hẳn so với NDWI và NDVI trong bóc tách bờ biển?

Chỉ số NDWI cổ điển sử dụng kênh NIR thường bị phản xạ mạnh bởi các công trình bê tông, khu dân cư và bề mặt cát khô, gây ra hiện tượng dương tính giả (False positive) khiến bờ biển bị nhận diện sai. Trong khi đó, nước hấp thụ gần như tuyệt đối bức xạ ở dải sóng hồng ngoại ngắn (SWIR), giúp MNDWI triệt tiêu độ phản xạ nền của đất và đô thị, tạo độ tương phản cực đại tại đường ranh giới mép nước.

4. Quy trình bảo trì và cập nhật dữ liệu định kỳ diễn ra như thế nào?

Quy trình có thể tự động hóa hoàn toàn theo chu kỳ 16 ngày (theo chu kỳ lặp của Landsat) hoặc 5 ngày (Sentinel-2). Script Python có thể thiết lập định kỳ kéo ảnh qua API của USGS/Copernicus, thực hiện tiền xử lý, bóc tách vector đường bờ và đẩy vào GeoDatabase để theo dõi cảnh báo biến động liên tục.

5. Chi phí triển khai hệ thống viễn thám so với đo đạc thực địa ước tính ra sao?

Dữ liệu Landsat và Sentinel được cung cấp hoàn toàn miễn phí từ USGS và ESA. Toàn bộ chi phí chủ yếu tập trung vào nhân lực phân tích và bản quyền phần mềm (có thể thay bằng QGIS/Python miễn phí). Chi phí tổng thể ước tính giảm hơn 90% so với việc thuê tàu khảo sát thủy đạc và trắc địa GPS trên toàn bộ dải bờ biển 103 km.


Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Vũ Tuấn Anh đã giải quyết thành công bài toán giám sát biến động đường bờ biển huyện Nghi Lan (Đông Bắc Đài Loan) giai đoạn 2000 – 2015 thông qua sự kết hợp hiệu quả giữa công nghệ viễn thám đa thời gian và GIS. Bằng việc ứng dụng chỉ số cải biên MNDWI trên nền dữ liệu vệ tinh Landsat 7 ETM+ và Landsat 8 OLI, nghiên cứu đã bóc tách chính xác vị trí đường bờ biển, xác lập bản đồ phân vùng biến động với diện tích bồi tụ đạt khoảng 1.04 km² tại khu vực hạ lưu sông Lanyang và làm rõ hiện tượng xói lở cục bộ nguy hiểm tại bãi biển Toucheng do ảnh hưởng của đê chắn sóng cảng cá Wushi.

Nghiên cứu khẳng định tính khả thi, độ tin cậy và hiệu quả kinh tế vượt trội của công nghệ viễn thám trong quản lý tổng hợp tài nguyên môi trường biển. Đây là tài liệu tham khảo có giá trị học thuật và thực tiễn cao, đóng góp nền tảng phương pháp luận vững chắc cho các dự án quy hoạch bảo vệ bờ biển và ứng phó với biến đổi khí hậu trong tương lai.