Từ kiến thức vật lý cho thấy các loại vật liệu khác nhau sẽ phản xạ khác nhau từ bức xạ cùng độ dài của sóng. Vì vậy máy quét sẽ ghi lại từng cường độ phản xạ tín hiệu trong độ dài sóng nhất định. Điều này chứng tỏ về các biên độ khác nhau và các trường sóng của các tín hiệu nhận được. Trong trường hợp phụ thuộc vào độ dốc của trường sóng, theo nguyên tắc có thể đặt thời điểm thời gian đầu tiên và cuối cùng xung của laser đến bộ thu, điều này sẽ ảnh hưởng trực tiếp tới sai số đo khoảng cách.
Theo kết quả quét laser mặt đất cho thấy việc khôi phục mô hình 3D bề mặt địa hình, để mô tả đúng nó thì mỗi địa vật quét chụp cần phản xạ càng nhiều xung laser, số lượng xung này phụ thuộc vào độ phân giải của quá trình quét.
2.2.3.4. Các sai số trong quá trình thực hiện chuyển đổi tham chiếu dữ liệu
Ở các tiểu mục trên đã nêu ra các nhóm sai số gây ra trong quá trình quét laser mặt đất. Tuy nhiên khi triển khai ngoài thực địa, ngoài việc cố định các trạm máy quét để giảm thiểu tối đa sự rung của máy trong quá trình quét, thì các sai số liên quan đến việc thu nhận, phản hồi dữ liệu của tia quét từ các địa vật ở dạng đám mây điểm để quy về hệ tham chiếu trắc địa. Các yếu tố ảnh hưởng này sẽ được mô tả khá chi
tiết đến chất lượng dữ liệu [25] bao gồm:
1. Độ chính xác xác định tọa độ điểm quét và tín hiệu phản hồi của tia quét: Giống như với toàn đạc điện tử được mô tả và xử lý qua ma trận hàm hiệp phương sai của máy quét Cmáyquét và của tín hiệu phản hồi CBS (Backsight). Ma trận
Cmáyquét được xử lý trên máy tính bằng hai cách, phụ thuộc vào vị trí của máy quét được
xác định như thế nào. Nếu trường hợp máy quét đặt trên tâm điểm đã có tọa độ
64
trắc địa và độ cao máy quét được đo bằng thước dây, ma trận Cmáyquét có thể được thể hiện ở dạng như sau:
Cmáyquét = C0 + CH
Ở đây C0 - là ma trận hàm hiệp phương sai của các tọa độ điểm trạm quét được định tâm, CH - là ma trận hàm hiệp phương sai của vị trí độ cao của tâm máy quét so với điểm gốc trạm quét
C = (0, 0, σ2 )
Ở đây: σH - là sai số xác định độ cao của tâm máy quét đo bằng thước.
Nếu vị trí của máy quét được xác định với độ chính xác tương đương như toàn đạc điện tử thì ma trận Cmáyquét tính toán sẽ đạt tương đương như sai số đo thông thường. Việc đo khoảng cách nghiêng, đo hướng ngang và đo độ dài thiên đỉnh tương tự như đo toàn đạc điện tử trong hệ tọa độ vuông góc.
2. Độ chính xác đo góc phương vị (K) từ trạm quét tới tiêu quét phản hồi:
σ = ±
Ở đây: σhor - là độ chính xác của tọa độ phẳng của máy quét và tọa độ điểm tiêu quét phản hồi là như nhau.
0− - là khoảng cách từ điểm trạm quét tới điểm tiêu quét phản hồi. Khoảng cách này được tính giá trị trung bình và chưa được hiệu chỉnh. Độ
chính xác đo góc phương vị tính theo đơn vị Radian. 3. Những sai số cài đặt thiết bị quét:
Việc cân bằng thiết bị quét theo chiều ngang và chiều đứng phải đảm bảo các sai số cho phép theo công thức sau:
Cân bằng máy quét:
σ = ±0.2 − sai số theo góc đứng
và
σ ℎ = ±σ tan − sai số theo góc ngang
Ở đây là bậc chia độ nhạy của ống bọt nước cân bằng máy quét hoặc bộ cân
bằng tự động; là giá trị đo góc đứng.
Độ chính xác định tâm quang học máy quét và tiêu quét
Ở đây: σ - là độ chính xác xác định tâm quang học (σ = 0,5mm/m giới hạn)
Độ chính xác tìm điểm tiêu quét phản hồi (với ống kính phóng quang học)
σ
Ở đây: - là độ phóng đại của ống kính.
Trường hợp máy quét không được trang bị ống kính phóng đại quang học trong quá trình tìm điểm quét thì việc xác định điểm quét sẽ hiện lên trong quá trình quét. Sai số này (tìm điểm quét) được bổ sung xử lý trong quá trình định hướng phẳng ngang trong chuyển đổi tham chiếu.
Ma trận hàm hiệp phương sai Cset các sai số cài đặt thiết bị máy quét được biểu thị như sau [25].
Ở đây các thành phần đường chéo chính là các giá trị phương sai của biên độ giá trị đo (coi như không có các sai số gốc cài đặt máy) cân bằng hướng, góc đứng, định tâm, tìm điểm tương ứng. Tổng độ chính xác khi trực tiếp thực hiện biến đổi tham chiếu trong các điểm thuộc đám mấy điểm được thực hiện mô tả theo ma trận hàm hiệp phương sai:
= C
sca nne r
Ma trận Rie( ) là ma trận quay giữa hai hệ thống, xác định theo công thức:
Trong đó: X máy quét là tọa độ gốc của trạm quét trong hệ tọa độ trắc địa; Xe là tọa độ của tiêu quét trong hệ tọa độ trắc địa; còn Xi là tọa độ của tiêu quét trong hệ tọa độ trạm quét. Và:
=
C inst
Trong đó: ; ; là các giá trị biên độ đo khoảng cách, góc ngang, góc đứng đối với điểm j tương ứng của đám mây điểm.
Tương ứng với giá trị σ2; σ2 ; σ2 là phương sai của các trị đo tương ứng. Các giá trị này đạt được tới hạn theo yêu cầu kĩ thuật
sản xuất của thiết bị. Theo các hãng và riêng σ là giá trị phương sai của trị đo chùm tia thu nhận được quy về sai số trong quá
trình thực hiện biến đổi tham chiếu dữ liệu.
Nếu định hướng tia quét được xác định rõ trong quá trình quét tiêu phản hồi thì vấn đề xác định tâm của tiêu phản hồi xảy ra, bởi tọa độ của tâm sẽ được ước tính từ số lượng của các tia laser phản xạ về từ bề mặt của tiêu quét.
Việc thu nhỏ tiêu quét phẳng đưa ra 3 hướng giải quyết xác định tâm của tiêu quét [39], đó là:
- Định vị nơi tín hiệu cực đại (max)/ khỏe (mạnh) nhất của tần số phản xạ. - Phương pháp radio xác định tâm của tất cả tia phản xạ.
- Phương pháp xác định 4 tia phản xạ mạnh nhất.
Tất cả các hướng giải quyết trên đều được dựa trên cơ sở là cường độ cực đại được ghi tại tâm tiêu phản xạ. Điều này đã được hãng Leica thiết kế và sản xuất giải quyết vấn đề trên. Ngoài ra có thể tự động hóa nhận biết tiêu quét bằng thuật toán phần mềm xử lý dữ liệu quét trên máy tính [39].
Các tọa độ của tiêu quét được ghi theo các hướng xác định bằng giá trị trung bình trọng số:
X =
Ở đây: XC - là véc tơ tọa độ của tâm tiêu quét
XJ - là véc tơ tọa độ của tia thứ j phản hồi
67
IJ - là cường độ của chính tia phản hồi thứ j
n - số lượng tia phản hồi (gần nhất n=4)
Đây là cơ sở để hoàn thiện và hiện đại hóa chế tạo thiết bị hệ thống quét laser mặt đất như Leica HDS_2500
Việc xác định các tọa độ tâm tiêu quét được xử lý bằng phần mềm và máy vi tính sử dụng phương pháp bình sai số bình phương nhỏ nhất.
Ma trận hiệp phương sai của tâm các tia phản chiếu thể hiện như sau:
Và tổng thể ma trận hàm hiệp thương sẽ có:
Việc lựa chọn ước tính độ chính xác phương sai σ cho xác định tâm tiêu quét phản chiếu [25] được tính như sau:
σ = ±
Ở đây ∆ là khoảng tái chia mẫu góc của tia quét, nó tương đương với định hướng cân bằng ngang và đứng máy quét. Ngoài giá trị
σ việc xác định vị trí tâm
tiêu quét còn phụ thuộc vào sai số xác định độ rộng của chùm tia quét σ và quan trắc tiếng ồn.
Để nâng cao độ chính xác xác định tâm tiêu quét cần tăng số lượng điểm tia quét phản hồi mà việc ước tính giá trị ∆ cho trước là nhỏ. Việc nghiên cứu khảo sát những yếu tố sai số của các tia phản chiếu trong đám mây điểm để đưa ra các giải pháp xử lý. Khi thực hiện chuyển đổi tọa độ dữ liệu sẽ góp phần nâng cao chất lượng và làm phong phú hơn cơ sở lý thuyết về quét laser mặt đất.
2.2.3.5. Sai số trong quá trình chuyển đổi khuôn dạng dữ liệu
Cũng như các quy trình xử lý số liệu đồ họa khác, việc chuyển đổi (export hoặc import) số liệu quét laser 3D đến hoặc đi từ các phần mềm khác nhau cũng có khả năng gây ra những sai số hoặc làm thay đổi chất lượng dữ liệu bởi mỗi phần mềm ứng dụng có những phương thức, định dạng và cơ chế hoạt động khác nhau.
Tới thời điểm hiện tại, số liệu đám mây điểm thường có các định dạng chuẩn như e57, ptx, pts, fls, lfc… đặc biệt mỗi hãng sản xuất máy quét lại tạo ra một định dạng chuẩn cho riêng mình, các phần mềm xử lý số liệu của nhà sản xuất thứ ba muốn sử dụng số liệu bắt buộc phải thực hiện quy trình chuyển đổi, theo đó sẽ có những ảnh hưởng nhất định tới chất lượng cũng như tính toàn vẹn của cả khối số liệu đám mây điểm tổng thể đã tạo ra. Đây cũng là điểm cần hết sức lưu ý trong quá trình xử lý và trao đổi dữ liệu giữa các phần mềm ứng dụng.
Tiểu kết Chương 2
Công nghệ quét laser mặt đất hoạt động trên nền tảng của khoa học kỹ thuật và công nghệ tin học phát triển, nó khẳng định sự vượt trội so với các công nghệ truyền thống khác và góp phần tăng cường, hiện đại hóa ngành đo đạc, bản đồ và thông tin địa lý;
Nội dung nghiên cứu của chương này góp phần bổ sung, hoàn thiện một số vấn đề về cơ sở lý thuyết của quét laser mặt đất phục vụ việc giảng dạy, nghiên cứu khoa học cũng như phát triển những ứng dụng vào thực tiễn sản xuất ở nước ta;
Thông qua việc phân tích, đánh giá các nguồn sai số ảnh hưởng đến kết quả quét laser mặt đất nhận thấy sai số trong quá trình thực hiện chuyển đổi tham số dữ liệu có ảnh hưởng lớn và mang tính hệ thống đến độ chính xác của đám mây điểm. Từ đó làm cơ sở đưa ra giải pháp xử lý yêu cầu về chất lượng dữ liệu thu thập và lựa chọn thiết bị phù hợp khi thực hiện những mục đích công việc khác nhau.
Dữ liệu thu được từ quá trình quét laser mặt đất lớn và đa dạng, để phục vụ cho mục đích địa hình và phi địa hình thì cần phải có chương trình phụ trợ cho phép tự động hóa quá trình xử lý dữ liệu đám mây điểm, đảm bảo chất lượng mô hình. Vấn đề này được giải quyết ở chương tiếp theo.
69
Chương 3.
XÂY DỰNG CHƯƠNG TRÌNH XỬ LÝ DỮ LIỆU QUÉT LASER 3D MẶT ĐẤT TRONG LĨNH VỰC ĐỊA HÌNH
Trong chương này, trình bày công tác thành lập mô hình số độ cao (DEM) bằng công nghệ quét laser 3D mặt đất theo quy định đã ban hành, ở đây tác giả chú trọng vào quá trình thực hiện định hướng ngoài các tia quét để đảm bảo độ chính xác xây dựng mô hình số địa hình. Đồng thời tiến hành xây dựng một số MO-DUL, chương trình phụ trợ cho phép tối ưu hoá và làm chủ việc tự động hóa quá trình xử lý dữ liệu đám mây điểm nhằm đẩy nhanh tiến độ, nâng cao chất lượng cũng như hiệu quả trong sản xuất.
3.1. Đặt vấn đề
Việc thành lập mô hình số địa hình (DTM - Digital Terrain Model) là phần cốt lõi trong lĩnh vực địa hình. DTM là mô hình mà ở đó người ta dùng số để mô tả bề mặt các khu đo bằng các điểm và các đường.
Thông thường thì sự thay đổi độ cao địa hình được thể hiện bằng một loạt đường đồng mức mà các điểm trên một đường đồng mức có cùng một giá trị độ cao. Các đường này là đường cong khép kín mà trong hệ thống thông tin địa lý (GIS -
Geographic Information Systems) người ta gọi là các polygon. Bằng phương pháp
này thì yếu tố địa hình cũng được thể hiện và lưu trữ trong GIS như trong các bản đồ số chuyên dùng khác [71].