Блог

Я ищу функцию R (или код), эквивалентную функции MATLAB scxsc

Это дает точки пересечения для пар маленьких окружностей на сфере.

Приложение выглядит так: «судно находится в азимуте X1 и на расстоянии d1 от точки 1, а также в азимуте X2 и на расстоянии d2 от точки 2. Какова его позиция в широтах, широтах?»

Я вижу множество примеров пересечения больших кругов, но не с маленькими кругами

Я нашел ответ на https://gis.stackexchange.com/questions/48937/calculating-intersection-of-two-circles/273855#273855

NB У меня есть подозрение, что в двух (неизмененных) строках ниже…

 lat1 = rad2deg(atan2(point1[2] ,point1[1]))
lon1= rad2deg(asin(point1[3]))
 

и

 lat2 = rad2deg(atan2(point2[2] ,point2[1]))
lon2 = rad2deg(asin(point2[3]))
 

… широта и долгота инвертированы !

С этими изменениями у меня получилось намного лучше!