How can I obtain the distance between two tracks of 3d GPS coordinates?

Question

Nicolas Bühler 2022-5-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1724435-how-can-i-obtain-the-distance-between-two-tracks-of-3d-gps-coordinates

回答： Adit Alware 2022-5-25

I have two matrices A and B with GPS track data points (longitude, latitude, altitude). However they do not have the same size (rows). I am looking for a way to obtain a matrix which contains the distance in 3d space between the two tracks as well as a way to obtain the distance component of each axis respectively (lon, lat, alt), namely the distance from each data point in A to it's closest data point in B.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Adit Alware 2022-5-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1724435-how-can-i-obtain-the-distance-between-two-tracks-of-3d-gps-coordinates#answer_971225

在 MATLAB Online 中打开

Hi Nicolas,

I understand that you have two matrices A and B with and 3 columns(lat, long, alt) and different number of rows.

To obtain the distance between coordinates with given latitude and longitude, you may use the "distance" function from Mapping Toolbox.

dist = distance(lat1,lon1,lat2,lon2)

The output distance will be in degrees. To convert it into km or miles you may use:

dist = deg2km(dist) % convert to kilo meters
dist = deg2nm(dist) % convert to nautical miles

To obtain the final distance accomodating altitude you may calculate the folllowing :

finaldist = sqrt(dist^2 + (alt1-alt2)^2)

Suppose that we have C number of rows in A and D number of rows in B.

Now you may run a nested for loop in a manner such that they compute a matrix R of size C X D where each entry (i,j) represents the distance between ith coordinate of A and jth coordinate of B.

Now for example, to calculate the distance component of a coordinate for latitude axis, you may find the distance between point (lat,long,alt) and (lat,0,0).

After that you may calculate the minimum distance of point in A from all the points in B using the "min" function as follows:

minDist = min(R')

The output will be a row vector where ith entry will represent minimum distance of ith coordinate in A from all the coordinates in B.

Thanks