Привязка мозаики из растров???

0 голосов
спросил 20 Дек, 05 от Vadim (28,120 баллов) в категории Программные продукты Esri
Есть мозаика из n-фрагментов растров, куторую нужно привязать. Есть ли идеи как оптимально быстро и просто сделать это? В итоге должны получить n-привязанных растров. Сшивать всю мозаику в единый файл, а потом привязывать, а потом снова разбивать не получится, так как идоговый файл будет очень большим.

76 Ответы

0 голосов
ответил 26 Дек, 05 от geologic (39,860 баллов)

Ну собстно при проработке геопривязываешься все равно по точкам (по одной, по крайней мере). Но простая логика позиционирует эту методику заведомо ниже по гамбургскому счету. Угадать проекцию куда более почетно. Впрочем, это идеи не мои а по сути Митрича - он у нас первый с проекциями разбирался, пока мы в кадах растры точковали, как в старину водилось :)

Результатами пока не могу, а соображениями запросто. Привязал растр в ГК, УТМ и меркаторе, используя вышеприведенную методику, т.е. начиная с угадыванья размера пикселя. Мысли были такие (слегка повторюсь, прошу пардону):

1. Пиксель при максимальном увеличении равен ровно 0,6 метра в перидоде, это читается элементарно из картинки. Это соответствует реальности на местности, и с таким мастшабом картинка, естественно, хорошо садится и на ГК и на УТМ. Однако и там и там видно, что нужен слабый поворот и небольшое трапецоидное искажение. Т.е. "НЕ ПОПАЛИ".

2. Программа скачивания дает пиксель 1,19, казалось бы странно... И не соответствует местности нифига - все будет растянуто в полтора раза. Это наводит на мысль, что проекция-таки НЕ равнопромежуточная. Ставим стандартный меркатор... Кладем растр... Улетает черт знает куда, наплевать, двигаем лапой "на место"... Садится очень красиво. Линии выдержаны ровно, "трапеции" нет, "поворота" нет, на первый взгляд - оно самое. НО - расстояние искажены примерно в два раза по долготе. Оно и логично - меркатор проекция равноугольная. Т.е. форма в целом похоже меркаторская, а на экране виртуально сжимается до нужного масштаба. UTM не канает, не советую.

Это, впрочем, результаты первой прмерки, на примере снимка самого крупного масштаба, где мало четких элементов привязки. Надо повторять то же на мелком масштабе с использованием крупных элементов - шоссе, ж/д и т.п, где видно схождение-расхождение. Так что по точкам или не по точкам, а советую вязать сначала в меркатора, не в UTM  и уж конешно не в ГК.

 

0 голосов
ответил 27 Дек, 05 от geologic (39,860 баллов)

проверил на примере севера МО - да, все сходится. Именно меркатор, и мало того, тот парень что написал софтину этот меркатор тоже угадал - спасиб ему, не надо вычислять размер пиксела. Волд-файл который генерит GoogleDown, именно меркаторский. Сдвиг правда офигенный, порядка километра откуда-то, зато масштаб тютелька в тютельку и - главное - ориентация. ну и форма растра, конечно. Далее все просто - уточняешь сдвиг по одной любой точке, вся ситуация "красиво" садится на место... Девятка позволяет всю плиту тут же перепроецировать в ГК, где масштаб карты становится близким к реальности, т.е. в два раза меньше. Посадку проверил, кстати, кроме всякой топографии собственным ЖПС-треками - точность с ходу до 50 м, и это при посадке по одной точке. Дальше можно уточняться точками, но уже на "плоскости", не на сферах.

В обчем, задачу для среднего масштаба можно считать решенной. Уточнять датум и прочая конешно нужно, тогда и 50 м дельта возможно уйдет.

0 голосов
ответил 27 Дек, 05 от The_Great_Zopuh (140 баллов)

Извиняюсь за тупость, я не совсем понимаю о чем тут идет речь, но возможно вам поможет статья:

"Как устроен Google Maps"
http://great-zopuh.chat.ru/Programs/MapBuilder/gm.htm

 

0 голосов
ответил 28 Дек, 05 от geologic (39,860 баллов)

Спасибо, четко все расписано. Однако где ж хоть какой картограф, пусть скажет, вот это вот

X = floor (BitmapOrigo[Level] + Lon * PixelsPerLonDegree[Level]);
Y = floor (BitmapOrigo[Level] - 0.5 * ln((1+z)/(1-z)) * PixelsPerLonRadian[Level]),

Дает Меркатор или нет? У нас-то в гисах формулой координаты не задашь, нужно использовать "проекцию", т.е. готовый алгоритм, ее уже настраивать.

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

0 голосов
ответил 29 Дек, 05 от Inna (300 баллов)

Может кто и кинет в меня камень, но у меня гугл на вост сибирь встал по транс меркатору на крассовском со скале фактором 1 без проблем.

0 голосов
ответил 29 Дек, 05 от geologic (39,860 баллов)
Эт дело нехитрое... А не пробовали вообще менять проекцию? Небось тоже подойдет и вас устроит :)))
0 голосов
ответил 29 Дек, 05 от valery (7,040 баллов)
А чего вы от этого хотите?
____________________

image
0 голосов
ответил 22 Фев, 06 от Vadim (28,120 баллов)

С проекцией был не прав, действительно Меркатора. В в волд-файле ставлю 0.6 м и все путем.

 

Волд-файл который генерит GoogleDown, именно меркаторский.
  Где он его генерит? ссылка
http://rock-et-al.webhost.ru/gdown/ не работаетimage
0 голосов
ответил 26 Фев, 06 от geologic (39,860 баллов)

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

- По расчетам 0.67 м получатся для Москвы. Это, кстати, подтверждаетсягугловской масштабной линейкой - на самом детальном уровне 50 метров содержат 75 пикселей. Но надо еще проверять, конечно - пока уверенности нет.

- для привязки в ГИС хорошо все-таки кроме разрешения уметь считать положение растра, в вышеупомянутой ссылке формулы для этого есть, но пока мозги не созрели разобраться. Какую проекцию имел в виду автор софта, неясно - его привязка в меркатор ложится лучше всего, но и она дает ошибку ~ 25 км.

0 голосов
ответил 27 Фев, 06 от Vadim (28,120 баллов)

А ссылка так и не работаетimage. Может кто куда-нить выгрузит программульки, или на мыло если небольшие jahvad@mail.ru  По последним подсчетам ложиться на Москву с размером пиксела 0.598432 метра. Вяжу к основе 1:10 000, составленной по планам 1:2000 (Мосгоргеоирест).

Приведенные формулы:

X = floor (BitmapOrigo[Level] + Lon * PixelsPerLonDegree[Level]);
Y = floor (BitmapOrigo[Level] - 0.5 * ln((1+z)/(1-z)) * PixelsPerLonRadian[Level]),

смутно напоминают уравнения проекции Меркатора. А форма картографического изображения не оставляет никаких сомнений. Будем разбираться

Добро пожаловать на сайт Вопросов и Ответов, где вы можете задавать вопросы по GIS тематике и получать ответы от других членов сообщества.
...