Страница 1 из 2
"сморщенные" карты
Добавлено: 20 июл 2016, 08:37
nadiopt
коллеги, мне тут дали задание, связанное с картами. это растр и к нему tifw
есть вполне прямой исходник (страничка из атласа)
в кугисе привязанный растр являет собой "Крик" Мунка:(
картинка тут.
https://yadi.sk/i/L8ZYjGO_tT7Q3
два вопроса: как он так сморщился и нельзя ли как-то эту привязку изничтожить?
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 08:42
Denis Rykov
Не понял, world файл же не может таким образом покорежить растр.
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 08:46
nadiopt
вот и я в непонятках. кладу исходники в личку
upd: привязка не при чем, сам растр сморщенный:(
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 09:07
Denis Rykov
Посмотрел. Файл привязки тут не причем, исходный растр уже вывернут.
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 10:06
nadiopt
в общем, я его обрезала, изничтожила на всякий случай внутреннюю привязку командой gdal_translate -co PROFILE=BASELINE и вяжу заново
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 16:21
Игорь Белов
nadiopt писал(а):есть вполне прямой исходник (страничка из атласа)
Разве это не решение вопроса?
Re: "сморщенные" карты
Добавлено: 20 июл 2016, 19:26
ericsson
nadiopt, ну вы же не вчера родились, сами знаете, что такие вопросы нужно не только скриншотами снабжать, но и выводом gdalinfo.
Re: "сморщенные" карты
Добавлено: 21 июл 2016, 23:53
Boris
Я присоединяюсь к ericsson , хотелось бы понять как так можно растр с помощью привязки уделать. Или все же его кто то таким отсканировал - обработал?
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 08:55
gimran
не район полюса?
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 11:53
Ariki
Должно быть, какая-нибудь трансформация по опорным точкам.
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 13:16
Игорь Белов
Boris писал(а):хотелось бы понять как так можно растр с помощью привязки уделать.
Есть варианты:
- трансформация из оригинальной проекции скана в проекцию, по какой-то причине нужную мне;
- не могу определить исходную проекцию, либо она отсутствует в матаппарате моего ПО, и я объявляю при геопривязке другую.
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 18:15
Boris
Игорь Белов писал(а):Boris писал(а):хотелось бы понять как так можно растр с помощью привязки уделать.
Есть варианты:
- трансформация из оригинальной проекции скана в проекцию, по какой-то причине нужную мне;
- не могу определить исходную проекцию, либо она отсутствует в матаппарате моего ПО, и я объявляю при геопривязке другую.
Уважаемая
nadiopt, не прояснила тему для нашего нездорового интереса, но я готов с вами не согласится. Или мы на разные картинки смотрим?
Нам не известно ПО, которое отображало растр(?) или вектор(?). На самом деле я перечитал, и что это растр не нашел 100% указаний.
Обычное ГИС ПО может так исказить вектор, но и это мало вероятно. Это же не полюс. Это крупный, но не очень кусок российско-китайской границы. Сделать трансформацией в стандартном ГИС ПО дырку в растре?
У меня такое получалось, но только в ScanEx Image Proc при "улучшении" плохой трансформации увеличением полинома до 6-10 степени. Обычное аффинное или проекционное преобразование рвать растр не может.
А, да! Еще такое было, в случае интерполяции сплайном, когда были известные внутренние точки, а опыта, что надо еще фиксировать внешние - не было. Иначе "резиновый лист" превращается в "резиновую чашу".
Но для обычных ГИС при самых косых привязках такого не удавалось.
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 19:58
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
да оно даже в паинте так открывает, с дырой. проблема в том, что "исходную" страницу атласа именно для этой карты мне не дали, в итоге я просто взяла этот крик и привязала заново по опорным точкам. Вопрос был скорее теоретический, что надо было сделать с картой, чтобы она стала выглядеть подобным образом?
Re: "сморщенные" карты
Добавлено: 22 июл 2016, 20:31
Игорь Белов
Вот это особенно интересно:
И я, в отличие от Бориса, ничего не смог разглядеть в слепой картинке. Видимо, материал секретный, раз можно получить его только в личку.
Re: "сморщенные" карты
Добавлено: 24 июл 2016, 03:30
Boris
Игорь Белов писал(а):Вот это особенно интересно:
И я, в отличие от Бориса, ничего не смог разглядеть в слепой картинке. Видимо, материал секретный, раз можно получить его только в личку.
Я то же ничего в личку не получал, только рассматривал картинку из начала темы по ссылке. В ней ясно написано "Россия", "КНР","Россия" и еще раз, при наличии воображения, "КНР". Такое место в мире одно - граница РФ и КНР правее Монголии. С учетом того, что центральная часть сохранила похожие очертания, очевидно, что:
- это трансформация полиномом высокого порядка или трансформация сплайнами. Сплайном скорее. Полиномы как правило сохраняют непрерывность плоскости.
- кто то не зная проекции, от "большого ума" решил использовать "резиной лист", задав только внутренние точки вокруг центральной части. Предполагаю, что в паре точек он перепутал широту с долготой, или как то кардинально неправильно посчитал их.
- все остальное согласно законам расчета сплайнов - "струна" зафиксированная, но противоречиво изогнутая, сохраняет гладкость внутри кольца расчета, а за пределами ее изгибает как угодно.
- пространство - рвется.
В свое время, мне пришлось привязывать план-карту на старую часть Серпухова. Поскольку привязка велась по номерам домов, а их "рисовали" со времен постройки монастыря: какой номер понравился тот и присваивали, то в какой то момент, после установки очередного четного дома на нечетной стороне, пространство свернулось в спираль внутри улицы с домами и превратилось во что то похожее на приведенную картинку на остальной части города.
После этого, было установлено правило - пользоваться сплайн-интерполяцией только с разрешения старших товарищей, обосновав невозможность использования других методов, и обязательно "закреплять" периметр растра двумя рядами точек, выставляемыми точка в точку(исходный - конечный), что бы сплайн не "изгибался" за пределами зоны расчета.
И сейчас часто приходят "картинки" от "не вру ... - это очень точный чертеж, его пан геодезист делал", в которых дома стоят немного не так как на снимке, размеры их слегка иные, углы или одна из граней слегка нафантазированы поскольку из-за забора пришлось часть додумать. А все вместе, если на эти фантазии начинаешь снимок "натягивать", дает такое искажение пространства, что только полином 6 степени уменьшает погрешность до 1-2 пикселей. И вот тогда результат, особенно, если точки есть только в части снимка, получается на загляденье.
Воронки, спирали, волны. НО, повторюсь, что бы порвать растр или сделать в нем дырки - это нужно особо постараться и с точками, и с софтом.