Lviv eXtreme club

 ДопомогаПравила   ДопомогаДопомога   ПошукПошук   Список учасниківСписок учасників   Мобільна версіяМобільна версія   ГрупиГрупи   ЗареєструватисьЗареєструватись 
 ПрофільПрофіль   Увійти, щоб переглянути приватні повідомленняУвійти, щоб переглянути приватні повідомлення   ВхідВхід 

Перетворення географічних координат

Нова тема   Відповісти    Форум Lviv eXtreme club -> Тощо... Сторінка 1 з 1
Ваш часовий пояс: GMT + 2 Години
Попередня тема :: Наступна тема  
Повідомленняnomeaning © 14.10.09 13:11:34    Перетворення географічних координат

Хто можне навчити переводити географічні координати з усіма нюансами?)) SOS

Наприклад маю координату в прямокутній системі
ш 49,28 059(33)
д 23,52349(33)
38  еліпсоїд

Таку ж (похибка в пару метрів) В системі unknown (Google Earth):
49° 16' 50.28" ПН
23° 31' 24.49" СХ

В географічній системі
ш ?
д ?

Треба всі нюанси як переводити з одної в іншу і назад))
ПовідомленняYdrozhdzhal © 14.10.09 19:03:56    

А ти в гугла або у вікіпедії не питав як то зробити?  Wink
_________________
Jedem das Seine («Кожному своє»).
Повідомленняnomeaning © 14.10.09 19:57:13    

Чом би було не запитати?)) От тобі лінки деякі:

http://astronet.ru:8100/db/msg/1169819/node3.html
http://uk.wikipedia.org/wiki/Географічні_координати/
http://tvsh2004.narod.ru/geo_koor.htm

Переведеш?) Думаєш я створював тему для того, щоб мене питали чи я гугла питав?)) Можу навіть книжку ще вислати)) Просто для когось це може заняти 5хв викладення нюансів ніж мені перелопачування десятків сторінок літератури...
Повідомлення7eMEN-off © 15.10.09 9:03:46    Перетворення географічних координат

WGS84 → UTM і UTM → WGS84

І. Зауваження №1
ГОСТ 51794
X = N * CosB * CosL
Y = N * CosB * SinL , (в 2D проекції)

ІІ. Зауваження №2

(ГОСТ 51794-2001).
Вариант в Візуал Бейсик
Функции рабочей книги для преобразования геодезических координат из координатной системы Пулково 1942 в координатную систему WGS84 и обратно
Все угловые значения передаются и возвращаются в десятичных градусах (dd.ddddd), высоты передаются и возвращаются в метрах
Const Pi As Double = 3.14159265358979 ' Число Пи
Const ro As Double = 206264.8062 ' Число угловых секунд в радиане

' Эллипсоид Красовского
Const aP As Double = 6378245 ' Большая полуось
Const alP As Double = 1 / 298.3 ' Сжатие
Const e2P As Double = 2 * alP - alP ^ 2 ' Квадрат эксцентриситета

' Эллипсоид WGS84 (GRS80, эти два эллипсоида сходны по большинству параметров)
Const aW As Double = 6378137 ' Большая полуось
Const alW As Double = 1 / 298.257223563 ' Сжатие
Const e2W As Double = 2 * alW - alW ^ 2 ' Квадрат эксцентриситета

' Вспомогательные значения для преобразования эллипсоидов
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

' Линейные элементы трансформирования, в метрах
Const dx As Double = 23.92
Const dy As Double = -141.27
Const dz As Double = -80.9
' Угловые элементы трансформирования, в секундах
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
WGS84_SK42_Lat = Bd - dB(Bd, Ld, H) / 3600
End Function

Function SK42_WGS84_Lat(Bd, Ld, H) As Double
SK42_WGS84_Lat = Bd + dB(Bd, Ld, H) / 3600
End Function

Function WGS84_SK42_Long(Bd, Ld, H) As Double
WGS84_SK42_Long = Ld - dL(Bd, Ld, H) / 3600
End Function

Function SK42_WGS84_Long(Bd, Ld, H) As Double
SK42_WGS84_Long = Ld + dL(Bd, Ld, H) / 3600
End Function

Function dB(Bd, Ld, H) As Double
Dim B, L, M, N As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
M = a * (1 - e2) / (1 - e2 * Sin(B) ^ 2) ^ 1.5
N = a * (1 - e2 * Sin(B) ^ 2) ^ -0.5
dB = ro / (M + H) * (N / a * e2 * Sin(B) * Cos(B) * da _ + (N ^ 2 / a ^ 2 + 1) * 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, Ld, H) As Double
Dim B, L, N As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
N = a * (1 - e2 * Sin(B) ^ 2) ^ -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
Dim B, L, N, dH As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
N = a * (1 - e2 * Sin(B) ^ 2) ^ -0.5
dH = -a / N * da + N * Sin(B) ^ 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 ^ 2 / N + H) * ms
WGS84Alt = H + dH
End Function

ІІІ. Зауваження №3
а) Офійійний сайт розробника ПО для растрових карт на ГПС
б) варіанти перетворення з СК-42 в ВЖС-84  і назад в прогр. коді (можна відстежити формулу):

ІV. Зауваження №4
Особливість системи ВЖС 84 те, що в ній використовується еліпсоїдна площина, яка відрізняється від абсолютної в географічних координатах. Географічні координати уявляють землю в площині, абсолютна площина спотворює реальну поверхню від 0.2 до 2.3 %, тобто в середньому 30-40 см, а по висоті то до кількох метрів.

Для початку потрібно вирішити  яку координатну сітку потрібно переводити:
ГПС сітки: ITRS, WGS 84, Г13-90, ETRS, Nad-83, Navd-4
Географічні сітки: СК-42,Ск-95, UTM.

V. Зауваження №5
Системи переводу "Гугл Мапс" (книжка з формулами на англ. мові) WGS84-UMT

_________________
И ты будешь волков на земле плодить,
И учить их вилять хвостом!
А то, что придется потом платить,
Так ведь это ж, пойми, - потом!
Повідомленняnomeaning © 15.10.09 18:20:54    

7eMEN-off, дуже дякую!!))
навіть не думав що їх так багато, тепер буду мати чим бавитись)
ПовідомленняНаталі © 15.10.09 19:45:09    

Слухайте, хлопці, а нашо Вам такі дані? Ви що зібралися робити, якусь точкову атаку на нашу бідну Землю?? (7eMEN-off, по моєму вже й детальний план розробив..)
Surprised

_________________
Чтобы сохранить ангельский характер, нужно дъявольское терпение!
Повідомлення7eMEN-off © 16.10.09 8:16:56    

Наталі написав(ла):
7eMEN-off, по моєму вже й детальний план розробив..)
шоб атаку то ні, ну хібашо так трошка.

nomeaning написав(ла):
тепер буду мати чим бавитись)
Все це, я знайшов за 20 хв в Гуглі
_________________
И ты будешь волков на земле плодить,
И учить их вилять хвостом!
А то, что придется потом платить,
Так ведь это ж, пойми, - потом!
ПовідомленняYdrozhdzhal © 17.10.09 8:47:42    

От бачите - на рахунок Гугла я ж писав, хіба ні?  Acute
_________________
Jedem das Seine («Кожному своє»).
Повідомленняnomeaning © 17.10.09 15:34:20    

Наталі написав(ла):
Ви що зібралися робити, якусь точкову атаку на нашу бідну Землю??

просто забавка для вільного часу)) а якщо ще якись результат з тої забавки вийде то бде добре))
ydrozhdzhal написав(ла):
Гугла я ж писав, хіба ні?

писав, писав...))
Повідомлення7eMEN-off © 20.05.10 17:07:22    

БУМАЖНЫЕ КАРТЫ И GPS

На сегодняшний день существует и практически доступен большой объем бывшей секретной картографической продукции, которую предлагают различные картографические предприятия. Прочитав вышеизложенные теоретические выкладки, можно подумать: бери приемник GPS, покупай топографическую карту и смело отправляйся в поход! Но, к сожалению, ориентироваться на местности с большинством таких карт, используя приемник GPS, практически невозможно. Дело в том, что с секретных карт до сих пор вместе с секретными объектами удаляется вся координатная оцифровка. Однако, используя вышеполученные знания эту оцифровку с достаточной для самостоятельных путешествий точностью, можно восстановить. В этой статье будет предложена методика адаптации украинской картографической продукции, которая может использоваться и в других регионах.

Территория Украины, согласно атласа, изданного главным управлением геодезии и картографии при совете министров СССР в 1988 г., расположена в зонах - по долготе: 4 (18о - 24о), 5 (24о - 30о), 6 (30о - 36о), 7 (36о - 42о) и небольшой части 8-й (42о - 48о); по широте: от 42о до 54о.

В настоящее время наиболее полная картографическая информация представлена Киевской военно-картографической фабрикой - практически вся территория Украины покрыта картами масштаба 1:100 000. Однако эти карты не имеют координатной оцифровки. Используя вышеприведенные несекретные сведения градусную координатную сетку легко можно восстановить.

Карты масштаба 1:100 000 с градусной оцифровкой выпустили УкрАГП и НПП "Картография" в 1998 г. Однако этими картами покрыта только небольшая часть Украины. На сегодняшний день выпуск этой продукции прекращен.

Наибольший интерес для путешественников представляют карты масштаба 1:50 000 и более,которые выпускаются теми же УкрАГП и НПП "Картография". Этими картами покрыта часть Крыма и часть Карпатского региона. Градусные координаты на этих картах отсутствуют, разграфка может оказаться нестандартной. Поскольку эта продукция представляет наибольший интерес для владельцев приемников GPS приведем несколько соображений по восстановлению рамки этих карт.

В идеальном случае необходимо иметь координаты (показания приемника GPS) двух точек, определенных на карте, и знать масштаб этой карты. Чем дальше эти точки будут удалены друг от друга, тем точнее будет результат вычислений. Решение этой задачи предлагается производить в электронных таблицах, например, Excel. Идея метода показана на рисунке, алгоритм приводится ниже.

➊ Находим длину 1 сек дуги меридиана на широте точки A в метрах
ДCДМA = (111,143-0,562COS(2ШA))/3,6
где ШA - географическая широта точки A.
➋ Находим разность координат точек A и B по широте в секундах - РКСШ.
➌ Находим разность координат точек A и B по широте в метрах
РКМШ = РКСШ х ДCДМA
➍ Находим разность координат точек A и B по широте в сантиметрах на карте
РККШ = РКМШ : М
где М - масштаб или количество метров в 1 см карты.
➎ Находим длину 1 сек дуги параллели на широте точки A в метрах
ДCДПA = ((111,321COS(ШA) - 0,094COS(3ШA))/3,6
где ШA - географическая широта точки A.
➏ Находим разность координат точек A и B по долготе в секундах - РКСД.
➐ Находим разность координат точек A и B по долготе в метрах
РКМД = РКСД х ДCДПA
➑ Находим разность координат точек A и B по долготе в сантиметрах на карте
РККД = РКМД : М
➒ При помощи циркуля и линейки строим прямоугольный треугольник ABC у которого катет AC должен быть параллелен верхней и нижней стороне рамки карты. В случае непараллельности можно сделать вывод о том, что данная карта не соответствует стандартной разграфке. В этом случае необходимо нанести новую рамку. Нижнюю и верхнюю стороны рамки можно начертить произвольно с соблюдением параллельности AC. Для территории Украины на картах масштаба 1:50 000 и более сближение меридианов можно не учитывать. Поэтому можно считать, что левая и правая стороны рамки карты параллельны катету BC. В случае непараллельности эти стороны также чертим самостоятельно.
➓ Используя свойства пропорций находим координаты углов рамки (существующей или вновь построенной) и наносим на рамку координатные деления.

_________________
И ты будешь волков на земле плодить,
И учить их вилять хвостом!
А то, что придется потом платить,
Так ведь это ж, пойми, - потом!
Повідомленняksenia © 20.05.10 20:36:53    

Якщо працюєте в ОзіЕкспорер, то там є такий момент, коли переводиться із системи в систему і є похибка, але це вирішується і про це окрема мова. Ще декілька закладок (по суті, шукаємо конвертор координат):

- Он-лайн конвертор координат GPS
- Програмуля (конвертер координат XYZ ←→ blh)
- Coordinate converter v1.00
- Тут багато цікавого
Повідомлення7eMEN-off © 20.05.10 20:52:21    

тодо жопи... нагуглити я того сам можу...
А вот як взяти звичайну карту і на ній насічки зробити в системі WGS84?
Отут вам ні один гугл не поможе... тут треба переміряти все циркулем, лінійкою, врахувати похибку проекції Гаусса-Крюгера та похибки радянської системи Пулково-1942...
Якщо все правильно зробите то все рівно похибка від 2 до 5 секунд, а це від 50 до 160 м на місцевості...

_________________
И ты будешь волков на земле плодить,
И учить их вилять хвостом!
А то, что придется потом платить,
Так ведь это ж, пойми, - потом!
Показувати:   
Нова тема   Відповісти    Форум Lviv eXtreme club -> Тощо... Ваш часовий пояс: GMT + 2 Години
Сторінка 1 з 1

 
Перейти до:  






© 2001-2017, eXtreme.lviv.ua     При використанні матеріалів посилання на джерело обов'язкове
СТАНЬ
ОДНИМ
З НАС




Экстремальный портал VVV.RU