"сморщенные" карты

[nadiopt]

коллеги, мне тут дали задание, связанное с картами. это растр и к нему tifw
есть вполне прямой исходник (страничка из атласа)
в кугисе привязанный растр являет собой "Крик" Мунка:(
картинка тут. https://yadi.sk/i/L8ZYjGO_tT7Q3
два вопроса: как он так сморщился и нельзя ли как-то эту привязку изничтожить?

[Denis Rykov]

Не понял, world файл же не может таким образом покорежить растр.

[nadiopt]

вот и я в непонятках. кладу исходники в личку
upd: привязка не при чем, сам растр сморщенный:(

[Denis Rykov]

Посмотрел. Файл привязки тут не причем, исходный растр уже вывернут.

[nadiopt]

в общем, я его обрезала, изничтожила на всякий случай внутреннюю привязку командой gdal_translate -co PROFILE=BASELINE и вяжу заново

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

есть вполне прямой исходник (страничка из атласа)

Разве это не решение вопроса?

[ericsson]

nadiopt, ну вы же не вчера родились, сами знаете, что такие вопросы нужно не только скриншотами снабжать, но и выводом gdalinfo.

[Boris]

Я присоединяюсь к ericsson , хотелось бы понять как так можно растр с помощью привязки уделать. Или все же его кто то таким отсканировал - обработал?

[gimran]

не район полюса?

[Ariki]

Должно быть, какая-нибудь трансформация по опорным точкам.

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

хотелось бы понять как так можно растр с помощью привязки уделать.

Есть варианты:

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

[Boris]

хотелось бы понять как так можно растр с помощью привязки уделать.

Есть варианты:

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

Уважаемая nadiopt, не прояснила тему для нашего нездорового интереса, но я готов с вами не согласится. Или мы на разные картинки смотрим?
Нам не известно ПО, которое отображало растр(?) или вектор(?). На самом деле я перечитал, и что это растр не нашел 100% указаний.
Обычное ГИС ПО может так исказить вектор, но и это мало вероятно. Это же не полюс. Это крупный, но не очень кусок российско-китайской границы. Сделать трансформацией в стандартном ГИС ПО дырку в растре?
У меня такое получалось, но только в ScanEx Image Proc при "улучшении" плохой трансформации увеличением полинома до 6-10 степени. Обычное аффинное или проекционное преобразование рвать растр не может.
А, да! Еще такое было, в случае интерполяции сплайном, когда были известные внутренние точки, а опыта, что надо еще фиксировать внешние - не было. Иначе "резиновый лист" превращается в "резиновую чашу".
Но для обычных ГИС при самых косых привязках такого не удавалось.

[nadiopt]

прошу прощения, было очень не до того
Борис, это растр в формате тифф
вот гдалинфо

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

Files: D:\far_east\arial_oak\sourse\china_1937.tif
       D:\far_east\arial_oak\sourse\china_1937.tif.ovr
       D:\far_east\arial_oak\sourse\china_1937.tif.aux.xml
Size is 1133, 1313
Coordinate System is `'
Metadata:
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=96.012001
  TIFFTAG_YRESOLUTION=96.012001
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 1313.0)
Upper Right ( 1133.0,    0.0)
Lower Right ( 1133.0, 1313.0)
Center      (  566.5,  656.5)
Band 1 Block=1133x1313 Type=Byte, ColorInterp=Red
  Overviews: 567x657, 284x329, 142x165
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 567x657, 284x329, 142x165
Band 2 Block=1133x1313 Type=Byte, ColorInterp=Green
  Overviews: 567x657, 284x329, 142x165
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 567x657, 284x329, 142x165
Band 3 Block=1133x1313 Type=Byte, ColorInterp=Blue
  Overviews: 567x657, 284x329, 142x165
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 567x657, 284x329, 142x165
Band 4 Block=1133x1313 Type=Byte, ColorInterp=Alpha
  Overviews: 567x657, 284x329, 142x165

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

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

Вот это особенно интересно:

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

Coordinate System is `'

И я, в отличие от Бориса, ничего не смог разглядеть в слепой картинке. Видимо, материал секретный, раз можно получить его только в личку.

[Boris]

Вот это особенно интересно:

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

Coordinate System is `'

И я, в отличие от Бориса, ничего не смог разглядеть в слепой картинке. Видимо, материал секретный, раз можно получить его только в личку.

Я то же ничего в личку не получал, только рассматривал картинку из начала темы по ссылке. В ней ясно написано "Россия", "КНР","Россия" и еще раз, при наличии воображения, "КНР". Такое место в мире одно - граница РФ и КНР правее Монголии. С учетом того, что центральная часть сохранила похожие очертания, очевидно, что:

1. это трансформация полиномом высокого порядка или трансформация сплайнами. Сплайном скорее. Полиномы как правило сохраняют непрерывность плоскости.
2. кто то не зная проекции, от "большого ума" решил использовать "резиной лист", задав только внутренние точки вокруг центральной части. Предполагаю, что в паре точек он перепутал широту с долготой, или как то кардинально неправильно посчитал их.
3. все остальное согласно законам расчета сплайнов - "струна" зафиксированная, но противоречиво изогнутая, сохраняет гладкость внутри кольца расчета, а за пределами ее изгибает как угодно.
4. пространство - рвется.

В свое время, мне пришлось привязывать план-карту на старую часть Серпухова. Поскольку привязка велась по номерам домов, а их "рисовали" со времен постройки монастыря: какой номер понравился тот и присваивали, то в какой то момент, после установки очередного четного дома на нечетной стороне, пространство свернулось в спираль внутри улицы с домами и превратилось во что то похожее на приведенную картинку на остальной части города.
После этого, было установлено правило - пользоваться сплайн-интерполяцией только с разрешения старших товарищей, обосновав невозможность использования других методов, и обязательно "закреплять" периметр растра двумя рядами точек, выставляемыми точка в точку(исходный - конечный), что бы сплайн не "изгибался" за пределами зоны расчета.
И сейчас часто приходят "картинки" от "не вру ... - это очень точный чертеж, его пан геодезист делал", в которых дома стоят немного не так как на снимке, размеры их слегка иные, углы или одна из граней слегка нафантазированы поскольку из-за забора пришлось часть додумать. А все вместе, если на эти фантазии начинаешь снимок "натягивать", дает такое искажение пространства, что только полином 6 степени уменьшает погрешность до 1-2 пикселей. И вот тогда результат, особенно, если точки есть только в части снимка, получается на загляденье.
Воронки, спирали, волны. НО, повторюсь, что бы порвать растр или сделать в нем дырки - это нужно особо постараться и с точками, и с софтом.