Tổng quan luận án

Luận án tiến sĩ Địa lý của tác giả Lại Anh Khôi có tiêu đề "Ứng dụng, nghiên cứu phát triển phương pháp xử lý ảnh số, theo dõi biến động tài nguyên thiên nhiên mặt đất vùng trung - hạ lưu sông Đà", thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội năm 2007 (Mã số chuyên ngành: 62.01). Luận án được trình bày trong 149 trang (chưa kể danh mục phụ lục và tài liệu tham khảo), tập trung giải quyết các vấn đề lý thuyết và thực tiễn trong kỹ thuật phân loại ảnh viễn thám số phục vụ quản lý tài nguyên môi trường.

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

Kỹ thuật xử lý ảnh số đã xuất hiện tại Việt Nam từ năm 1983 với hệ thống Robotron (Cộng hòa Dân chủ Đức) trong khuôn khổ dự án VIE/79/001 và VIE/83/04 do UNDP/FAO tài trợ, tiếp theo là phần mềm nội địa MDASER (1989). Dù các hệ thống phần cứng và phần mềm chuyên dụng (như ERDAS, PCI, ENVI, ILWIS) đã được trang bị ở nhiều cơ quan chuyên ngành, việc ứng dụng xử lý ảnh số vẫn chưa đạt hiệu quả tương xứng với vốn đầu tư. Luận án chỉ rõ hai nguyên nhân chính dẫn đến hạn chế này:

  1. Thiếu đội ngũ cán bộ chuyên sâu có hiểu biết bản chất về nguyên lý, khả năng và giới hạn của kỹ thuật xử lý ảnh số để làm chủ công nghệ.
  2. Còn thiếu các nghiên cứu mang tính học thuật chuyên sâu về các phương pháp xử lý ảnh số; phân loại ảnh số vốn là khâu được kỳ vọng nhiều nhất nhưng trên thực tế lại là khâu yếu nhất về độ tin cậy.

Về mặt thực tiễn, vùng Trung - Hạ lưu sông Đà là địa bàn có vị trí quan trọng về kinh tế, xã hội và quốc phòng. Sự vận hành của nhà máy thủy điện Hòa Bình và quá trình xây dựng nhà máy thủy điện Sơn La, cùng với các định hướng tái cơ cấu không gian lãnh thổ Tây Bắc, đã tạo ra những biến động sâu sắc về tài nguyên đất, rừng và nước mặt. Do đó, việc nghiên cứu một quy trình phân loại ảnh số có độ tin cậy cao để theo dõi biến động tài nguyên khu vực này là yêu cầu cấp bách.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu tổng quát của luận án là đề xuất các giải pháp nâng cao độ tin cậy của phương pháp phân loại ảnh số và vận dụng để nghiên cứu theo dõi biến động tài nguyên mặt đất vùng Trung - Hạ lưu sông Đà.

Luận án xác định 3 nhiệm vụ cụ thể:

  1. Nghiên cứu bản chất, đặc điểm, phân tích điểm mạnh, yếu của các phương pháp phân loại ảnh số trong viễn thám;
  2. Đề xuất các giải pháp nâng cao độ tin cậy của kết quả phân loại ảnh số và các kỹ thuật sau phân loại (post-classification) nhằm cải thiện chất lượng ảnh phân loại;
  3. Ứng dụng để xây dựng các bản đồ hiện trạng làm căn cứ xác định biến động tài nguyên thiên nhiên mặt đất vùng Trung - Hạ lưu sông Đà giữa hai thời kỳ: đầu thập kỷ 1990 (khi nhà máy thủy điện Hòa Bình mới hoạt động) và giai đoạn 2000 (trong quá trình thi công xây dựng thủy điện Sơn La).

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

  • Đối tượng nghiên cứu: Các phương pháp phân loại ảnh số viễn thám (có giám sát, không giám sát, mạng nơ-ron), các nguyên tắc và giải pháp nâng cao độ tin cậy của kết quả phân loại, cùng các thuật toán lọc và hiệu chỉnh sau phân loại.
  • Phạm vi không gian: Lưu vực sông Đà đoạn từ sau vị trí đập thủy điện Sơn La đến cuối lưu vực (thuộc vùng Trung và Hạ lưu sông Đà).
  • Phạm vi thời gian: Giai đoạn 1993 – 2000, tương ứng với hai thời điểm chụp của ảnh vệ tinh Landsat TM (27/12/1993) và Landsat ETM+ (04/11/2000).
  • Phạm vi nội dung: Tập trung vào hiện trạng và biến động lớp phủ bề mặt, bao gồm tài nguyên rừng, tài nguyên nước mặt và hiện trạng sử dụng tài nguyên đất (đất nông nghiệp, đất thổ cư, đất chưa sử dụng).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu viễn thám và xử lý ảnh số trên thế giới

Luận án điểm lại lịch sử phát triển của viễn thám từ các mốc thử nghiệm chụp ảnh hàng không và đo đạc địa hình của Arago (1840), Aimé Laussedat (1849), khinh khí cầu (1858), khí tài diều (1880) và máy bay của Wilbur Wright (1909). Các nghiên cứu ứng dụng phản xạ phổ thực vật của Colwell (1956), cùng sự ra đời của các vệ tinh quan trắc Trái Đất như Tiros-1 (1960), Landsat-1/ERTS-1 (1972), các thế hệ SPOT (Pháp), ERS và Envisat (Châu Âu), IRS (Ấn Độ), Radarsat (Canada), IKONOS, QuickBird, Terra và Aqua mang các bộ cảm ASTER, MODIS (Mỹ) đã định hình sự phát triển của viễn thám hiện đại.

Về xử lý ảnh số, luận án dẫn chứng các nền tảng toán học và nhận dạng mẫu hình thành từ cuối thập niên 1950. Các phương pháp phân loại thống kê (như Maximum Likelihood, Minimum Distance, Mahalanobis) và phương pháp phi tham số (như Mạng nơ-ron nhân tạo) đã được chuẩn hóa trong nhiều tài liệu chuyên khảo quốc tế.

Tình hình nghiên cứu và ứng dụng viễn thám tại Việt Nam

Tại Việt Nam, sự phát triển viễn thám gắn liền với chương trình nghiên cứu 3 tầng (vệ tinh, máy bay, mặt đất) năm 1980 nhân chuyến bay của Phạm Tuân mang máy ảnh MKF-6. Luận án tổng hợp các công trình tiêu biểu của các tác giả trong nước qua các giai đoạn:

  • Sử dụng ảnh Landsat MSS kiểm kê và theo dõi biến động rừng toàn quốc giai đoạn 1975–1983 của Viện Điều tra Quy hoạch Rừng (1983).
  • Các công trình ứng dụng viễn thám đa ngành được báo cáo tại Hội nghị Viễn thám Khu vực (Hà Nội, 1989): nghiên cứu địa mạo của Lê Đức An và Võ Thịnh; địa chất của Nguyễn Thanh Giang; địa chất thủy văn của Nguyễn Thượng Hùng và Kim Minh Công; thổ nhưỡng của Nguyễn Bá Nhuận và Nguyễn Thị Choắt; thảm thực vật của Phan Kế Lộc và Phan Phú Bồng; hiện trạng sử dụng đất của Phạm Trung Lương, Trần Minh Ý, Lê Văn Thành, Trương Thị Hoà Bình; sinh thái ven biển của Trịnh Đình Cương; biến động đường bờ của Nguyễn Đình Dương, Võ Đức Tuyến, Nguyễn Hồng Châu, Phạm Việt Cường; cửa sông của Phạm Việt Cường, Nguyễn Hồng Châu, Trần Minh Hiền; cây dừa của Nguyễn Viết Chiến, Trần Trọng Tĩnh, Nguyễn Thị Thu Vân; bản đồ ngập lụt của Nguyễn Phương Thảo và Nguyễn Thị Kỳ Nam; độ đục nước biển của Bùi Doãn Trọng, Hoàng Viết Giao, Lại Anh Khôi, Trần Minh Hồng, Nguyễn Thu Nga; nhiệt độ mặt biển của Bùi Doãn Trọng, Đỗ Minh Phú, Bùi Trọng Tuyên, Lại Anh Khôi, Nguyễn Thành Long, Phạm Quốc Trung, Nguyễn Thế Tường, Trương Trọng Xuân.
  • Xây dựng phần mềm xử lý ảnh số nội địa: chương trình MDASER (1989) trên MS-DOS của Bùi Doãn Trọng, Lại Anh Khôi, Hoàng Viết Giao, Trần Minh Hồng, Đinh Thị Minh; phần mềm VM ImagePro và WinASEAN.
  • Các nghiên cứu sau năm 1993: ứng dụng ảnh radar trong nghiên cứu ngập lụt và sử dụng đất của Viện Khoa học và Công nghệ Việt Nam; nghiên cứu trượt lở, xói mòn, biến động lòng sông của Đại học Quốc gia Hà Nội và Đại học Mỏ - Địa chất; các luận án tiến sĩ của Trương Thị Hoà Bình, Phạm Văn Cự, Nguyễn Đình Minh, Nguyễn Ngọc Thạch, Vũ Anh Tuân.

Khoảng trống luận án lựa chọn giải quyết

Luận án chỉ ra rằng dù ứng dụng viễn thám tại Việt Nam phát triển rộng rãi, phần lớn công trình chỉ dừng lại ở mức vận dụng các quy trình có sẵn mà thiếu các nghiên cứu cơ bản về phương pháp luận xử lý ảnh. Đặc biệt, các phân tích về sự sai lệch giữa lớp phổ và lớp thông tin, việc thiết lập ngưỡng xác suất thống kê và sự phối hợp giữa các thuật toán tự động với giải đoán chuyên gia chưa được giải quyết một cách có hệ thống. Luận án tập trung vào khoảng trống này nhằm nâng cao độ tin cậy toán học cho kết quả phân loại.

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

Cơ sở lý thuyết và khung phân tích

Luận án dựa trên nền tảng của lý thuyết nhận dạng mẫu và không gian thuộc tính (feature space / spectral space) $N$ chiều, trong đó mỗi điểm ảnh được biểu diễn bằng một véc-tơ độ xám $\mathbf{x} = [x_1, x_2, \dots, x_N]^t$.

  1. Phương pháp phân loại có giám sát:

    • Phân loại hợp lý tối đa (Maximum Likelihood): Sử dụng nguyên lý Bayes để xác định hàm mật độ xác suất có điều kiện $p(\mathbf{x}/\omega_i)$ dưới giả định phân bố chuẩn nhiều chiều. Hàm phân tách được xác định theo công thức: $$g_i(\mathbf{x}) = \ln p(\omega_i) - \frac{1}{2} \ln |\mathbf{\Sigma}_i| - \frac{1}{2} (\mathbf{x} - \mathbf{m}_i)^t \mathbf{\Sigma}_i^{-1} (\mathbf{x} - \mathbf{m}_i)$$ Trong đó $\mathbf{m}_i$ là véc-tơ trung bình, $\mathbf{\Sigma}_i$ là ma trận hiệp biến của lớp $\omega_i$, và $p(\omega_i)$ là xác suất tiên nghiệm. Để kiểm soát điểm ảnh ngoại lai, luận án ứng dụng phân bố $\chi^2$ để lập ngưỡng phân loại $T_i$.
    • Phương pháp khoảng cách tối thiểu (Minimum Distance): Phân loại dựa trên khoảng cách Euclid hoặc khoảng cách theo cạnh khối (City-block Distance) từ điểm ảnh đến tâm $\mathbf{m}_i$. Luận án chứng minh phương pháp khoảng cách Euclid là trường hợp riêng của Maximum Likelihood khi ma trận hiệp biến là ma trận đường chéo đồng nhất $\mathbf{\Sigma}_i = \sigma^2 \mathbf{I}$.
    • Phương pháp khoảng cách Mahalanobis: Rút gọn từ Maximum Likelihood khi giả định các lớp có chung ma trận hiệp biến $\mathbf{\Sigma}$, loại bỏ ảnh hưởng của tương quan giữa các kênh phổ.
    • Phương pháp hình hộp (Parallelepiped): Xác định biên phân loại theo khoảng biến thiên độ lệch chuẩn trên từng kênh; khắc phục hiện tượng trùng phủ bằng phép biến đổi thành phần chính (PCA).
    • Mạng nơ-ron nhân tạo (Artificial Neural Network): Sử dụng mạng nhiều lớp lan truyền ngược (Backpropagation) với các hàm truyền dạng chữ S: hàm Logistic $g(u) = \frac{1}{1 + e^{-u}}$, hàm Hyperbolic $h(u) = \frac{1 - e^{-u}}{1 + e^{-u}}$, và hàm Tang-hyperbolic $\tanh(u) = \frac{e^u - e^{-u}}{e^u + e^{-u}}$.

  2. Phương pháp phân loại không giám sát: Thuật toán K-means và ISODATA (Interactive Self-Organizing Data Analysis Technique) tự động nhóm cụm các điểm ảnh dựa trên tiêu chuẩn khoảng cách phổ tối thiểu.

Phương pháp nghiên cứu thực tế

  • Phương pháp toán học và thống kê: Đánh giá cấu trúc ma trận hiệp biến, tính toán hàm mật độ xác suất, phân tích ma trận nhầm lẫn (confusion matrix).
  • Phương pháp viễn thám và xử lý ảnh số: Tiền xử lý ảnh (nắn chỉnh hình học, hiệu chỉnh bức xạ, cân bằng biểu đồ tần số độ xám), phân loại đa bước, lọc không gian sau phân loại (Majority filter, Sieve filter, lọc hình thái học Dilation và Erosion).
  • Phương pháp khảo sát thực địa và bản đồ: Kiểm chứng các vị trí mẫu giải đoán ngoài thực địa; kết hợp bản đồ địa hình và bản đồ hiện trạng sử dụng đất.
  • Phương pháp GIS: Chồng xếp lớp dữ liệu không gian, phân tích tương quan giữa biến động lớp phủ với các yếu tố địa hình (độ dốc).

Nguồn tư liệu và dữ liệu

Nhóm dữ liệu Chi tiết thông số Nguồn cung cấp / Mục đích sử dụng
Ảnh vệ tinh Landsat TM Thu nhận ngày 27/12/1993; Path 127, Row 45 và 46; độ phân giải không gian 30m Dữ liệu nền phân loại hiện trạng lớp phủ năm 1993
Ảnh vệ tinh Landsat ETM+ Thu nhận ngày 04/11/2000; Path 127, Row 45 và 46; độ phân giải không gian 30m Dữ liệu nền phân loại hiện trạng lớp phủ năm 2000
Bản đồ địa hình Tỷ lệ 1:100.000 Nắn chỉnh tọa độ hình học và xây dựng mô hình số độ cao (DEM), bản đồ độ dốc
Bản đồ chuyên đề Bản đồ hiện trạng sử dụng đất năm 2000; Atlas Quốc gia Cung cấp thông tin kiểm chứng và đối chiếu ranh giới
Tài liệu kinh tế - xã hội Báo cáo điều tra tổng hợp Tây Bắc, báo cáo ĐTM thủy điện Hòa Bình và Sơn La Phân tích nguyên nhân và động lực biến động tài nguyên
Phần mềm xử lý ENVI Thực hiện các thuật toán phân loại và biên tập sau phân loại

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

Chương 1: Viễn thám và xử lý ảnh số trong viễn thám

Chương 1 hệ thống hóa lịch sử phát triển của công nghệ viễn thám từ chụp ảnh trên khinh khí cầu, máy bay đến hệ thống vệ tinh quan trắc hiện đại. Tác giả phân tích 6 chức năng cơ bản của kỹ thuật xử lý ảnh số:

  1. Hiệu chỉnh ảnh: Đưa ảnh về trạng thái phản ánh trung thực bề mặt (khử nhiễu, nắn chỉnh hình học, định chuẩn bức xạ).
  2. Điều chế ảnh: Làm nổi bật các thông tin quan tâm (tăng cường độ tương phản, lọc thông gian, phân ngưỡng, tính chỉ số thực vật).
  3. Hiển thị thông tin: Chuyển đổi dữ liệu số thành biểu diễn trực quan.
  4. Chiết xuất thông tin: Nhận dạng và phân loại đối tượng tự động.
  5. Mã hóa thông tin: Nén và lưu trữ dữ liệu truyền thông.
  6. Mô hình hóa hệ thống thu ảnh: Phục vụ thiết kế thiết bị quang điện tử mới.

Đồng thời, chương này tổng kết toàn diện quá trình du nhập và phát triển viễn thám tại Việt Nam từ năm 1980 đến 2007, chỉ rõ sự mất cân đối giữa đầu tư trang thiết bị với nghiên cứu học thuật nền tảng.

Chương 2: Đánh giá về các phương pháp phân loại ảnh số

Chương 2 đi sâu phân tích bản chất toán học của bài toán nhận dạng mẫu trong không gian phổ $N$ chiều. Tác giả chỉ ra sự khác biệt bản chất giữa:

  • Lớp thông tin (Informational Classes): Các đối tượng thực tế trên mặt đất cần xác định (rừng, đất nông nghiệp, mặt nước).
  • Lớp phổ (Spectral Classes): Các cụm điểm ảnh có đặc tính phản xạ phổ tương đồng trong không gian dữ liệu. Một lớp thông tin trong thực tế thường phân hóa thành nhiều lớp phổ khác nhau do ảnh hưởng của độ ẩm đất, hướng phơi của sườn địa hình hay mật độ thảm thực vật.

Tác giả so sánh chi tiết ưu và nhược điểm của các thuật toán:

  • Maximum Likelihood: Chặt chẽ về mặt xác suất nhưng phụ thuộc tuyệt đối vào giả thiết phân bố chuẩn và chất lượng mẫu huấn luyện.
  • Minimum Distance: Tốc độ tính toán nhanh, phù hợp khi mẫu huấn luyện phân tán hoặc không đủ đại diện để tính ma trận hiệp biến, nhưng bỏ qua mối tương quan giữa các kênh phổ.
  • Mahalanobis: Giải pháp trung gian hữu hiệu khi các kênh ảnh có sự khác biệt về độ phân tán bức xạ nhưng không đủ điều kiện lập ma trận hiệp biến riêng cho từng lớp.
  • Phương pháp hình hộp: Ranh giới đơn giản, tốc độ xử lý nhanh, song tạo ra tỷ lệ điểm ảnh không được phân loại lớn khi các lớp có sự tương quan phổ.
  • Mạng nơ-ron: Khả năng xấp xỉ phi tuyến mạnh, không đòi hỏi giả thiết về phân bố chuẩn của dữ liệu, nhưng đòi hỏi thời gian huấn luyện dài và khó kiểm soát cơ chế trọng số nội tại.

Chương 3: Nâng cao độ tin cậy của kết quả phân tích ảnh

Dựa trên các hạn chế đã phân tích, tác giả đề xuất 5 nguyên tắc cơ bản nhằm nâng cao độ tin cậy của quá trình phân loại ảnh:

  1. Hệ thống phân loại phải phù hợp: Thiết kế dựa trên các lớp phổ thực tế trên ảnh, sau đó mới nhóm gộp thành các lớp thông tin.
  2. Số liệu mẫu phải thực sự đại diện cho mỗi lớp: Đảm bảo tính toán chính xác véc-tơ trung bình và ma trận hiệp biến.
  3. Tận dụng chiều của không gian phổ: Giữ đủ số kênh phổ cần thiết để tránh làm suy giảm khả năng phân tách giữa các lớp.
  4. Tích hợp các thông tin bổ trợ trong quá trình phân loại: Sử dụng tư liệu địa hình, thổ nhưỡng và bản đồ chuyên đề để hỗ trợ giải đoán.
  5. Sử dụng các giá trị ngưỡng để kiểm soát quá trình phân loại: Loại trừ các điểm ảnh nằm ngoài vùng phân bố thống kê tin cậy ($\chi^2$) đưa vào nhóm chưa phân loại để xử lý riêng.

Chương này xây dựng Quy trình phân tích ảnh tích hợp gồm các bước tuần tự:

Thực nghiệm kiểm chứng trên tập dữ liệu thực tế thông qua các ma trận nhầm lẫn (Bảng 3.1, 3.2, 3.3) chứng minh: việc phân loại trực tiếp theo lớp thông tin hoặc cắt giảm số kênh phổ làm tăng mạnh tỷ lệ sai số nhầm lẫn giữa các đối tượng; ngược lại, quy trình phân loại kết hợp theo lớp phổ có sử dụng ngưỡng kiểm soát mang lại độ chính xác cao nhất.

Chương 4: Ứng dụng kỹ thuật xử lý ảnh số trong nghiên cứu biến động tài nguyên vùng Trung-Hạ lưu sông Đà

Chương 4 áp dụng toàn bộ quy trình đề xuất vào việc phân loại ảnh vệ tinh Landsat TM 1993 và Landsat ETM+ 2000 trên lưu vực sông Đà từ sau đập Sơn La đến Hòa Bình và hạ lưu.

Quy trình thực nghiệm bao gồm: nắn chỉnh tọa độ hình học theo bản đồ địa hình 1:100.000, cân bằng biểu đồ độ xám giữa hai thời kỳ, phân loại sơ bộ bằng ISODATA, phân loại chi tiết bằng Maximum Likelihood có ngưỡng, lọc nhiễu không gian bằng bộ lọc đa số (Majority) và bộ lọc Sieve, kết hợp giải đoán mắt trên tổ hợp màu (R:4, G:2, B:1).

Kết quả xây dựng thành công bản đồ hiện trạng lớp phủ năm 1993 và năm 2000, từ đó tính toán biến động diện tích và cơ cấu sử dụng đất:

  • Đất có rừng: Có xu hướng phục hồi và gia tăng diện tích rõ rệt nhờ các chương trình trồng rừng và tái sinh tự nhiên.
  • Đất nông nghiệp và đất thổ cư: Gia tăng ở những mức độ khác nhau, tập trung dọc theo các thung lũng, bậc thềm sông và các khu vực tái định cư.
  • Diện tích mặt nước: Biến đổi cục bộ liên quan đến chế độ điều tiết nước của hồ chứa thủy điện Hòa Bình.
  • Đất chưa sử dụng: Thu hẹp đáng kể do được khai thác chuyển đổi sang đất lâm nghiệp và đất canh tác nông nghiệp.
  • Phân bố theo độ dốc: Tác giả tích hợp mô hình số độ cao để phân tích sự dịch chuyển của các loại hình lớp phủ trên các cấp độ dốc khác nhau, cung cấp cơ sở định lượng cho công tác quy hoạch sử dụng đất lưu vực.

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

Những đóng góp mới về mặt lý luận và phương pháp

  1. Đề xuất 5 nguyên tắc nâng cao độ tin cậy phân loại: Luận án đã hệ thống hóa các yêu cầu bắt buộc đối với việc lựa chọn lớp phổ, tính đại diện của mẫu, bảo toàn không gian chiều phổ, tích hợp dữ liệu bổ trợ và sử dụng ngưỡng thống kê $\chi^2$ trong xử lý ảnh viễn thám.
  2. Xây dựng quy trình phân tích ảnh tích hợp: Kết hợp chặt chẽ giữa thuật toán phân loại không giám sát (ISODATA) để phát hiện cấu trúc phổ tự nhiên, phân loại có giám sát (Maximum Likelihood) có kiểm soát ngưỡng xác suất, và kỹ thuật giải đoán bằng mắt kết hợp xử lý lọc sau phân loại. Quy trình này dung hòa được mối quan hệ giữa yêu cầu nhiệm vụ và khả năng phản xạ phổ thực tế của tư liệu ảnh.

Những đóng góp mới về mặt thực tiễn

  1. Bộ bản đồ và dữ liệu biến động tài nguyên vùng Trung - Hạ lưu sông Đà: Luận án đã cung cấp hệ thống số liệu định lượng và bản đồ hiện trạng lớp phủ chi tiết các năm 1993 và 2000 cho khu vực từ sau thủy điện Sơn La đến hạ lưu sông Đà.
  2. Đánh giá khách quan thực trạng sử dụng đất: Chứng minh rằng trong giai đoạn 1993 – 2000, diện tích đất chưa sử dụng tại khu vực nghiên cứu đã được thu hẹp đáng kể; sự gia tăng của đất lâm nghiệp, đất nông nghiệp và đất thổ cư phản ánh quá trình khai thác tài nguyên đất tích cực hơn và môi trường tự nhiên của khu vực có những chuyển biến rõ rệt.

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

  • Về công tác nghiên cứu: Cần đầu tư thích đáng cho các nghiên cứu cơ bản về cơ sở lý thuyết viễn thám và xử lý ảnh số, tránh việc áp dụng máy móc các phần mềm thương mại khi chưa nắm vững bản chất toán học của các thuật toán.
  • Về ứng dụng thực tiễn: Áp dụng quy trình phân tích tích hợp nhiều bước có kiểm soát ngưỡng cho các chương trình điều tra, kiểm kê tài nguyên định kỳ bằng ảnh viễn thám đa thời gian tại các vùng sinh thái nhạy cảm.

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

Hạn chế được xác định từ phạm vi và tư liệu

  • Độ phân giải không gian và ảnh hưởng thời tiết: Tư liệu sử dụng trong luận án là ảnh quang học Landsat TM và ETM+ với độ phân giải không gian 30m. Do đặc thù khu vực miền núi Tây Bắc thường xuyên có mây mù bao phủ, việc lựa chọn các cảnh ảnh hoàn toàn không có mây trùng khớp vào các thời điểm nghiên cứu gặp nhiều khó khăn.
  • Nguồn dữ liệu chuyên biệt: Tại thời điểm thực hiện luận án (2007), việc khai thác các tính năng lập thể của ảnh SPOT, kỹ thuật giao thoa của ảnh radar hay xử lý chuỗi phổ đa kênh MODIS chưa được triển khai đầy đủ do những hạn chế về phương pháp xử lý chuyên biệt tại Việt Nam.

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

  • Mở rộng ứng dụng quy trình phân loại tích hợp trên các nguồn ảnh viễn thám có độ phân giải không gian cao hơn (như SPOT-5, IKONOS, QuickBird) và ảnh đa phổ nhiều kênh (MODIS).
  • Nghiên cứu kết hợp dữ liệu ảnh quang học với dữ liệu radar viễn thám nhằm khắc phục triệt để ảnh hưởng của điều kiện mây mù nhiệt đới.
  • Tự động hóa sâu hơn các bước can thiệp của chuyên gia trong quá trình chỉnh sửa, biên tập và lọc không gian sau phân loại.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị trong các lĩnh vực Địa lý, Bản đồ học, Viễn thám và Hệ thông tin địa lý (GIS):

  • Đối với nghiên cứu sinh và học viên cao học: Luận án cung cấp cơ sở toán học chi tiết về các hàm phân tách xác suất, ma trận hiệp biến, nguyên lý mạng nơ-ron và các thuật toán phân cụm điểm ảnh (Chương 2 và Chương 3).
  • Đối với giảng viên và nhà nghiên cứu viễn thám: Quy trình tích hợp phân loại kết hợp với các kỹ thuật lọc hình thái học (Dilation, Erosion, Sieve, Majority) là tài liệu tham khảo chuẩn mực cho các bài thực hành xử lý ảnh chuyên sâu trên phần mềm ENVI.
  • Đối với các nhà quản lý và quy hoạch tài nguyên: Phương pháp luận và chuỗi số liệu biến động lớp phủ giai đoạn 1993 – 2000 tại vùng Trung - Hạ lưu sông Đà (Chương 4) cung cấp tài liệu cơ sở phục vụ đánh giá tác động môi trường và quy hoạch không gian lãnh thổ vùng Tây Bắc.

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

1. Bản chất sự khác biệt giữa "lớp thông tin" và "lớp phổ" trong luận án được giải thích như thế nào?

Trả lời: Lớp thông tin (Informational Class) là các hạng mục đối tượng thực tế trên mặt đất mà người sử dụng quan tâm (như đất rừng, đất nông nghiệp, mặt nước). Lớp phổ (Spectral Class) là nhóm các điểm ảnh có giá trị độ xám tương đồng nhau trong không gian phổ đa chiều. Trong thực tế, một lớp thông tin thường phân hóa thành nhiều lớp phổ khác nhau do sự khác biệt về độ ẩm đất, hướng phơi sườn dốc hoặc mật độ thực vật. Nếu ép buộc thuật toán phân loại trực tiếp theo lớp thông tin mà không tách thành các lớp phổ trung gian sẽ dẫn đến sai số nhầm lẫn lớn.

2. Ngưỡng phân loại trong phương pháp Hợp lý tối đa (Maximum Likelihood) được xác định bằng công cụ toán học nào?

Trả lời: Luận án xác định ngưỡng phân loại dựa trên biểu thức khoảng cách $(\mathbf{x} - \mathbf{m}_i)^t \mathbf{\Sigma}_i^{-1} (\mathbf{x} - \mathbf{m}_i)$, biểu thức này có quy luật phân bố $\chi^2$ (Chi-square). Bằng cách tra bảng phân bố $\chi^2$ với số bậc tự do tương ứng với số kênh ảnh $N$, tác giả xác định được giá trị ngưỡng $T_i$ để giữ lại tỷ lệ phần trăm điểm ảnh mong muốn gần tâm phân bố (ví dụ 95%) và đưa các điểm ảnh sai lệch ngoài rìa vào nhóm chưa phân loại.

3. Tại sao phương pháp Khoảng cách tối thiểu (Minimum Distance) được coi là trường hợp riêng của Maximum Likelihood?

Trả lời: Luận án chứng minh bằng toán học rằng khi ma trận hiệp biến của tất cả các lớp đều là ma trận đường chéo đồng nhất ($\mathbf{\Sigma}_i = \sigma^2 \mathbf{I}$, tức không có tương quan giữa các kênh ảnh và phương sai bằng nhau trên mọi kênh) đồng thời xác suất tiên nghiệm của các lớp bằng nhau, thì hàm phân tách của phương pháp Maximum Likelihood triệt tiêu các thành phần logarit và trở về dạng hàm tuyến tính tỉ lệ nghịch với bình phương khoảng cách Euclid của phương pháp Minimum Distance.

4. Ba nhóm kỹ thuật lọc sau phân loại (post-classification) được luận án đề cập gồm những gì?

Trả lời: Ba nhóm kỹ thuật gồm: (1) Lọc đa số (Majority filter) nhằm thay thế các điểm ảnh dị biệt đơn lẻ bằng giá trị chiếm đa số trong cửa sổ lân cận; (2) Lọc Sieve nhằm loại bỏ các cụm điểm ảnh nhỏ hơn một kích thước diện tích ngưỡng định trước; (3) Lọc hình thái học toán học (Morphological filters) bao gồm phép giãn mở (Dilation) và phép xói mòn (Erosion) với các phần tử cấu trúc khác nhau để làm mịn ranh giới các mảng phân loại.

5. Xu hướng biến động tài nguyên chính tại vùng Trung - Hạ lưu sông Đà giai đoạn 1993 – 2000 là gì?

Trả lời: Kết quả phân tích đa thời gian chỉ ra rằng diện tích đất chưa sử dụng giảm mạnh do được khai thác chuyển đổi. Diện tích đất có rừng gia tăng đáng kể nhờ chính sách phục hồi và trồng rừng; đất canh tác nông nghiệp và đất thổ cư mở rộng ở các khu vực thung lũng và ven trục giao thông, thể hiện sự gia tăng cường độ hoạt động kinh tế - xã hội trong khu vực lưu vực sông Đà.

Kết luận

Luận án tiến sĩ của tác giả Lại Anh Khôi đã giải quyết trọn vẹn cả hai phương diện lý thuyết nhận dạng mẫu trong viễn thám và ứng dụng thực tiễn theo dõi tài nguyên thiên nhiên. Bằng việc xây dựng 5 nguyên tắc cốt lõi và quy trình phân tích ảnh tích hợp đa bước, luận án đã nâng cao độ tin cậy của phương pháp phân loại ảnh số viễn thám. Kết quả thực nghiệm giai đoạn 1993 – 2000 tại vùng Trung - Hạ lưu sông Đà cung cấp luận cứ khoa học và bộ số liệu xác thực về diễn biến tích cực của lớp phủ bề mặt, phục vụ trực tiếp cho công tác quản lý tài nguyên và quy hoạch lãnh thổ khu vực Tây Bắc.