Nói tóm lại, khoảng cách có thể bị lỗi lên tới khoảng 22km hoặc 0,3%, tùy thuộc vào các điểm trong câu hỏi. Đó là:
Lỗi có thể được biểu thị theo nhiều cách tự nhiên, hữu ích , chẳng hạn như lỗi (i) (dư), bằng với chênh lệch giữa hai khoảng cách tính toán (tính bằng km) và (ii) lỗi tương đối, bằng với chênh lệch chia cho Giá trị "chính xác" (ellipsoidal). Để tạo ra các số thuận tiện để làm việc, tôi nhân các tỷ lệ này với 1000 để biểu thị sai số tương đối tính theo phần nghìn .
Các lỗi phụ thuộc vào các điểm cuối. Do sự đối xứng quay của ellipsoid và hình cầu và đối xứng hai bên (bắc-nam và đông-tây) của chúng, chúng ta có thể đặt một trong những điểm cuối ở đâu đó dọc theo kinh tuyến gốc (kinh độ 0) ở bán cầu bắc (vĩ độ từ 0 đến 90 ) và điểm cuối khác ở bán cầu đông (kinh độ từ 0 đến 180).
Để khám phá những phụ thuộc này, tôi đã vẽ các lỗi giữa các điểm cuối tại (lat, lon) = (mu, 0) và (x, lambda) là một hàm của vĩ độ x trong khoảng từ -90 đến 90 độ. (Tất cả các điểm trên danh nghĩa ở độ cao ellipsoid bằng 0.) Trong các hình, các hàng tương ứng với các giá trị của mu ở {0, 22,5, 45, 67,5} độ và các cột đối với các giá trị của lambda ở {0, 45, 90, 180} độ. Điều này cho chúng ta một cái nhìn tốt về phổ của các khả năng. Theo dự kiến, kích thước tối đa của chúng xấp xỉ bằng phẳng (khoảng 1/300) lần so với trục chính (khoảng 6700 km), hoặc khoảng 22 km.
Lỗi
Lỗi tương đối
Đường viền
Một cách khác để hình dung các lỗi là sửa một điểm cuối và để cho điểm khác thay đổi, tạo đường viền cho các lỗi phát sinh. Ở đây, ví dụ, là một đường đồng mức trong đó điểm cuối đầu tiên ở 45 độ vĩ bắc, 0 độ kinh độ. Như trước đây, các giá trị lỗi được tính bằng km và lỗi dương có nghĩa là phép tính hình cầu quá lớn:
Nó có thể dễ đọc hơn khi được bao quanh toàn cầu:
Dấu chấm màu đỏ ở phía nam nước Pháp cho thấy vị trí của điểm cuối đầu tiên.
Đối với bản ghi, đây là mã Mathicala 8 được sử dụng để tính toán:
WGS84[x_, y_] := GeoDistance @@ (GeoPosition[Append[#, 0], "WGS84"] & /@ {x, y});
sphere[x_, y_] := GeoDistance @@
(GeoPosition[{GeodesyData["WGS84", {"ReducedLatitude", #[[1]]}], #[[2]], 0}, "WGS84"] & /@ {x, y});
Và một trong các lệnh vẽ:
With[{mu = 45}, ContourPlot[(sphere[{mu, 0}, {x, y}] - WGS84[{mu, 0}, {x, y}]) / 1000,
{y, 0, 180}, {x, -90, 90}, ContourLabels -> True]]