![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
お世話になっております。ExcelVBAを使って2か所の
緯度と経度から距離を測定したいのですが上手くいきません。
データ(緯度と経度)が入っているセルはデータ型で小数点7ケタまで入っています。
これをVBAを使用し数値データからシリアル値に変換し計算をしたいのです。
緯度と経度から距離を求める式はヒュベニの計算式を使用しています。
公式は距離=sqrt((M*dP)*(M*dP)+(N*cos(P)*dR)*(N*cos(P)*dR))です。
P=2点の平均緯度 radians((緯度01+緯度02))*24
dP=2点の緯度差 (緯度01-緯度02)*24
dR=2点の経度差 (経度01-経度02)*24
M=子午線曲率半径 6334834/sqrt((1-0.006674*sin(P)*sin(P))^3)
N=卯酉線曲率半径 6377397/sqrt(1-0.006674*sin(P)*sin(P))
作成した関数は
Function kyori(X1,Y1,X2,Y2)
Dim Latit01 As Double
Dim Longi01 As Double
Dim Latit02 As Double
Dim Longi01 As Double
Dim LatitXX1 As Date
Dim LongiYY1 As Date
Dim LatitXX2 As Date
Dim LongiYY2 As Date
Dim D As Double
Dim P As Double
Dim dP As Double
Dim dR As Double
Dim M As Double
Dim N As Double
'値の代入
Latit01 = X1
Latit02 = X2
Longi01 = Y1
Longi02 = Y2
'シリアル値に変換
LatitXX1 = format((Latit01 /24), "[h]:mm:ss.000")
LongiYY1 = format((Longi01 /24), "[h]:mm:ss.000")
LatitXX2 = format((Latit02 /24), "[h]:mm:ss.000")
LongiYY2 = format((Longi02 /24), "[h]:mm:ss.000")
P = (Application.WorksheetFunction.Radians((LatitXX1 + LatitXX2) / 2)) * 24
dP = (LatitXX1 - LatitXX2) * 24
dR = (LongiYY1 - LongiYY2) *24
M = 6334834 / Application.WorksheetFunction.sqrt((1-0.006674 * Sin(P) * Sin(P))^3)
N = 6377397 / Application.WorksheetFunction.sqrt((1-0.006674 * Sin(P) * Sin(P))
kyori = Application.WorksheetFunction.sqrt((M*dP*Application.WorksheetFunction.Pi()/180)*(M*dP*Application.WorksheetFunction.Pi()/180)+(N*sos(P)*dR*Application.WorksheetFunction.Pi()/180)*(N*cos(P)*dR*Application.WorksheetFunction.Pi/180))
これを実行するとシリアル値に変換するときに止まってしまいます。
エラーは出ず、Excekのセルには#VALUEと返されます。
どなたか教えてください。
No.3ベストアンサー
- 回答日時:
ANo.2です。
> Excelでの結果は#NAME?が返されました。
>
> またエディターでブレイクポイントを打っても止まりませんでした。
こちらの環境(Excel2010)では、正常に値を返しますしブレイクポイントも効きます。
コードは標準モジュールに入れていますか?あと、Excelのバージョンは何でしょう?
![「ExcelVBAでi緯度と経度から距離を」の回答画像3](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/4/1241782_5497ea2406f72/M.jpg)
mt2008さま
ありがとうございます。
出来ました。
私の環境はExcel2010で標準モジュールにも書いておりましたが
先ほどまでできませんでしたが今、できました。
お騒がせしてすいませんでした。
ありがとうございます。
No.2
- 回答日時:
Date型変数に文字列を入れようとしてたり、かっこの数が合わなかったりして解りにくかったので、ざっくりと作り直してみました。
ちょっと自信がないですが大体こんな感じで良いんじゃないでしょうか。
Function kyori(X1, Y1, X2, Y2)
Dim D As Double
Dim P As Double
Dim dP As Double
Dim dR As Double
Dim M As Double
Dim N As Double
With WorksheetFunction
LatitXX1 = (X1 / 180) * .Pi
LongiYY1 = (Y1 / 180) * .Pi
LatitXX2 = (X2 / 180) * .Pi
LongiYY2 = (Y2 / 180) * .Pi
P = (LatitXX1 + LatitXX2) / 2
dP = (LatitXX1 - LatitXX2)
dR = (LongiYY1 - LongiYY2)
M = 6334834 / Sqr((1 - 0.006674 * Sin(P) * Sin(P)) ^ 3)
N = 6377397 / Sqr((1 - 0.006674 * Sin(P) * Sin(P)))
kyori = Sqr((M * dP) * (M * dP) + (N * Cos(P) * dR) * (N * Cos(P) * dR))
End With
End Function
この回答への補足
mt2008さま
ありがとうございます。
教えていただいた関数を実行してみました。
Excelでの結果は#NAME?が返されました。
またエディターでブレイクポイントを打っても止まりませんでした。
すいません。
No.1
- 回答日時:
テストプログラムとデータのサンプルを!
この回答への補足
JazzCorp様
ありがとうございます。
ExcelにはX1(緯度1)が35.658599、Y1(経度1)が139.745443。
X2(緯度2)35.71007,、Y2(経度2)139.80948となっております。
ちなみにテストデータで東京タワーとスカイツリーの緯度経度です。
よろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Visual Basic(VBA) VBAプログラミング 2 2022/11/27 12:13
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
- Visual Basic(VBA) VBAプログラミング 2 2022/11/27 12:07
- Excel(エクセル) VBA オリジナル関数で選択セルの合計を作成したい 3 2023/03/19 19:45
- Visual Basic(VBA) エクセル VBAについて教えてください 2 2023/04/26 13:25
- Visual Basic(VBA) VBAで質問ですが、皆さんはどの様に導き出しているのでしょうか? 6 2022/05/03 21:53
- Visual Basic(VBA) いつもお世話になっております、VBAで教えて頂きたいのですが 2 2022/05/05 22:20
- Excel(エクセル) excel2013 色付きセルの値合計 4 2023/02/28 12:53
- Visual Basic(VBA) ExcelのVBAコードについて教えてください。 1 2023/08/08 15:45
- Visual Basic(VBA) Sheet1のA列にコードB列にメアド、Sheet2のB列にコード一覧とD列にメアド一覧があり、Sh 3 2022/10/19 11:57
このQ&Aを見た人はこんなQ&Aも見ています
関連するカテゴリからQ&Aを探す
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
緯度、経度 1分は何km?
-
Excelの自動計算で緯度・経度を...
-
EXCELで、60進法の緯度経度を10...
-
Excelでの緯度経度入力方法
-
Excelについて質問です。 画像...
-
ExcelVBAでi緯度と経度から距離...
-
1秒(角度の)辺りの距離は何...
-
緯度経度の読み方教えてください。
-
札幌の対蹠点の経度と緯度は?
-
緯度と経度や赤道線を習うのは...
-
緯度経度1秒=?m
-
北緯・緯度の読み方について
-
エクセルVBA IEを操作して自動...
-
GIS、GPS,WebGIS...
-
緯度と経度の秒とその下の単位...
-
北緯30から北緯40度の距離
-
北九州市は、思いのほか寒いと...
-
中緯度、高緯度、低緯度って具...
-
金髪で青い目の人が多い国はど...
-
緯度経度の2桁秒を3桁に変換...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
緯度、経度 1分は何km?
-
Excelでの緯度経度入力方法
-
北緯・緯度の読み方について
-
EXCELで、60進法の緯度経度を10...
-
東京(北緯34度41分、東経139度...
-
Excelについて質問です。 画像...
-
Excelの自動計算で緯度・経度を...
-
緯度と経度や赤道線を習うのは...
-
GIS、GPS,WebGIS...
-
金髪で青い目の人が多い国はど...
-
緯度経度1秒=?m
-
インドのニューデリーは鹿児島...
-
1秒(角度の)辺りの距離は何...
-
経度緯度の小数点表記への変換...
-
緯度と経度の秒とその下の単位...
-
緯度経度の読み方教えてください。
-
秒数が3ケタの緯度経度の意味は?
-
住所から直線距離を算出したい...
-
札幌の対蹠点の経度と緯度は?
-
エクセルVBA IEを操作して自動...
おすすめ情報