Tổng quan nghiên cứu

Công nghệ quét laser viễn thám LiDAR (Light Detection And Ranging) là bước tiến đột phá trong việc thu nhận dữ liệu không gian ba chiều với độ chính xác mặt bằng đạt dưới 25cm và sai số độ cao dưới 20cm. Tốc độ phát xung của các hệ thống hiện đại đạt từ 150.000 đến 240.000 xung mỗi giây, cho phép đo quét diện tích 1.000 km2 chỉ trong khoảng 25 đến 30 giờ bay thực địa. Tuy nhiên, việc khai thác hiệu quả nguồn dữ liệu này gặp rào cản lớn do kích thước dữ liệu đám mây điểm quá lớn, thường dao động từ hàng triệu đến hàng trăm triệu điểm không có cấu trúc. Các giải pháp lọc tự động truyền thống hiện nay mới chỉ đạt hiệu suất xử lý khoảng 90% đến 95%, vẫn đòi hỏi can thiệp thủ công từ 5% đến 10% khối lượng công việc, làm tăng chi phí và kéo dài thời gian xử lý dữ liệu.

Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết bài toán phân loại tự động đám mây điểm LiDAR, nhằm phân tách chính xác các lớp điểm bề mặt mặt đất (Ground) và phi mặt đất (Non-Ground) phục vụ xây dựng mô hình số độ cao (DEM) và mô hình số bề mặt (DSM). Mục tiêu cụ thể là nghiên cứu cơ chế hoạt động, triển khai thực nghiệm và so sánh toàn diện hiệu năng giữa thuật toán phân cụm K-Means và mạng nơ-ron học sâu PointNet.

Phạm vi nghiên cứu được thực hiện trên các bộ dữ liệu bay quét thực tế tại khu vực Nghệ An và Quảng Ninh, kết hợp bộ dữ liệu chuẩn quốc tế Sydney Urban Objects trong giai đoạn năm 2021. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật to lớn, giúp rút ngắn thời gian thành lập bản đồ địa hình từ 10 ngày xuống còn vài giờ, đồng thời giảm thiểu sự phụ thuộc vào các gói phần mềm thương mại đắt đỏ.

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 dựa trên sự tích hợp của ba hệ thống thành phần trong công nghệ LiDAR: thiết bị phát xung laser dải sóng cận hồng ngoại, hệ thống định vị vệ tinh vi phân toàn cầu GNSS và hệ thống dẫn đường quán tính INS. Cấu trúc đồng bộ này xác định chính xác tọa độ không gian ba chiều (X, Y, Z), góc quét từ âm 90 độ đến dương 90 độ, cùng cường độ phản xạ ánh sáng với khả năng ghi nhận tối đa 4 lần xung phản hồi trên mỗi tia phát.

Mô hình lý thuyết thứ nhất là thuật toán phân cụm không giám sát K-Means, vận hành theo nguyên lý tối thiểu hóa tổng bình phương khoảng cách Euclid giữa các điểm dữ liệu tới K tâm cụm khởi tạo qua nhiều vòng lặp. Mô hình lý thuyết thứ hai là kiến trúc học sâu PointNet, giải quyết trực tiếp dữ liệu đám mây điểm thưa thớt bằng cách kết hợp các tầng mạng perceptron đa lớp (Shared MLP), mạng biến đổi không gian T-Net và hàm đối xứng Max Pooling nhằm đảm bảo tính bất biến đối với thứ tự điểm và hướng xoay không gian.

Các khái niệm chuyên ngành then chốt bao gồm: định dạng dữ liệu nhị phân chuẩn LAS của Hiệp hội Đo ảnh và Viễn thám Hoa Kỳ (ASPRS), mô hình số độ cao (DEM), mô hình số bề mặt (DSM), và kỹ thuật Voxelization chuyển đổi không gian điểm thành lưới thể tích ba chiều.

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

Dữ liệu nghiên cứu bao gồm bộ dữ liệu thực địa 485 điểm đo tại Nghệ An, tệp dữ liệu bay quét laser hàng không (ALS) định dạng LAS tại Quảng Ninh và tập dữ liệu chuẩn Sydney Urban Objects gồm 100 đối tượng đô thị phân lớp chi tiết.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu áp dụng phương pháp lấy mẫu phân tầng ngẫu nhiên với 485 điểm thử nghiệm cho thuật toán K-Means nhằm đại diện cho nhiều dạng địa hình biến thiên. Đối với mạng PointNet, mỗi đám mây điểm đối tượng được chuẩn hóa về cỡ mẫu cố định 1024 điểm nhằm tối ưu hóa độ phân giải hình học và cân bằng tải tính toán của bộ nhớ GPU.

Lý do lựa chọn phương pháp phân tích bắt nguồn từ nhu cầu so sánh giữa thuật toán gom cụm hình học truyền thống và kiến trúc mạng học sâu hiện đại. Thuật toán K-Means được xử lý trên nền tảng phần mềm IBM SPSS 20 dựa trên thuộc tính độ cao Z. Mạng PointNet được lập trình trên nền tảng TensorFlow, kết hợp ngôn ngữ C# và bộ công cụ LASTools trên môi trường Visual Studio 2015 để đọc, ghi và trực quan hóa tệp LAS. Toàn bộ lộ trình nghiên cứu và thực nghiệm được tiến hành liên tục trong thời gian 24 tháng.

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

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

Thứ nhất, thực nghiệm thuật toán K-Means trên 485 điểm mẫu địa hình tại Nghệ An cho thấy tỷ lệ điểm không gán được vào cụm (điểm lỗi) chiếm 31,55% (tương đương 153 điểm). Khi thiết lập K=5 qua 10 vòng lặp, cụm 3 chiếm tỷ trọng lớn nhất với 208 điểm (42,88%), cụm 4 có 110 điểm (22,68%), các cụm 1, 2 và 5 chỉ chiếm từ 2 đến 10 điểm. Khi tăng lên K=7, số lượng điểm lỗi vẫn giữ nguyên ở mức 153 điểm (31,55%), trong đó cụm 3 chiếm 179 điểm (36,90%) và cụm 4 chiếm 99 điểm (20,41%).

Thứ hai, mạng học sâu PointNet thể hiện hiệu năng xử lý vượt bậc khi đạt tốc độ phân loại hơn 1.000.000 điểm mỗi giây (tương đương xấp xỉ 1.000 đối tượng mỗi giây) trên phần cứng đồ họa GPU GTX 1080X, cho phép trích xuất đặc trưng hình học trực tiếp mà không gây suy giảm chất lượng dữ liệu như kỹ thuật tạo lưới Voxel.

Thứ ba, trên tập dữ liệu Sydney Urban Objects với số lượng mẫu hạn chế, mô hình phân loại PointNet đạt độ chính xác xác thực xấp xỉ 60% sau khi áp dụng kỹ thuật tối ưu hóa Adam kết hợp chuẩn hóa trọng số L2.

Thứ tư, chương trình ứng dụng xây dựng trên ngôn ngữ C# kết hợp LASTools trên dữ liệu Quảng Ninh đã phân tách thành công đám mây điểm thành hai lớp Ground và Non-Ground, thành lập mô hình DEM và DSM đạt độ phân giải không gian 1m x 1m.

Thảo luận kết quả

Tỷ lệ sai số 31,55% của K-Means bắt nguồn từ việc thuật toán phụ thuộc tuyệt đối vào khoảng cách hình học của trục độ cao Z và độ nhạy cảm với vị trí tâm cụm khởi tạo. Khi địa hình có độ dốc phức tạp, K-Means dễ rơi vào cực trị địa phương. Ngược lại, PointNet khắc phục hoàn toàn nhược điểm này nhờ cơ chế ánh xạ đặc trưng đa chiều từ 3 chiều lên 1024 chiều, kết hợp lớp gộp Max Pooling giúp trích xuất thông tin toàn cục ổn định.

So với các thuật toán lọc cổ điển như mô hình góc nghiêng di động của Vosselman hay mô hình lưới tam giác di động TIN của Axelsson, phương pháp học sâu PointNet thể hiện khả năng khái quát hóa vượt trội trên các cấu trúc hình học phức tạp. Dữ liệu thực nghiệm được trực quan hóa sinh động qua biểu đồ Histogram thể hiện sự phân bổ tách biệt rõ ràng giữa hai lớp điểm mặt đất và vật thể, cùng mô hình 3D địa hình phân giải 1m phản ánh chân thực hiện trạng khu vực khảo sát.

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

Thứ nhất, tăng cường kích thước và độ đa dạng của tập dữ liệu huấn luyện. Các đơn vị nghiên cứu cần áp dụng kỹ thuật Data Augmentation (xoay góc, co giãn, thêm nhiễu ngẫu nhiên) để mở rộng tập mẫu từ 100 đối tượng lên hơn 10.000 đối tượng trong vòng 6 tháng tới. Mục tiêu cụ thể là nâng độ chính xác xác thực của mô hình PointNet từ 60% lên trên 85% đối với các lớp đối tượng phức tạp như thảm thực vật và xe cộ.

Thứ hai, thiết kế mô hình phân loại lai kết hợp giữa lọc hình học và mạng nơ-ron học sâu. Nhóm phát triển phần mềm nên ứng dụng thuật toán lọc thô của LASTools để loại bỏ 95% điểm nhiễu ban đầu, sau đó chuyển giao các điểm biên phức tạp cho mạng C-PointNet xử lý phân lớp ngữ nghĩa chi tiết. Thời gian hoàn thiện giải pháp dự kiến trong 12 tháng do các kỹ sư hệ thống thông tin đảm nhiệm.

Thứ ba, xây dựng phần mềm xử lý dữ liệu LiDAR nội địa trên nền tảng .NET và điện toán đám mây. Doanh nghiệp công nghệ địa không gian cần tích hợp các mô hình học máy vào ứng dụng máy tính độc lập để thay thế các phần mềm thương mại đóng đắt đỏ. Giải pháp này giúp tiết kiệm khoảng 70% chi phí đầu tư công nghệ cho các cơ quan khảo sát trắc địa trong lộ trình 18 tháng.

Thứ tư, chuẩn hóa quy trình phân loại dữ liệu theo tiêu chuẩn quốc tế ASPRS. Các cơ quan quản lý nhà nước về tài nguyên môi trường cần ban hành khung phân loại tối thiểu 20 lớp đối tượng địa vật, thực hiện kiểm chuẩn chất lượng dữ liệu định kỳ 6 tháng một lần nhằm giữ tỷ lệ lỗi dữ liệu dưới 5%.

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

Thứ nhất, học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Hệ thống thông tin, Khoa học máy tính và Viễn thám. Luận văn cung cấp cơ sở toán học chặt chẽ về mạng PointNet và thuật toán K-Means, hỗ trợ trực tiếp cho việc phát triển các thuật toán phân tích dữ liệu không gian 3D.

Thứ hai, kỹ sư trắc địa, chuyên viên GIS và bản đồ số tại các cơ quan đo đạc địa chính. Nghiên cứu cung cấp tài liệu thực hành chi tiết về quy trình chuyển đổi tệp tin LAS, trích xuất điểm mặt đất và thành lập mô hình DEM/DSM độ phân giải 1m phục vụ quy hoạch hạ tầng.

Thứ ba, chuyên gia phát triển phần mềm và kiến trúc sư hệ thống tại các doanh nghiệp công nghệ địa không gian. Nhóm này có thể ứng dụng mã nguồn, giải pháp kết nối C# với LASTools để xây dựng các sản phẩm thương mại xử lý dữ liệu đám mây điểm chi phí thấp.

Thứ tư, cán bộ quản lý quy hoạch đô thị, giao thông vận tải và lâm nghiệp. Luận văn hỗ trợ ra quyết định ứng dụng công nghệ quét LiDAR trong việc giám sát biến động rừng, quản lý hành lang an toàn lưới điện và mô phỏng nguy cơ ngập lụt đô thị.

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

Công nghệ LiDAR có ưu điểm gì vượt trội so với phương pháp đo ảnh lập thể truyền thống? Công nghệ LiDAR thu nhận trực tiếp tọa độ không gian 3D với độ chính xác mặt bằng dưới 25cm và độ cao dưới 20cm mà không phụ thuộc vào điều kiện chiếu sáng. Xung laser có khả năng xuyên qua tán cây với tối đa 4 lần phản xạ để ghi nhận chính xác bề mặt đất bên dưới.

Tại sao thuật toán K-Means lại có tỷ lệ phân loại lỗi lên tới 31,55% trong thực nghiệm? Trong thực nghiệm với 485 điểm tại Nghệ An, K-Means chỉ sử dụng duy nhất giá trị độ cao Z để phân cụm. Thuật toán phân nhóm thuần túy theo khoảng cách hình học nên khi gặp các điểm dị biệt hoặc địa hình biến thiên mạnh, 153 điểm dữ liệu đã không thể gán chính xác vào các tâm cụm đại diện.

Kiến trúc mạng PointNet giải quyết tính chất không có thứ tự của đám mây điểm như thế nào? PointNet sử dụng các mạng nơ-ron đa tầng perceptron kết hợp với lớp biến đổi không gian T-Net để chuẩn hóa dữ liệu đầu vào. Sau đó, mạng áp dụng hàm đối xứng Max Pooling để tổng hợp toàn bộ các vector đặc trưng, đảm bảo kết quả đầu ra hoàn toàn bất biến đối với thứ tự sắp xếp của tập điểm.

Tệp định dạng LAS đóng vai trò như thế nào trong bài toán phân loại dữ liệu LiDAR? Tệp LAS là định dạng nhị phân tiêu chuẩn được ASPRS thiết kế riêng để lưu trữ dữ liệu đám mây điểm. Tệp quản lý đồng thời tọa độ X, Y, Z, cường độ phản xạ, số lượng xung hồi đáp và các trường mã phân loại, giúp phần mềm truy xuất và xử lý hàng triệu điểm mà không làm giảm tốc độ hệ thống.

Mô hình số độ cao DEM và mô hình số bề mặt DSM thành lập từ dữ liệu LiDAR ứng dụng trong thực tế ra sao? Mô hình DEM phản ánh địa hình bề mặt đất tự nhiên không chứa vật thể, ứng dụng chính trong phân tích dòng chảy, quy hoạch thoát nước và cảnh báo trượt lở đất. Mô hình DSM bao gồm cả thảm thực vật và công trình xây dựng, phục vụ đắc lực cho quản lý chiều cao đô thị và hành lang truyền tải điện.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết công nghệ LiDAR và cấu trúc lưu trữ dữ liệu đám mây điểm theo định dạng tiêu chuẩn LAS của ASPRS.
  • Thực nghiệm chứng minh rõ giới hạn của thuật toán K-Means khi phân loại theo trục độ cao đơn lẻ với tỷ lệ điểm lỗi chiếm 31,55% trên 485 điểm mẫu.
  • Mô hình học sâu PointNet khẳng định tính ưu việt với tốc độ xử lý đạt trên 1.000.000 điểm mỗi giây trên môi trường GPU, mở ra tiềm năng ứng dụng phân loại thời gian thực.
  • Xây dựng thành công ứng dụng C# tích hợp thư viện LASTools, tự động hóa quy trình phân tách lớp điểm và thành lập mô hình DEM, DSM độ phân giải 1m tại Quảng Ninh.
  • Đề tài đóng góp giải pháp công nghệ quan trọng, giúp tự động hóa trên 90% công đoạn phân loại dữ liệu địa không gian và cắt giảm đáng kể chi phí sản xuất bản đồ.

Đóng góp lớn nhất của công trình là làm chủ quy trình phân loại dữ liệu LiDAR bằng trí tuệ nhân tạo tại Việt Nam. Lộ trình phát triển trong 12 tháng tới sẽ tập trung tối ưu hóa kiến trúc C-PointNet trên các bộ dữ liệu quy mô lớn. Hãy tải toàn văn tài liệu nghiên cứu và áp dụng ngay các giải pháp học máy tiên tiến này vào dự án xử lý dữ liệu không gian của bạn.