GDAL, SXF и наличие одновременно несколько проекций

[glax2020]

Всех приветствую.

Столкнулся с проблемой и не знаю как к ней подступиться.

Проблема в следующем. При обработке номенклатурного листа O_36_82.sxf (Тверская область) обнаружил что присутствуют метрические данные в разных проекциях.

Например (извлекаю с помощью GDAL):

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

Объект (72310000,"13837_0_БОЛОТА ПРОХОДИМЫЕ")

Region 1
  144
6592329.84 6431779.42
6431877.25 -111111
-111111 6592322.51
6592217.06 6431997.37
6431882.17 -111111
-111111 6592119.02
6592097.77 6431900.6
...

Как можно выше видеть - это все неподходящие значения... -- похоже на метрическую проекцию?

Следующий объект "нормальный" (извлекаю с помощью GDAL):

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

72310000,"13838_0_БОЛОТА ПРОХОДИМЫЕ"

Region 1
  84
34.5724174047872 57.980082614414
34.572516620711 57.9798662282021
34.5726501415347 57.9796843497874
34.573535938479 57.9792527738223
34.5743593617355 57.978686815157
...

Но Панорама умеет и корректно обрабатывает Объект (72310000,"13837_0_БОЛОТА ПРОХОДИМЫЕ"), который я пока не умею обрабатывать:

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

Region  1
  144
34.561189 57.993872
34.561103 57.994752
34.559368 57.995852
34.557666 57.994838
34.557314 57.995008
34.556917 57.995295
34.556519 57.995791
34.555384 57.996011
...

Получается на одном номенклатурном листе O_36_82.sxf метрические объекты в разной проекции? Я имея опыт работы пока только с 10 номенклатурными листами, и только на этом одном O_36_82.sxf столкнулся с этой проблемой, и насчитал подобных объектов порядка 452 (тех которым я не могу присвоить корректные координаты) из 11698 на этом номенклатурном листе O_36_82.sxf.

Обратился естественно к документации ( http://gistoolkit.ru/download/doc/sxf4bin.pdf ), но не смог найти в описании метрики к записи объекта ссылку на проекцию. Допуская что подобная ссылка могла бы прийти через семантику -- но там нет такой ссылки (приходит в семантике одна единственная ссылка на "63 ПРОХОДИМОСТЬ").

Может кто нибудь из Вас сможет мне объяснить или посоветовать что почитать или на что обратить внимание, чтобы разрешить это проблему. Буду чрезмерно благодарен.

С уважением, Александр Глушко
(ИБРАЭ РАН).

[gamm]

похоже, что -111111 это не координаты, а признак отсутствия данных. Детали нужно спрашивать на форуме Панорамы.

[glax2020]

похоже, что -111111 это не координаты, а признак отсутствия данных. Детали нужно спрашивать на форуме Панорамы.

Спасибо, очень ценный совет. Открыл тему на форуме Панорама GDAL, SXF и наличие...

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

[Игорь Белов]

Например (извлекаю с помощью GDAL):

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

Объект (72310000,"13837_0_БОЛОТА ПРОХОДИМЫЕ")

Region 1
  144
6592329.84 6431779.42
6431877.25 -111111
-111111 6592322.51
6592217.06 6431997.37
6431882.17 -111111
-111111 6592119.02
6592097.77 6431900.6
...

Здесь координаты в СК-42 зона 6. Записи с единицами должны интерпретироваться так:

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

6431877.25 -111111
-111111 6592322.51
...
6431882.17 -111111
-111111 6592119.02

надо преобразовать в

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

6592322.51 6431877.25
...
6592119.02 6431882.17

P. S. Красные точки построены по исправленным таким образом координатам. Синие по широтам/долготам.

[Franklin1967]

похоже, что -111111 это не координаты, а признак отсутствия данных. Детали нужно спрашивать на форуме Панорамы.

-11111 это признак отсуствия высоты в точке. Выгружен объект с в трёхмерном виде без наличия высоты в точках

[glax2020]

P. S. Красные точки построены по исправленным таким образом координатам. Синие по широтам/долготам.

Большое спасибо. Буду осмысливать...

[glax2020]

Александр Глушко пишет:
Получается на одном номенклатурном листеO_36_82.sxf метрические объекты в разной проекции?

Да, я уже начал догадываться что здесь проблема в другом. Читая тему Проблема при импорте из файла SXF, Расхождение координат, полученных из sxf, с теми, которые отображаются в ГИС Панорама, возникает представление и здесь суть в том, что присутствуют данные в "устаревшем" формате, с опорой на рамку - т.е. с координатами условными?, в отличие от тех как это уже принято на протяжении последних лет - карт подготовленными с координатами реальными.

Обратно же, для того чтобы хотя диагностировать наличие метрических данных в записи в формате "условных", а не реальных координат - пока мне по документации и по формату записи и метрики не удается найти признак какого рода данные присутствуют в метрике - реальные или условные

[Дмитрий Барышников]

Признак есть:

3.
Флаг наличия реальных координат (2 бита):
ххх00xxx - вся метрика объектов представлена в условной системе координат ( в дискретах );
ххх11xxx - вся метрика объектов представлена в реальных координатах на местности в соответствии с паспортными данными листа (проекция, система координат, единица измерения), данные о масштабе и
дискретности цифрования носят справочный характер.

Но этот признак для всего листа, а не для записи.

Я посмотрел исходники - в текущей версии драйвера информации о координате Z не считывается - код есть, но он не используется при создании OGRGeometry.
А есть информации о версии sxf карты? Если меньше 3 - то поддержки такой версии в текущем драйвере нет.

[Дмитрий Барышников]

Исправил в r27998

[glax2020]

Исправил в r27998

Большое спасибо, Дмитрий!

Теперь все генерируется без замечаний.

[glax2020]

А есть информации о версии sxf карты? Если меньше 3 - то поддержки такой версии в текущем драйвере нет.

Версия sxf карты для бывшего проблемного sxf: 4 ( и в нем как выяснилось были объекты с координатой Z)