Страница 1 из 2

Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 04 дек 2013, 04:31
wowka1319
Имеется GeoTiff. Открываю его в QGIS и вижу декартовы координаты (внизу в строке состояния), причем нулевой или близкой к ней точки нет. Полагаю это спроецированные и привязанные координаты.
Необходимо получить/просмотреть географические координаты, но как я не крутил настройки QGIS (задавал эллипсоид где только можно) - ничего не меняется.

И еще: я понимаю что есть WGS84 и Pulkovo1942 (и еще куча других) и что это параметры эллипсоидов, но что значат другие слова/значения (какие-то зоны) в названиях, типа:
"WGS84 UTM N20",
"Pulkovo 1942 / Caspian Sea Mercator",
"Pulkovo 1942(58)",
"Pulkovo 1942(83)",
"WGS 84 / Mercator 41",
"WGS 84 / Pseudo Mercator",
???
и что значат идентификаторы EPSG:#### ???


//Добавлено:
Уточню исходные данные (то что удалось вытащить):
GeoTiff, проекция Гаусса-Крюгера, вродебы Pulkovo1942, какой-то Красовский 1940 и еще есть ГК углов проекции с точностью до минуты. В общем теряюсь что куда и почему.

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 04 дек 2013, 06:32
dab
wowka1319 писал(а):Необходимо получить/просмотреть географические координаты, но как я не крутил настройки QGIS (задавал эллипсоид где только можно) - ничего не меняется.
Исходя из моего предположения, что вы используете QGIS 2:
"Проект" / "Свойства проекта..." / "Система координат" / "Включить автоматическое перепроецирование координат" / WGS 84

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 04 дек 2013, 08:58
wowka1319
Спасибо! А как насчет остальных вопросов? что это за деление на зоны - какой в этом смысл? разве та же WGS84 не универсальна для всего земного шара?

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 04 дек 2013, 09:15
Максим Дубинин

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 05 дек 2013, 07:35
wowka1319
Я читал FAQ, но там ни слова про зоны, зачем они нужны и EPSG.
Так что вопрос открыт.

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 05 дек 2013, 07:57
Petruxin
Про зоны
miltop.narod.ru/Nomenclature/projections.htm
http://miltop.narod.ru/Flat/concept.htm

[ Сообщение с мобильного устройства ]

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 05 дек 2013, 07:57
Petruxin
http://miltop.narod.ru/Nomenclature/projections.htm

[ Сообщение с мобильного устройства ]

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 05 дек 2013, 08:00
trir
Плохо читал.
Зоны нужны потому, что в поперечно-цилиндрических проекциях набегает ошибка по длине, при удалении от осевого меридиана - зоны отличаются положением осевого меридиана (например на 6°)
EPSG - это просто европейские нефтяники ведут базу со всеми СК и это ключ в этой базе

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 05 дек 2013, 13:13
dab
wowka1319 писал(а):Я читал FAQ, но там ни слова про зоны, зачем они нужны и EPSG.
Так что вопрос открыт.
Вы бы пояснили, с какой целью интересуетесь - с чисто образовательной или практическая задача имеется?

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 01:07
wowka1319
Есть практическая задача - нписать програмку для нужд на работе:
На работе принимаем снимки со спутников и некоторые обрабатываем.
Суть обработки:
1 - вырезать интересные области
2 - наложить на них сетку широта/долгота
3 - оформить подложку (белый фон с надписями)
4 - подписать некоторые объекты (реки, озера, заливы, моря, города и т.д.)

1 пункт выполняется в MeteorSatN и проблем не встречается
2 и 3 пункты выполняются в ArcMap
- вот здесь я и вижу проблему, а вернее возможность автоматизировать процесс. Подробнее:
Снимок приходящий в MeteorSatN (1 пункт) несет в себе информацию о координатах, но все форматы в которые он сохраняет (в т.ч. GeoTiff) НЕ понимает ArcMap. Поэтому в ArcMap мы загружаем снимок в BMP, теряя координаты, и восстанавливаем координаты по подложенному под снимок атласу при помощи ключевых точек (снимок при этом сжимается/деформируется). А это занятие нудное и неблагодарное, да еще и ArcMap глючит и тормозит безбожно.
Ну а в фотошопе (п.4) мы подписываем объекты, причем каждый раз вручную, несмотря на то, что работаем мы каждый день со снимками одних и тех же областей (просто немного сдвинутых при каждом витке). Здесь я хочу сделать некоторую БД вида "надпись объекта - ГК", чтобы объекты подписывались автоматом (ну с возможностью корректировки).

Вот что бы написать программу автоматизации этих "задротских" процессов я и пытаюсь вникнуть в это.
Например я проверил ГК одного приходящего снимка и увидел погрешность около Байкала по широте/долготе - 7'/2' что очень много и следовательно неприемлемо (проверял по атласу OziExplorer). И я сомневаюсь что эта погрешность имеет линейную закономерность и поэтому необходимо знать всю "химию" преобразований. Может дело как раз в этих самых зонах (к нам приходят снимки в границах по долготе от Байкала до Японии и по широте от Якутии-Камчатка до границ с Монголией и Китаем).

На всякий случай приложу gdalinfo геотифа (планирую работать с этим форматом):

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

Driver: GTiff/GeoTIFF
Files: TEST_01.tiff
Size is 4483, 2494
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Pulkovo 1942",
        DATUM["Pulkovo_1942",
            SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
                AUTHORITY["EPSG","7024"]],
            AUTHORITY["EPSG","6284"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4284"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",111],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (73671.999948490688000,5838215.177363829700000)
Pixel Size = (50.000000000000057,-50.040096230954290)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (   73672.000, 5838215.177) (112d 5'20.30"E, 52d40'0.97"N)
Lower Left  (   73672.000, 5713415.177) (112d 3'43.17"E, 51d32'43.98"N)
Upper Right (  297822.000, 5838215.177) (115d23'44.13"E, 52d35'25.17"N)
Lower Right (  297822.000, 5713415.177) (115d17'13.46"E, 51d28'19.01"N)
Center      (  185747.000, 5775815.177) (113d42'34.71"E, 52d 4'48.09"N)
Band 1 Block=4483x1 Type=UInt16, ColorInterp=Red
Band 2 Block=4483x1 Type=UInt16, ColorInterp=Green
Band 3 Block=4483x1 Type=UInt16, ColorInterp=Blue

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 03:53
trir
Может для начала взять прогу которая умеет работать с Geotiff? QGIS например...

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 06:40
wowka1319
Зачем? (Qgis давно установлен)

[ Сообщение с мобильного устройства ]

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 07:27
Игорь Лебедь
wowka1319 писал(а):Снимок приходящий в MeteorSatN (1 пункт) несет в себе информацию о координатах, но все форматы в которые он сохраняет (в т.ч. GeoTiff) НЕ понимает ArcMap.

На всякий случай приложу gdalinfo геотифа (планирую работать с этим форматом):
Driver: GTiff/GeoTIFF
Files: TEST_01.tiff
***
AUTHORITY["EPSG","7024"]],
AUTHORITY["EPSG","6284"]],
***
AUTHORITY["EPSG","4284"]],
***
AUTHORITY["EPSG","9001"]]]
Звучит по меньшей мере странно. У Вас ArcMap какой версии? Ещё, судя по тому, что EPSG упоминается 4 раза, то файлу вашему назначили (или пытались назначить) как минимум 4 разных системы координат, в стандартных настройках ArcGIS есть только один из этих номеров - EPSG: 4284
prc.jpg
prc.jpg (36.26 КБ) 10237 просмотров

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 07:51
trir
"4 разных системы координат"
Нет, просто даны EPSG для : SPHEROID, DATUM, GEOGCS, PROJECTION

Re: Проблема перевода GeoTiff Pulkovo42: ПСК -> ГСК

Добавлено: 06 дек 2013, 09:11
ericsson
wowka1319 писал(а):Зачем? (Qgis давно установлен)
Ну как зачем, чтобы хранить все эти ваши подписи тоже в geotiff или в векторе, грузить их вместе со снимком, получать точное наложение одного на другое, а потом обрезать и экспортировать, как захочется.
Использовать какой-то там атлас из OZIExplorer для проверки координат в профессиональной задаче - вот тут моя очередь спросить "зачем", учитывая что программку эту писал человек, который имеет весьма странные представления об электронной картографии.