Помогите с расчетом направления отрезка

[klerk184]

Точно:
φ = (φ₁ + φ₂) / 2

Спасибо огромное! Буду сейчас пробовать делать!) А то уже чего только не пытался сделать)))

[klerk184]

Чтото не совсем то получается. Разве можно при нахождении делта Х например, от градусов отнимать градусы?

[Игорь Белов]

Разве можно при нахождении делта Х например, от градусов отнимать градусы?

В первом случае - с длинами - должны быть радианы. Во втором случае - только с азимутами - можно использовать и градусы, но Cos() требует аргумента в радианах, а Atn() в радианах же результат выдаёт.

[klerk184]

Разве можно при нахождении делта Х например, от градусов отнимать градусы?

В первом случае - с длинами - должны быть радианы. Во втором случае - только с азимутами - можно использовать и градусы, но Cos() требует аргумента в радианах, а Atn() в радианах же результат выдаёт.

т.е. вот так:

Код: Выделить всё

Dim dX, dY, rumb, alfa, Radius, secondX As Float
						
						Radius = 6370997
						For next_one = 1 To end_for
							If next_one < end_for Then
								start_x = massX(next_one)
								start_y = massY(next_one)
								next_two = next_one + 1
								If next_two <= end_for Then							
									end_x = massX(next_two)
									end_y = massY(next_two)
								End If
								secondX = (start_x + end_x) / 2							
								dX = Radius * (end_x - start_x) * (cos(secondX) * RAD_2_DEG)
								dY = Radius * (end_y - start_y)
								'Note dX
									If end_x > start_x And end_y > start_y Then
									' 1 четверть
										direction = "СВ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = rumb
									End If
									If end_x > start_x And end_y < start_y Then
									' 2 четверть
										direction = "ЮВ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 180 - rumb
									End If
									If end_x < start_x And end_y < start_y Then
									' 3 четверть
										direction = "ЮЗ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 180 + rumb
									End If
									If end_x < start_x And end_y > start_y Then
									' 4 четверть
										direction = "СЗ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 360 - rumb
									End If

[Игорь Белов]

1. Плохо представляете, что происходит с индексами. В полигонах последняя точка совпадает с первой. Пишем "For i = 1 To n - 1", и никаких "If next_one", "If next_two".
2. Определитесь, что из φ, λ принимать за x, y. В Вашем коде сразу так и эдак в одной формуле.
3. Если нужны формулы с румбами, используйте абсолютные значения ∆x, ∆y: Abs(). Вообще-то румбы из эпохи таблиц, когда не только компьютеров, но и калькуляторов с функциями не было.

[klerk184]

1. Плохо представляете, что происходит с индексами. В полигонах последняя точка совпадает с первой. Пишем "For i = 1 To n - 1", и никаких "If next_one", "If next_two".
2. Определитесь, что из φ, λ принимать за x, y. В Вашем коде сразу так и эдак в одной формуле.
3. Если нужны формулы с румбами, используйте абсолютные значения ∆x, ∆y: Abs(). Вообще-то румбы из эпохи таблиц, когда не только компьютеров, но и калькуляторов с функциями не было.

* тогда получается вот так, так и правда лучше, меньше циклов лишних:

Код: Выделить всё

For next_one = 1 To end_for - 1
							'If next_one < end_for Then
								start_x = massX(next_one)
								start_y = massY(next_one)
								next_two = next_one + 1
								'If next_two <= end_for Then							
									end_x = massX(next_two)
									end_y = massY(next_two)
								'End If

* φ = У , λ = Х. Заметил свою ошибку

* получилось вот что:

Код: Выделить всё

Dim dX, dY, rumb, alfa, Radius, secondY As Float
						
						Radius = 6370997
						For next_one = 1 To end_for - 1
							'If next_one < end_for Then
								start_x = massX(next_one)
								start_y = massY(next_one)
								next_two = next_one + 1
								'If next_two <= end_for Then							
									end_x = massX(next_two)
									end_y = massY(next_two)
								'End If
								secondY = ((start_y + end_y) / 2) * DEG_2_RAD						
								dX = Abs(Radius * (end_x - start_x) * (cos(secondY) * RAD_2_DEG))
								dY = Abs(Radius * (end_y - start_y))
									If end_x > start_x And end_y > start_y Then
									' 1 четверть
										direction = "СВ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = rumb
									End If
									If end_x > start_x And end_y < start_y Then
									' 2 четверть
										direction = "ЮВ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 180 - rumb
									End If
									If end_x < start_x And end_y < start_y Then
									' 3 четверть
										direction = "ЮЗ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 180 + rumb
									End If
									If end_x < start_x And end_y > start_y Then
									' 4 четверть
										direction = "СЗ"
										rumb = Atn(dX / dY) * RAD_2_DEG
										alfa = 360 - rumb
									End If
									If dX = 0 And dY > 0 Then
										direction = "С"
									End If
									If dX = 0 And dY < 0 Then
										direction = "Ю"
									End If
									If dX > 0 And dY = 0 Then
										direction = "В"
									End If
									If dX < 0 And dY = 0 Then
										direction = "З"
									End If

Вроде что-то ближе к правде

[Игорь Белов]

Код: Выделить всё

dX = Abs(Radius * (end_x - start_x) * (cos(secondY) * RAD_2_DEG))

Зачем же косинус в градусы переводить? Это не угол, а величина безразмерная.

[gornak]

я с mapbasic-ом не знаком,
но догадываюсь какую задачу вы решаете, по-видимому вам необходимо построить геоход по
полигональным объектам, скорее всего это делянки(т.е. небольшие объекты)
поэтому советую:
1)не использовать географическую ск
тогда для нахождения угла вам необходимо
2)вычислить dx и dy
3)tan (a)=|dx|/|dy|
4)a=atan(a) угол
5)не забыть радианы перевести в градусы(десятичные)
6)если надо градусы-в градусы-(минуты-секунды)?
7)румб определяете по знакам dx и dy
l=(dx^2+dy^2)^0.5 - длина

[klerk184]

Код: Выделить всё

dX = Abs(Radius * (end_x - start_x) * (cos(secondY) * RAD_2_DEG))

Зачем же косинус в градусы переводить? Это не угол, а величина безразмерная.

Убрал. Что-то и не подумал)))

[klerk184]

я с mapbasic-ом не знаком,
но догадываюсь какую задачу вы решаете, по-видимому вам необходимо построить геоход по
полигональным объектам, скорее всего это делянки(т.е. небольшие объекты)
поэтому советую:
1)не использовать географическую ск
тогда для нахождения угла вам необходимо
2)вычислить dx и dy
3)tan (a)=|dx|/|dy|
4)a=atan(a) угол
5)не забыть радианы перевести в градусы(десятичные)
6)если надо градусы-в градусы-(минуты-секунды)?
7)румб определяете по знакам dx и dy
l=(dx^2+dy^2)^0.5 - длина

Да, сейчас буду переводить градусы в минуты и секунды