0 
%Great Circle Distnace -- Simplified
%% 12.18993,133.45898 %% point 1 (lat/long) 
%% 14.34243,65.12750 %% point 2 (lat/long)

%%VARIABLES%%
phi_1=12.18993; %lat_1
phi_2=14.34243; %lat_2
gam_1=133.45898; %long_1
gam_2=65.12750;  %long_2
delt_gam = abs(gam_1 - gam_2); %absoulte difference in longitudes
R_Earth = 6371000; %mean radius of the earth in meters, change to FT to get distance accordingly

%%Unsimplified Great-Circle Equation -- Breaking it up into numerator and
%%denominator sections to avoid more problems -- Spherical Case of the
%%Vincenty Formula
Numer_sec1= ((cos(phi_2))*(sin(delt_gam))^2);
Numer_sec2=(((((cos(phi_1))*(sin(phi_2)))+((sin(phi_1))*(cos(phi_2))*(delt_gam))))^2);
Denom_1= (((sin(phi_1))*(sin(phi_2)))+((cos(phi_1))*(cos(phi_2))*delt_gam));

delt_sig2=atan((sqrt(Numer_sec1+Numer_sec2))/(Denom_1));

delt_GC2=R_Earth*delt_sig2;

disp(delt_GC2)

こんにちは、現在、MatLab の Vincenty 式の球面ケースを使用して、2 つの緯度/経度ポイント間の距離を計算しようとしています。http://en.wikipedia.org/wiki/Great-circle_distanceを参照しています

そこから、上記の MatLab コードを作成しました。私は与えられた最初の方程式を試しました (より単純化されたバージョンですが、どちらも役に立たなかった) ので、Vincenty の場合を使用します。コードの先頭にリストされている 2 つの緯度/経度ポイント (10 進形式) を考えると、プログラムで 2 つのポイント間の正しい距離をまだ計算していません。何が起こっているのか分からないので、何か解決方法があれば教えてください。

よろしくお願いします。投稿を頻繁に確認して、これまでのコードについての質問にお答えできるようにします。

このウェブサイトに基づいて: http://williams.best.vwh.net/gccalc.htm距離は 7381.56km である必要があります。

以下の最初の回答は、マッピング ツールボックスを持っていることを思い出させてくれましたが、取得した結果を解釈する方法がわかりません。以下に投稿したコメントを確認してください。 [ARCLEN, AZ] = distance(LAT1,LON1,LAT2,LON2)

これは実際には機能しますが、生成された弧の長さや方位角をどうするかわかりません。

明けましておめでとうございます。

4

2 に答える 2

1

アルゴリズムをプログラミングせず、Mapping Toolbox にお金を払わずに WGS84 の答えが必要な場合は、Matlab パッケージGeodesics on an ellipsoid of Revolution をダウンロードしてください。これには、geoddistance と呼ばれる Mapping Toolbox 機能の改善が含まれます。問題を解決するには

format long;
geoddistance(12.18993,133.45898,14.34243,65.12750)
->
7381566.23351761

geoddistance の引数は度単位で、結果はメートル単位です。これは、WGS84 楕円体の計算を行います。差分楕円体を使用する場合は、5 番目の引数 [a,e] (赤道半径、離心率) を指定します。(球の場合は e = 0 に設定します。長楕円体を指定する場合は、e を純粋な虚数に設定します。|e| < 0.2 の場合、正確な答えが返されます。)

ちなみに、楕円体測地線に関するウィキペディアの記事に示されている測地線の写真の多くは、このパッケージで描かれています。

于 2014-01-03T17:13:25.220 に答える