Làm thế nào để bạn tính toán bán kính trái đất ở một vĩ độ trắc địa nhất định?


19

(Tôi thấy có một phương trình trên wikipedia thực hiện chính xác những gì tôi đang hỏi nhưng không có tài liệu tham khảo. Tôi không có cách nào để xác nhận tính hợp lệ của phương trình này!)

Tôi đã hiểu sự khác biệt giữa Vĩ độ địa tâm so với Vĩ độ trắc địa.

Giả sử bán kính đã biết chính, avà bán phụ, bđược đưa ra. Làm thế nào để bạn tính toán bán kính ở một vĩ độ trắc địa nhất định?

Tôi cần một số loại xác nhận chuyên gia (phái sinh, liên kết đến phái sinh, xác nhận từ chuyên gia, giải thích, vv).

Câu trả lời:


25

Câu hỏi này giả định một mô hình ellipsoidal của trái đất. Bề mặt tham chiếu của nó thu được bằng cách xoay một hình elip xung quanh trục nhỏ của nó (được vẽ theo chiều dọc theo quy ước). Một hình elip như vậy chỉ là một vòng tròn được kéo dài theo chiều ngang bởi một yếu tố a và theo chiều dọc bởi một yếu tố b . Sử dụng tham số tiêu chuẩn của vòng tròn đơn vị,

t --> (cos(t), sin(t))

( định nghĩa cosine và sin), chúng ta có được một tham số hóa

t --> (a cos(t), b sin(t)).

(Hai thành phần của tham số hóa này mô tả một chuyến đi quanh đường cong: chúng chỉ định, theo tọa độ của Cartesian, vị trí của chúng tôi tại "thời gian" t .)

Các vĩ độ trắc địa , f , của bất kỳ điểm nào là góc mà "lên" làm cho với mặt phẳng xích đạo. Khi a khác với b , giá trị của f khác với giá trị của t (ngoại trừ dọc theo đường xích đạo và ở hai cực).

Nhân vật

Trong ảnh này, đường cong màu xanh là một góc phần tư của hình elip như vậy (được phóng đại rất nhiều so với độ lệch tâm của trái đất). Dấu chấm màu đỏ ở góc dưới bên trái là trung tâm của nó. Đường đứt nét chỉ định bán kính đến một điểm trên bề mặt. Hướng "lên" của nó được hiển thị với một đoạn màu đen: theo định nghĩa, nó vuông góc với hình elip tại điểm đó. Do độ lệch tâm phóng đại, dễ dàng thấy rằng "lên" không song song với bán kính.

Theo thuật ngữ của chúng tôi, t liên quan đến góc tạo bởi bán kính so với phương ngang và f là góc tạo bởi đoạn màu đen đó. (Lưu ý rằng bất kỳ điểm nào trên bề mặt đều có thể được xem từ phối cảnh này. Điều này cho phép chúng tôi giới hạn cả tf nằm trong khoảng từ 0 đến 90 độ; cosin và sin của chúng sẽ dương, vì vậy chúng tôi không phải lo lắng về tiêu cực căn bậc hai trong công thức.)

Bí quyết là để chuyển đổi từ các t -parameterization đến một trong các điều khoản của f , bởi vì trong điều kiện của t bán kính R là dễ dàng để tính toán (thông qua định lý Pythagore). Hình vuông của nó là tổng bình phương của các thành phần của điểm,

R(t)^2 = a^2 cos(t)^2 + b^2 sin(t)^2.

Để thực hiện chuyển đổi này, chúng ta cần liên kết hướng "lên" f với tham số t . Hướng này vuông góc với tiếp tuyến của hình elip. Theo định nghĩa, một tiếp tuyến của một đường cong (được biểu thị dưới dạng vectơ) có được bằng cách phân biệt tham số hóa của nó:

Tangent(t) = d/dt (a cos(t), b sin(t)) = (-a sin(t), b cos(t)).

(Sự khác biệt tính toán tốc độ thay đổi. Tốc độ thay đổi vị trí của chúng ta khi chúng ta di chuyển quanh đường cong dĩ nhiên là vận tốc của chúng ta và luôn luôn chỉ dọc theo đường cong.)

Xoay theo chiều kim đồng hồ 90 độ để có được đường vuông góc, được gọi là vectơ "bình thường":

Normal(t) = (b cos(t), a sin(t)).

Độ dốc của vectơ bình thường này, bằng (a sin (t)) / (b cos (t)) ("tăng lên khi chạy"), cũng là tiếp tuyến của góc mà nó tạo ra theo phương ngang, từ đó

tan(f) = (a sin(t)) / (b cos(t)).

Tương đương,

(b/a) tan(f) = sin(t) / cos(t) = tan(t).

(Nếu bạn có cái nhìn sâu sắc tốt vào Euclid, bạn có thể có được mối quan hệ này trực tiếp từ định nghĩa của một hình elip mà không cần trải qua bất kỳ trang điểm hay tính toán, đơn giản bằng cách công nhận rằng mở rộng theo chiều ngang và thẳng đứng kết hợp bởi mộtb tương ứng có tác dụng thay đổi tất cả các độ dốc theo yếu tố này b / a .)

Nhìn lại công thức của R (t) ^ 2: chúng ta biết ab - họ xác định hình dạng và kích thước của hình elip - vì vậy chúng ta chỉ cần tìm cos (t) ^ 2 và sin (t) ^ 2 về mặt f , phương trình trước cho phép chúng ta thực hiện dễ dàng:

cos(t)^2 = 1/(1 + tan(t)^2) 
         = 1 / (1 + (b/a)^2 tan(f)^2) 
         = a^2 / (a^2 + b^2 tan(f)^2);
sin(t)^2 = 1 - cos(t)^2 
         = b^2 tan(f)^2 / (a^2 + b^2 tan(f)^2).

(Khi tan (f) là vô hạn, chúng ta ở cực, vì vậy chỉ cần đặt f = t trong trường hợp đó.)

Đây là kết nối chúng ta cần. Thay thế các giá trị này cho cos (t) ^ 2 và sin (t) ^ 2 vào biểu thức cho R (t) ^ 2 và đơn giản hóa để có được

R(f)^2 = ( a^4 cos(f)^2 + b^4 sin(f)^2 ) / ( a^2 cos(f)^2 + b^2 sin(f)^2 ).

Một phép biến đổi đơn giản cho thấy phương trình này giống như phương trình được tìm thấy trên Wikipedia. Bởi vì a ^ 2 b ^ 2 = (ab) ^ 2 và (a ^ 2) ^ 2 = a ^ 4,

R(f)^2 = ( (a^2 cos(f))^2 + (b^2 sin(f))^2 ) / ( (a cos(f))^2 + (b sin(f))^2 )

+1 .. ngoại trừ tôi nghĩ rằng công thức cuối cùng có một vị trí không phù hợp ... không (b^4 sin(f))^2nên thay đổi thành (b^4 sin(f)^2)?
Kirk Kuykendall

thực sự vui mừng vì có một số chuyên gia xung quanh về chủ đề này =).
Trevor Boyd Smith

Tập tin Geogebra (html) có thể được đăng lên trang web này không? Tôi có bán kính của phương thẳng đứng có thể chứng minh trực quan những gì đang diễn ra.

Bạn có thể xuất bản gốc ở định dạng .png, @Dan: sử dụng hộp thoại Tệp | Xuất. Tôi khuyên bạn nên sử dụng phông chữ lớn (16 hoặc 18 điểm dường như hoạt động tốt) và phóng to đến tận hình ảnh như bạn có thể trước tiên.
whuber

Tôi cho rằng sự tương tác sẽ bị mất sau đó. Bản demo cho thấy mức độ thay đổi bán kính và vĩ độ thay đổi của các thuộc tính.

3

Thật thú vị khi thấy rằng giải pháp toán học không biết chữ của tôi đã thực hiện công việc với 5 phút suy nghĩ và mã hóa, liệu yếu tố làm phẳng có phải được xem xét chứ không phải là một mô hình elip hoàn hảo?

        double pRad = 6356.7523142;
        double EqRad = 6378.137;                      
        return pRad + (90 - Math.Abs(siteLatitude)) / 90 * (EqRad - pRad); 

1
Trong đó pRad là bán kính cực và EqRad là bán kính xích đạo.
Stefan Steiger

đây là câu trả lời duy nhất tôi có thể đọc Nó dường như làm việc cho tôi.
Sean Bradley

1
Tôi thấy rằng bạn đang thực hiện một phép nội suy tuyến tính bán kính giữa cực và xích đạo. Mặc dù không có lý do nào để tin rằng phép nội suy tuyến tínhchính xác , tôi sẽ sử dụng điều này là "đủ tốt" cho Trái đất, do yếu tố làm phẳng nhẹ của nó. BTW Tôi nghĩ rằng việc đọc tương đương dễ dàng hơn một chút: return E + (P - E) * Abs(Lat) / 90vì vậy không cần phải có 90 - ...trong công thức.
ToolmakerSteve

2

nhập mô tả hình ảnh ở đây

Ít nhất đó là công thức tôi tìm thấy ở Phân tích và đánh giá Trung tâm Dữ liệu Mỹ (DAAC) cho Bộ Quốc phòng Mỹ (DoD) High Performance Computing hiện đại hóa chương trình (HPCMP) wiki . Nó nói rằng họ đã vay rất nhiều từ mục nhập của Wikipedia . Tuy nhiên, thực tế là họ giữ lại công thức đó nên được tính cho một cái gì đó.


Bạn có thể cung cấp một liên kết đến nội dung?
Trevor Boyd Smith

trong đó φ là vĩ độ trắc địa, và một (trục bán chính) và b (trục bán phụ) lần lượt là bán kính xích đạo và bán kính cực. var a = 6378137; // m var b = 6356752.3142; // m en.wikipedia.org/wiki/Semi-major_and_semi-minor_axes en.wikipedia.org/wiki/World_Geodetic_System
Stefan Steiger
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.