ポイントABとCの位置に基づいて、ポイントDの位置を計算したいと考えています.Dに対するポイントAの角度と、Bに対するDとDに対するCの角度を知っています.
実際には、ポイント AB と C は、私が GPS でマークした 3 つの場所であり、ポイント D は、GPS の場所を取得しようとしている電波望遠鏡を装着した動物の場所です。無線首輪を付けた動物が北に対してどの方向を向いているかを知ることで得られる角度。
私はアルゴリズムを書きましたが、GPS 座標をそのままアルゴリズムに入れることはできず、それらを変換してから再度変換する必要があることはわかっています。私はグーグルで調べてきましたが、少し混乱しています。これには、デカルトまたはUTMの使用がより適切ですか?
GPS を UTM に変換するにはどうすればよいですか? 検索しましたが、少し混乱しています。一部の変換では、度分と秒が表示されます。私の GPS では、これに追加の数値が表示されるようです。つまり、N 68.21.446 と `w 12.14.284 です。
関連する場合は、計算を少し簡単にするために、面積が2次元であると想定しました。
必要かどうかはわかりませんが、コードは次のとおりです。
#10/09/2013
#Enter your points for locations A B and C
#AN and AW is your first GPS points AA is the angle
AN<-10
AW<-0
AA<-45
#BN and BW are your second
BN<-10
BW<-0
BA<-0
#CN and CW are your third
CN<-0
CW<-10
CA<-90
#Convert these to ?
#work out distance
#For each co ordinate and angle, you need to calculate y=mx+c to make a line
#From these 3 lines, you can work out where they intersect
#If the angle is 0 it wont work, so make it very close to 0.
if(AA==0) {AA<-0.00001}
if(BA==0) {BA<-0.00001}
if(CA==0) {CA<-0.00001}
#Convert all angles to radians
AAr<-(AA*pi)/180
BAr<-(BA*pi)/180
CAr<-(CA*pi)/180
#Calculate M which is 1/tan(b)
AM<-1/tan(AAr)
BM<-1/tan(BAr)
CM<-1/tan(CAr)
#Calculate C the equation constant
#c=y-m*x
AC<-AW-AM*AN
BC<-BW-BM*BN
CC<-CW-CM*CN
#Caclulate intersections
#A and B
XAB<-(AC-BC)/(BM-AM)
YAB<-(AM*XAB+AC)
#B and C
XBC<-(BC-CC)/(CM-BM)
YBC<-(BM*XBC+BC)
#C and A
XAC<-(CC-AC)/(AM-CM)
YAC<-(CM*XAC+CC)
#Work out average of these 3 points
(XofABC<-(XAB+XBC+XAC)/(3))
(YofABC<-(YAB+YBC+YAC)/(3))
#Convert this back into GPS coordinate
`