UR5FCM Logger ver. 3b

[UR5FCM]

Есть 2 точки: RX4HX и W7AA. Начальное направление 0 градусов. Рисуем линию длинной 100 км от RX4HX в сторону W7AA. Теперь считаем направление от этой точки до W7AA - оно будет уде не 0 градусов, а 359,9. Рисуем опять линию 100 км. И т.д. Принцип ясен?

Принцип то ясен, но как вычислить зная 100 км координаты конца отрезка?
Начальные угол азимута всегда 0?

[RX4HX]

Принцип то ясен, но как вычислить зная 100 км координаты конца отрезка?
Начальные угол азимута всегда 0?

Нет конечно, не всегда 0. 0 в данном случае - это начальный азимут от RX4HX до W7AA. Делается такой цикл:

Точка1 = RX4HX; // Координаты X,Y на карте
Точка2 = W7AA;

пока Точка1 <> Точка2 делаем
определяем азимут от Точка1 до Точка2.
Точка3 = от Точка1 определенным выше азимутом длинной на корте равной 100 км.
Нарисовать линию(Точка1, Точка3);
Точка1 = Точка3;
конец цикла пока.

Все! Очень просто.

[UR5FCM]

Точка1 = RX4HX; // Координаты X,Y на карте
Точка2 = W7AA;

пока Точка1 <> Точка2 делаем
определяем азимут от Точка1 до Точка2.
Точка3 = от Точка1 определенным выше азимутом длинной на корте равной 100 км.
Нарисовать линию(Точка1, Точка3);
Точка1 = Точка3;
конец цикла пока.

Азимут от Точки1 до Точки 2 я знаю
Где взять координаты Точки 3?
Точка2 это W7AA?
Да так там еще же x,y

[RX4HX]

Азимут от Точки1 до Точки 2 я знаю
Где взять координаты Точки 3?
Точка2 это W7AA?
Да так там еще же x,y

Для того, чтоб решить задачу с картотой Вам прежде всего нужно написать следущие ф-ции:
1. Перевод географических координат в X,Y координаты Вашей карты.
2. Зная координаты начальной точки и азимут и дистанцию до искомой точки = получить координаты искомой точки. Сами эту ф-ци. Вы не напишите - она очень сложная. Вот ее реализация:

procedure AzimuthDistToLatLon(const LatIn, LonIn, Azimuth, Distance: extended;
var LatOut, LonOut: extended);
const
EARTHRADIUS = 6366.71;
var
dlon, sinlat1, coslat1, sind, cosd, sintc, costc: extended;

begin
SinCos(LatIn, sinlat1, coslat1);
SinCos(Distance / EARTHRADIUS, sind, cosd);
SinCos(Azimuth, sintc, costc);
LatOut := arcsin(sinlat1 * cosd + coslat1 * sind * costc);
dlon := arctan2(sintc * sind * coslat1, cosd - sinlat1 * sin(LatOut));
LonOut := SphericalMod(LonIn - dlon + LocalPI) - LocalPI;
LatOut := LatOut * R2D;
LonOut := LonOut * R2D;
if LonOut < -180 then
LonOut := 360 + LonOut;
end;

Ну думаю теперь все понятно.

[UR5FCM]

Спасибо!
Что такое SinCos, R2D, SphericalMod
atan2 -это арктангенс?

Не понятно где брать sinlat и т.д и что такое extended?

Куда теперь эту функцию
Еще dll

Координаты привязаны к X,Y c самого начала

Ух не люблю я эту математику :russian_roulette:

Давай по e-mail

Добавлено через 1 час 34 минуты

Я надеюсь, что кроме своего лога на вашем компе есть еще другие логи ? Александра или Алексея, например, или Logger, и повводите для примера позывные штатов, японии, VK или ZL, посмотрите, какие направления на корреспондентов покажут вам программы...

Да азимут считает у меня правильно
Тут вопрос в отображении на карте. Понимаете

[UA3GDW]

Что такое SinCos, R2D, SphericalMod
atan2 -это арктангенс?

SinCos, ArcTan2 - это из модуля Math, а вот, что такое R2D и SphericalMod - действительно не понятно.

[RX1AL]

Игорь, я же тебе линки послал... на мыло... Разве не хватает? Тем более там готовый и практический пример есть. С полным кодом... а не с одной функцией, непонятно откуда взятой... и к тому же не полной, поскольку нет некоторых функций. Если уж приводится фрагмент кода, то неплохо иметь и первоисточник или ссылку на него... нормальная практика...

[UR5FCM]

Игорь, я же тебе линки послал... на мыло... Разве не хватает? Тем более там готовый и практический пример есть. С полным кодом... а не с одной функцией, непонятно откуда взятой... и к тому же не полной, поскольку нет некоторых функций. Если уж приводится фрагмент кода, то неплохо иметь и первоисточник или ссылку на него... нормальная практика...

Да это было до того как получил линки
В принципе уже получается дуга надо еще потестироватьна направлениях
Пришлю скриншот на мыло попозже
Ну да некоторые функции отстутствуют
Ладно пришлось ArcSin ATan2 поискать и с радианами повозиться

Спасибо

[RX4HX]

procedure SinCos(Theta: extended; var eSin, eCos: extended);
asm
FLD Theta
FSINCOS
FSTP tbyte ptr [EDX] // Cos
FSTP tbyte ptr [EAX] // Sin
FWAIT
end;

R2D: double = 57.295781; // Константа для преобразования радиан в градусы

//********************************************************************** *
// X значение в радиантах по модулю
// Возвращает X в диапазоне -LocalPi to +LocalPi
//********************************************************************** *

function SphericalMod(X: extended): extended;
begin
Result := X;
if X < -LocalPi then
Result := X + DoublePi
else if X > LocalPi then
Result := X - DoublePi;
end;

[UR5FCM]

26.01.2009

Build ver. 200

+ Добавлено возможность скрывать / показывать колонки в логе
+ Изменять размеры колонок, шрифт, цвет шрифта
* Исправлено отображение линии между корреспондентами на карте QSO
+ Опция в поиске. Можно теперь выбрать связи по новым критериям (QSL и eQSL получена/отправлена)
* Исправлены некоторые ошибки в коде
+ новая библиотека для CW

Установка распакуйте архив в каталог с программой
Файлы в архиве
- Cweng.dll - новая библиотека CW
- config.mdb - база данных для сохранения настроек
-UR5FCMLogger.exe - новый исполняемый файл программы

Ссылка для скачивания
ur5fcm-new_200.exe

[UR5FCM]

И снова новая тестовая версия 201
Установка распакуйте архив в каталог с программой
Файлы в архиве
- Cweng.dll - новая библиотека CW
- config.mdb - база данных для сохранения настроек (должен быть в
папке с программой проверьте)
-UR5FCMLogger.exe - новый исполняемый файл программы

Вид программы
1. Настройки колонок лога
2. Вид карты (исправленная)

[RX1AL]

UR5FCM:

Вот... уже не слово мальчика, а мужа! Так и продолжать...

Из замечаний... Я одного не понял: ты считаешь дистанцию от себя до корреспондента всегда по длинному пути? Если так, то это не совсем верно. Нужно в обе стороны считать. У нас же есть Short Path и Long Path. И отображать их разным цветом в таком случае. Внизу можно дать легенду по цветам, или добавить тултип для линии. Ты же наверху данные отображаешь для направления антенны, ну вот...

Чтобы не думать, о том какую формулу использовать для нахождения кратчайшей дистанции с учетом кривизны поверхности земли, сразу скажу, что это формула Haversine.

ПС То есть завтра у нас будет уже 202 билд...

[UR5FCM]

Из замечаний... Я одного не понял: ты считаешь дистанцию от себя до корреспондента всегда по длинному пути? Если так, то это не совсем верно. Нужно в обе стороны считать. У нас же есть Short Path и Long Path. И отображать их разным цветом в таком случае. Внизу можно дать легенду по цветам, или добавить тултип для линии. Ты же наверху данные отображаешь для направления антенны, ну вот...

Да вроде этого бы хватило
Я проверял в ДХАтласе такую же линию рисует
Haversin есть в VB6 или нет
Я думаю нет пока нужды рисовать обратный азимут

[RX1AL]

UR5FCM:

Обратный азимут тоже надо рисовать... а как же... У тебя же направление зависит от нескольких факторов: прохождения и направления антенны.
Haversinе прамой реализации нет, но написать ее дело 10-15 минут. Тем более, у меня и код готовый есть... Когда-то писал... На мыло могу заслать, если скажешь, что надо...

Из добавлений... Если уж ты имеешь карту, то неплохо было бы на ней день и ночь отображать, в зависимости от солнца... И не в статике, а динамике по изменению времени. Очень полезная функциональность. А если еще связать с MUF - так вообще будет то, что доктор прописал.

Ниже код функции на VB 6.0:

Public Function Haversine(ByVal Latitude1 As Double, _
ByVal Latitude2 As Double, _
ByVal Longitude1 As Double, _
ByVal Longitude2 As Double) As Double

Const EarthRadius As Integer = 6371 'Earth radius in kilometer

Dim DeltaLatitude As Double, DeltaLong As Double
Dim a As Double, c As Double
Dim Pi As Double

On Error GoTo ErrorExit

Pi = 4 * Atn(1)

'Convert Latitude1, Latitude2, Longitude1, Longitude2 from decimal degrees into radians
Latitude1 = Latitude1 * Pi / 180
Latitude2 = Latitude2 * Pi / 180
Longitude1 = Longitude1 * Pi / 180
Longitude2 = Longitude2 * Pi / 180

'Calculate change in Latitude and Longitude
DeltaLatitude = Abs(Latitude2 - Latitude1)
DeltaLongitude = Abs(Longitude2 - Longitude1)

a = ((Sin(DeltaLatitude / 2)) ^ 2) + (Cos(Latitude1) * Cos(Latitude2) * ((Sin(DeltaLongitude / 2)) ^ 2))
c = 2 * Atan((Sqr(1 - a)) / (Sqr(a)))
Haversine = R * c
ErrorExit:
End Function

[UR5FCM]

Ниже код функции на VB 6.0:

Спасибо
и куда это
я не понял
Да хватит отображения одного азимута
Я же рисую то считает программа по координатам без учета прохода и т.д.