данные в двух СК

[Александр Мурый]

Уваж. знатоки! Сразу скажу - форум читал, мануалы cs2cs, proj, ogr2ogr смотрел.

Есть отдельные наборы векторных / растровых карт в двух СК. Их параметры в формате proj:

1. +proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=500000 +y_0=0 +no_defs +a=6378245 +rf=298.3 +to_meter=1
В QGIS проекция определилась как Pulkovo 1995 / Gauss-Kruger CM 39E

2. +proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=7500000 +y_0=0 +no_defs +a=6378245 +rf=298.3 +towgs84=28.000,-130.000,-95.000 +to_meter=1
проекция Pulkovo 42 / GK zone 7

Как видно, и то, и другое -- Gauss-Kruger, только СК95 / СК42.

Мне было сказано, что №1 - это местная СК.
Вопрос: Похоже на местную СК? Я не имел с ними дело..
Задача: совместить данные

Пробовал переводить координаты в ogr2ogr, затем "совмещать" данные в QGIS поиграться -- запутался

[nadiopt]

интересное совпадение, как раз сижу бьюсь над переводом из Пулково-95 в Пулково-42
в аркгис пишет "Трансформация не определена"
в Кугис молча преобразует, но явно неправильно, без учета эллипсоидов. Там у обеих Красовский или у 95-ой другой какой-то?
Гуру, помогите:)

[Александр Мурый]

гуру пока нет. напишу я)
как я прочитал где-то, у СК95 тоже эллипсоид Красовского, но "подправленный"

[nadiopt]

пока что мы сделали следующим образом
была СК-95 с трехградусными зонами Гаусса-Крюгера
перевели ее в СК-95 Гаусс-Крюгер зона 7Н
изничтожили СК, объявили СК Пулково-42 зона 7Н
перевели в ВГС
на гугль садится прилично:)

[Александр Мурый]

Да, надо и мне пойти через WGS84. Еще поучу матчасть..

[Mavka]

В обоих вариантах эллипсоид вообще не указан, должно быть типа "+ellps=krass".
Pulkovo 1995 / Gauss-Kruger CM 39E

[Александр Мурый]

Спасибо, эллипсоид-то пропустил.. Задача немного изменилась:

Pulkovo 1995 / Gauss-Kruger CM 39E ---> Pulkovo 42 / GK zone 7 <--- WGS84

В общем, делаю:
1. Pulkovo 1995 / Gauss-Kruger CM 39E --> Pulkovo 42 / GK zone 7

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

ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=500000 +y_0=0 +ellps=krass +units=m +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=7500000 +y_0=0 +ellps=krass +units=m +no_defs" -f "ESRI Shapefile" example_28407.shp example.shp

2. WGS84 --> Pulkovo 42 / GK zone 7

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

ogr2ogr -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" -t_srs "+proj=tmerc +lat_0=0 +lon_0=39 +k=1 +x_0=7500000 +y_0=0 +ellps=krass +units=m +no_defs" -f "ESRI Shapefile" graniza_28407.shp  graniza.shp

3. Ставлю в QGIS в свойствах проекта Pulkovo 42 / GK zone 7 и "преобразование координат на лету"

В результате по X порядок координат ~ 7 000 000 и 7 500 000 (различаются на 500 000)
по Y ~ 50 000 и 4 850 000 (различаются на 4 800 000)

Где я ступил, подскажите, а? Какие-то коэффициенты пересчета надо? Или неправильно указаны +x_0 и +y_0?

[SS_Rebelious]

В результате по X порядок координат ~ 7 000 000 и 7 500 000 (различаются на 500 000)
по Y ~ 50 000 и 4 850 000 (различаются на 4 800 000)

А нельзя ли пример начальных и конечных координат?

[Александр Мурый]

Пример:
1. Pulkovo 1995 / Gauss-Kruger CM 39E
X: 11 936,73
Y: 103 153,8
Pulkovo 42 / GK zone 7, та же точка объекта
X: 7 011 936,73
Y: 103 153,8

При включенном "преобразовании координат на лету" они совмещаются.

2. WGS84 (градусы)
X: 40,661097
Y: 43,561907

Pulkovo 42 / GK zone 7, та же точка

X: 7 634 209,886
Y: 4 826 572,81

То есть в Pulkovo 42 / GK zone 7 разница порядка 4 800 000 по Y

[SS_Rebelious]

По направлению север-юг запад-восток для X и Y для каждой из СК уточнить тоже не мешает.

Pulkovo 42 / GK zone 7, та же точка объекта
X: 7 011 936,73
Y: 103 153,8

здесь у координаты Y (на север или на восток?) не хватает миллионного значения (объект же не на экваторе?).

Прежде, чем продолжим дальше, следует выяснить: точно ли известны все необходимые параметры систем координат в которых находятся точки?

если широта 4103153.8, а долгота 7 011 936,73 в СК-42 зона 7 (шестиградусная), то в геокалькулятор фотомода даёт в WGS значения 36.932374136 и 33.522832811 соответственно, что далеко не 40 или 43.

Вообще, я сейчас в ирке, можем там обсудить.

[trir]

В СК-95 datum - ПЗ-90, а проекция кажется та же, что и в СК-42.
http://ru.wikipedia.org/wiki/%D0%9F%D0%97-90
ПЗ-90 - почти то же самое, что и WGS-84. Расхождение осей - метры, ориентация та же.
Конкретней смотри ссылки в википедии.
А у тебя не местная, а кажется просто обрезан номер зоны, типа для удобства. У нас так же делают,
только с СК-63 и называют это "условной системой координат" - чтоб враги не догадались ;=)

[Александр Мурый]

И как вернуть в "нормальный" вид? То есть "ключики" перехода от МСК к ГСК не нужны, так?

просто СК95 --> СК42 не получается

[trir]

Попробуй вариант - datum wgs84, а проекция GK zone 7, но без первой цифры

[trir]

не, гоню. смотри здесь:
viewtopic.php?f=34&t=5061

[Александр Мурый]

понял.
посмотрю параметры на http://gis-lab.info/qa/datum-transform-sets.html