Форумы GIS-Lab.info

Коллеги, помогите с датумом для 9 зоны ГК. Получил из-под Арки растр с параметрами

PROJCS["Pulkovo_1942_GK_Zone_9",
    GEOGCS["GCS_Pulkovo_1942",
        DATUM["unknown",
            SPHEROID["Krasovsky_1940",6378245,298.3,
                AUTHORITY["EPSG","0"]],
            AUTHORITY["EPSG","6284"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",51],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",9500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]

у меня программа распознает датум как нестандартный, а стандартным считает датум с параметрами

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

PROJCS["unnamed",
    GEOGCS["Krassovsky, 1942",
        DATUM["unknown",
            SPHEROID["krass",6378245,298.3],
            TOWGS84[24,-123,-94,0.02,-0.25,-0.13,1.1]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",51],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",9500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]

это вывод gdalinfo, судя по обсуждению на форуме тут, это то, что принято в Мапинфо, и по словам уважаемого ErnieBoyd в этом обсуждении, возможно имеет официальные корни.

Вопрос, собственно, какие параметры TOWGS84 предполагает первый датум? и почему имена сфероидов разные, если параметры одинаковые?

EPSG:6284 — так называемый геодезический датум. Какую следует использовать трансформацию для перехода от него к WGS 84? Надо смотреть на происхождение данных.

gamm писал(а):Вопрос, собственно, какие параметры TOWGS84 предполагает первый датум? и почему имена сфероидов разные, если параметры одинаковые?

ArcGIS не предполагает датума: выбор между трансформациями за пользователем. А имена не имеют значения. Сам GDAL вон пишет "Krassowsky".

ErnieBoyd писал(а):EPSG:6284 — так называемый геодезический датум. Какую следует использовать трансформацию для перехода от него к WGS 84? Надо смотреть на происхождение данных.

данные из АркГИС. Что нужно спросить у поставщика (чтобы он посмотрел в Арке), чтобы узнать все параметры проекции для PROJ4? Может, они выставили в Арке "Pulkovo_1942_GK_Zone_9", и для этого варианта какие-то параметры Арка сама использует автоматом? как было хорошо, пока вместо нормальных параметров проекции не стали писать в TIFF эту белиберду EPSG ...

gamm писал(а):Что нужно спросить у поставщика (чтобы он посмотрел в Арке), чтобы узнать все параметры проекции для PROJ4?

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

gamm писал(а):данные из АркГИС.

Я не об этом. Предположим, речь идёт о сантиметрвой точности. Геодезисты заложили пункты, увязали с ГГС, передали координаты на опознаки, обработали фотограмметрию, выдали крупномасштабные ортофотопланы, выполнили полевое дешифрирование. Поставщик всё векторизовал и наделал векторных слоёв. Заметьте, всё началось с координат пунктов ГГС и закончилось векторами в одной и той же системе координат. Какие там были параметры трансформации датума? Не имеет значения, всё равно все результаты получены из относительных измерений от исходных пунктов.

В общем, не в программе дело. Важны данные и задачи.

ErnieBoyd писал(а):Не имеет значения, всё равно все результаты получены из относительных измерений от исходных пунктов.

переформулирую слегка - а какие параметры датума у генштабовских карт? скорее всего EPSG:6284 был просто назначен, данные цифровались с карт. Нужно попасть в снимки, которые в UTM, точность низкая (порядка 10-15 метров), но разница для разных датумов больше.

gamm писал(а):а какие параметры датума у генштабовских карт? … Нужно попасть в снимки, которые в UTM, точность низкая (порядка 10-15 метров), но разница для разных датумов больше.

Ясно. Параметры датума в этом вопросе вторая проблема.

Первая проблема — разное происхождение данных. Недавно на форуме всплывал факт, что карты разных лет съёмки различаются местами весьма заметно. И это естественно: точность привязки опознаков для карт масштабов 1:10000–1:25000 составляла 4 метра, а ведь они привязывались в поле далеко не для каждой пары снимков. Дальше накопление ошибок измерений на аналоговых фотограмметрических приборах, ошибки черчения/гравировки пластика. Со своей стороны снимки, которые в UTM, обычно используют слишком грубые модели рельефа для орторектификации. Короче, точность попадания в 10–15 метров не низкая, а недостижимо высокая. Никакой датум не спасёт.

Теперь о датуме. Предположим, снимки в UTM трансформировались по точкам с генштабовских карт. Однако они на WGS 84, значит, должно было использоваться какое-то трансформирование "Pulkovo 1942 to WGS 84". По моему скромному мнению, для территории России мог использоваться один из трёх вариантов. Зарубежные партнёры, вероятно, использовали бы EPSG:15865 "Pulkovo 1942 to WGS 84 (16)" либо EPSG:1267 "Pulkovo 1942 to WGS 84 (17)". Первый был действительно первым, официально созданным ЦНИИГАиК для СССР, известен по публикации в журнале «Геодезия и картография» 1993 года. Второй рекомендован ГОСТ Р 51794—2001 и, в частности, привязан к "Pulkovo 1942" в PROJ.4. Российские специалисты могли бы использовать EPSG:5044 "Pulkovo 1942 to WGS 84 (20)", который появился в ГОСТ Р 51794—2008 и перешёл в ГОСТ 32453—2013. С практической точки зрения все три варианта равнозначны, разница получается в первых метрах.

я так понимаю, смотреть параметры нужно здесь:

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

 <gml:sourceCRS xlink:href="urn:ogc:def:crs:EPSG::4284" />
 <gml:targetCRS xlink:href="urn:ogc:def:crs:EPSG::4326" />

Можно так. Смотрим EPSG:15865 "Pulkovo 1942 to WGS 84 (16)". В самом начале читаем:

Method: Coordinate Frame Rotation (geog2D domain)
…
Description: (25.0, -141.0, -78.5, 0.0, -0.35, -0.736, 0.0)

Поскольку PROJ.4 использует метод "Position Vector transformation", меняем знаки параметров вращения и пишем:

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

+towgs84=25,-141,-78.5,0,0.35,0.736,0

Аналогично получаем для EPSG:1267 "Pulkovo 1942 to WGS 84 (17)"

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

+towgs84=23.92,-141.27,-80.9,0,0.35,0.82,-0.12

и для EPSG:5044 "Pulkovo 1942 to WGS 84 (20)"

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

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22

ErnieBoyd писал(а):Аналогично получаем для EPSG:1267 "Pulkovo 1942 to WGS 84 (17)"
Код: Выделить всё
+towgs84=23.92,-141.27,-80.9,0,0.35,0.82,-0.12

Спасибо!

у меня программа выдала именно этот вариант (обозначив его Custom), просто возникли сомнения - увидев EPSG:6284 в самом файле, стал искать, что про него известно.

Форумы GIS-Lab.info

Не смог найти параметры EPSG:6284

Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284

Re: Не смог найти параметры EPSG:6284