Giải quyết sự khơng đồng nhất có thể trong điều kiện nếu như xác định được vị trí hai hàm cực đại tín hiệu đầu vào. Chất lượng giải quyết bài toán này phụ thuộc vào sự kéo dài của xung laser và khoảng cách giữa hai đối tượng địa vật A và B (xem hình 2.15). Khoảng cách giữa hai đối tượng cực tiểu, trong đó có thể chia tách một hàm cực đại với hàm kia có thể gọi là phương pháp giải quyết tách đôi hai đối tượng, phụ thuộc vào sự kéo dài xung của nguồn laser.
Cấu tạo và màu của đối tượng quét trực tiếp ảnh hưởng đến độ mạnh của tín hiệu phản xạ, từ đó dẫn đến thay đổi hình dáng của hàm tín hiệu phản xạ. Ví dụ hình 2.16 chỉ rõ, khi quét đối tượng phẳng, phần có màu trắng, phần khác có màu đen, kết quả chụp quét sẽ có thành phần nhiễu ồn lớn tại ranh giới hai màu.
Quá trình truyền xung laser vào khơng gian rất phức tạp. Hàm tín hiệu vào có dạng chung như sau:
𝐹𝑣à𝑜(t) =∫−∆𝜑∆𝜑 ∫−∆𝜐∆𝜐 [𝐹𝑟𝑎(𝑡 −2𝐷(𝜑,𝜐)
𝐶 ) 𝐾𝑝ℎ.𝑥𝑎(𝜑, 𝜐) 𝜓 (𝜑, 𝜐)]d𝜑. D𝜐 + η(t) (2.23) Ở đây: Fvào (t) – là hàm tín hiệu vào;
D (φ, υ) – là hàm mô tả thay đổi từ nguồn phát laser đến các điểm khác nhau của đối tượng địa vật (phản ảnh hình dáng địa vật);
Kph.xa (φ, υ) – là hàm mang đặc tính hệ số phản xạ tín hiệu từ các điểm khác nhau của địa vật;
ψ (φ, υ) - là hàm mô tả phân bố tín hiệu trong khơng gian, đại lượng phân tán tia laser;
ɳ(t) - là hàm mang đặc tính thành phần nhiễu ồn của kết quả quét. Sử dụng cơng thức (2.23) có thể thực hiện mơ hình hóa tốn học q trình chụp qt laser và nghiên cứu ảnh hưởng của các đặc tính đối tượng địa vật như hình dáng, cấu trúc, màu tới độ chính xác nhận được những đặc tính đó.
Kinh nghiệm thực tiễn chụp quét laser cho thấy, sai số đo khoảng cách liên quan đến các tính chất đo của đối tượng địa vật, có thể tới hàng mét (nếu sự phân tán tia laser lớn hơn 3,5 phút), thường trung bình đạt khoảng 1–3 cm.
Như vậy phần lớn các vấn đề liên quan đến sự ảnh hưởng của cấu trúc, hình dáng, định vị và màu sắc của đối tượng địa vật đều ảnh hưởng đến kết quả quét. Có thể tránh nếu chúng ta giảm độ phân tán tia laser, tức là làm tăng mức độ chuẩn hội tụ, tốt nhất là tăng chất lượng đặc tính đo của quét laser.
Việc tăng độ chính xác quét laser mặt đất đối với các đối tượng địa vật giải quyết tốt nhất bằng các phương pháp tính tốn mơ hình hóa q trình chụp qt laser đối với đối tượng cụ thể với việc sử dụng công thức (2.23).
2.2.3.3. Ảnh hưởng các tham số quét và các đặc tính máy quét đến độ chính xác kết quả đo vẽ
Khi thực hiện quét laser mặt đất cần thiết phải lựa chọn thiết bị, các thông số kỹ thuật và quy trình cơng nghệ để cho phép xây dựng mơ hình số địa hình và mơ hình bề mặt theo yêu cầu về độ chính xác cũng như mức độ chi tiết.
Xác định các đặc tính của trạm quét laser mặt đất bao gồm: - Độ phân giải quét đứng và quét ngang cực đại;
- Độ dài sóng bức xạ laser;
- Độ chính xác đo cạnh, đo các góc đứng và ngang; - Độ phân tán của tia laser.
Nếu tiến hành với các hệ thống quét chụp laser mặt đất thì ba đặc tính đầu tiên nêu trên sẽ xác định độ chính xác hình học của dịng qt, đặc tính sau sẽ là độ phân giải của hệ thống, chính đại lượng này cho thấy kích thước nhỏ nhất của địa vật khi ở khoảng cách xa định trước từ máy qt để tìm được hàm tín hiệu ra cực đại duy nhất. Độ phân giải của hệ thống ở mức độ thấp phụ thuộc vào độ phân giải quét, càng tránh phân tán tia laser, hơn nữa trong khi giảm giá trị góc của bước qt thì vùng phủ giữa hai đối tượng quét sẽ được tăng lên theo hình 2.17 [9].
Như vậy độ phân giải máy quét laser cơ bản phụ thuộc vào độ phân tán tia laser. Từ hình 2.17 có thể thấy rằng độ phân giải góc có giá trị nhỏ hơn là độ phân tán tia laser. Điều này có nghĩa là khi định trước độ phân tán tia laser trong máy qt cụ thể khơng thể nhận được trên dịng qt những chi tiết nhỏ của địa vật chỉ nhờ tính đến tăng độ phân giải quét.
a). Độ phân giải cực đại khi quét b) Độ phân giải cực tiểu khi qt Hình 2.17. Phân tích độ phân giải trạm quét laser mặt đất [9]
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 qt 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 qt 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 q 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 độ
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 (2.24)
Ở đâ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) (2.25)
Ở đâ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ư tồn đạc điện tử thì ma trận Cmáyqt tính tố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 độ vng góc.
2. Độ chính xác đo góc phương vị (K) từ trạm qt tới tiêu quét phản hồi:
σ𝑘 = ±√2σℎ𝑜𝑟
𝑑0−𝐵𝑆 [𝑟𝑎𝑑] (2.26)
Ở đâ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 (2.27) và
σ𝑙ℎ = ±σ𝑙𝑣tan𝜃 − sai số theo góc ngang (2.28)
𝑀
Ở đây: 𝑀 - là độ phóng đại của ống kính.
Trường hợp máy qt khơng được trang bị ống kính phóng đại quang học trong q trình tìm điểm qt 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]. C𝑠𝑒𝑡 = [ 0 0 0 0 σ𝑙ℎ2 + σ𝑐𝑒𝑛𝑡2 + σ𝑝𝑜𝑖𝑛𝑡2 0 0 0 σ𝑙𝑣2 ] (2.31)
Ở đâ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:
𝐶𝑋𝑒 = Cscanner+ [Rie(𝑘)J(Cinst+ Cset)JT[Rie(𝑘)]T+ σ𝑘2 ∂Xe
∂𝑘 (∂Xe
∂𝑘)T] (2.32) Ma trận Rie(𝑘) là ma trận quay giữa hai hệ thống, xác định theo công thức:
Xe = X máy quét + Rie(𝑘)Xi (2.33)
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à:
𝐽 = [ ∂𝑥𝑗 ∂𝑟𝑗 ∂𝑥𝑗 ∂𝜑𝑗 ∂𝑥𝑗 ∂𝜃𝑗 ∂𝑦𝑗 ∂𝑟𝑗 ∂𝑦𝑗 ∂𝜑𝑗 ∂𝑦𝑗 ∂𝜃𝑗 ∂𝑧𝑗 ∂𝑟𝑗 ∂𝑧𝑗 ∂𝜑𝑗 ∂𝑧𝑗 ∂𝜃𝑗] (2.34) Cinst = 𝑑𝑖𝑎𝑔(σ𝑟2, σ𝜑2 + σ𝑏𝑒𝑎𝑚2 , σ𝜃2 + σ𝑏𝑒𝑎𝑚2 ) (2.35)
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. Ngồi ra có thể tự động hóa nhận biết tiêu qt bằng thuật tố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𝑐 =∑ 𝐼𝑗X𝑗
𝑛 𝑗=1
∑𝑛𝑗=1𝐼𝑗 (2.36) Ở đây: XC - là véc tơ tọa độ của tâm tiêu qué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: 𝐶𝑋𝑗 = 𝐽Cinst𝐽𝑇 = [ σ𝑥2𝑗 σ𝑥2𝑗𝑦𝑗 σ𝑥2𝑗𝑧𝑗 σ𝑦2𝑗𝑥𝑗 σ𝑦2𝑗 σ𝑦2𝑗𝑧𝑗 σ𝑧𝑗𝑥𝑗2 σ𝑧𝑗𝑦𝑗2 σ𝑧𝑗2 ] (2.37) Và tổng thể ma trận hàm hiệp thương 𝐶𝑋𝑐 sẽ có: 𝐶𝑋𝑐 =∑ 𝐼𝑗 2C𝑋𝑗 𝑛 𝑗=1 (∑𝑛𝑗=1𝐼𝑗)2 (2.38)
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: σ𝑐 = ± ∆
2√3 (2.39)
Ở đâ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 qt 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.
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 tồ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, hồ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 q 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.
chương trình phụ trợ cho phép tối ưu hố và làm chủ việc tự động hóa q 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 đồ