Это умная програмка перевода с координат системы GPS в наши метровые координаты. Очень полезная вещь. Только написана она не на делфи.




1 Функции рабочей книги для преобразования
1 геодезических координат из координатной системы Пулково 1942
' в координатную систему WGS84
' Все угловые значения передаются и возвращаются в градусах, ' высоты передаются и возвращаются в метрах
1 Сделано на основании ГОСТ Р 51794-2001 Российской Федерации
Const Pi As Double = 3.14159265358979 ' Число, Пи
Const ro As Double = 206264.8062 ' Число угловых секунд в радиане
1 Эллипсоид Красовского (Пулково 1942)
Const aP As Double = 6378245 ' Большая полуось
Const alP As Double = 1 / 298.3 ' Сжатие
Const e2P As Double = 2 * alP - alP A 2 ' Квадрат эксцентриситета
1 Эллипсоид GRS80 (WGS84)
Const aW As Double = 6378137 ' Большая полуось
Const alW As Double = 1 / 298.257223563 ' Сжатие
Const e2W As Double = 2 * alW - alW x 2 ' Квадрат эксцентриситета
1 Вспомогательные значения для преобразования эллипсоидов Const a As Double = (aP + aW) /2 Const e2 As Double = (e2P + e2W) / 2 Const da As Double = aW - aP Const de2 As Double = e2W - e2P
1 Линейные элементы трансформирования, в метрах
Const dx As Double =23.92
Const dy As Double = -141.27
Const dz As Double = -80.9
1 Угловые элементы трансформирования, в секундах
Const wx As Double = 0
Const wy As Double = 0
Const wz AS Double = 0
' Дифференциальное различие масштабов
Const ms As Double = 0
Function WGS84_SK42_Lat(Bd, Ld, H As Double) As Double
WGS84_SK42_Lat = Bd - dB(Bd, Ld, H) / 3600
End Function
Function SK42_WGS84_Lat(Bd, Ld, H As Double) As Double
SK42__WGS84__Lat = Bd + dB(Bd, Ld, H) / 3600
End Function
Function WGS84_SK42__Long(Bd, Ld, H As Double) As Double
WGS84_SK42__Long = Ld - dL(Bd, Ld, H) / 3600
End Function
Function SK42_WGS84_Long(Bd, Ld, H As Double) As Double
SK42_WGS84_Long = Ld + dL(Bd, Ld, H) / 3600
End Function
Function dB(Bd, Ld, H As Double) As Double
Dim B, L, M, N- As Double
В = Bd * Pi / 180
L * Ld * Pi / 180
M = a * (1 - e2) / (1 - e2 * Sin(B) A 2) A 1.5
N = a * (1 - e2 * Sin(B) A 2) л -0.5
dB = ro / (M + H) * (N / a * e2 * Sin(B) * Cos(B) * da _ + (NA2/aA2+l) *N* Sin(B) * Cos(B) * de2 / 2
- (dx * Cos(L) + dy * Sin(L)) * Sin(B) + dz * Cos(B)) _
- wx * Sin(L) * (1 + e2 * Cos ,(2 * B) ) _
+ wy * Cos(L) * (1 + e2 * Cos(2 * B)) _
- ro * ms * e2 * Sin(B) * Cos(B)




End Function
Function dL(Bd, Lid, H As Double) As Double
Dim B, L, N As Double
В = Bd * Pi / 180 L = Ld * Pi / 180 N = a * (1 - e2 * Sin(B) л 2) A -0.5
dL = ro / ((N + H) * Cos(B)) * (-dx * Sin(L) + dy * Cos(L))
+ Tan(B) * (1 - e2) * (wx * Cos(L) + wy * Sin(L)) - wz
End Function
Function WGS84Alt(Bd, Ld, H As Double) As Double
Dim B, L, N, dH As Double
В = Bd * Pi / 180 L = Ld * Pi / 180 N = a * (1 - e2 * Sin(B) A 2) x -0.5
dH = -a / N * da + N * Sin(B) A 2 * de2 / 2
+ (dx * Cos(L) + dy * Sin(L)) * Cos(B) + dz * Sin(B) _
- N * e2 * Sin(B) * Cos(B) * (wx / ro * Sin(L) - wy / ro * Cos(L)) _
+ (a A 2 / N + H) * ms
WGS84Alt = H + dH End Function

ПОПРАВОЧКА


Вместо л и А ставится ^ это в квадрате просто фаин ридер плохо раскидал символы

Модератор: друзья! Если у вас назрел вопрос по Бэйсику или по Делфи или по этим обоим языкам вместе взятым, то зачем же постить вопрос в раздел Паскаля? В общем, перенесено из Паскаля в Центр помощи.

Спасибо что пернесли туда куда надо. но мне надо решить этот вопрос. с переводом с системы измерения GPS в нашё метровую систему координат

Ну и в чем собсно проблема? Инструкции вида

Цитата(VirtualHouse @ 20.8.2005, 03:04)

Function WGS84_SK42_Lat(Bd, Ld, H As Double) As Double

Заменяешь на

Код

Function WGS84_SK42_Lat(Bd, Ld, H: Double) : Double; begin

Цитата(VirtualHouse @ 20.8.2005, 03:04)

WGS84_SK42_Lat = Bd - dB(Bd, Ld, H) / 3600

на

Код

WGS84_SK42_Lat := Bd - dB(Bd, Ld, H) / 3600;

Цитата(VirtualHouse @ 20.8.2005, 03:04)

Const aW As Double = 6378137 ' Большая полуось

на

Код

Const aW : Double = 6378137; //Большая полуось

(или как там в пасе консты объявляются?)

Цитата(Mayk @ 20.8.2005, 09:31)

или как там в пасе консты объявляются?)

Верно.. так...

Спасибо но он мне тут даёт ошибку В = Bd * Pi / 180

Цитата(VirtualHouse @ 20.8.2005, 11:37)

В = Bd * Pi / 180

Код

В := Bd * Pi / 180

Вроде накатал нечто похожее на этот код. Только перед работой нужно будет вызвать процедуру InitVariables, т.к. Дельфи не понимает объявлений вида "da As Double = aW - aP" ни для констант, ни для переменных.

Код

unit gps; interface uses Math;    procedure InitVariables;    function WGS84_SK42_Lat(Bd, Ld, H :Double) :Double;    function SK42_WGS84_Lat(Bd, Ld, H :Double) :Double;    function WGS84_SK42__Long(Bd, Ld, H :Double) :Double;    function SK42_WGS84_Long(Bd, Ld, H :Double) :Double;    function dB(Bd, Ld, H :Double) :Double;    function dL(Bd, Ld, H :Double) :Double;    function WGS84Alt(Bd, Ld, H :Double) :Double; const    Pi   :Double = 3.14159265358979;    ro   :Double = 206264.8062;    aP   :Double = 6378245;    alP  :Double = 1 / 298.3;    aW   :Double = 6378137;    alW  :Double = 1 / 298.257223563;    dx   :Double = 23.92;    dy   :Double = -141.27;    dz   :Double = -80.9;    wx   :Double = 0;    wy   :Double = 0;    wz   :Double = 0;    ms   :Double = 0; var    e2W, a, e2P, de2, e2, da :Double; implementation procedure InitVariables; begin    e2W:= 2 * alW - Sqr(alW);    a:= (aP + aW) / 2;    e2P:= 2 * alP - Sqr(alP);    de2:= e2W - e2P;    e2:= (e2P + e2W) / 2;    da:= aW - aP; end; function WGS84_SK42_Lat(Bd, Ld, H :Double) :Double; begin    Result:= Bd - dB(Bd, Ld, H) / 3600; end; function SK42_WGS84_Lat(Bd, Ld, H :Double) :Double; begin    Result:= Bd + dB(Bd, Ld, H) / 3600; end; function WGS84_SK42__Long(Bd, Ld, H :Double) :Double; begin    Result:= Ld - dL(Bd, Ld, H) / 3600; end; function SK42_WGS84_Long(Bd, Ld, H :Double) :Double; begin    Result:= Ld + dL(Bd, Ld, H) / 3600; end; function dB(Bd, Ld, H :Double) :Double; var    B, L, M, N :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    M:= a * (1 - e2) / Power(1 - e2 * Sqr(Sin(B)), 1.5);    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    Result:= ro / (M + H) * (N / a * e2 * Sin(B) * Cos(B) * da +       (Sqr(N) / Sqr(a)+l) * N * Sin(B) * Cos(B) * de2 / 2 -       (dx * Cos(L) + dy * Sin(L)) * Sin(B) + dz * Cos(B)) -       wx * Sin(L) * (1 + e2 * Cos(2 * B)) + wy * Cos(L) *       (1 + e2 * Cos(2 * B)) -  ro * ms * e2 * Sin(B) * Cos(B); end; function dL(Bd, Ld, H :Double) :Double; var    B, L, N :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    Result:= ro / ((N + H) * Cos(B)) * (-dx * Sin(L) + dy * Cos(L)) +       Tan(B) * (1 - e2) * (wx * Cos(L) + wy * Sin(L)) - wz; end; function WGS84Alt(Bd, Ld, H :Double) :Double; var    B, L, N, dH :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    dH:= -a / N * da + N * Sqr(Sin(B)) * de2 / 2 +       (dx * Cos(L) + dy * Sin(L)) * Cos(B) + dz * Sin(B) -       N * e2 * Sin(B) * Cos(B) * (wx / ro * Sin(L) - wy / ro * Cos(L)) +       (Sqr(a) / N + H) * ms;    Result:= H + dH; end; end.

Это сообщение отредактировал(а) rsm - 20.8.2005, 12:29

Цитата(rsm @ 20.8.2005, 13:15)

Только перед работой нужно будет вызвать процедуру InitVariables

Не нужно:

Код

unit gps; interface uses Math;    function WGS84_SK42_Lat(Bd, Ld, H :Double) :Double;    function SK42_WGS84_Lat(Bd, Ld, H :Double) :Double;    function WGS84_SK42__Long(Bd, Ld, H :Double) :Double;    function SK42_WGS84_Long(Bd, Ld, H :Double) :Double;    function dB(Bd, Ld, H :Double) :Double;    function dL(Bd, Ld, H :Double) :Double;    function WGS84Alt(Bd, Ld, H :Double) :Double; const    Pi   :Double = 3.14159265358979;    ro   :Double = 206264.8062;    aP   :Double = 6378245;    alP  :Double = 1 / 298.3;    aW   :Double = 6378137;    alW  :Double = 1 / 298.257223563;    dx   :Double = 23.92;    dy   :Double = -141.27;    dz   :Double = -80.9;    wx   :Double = 0;    wy   :Double = 0;    wz   :Double = 0;    ms   :Double = 0; var    e2W, a, e2P, de2, e2, da :Double; implementation procedure InitVariables; begin    e2W:= 2 * alW - Sqr(alW);    a:= (aP + aW) / 2;    e2P:= 2 * alP - Sqr(alP);    de2:= e2W - e2P;    e2:= (e2P + e2W) / 2;    da:= aW - aP; end; function WGS84_SK42_Lat(Bd, Ld, H :Double) :Double; begin    Result:= Bd - dB(Bd, Ld, H) / 3600; end; function SK42_WGS84_Lat(Bd, Ld, H :Double) :Double; begin    Result:= Bd + dB(Bd, Ld, H) / 3600; end; function WGS84_SK42__Long(Bd, Ld, H :Double) :Double; begin    Result:= Ld - dL(Bd, Ld, H) / 3600; end; function SK42_WGS84_Long(Bd, Ld, H :Double) :Double; begin    Result:= Ld + dL(Bd, Ld, H) / 3600; end; function dB(Bd, Ld, H :Double) :Double; var    B, L, M, N :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    M:= a * (1 - e2) / Power(1 - e2 * Sqr(Sin(B)), 1.5);    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    Result:= ro / (M + H) * (N / a * e2 * Sin(B) * Cos(B) * da +       (Sqr(N) / Sqr(a)+l) * N * Sin(B) * Cos(B) * de2 / 2 -       (dx * Cos(L) + dy * Sin(L)) * Sin(B) + dz * Cos(B)) -       wx * Sin(L) * (1 + e2 * Cos(2 * B)) + wy * Cos(L) *       (1 + e2 * Cos(2 * B)) -  ro * ms * e2 * Sin(B) * Cos(B); end; function dL(Bd, Ld, H :Double) :Double; var    B, L, N :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    Result:= ro / ((N + H) * Cos(B)) * (-dx * Sin(L) + dy * Cos(L)) +       Tan(B) * (1 - e2) * (wx * Cos(L) + wy * Sin(L)) - wz; end; function WGS84Alt(Bd, Ld, H :Double) :Double; var    B, L, N, dH :Double; begin    B:= Bd * Pi / 180;    L:= Ld * Pi / 180;    N:= a * Power(1 - e2 * Sqr(Sin(B)), -0.5);    dH:= -a / N * da + N * Sqr(Sin(B)) * de2 / 2 +       (dx * Cos(L) + dy * Sin(L)) * Cos(B) + dz * Sin(B) -       N * e2 * Sin(B) * Cos(B) * (wx / ro * Sin(L) - wy / ro * Cos(L)) +       (Sqr(a) / N + H) * ms;    Result:= H + dH; end; initialization   InitVariables; end.

Спрасибо большое. Кстати эта нужная прога для всех кто пользуется GPS