Объединение космоснимка, кварт. сетки, поквартальных итогов

[bolotoved]

1. Давайте пока классификацию отложим, а то у вас и так задача многосоставная. (В Envi можно делать классификацию, но для этого нужно знать местность, где у вас что, сама енви не знает, но зная как на снимке выглядят те или иные объекты можно автомаически выявить их на всем снимке, что иногда облегчает работу, по сравнению с ручной обводкой, а иногда -- только усложняет).
2. На счет смещения, нужно выяснить, систематическая ли ошибка, т.е. по всей площади вектора все сдвинуто на одинаковое расстояние и в одинаковом направлении, или наблюдается разнобой.
Вообще какое происхождение вектора? Может он изначально делался не по космосъемке и привязывался к какой-нибудь не очень качественной топооснове?
3. Снимок с перепутанными слоями (тем более если их у вас 5) лучше сшить заново, замучаетесь комбинации перебирать, да и в мозаику сшить вы его уже не сможете адекватно: сшивка идет послойно.

[bolotoved]

Кстати, снимки тоже могут быть привязаны не очень точно. Ориентироваться надо в первую очередь на наземные данные с GPS. Если у вас есть треки на эту территорию, то в первую очередь нанесите на снимок эти треки и точки перекрестий дорог, профилей и т.п. и посмотрите, что привязано лучше, снимок или вектор.

[nadiopt]

у него лесоустройство, делалось соответственно по АФС и топухам. Это Новгородская область, и устроена она плохо, смещения 100-150 м - обычное дело
я все пытаюсь добиться, чтобы Ночной Лесник соизволил уже проверить точность привязки Рапида хотя бы по Ландсату и решить, что сидит лучше - рапид или вектор. Но сижу глубоко в игноре с этой мыслью

[Ночной Лесник]

Nadiopt, вечером постараюсь,а какой заказать снимок? Lansat legacy ETM+ Pan Mosaics?

[nadiopt]

можно просто 5-ый ландсат свежий, чтоб мозаику не качать:)

[Ночной Лесник]

Ну вот взял лансат для этого участка....расстояние между снимками 36 км...
видимо делаю что то не так

[bolotoved]

У вас есть GPS треки на данную территорию? если нет, о можно взять дороги с OSM -- чаще всего они делаются именно по GPS-трекам. По субъектам РФ, сконвертированные уже Shape-формат лежат тут:http://gis-lab.info/projects/osm-export.html
Скачайте шейп на вашу территорию, наложите дорожную сеть на снимок и посмотрите, насколько адекватно привязан сам ландсат, насколько -- ваши квартальные сетки. (ландсат тоже может быть плохо привязан).

[nadiopt]

нужно понять, что куда отлетает. У вас снимки и квартальная сеть хоть как-то коррелируют или полностью разлетаются?

[ericsson]

Вывод gdalinfo для растров и ogrinfo для векторных данных - в студию.
Так будет проще.

[Ночной Лесник]

нужно понять, что куда отлетает. У вас снимки и квартальная сеть хоть как-то коррелируют или полностью разлетаются?

если условно принять мозаику за эталон то получается следующее:
1. Квартальная сеть - вектор и собранная мозаика из 3х снимков рапид похожи, (просеки видны на растре, вроде как к делу), кварталка "уплывает" метров на 100 северо восток 45 градусов вот так :

2.Снимок Lansat legacy ETM+ Pan Mosaics лежит на востоке от всего вышеперечисленного в 36 километрах.
По совету bolotoved скачал шейп региона, открыл кугисом файл novgorod.qgs
открылся проект с кучей слоев, выключил все кроме дорог.
Затем, Добавить растровый слой, открываю свежескачанный Lansat legacy ETM+ Pan Mosaics а в ответ тишина.(просто никаких движений, выключаю отображение всех слоев кроме Lansat указываю "Увеличить до границ слоя, ноль результата..
То же самое с Рапидовской мозаикой.

[Ночной Лесник]

Вывод gdalinfo для растров и ogrinfo для векторных данных - в студию.
Так будет проще.

вот что выдает...
Кугис, свежий устанавливал по стандарту

[bolotoved]

Скорее всего у вас не стоит крыжик "включить перепроецирование на лету", в таком режиме нужно вначале открывать растр, а потом вектор. В противном случае, растры не перепроецируются. Поставьте крыжик Установки/свойства/система координат.
Не загружайте весь проект OSM, а загрузите слой дорог из него (-roads-l.shp) (Слой/Добавить векторный слой) в уже созданный вами проект со снимками и сеткой.

[Ночной Лесник]

Скорее всего у вас не стоит крыжик "включить перепроецирование на лету", в таком режиме нужно вначале открывать растр, а потом вектор. В противном случае, растры не перепроецируются. Поставьте крыжик Установки/свойства/система координат.

Крыжик ставлю в свойствах проекта сразу же до "открывания" слоев.

Не загружайте весь проект OSM, а загрузите слой дорог из него (-roads-l.shp) (Слой/Добавить векторный слой) в уже созданный вами проект со снимками и сеткой.

делал так, пишет дословно:
-roads-l.shp Не является действительным источником данных
После некоторых попыток:

[Ночной Лесник]

вот так получается, примерно 88 метров смещение (мерил линейкой Qgis)

[Ночной Лесник]

На самом деле векторных карту меня несколько, подскажите как оставить только контуры кварталов, без заливки.
Выглядит примерно вот так: