Расчет перекрытий полигонов по заданным параметрам

[geo-aleksey]

Аэросъемка. Из рамок снимков сделаны полигоны.
Задача. Рассчитать перекрытие каждых соседних снимков в процентном содержании. Т.е. первого со вторым, второго с третьим и т.д., и выдать в одном отчете результаты на весь блок.

Есть ли программа, в которую можно импортировать полигоны с привязкой и автоматом получить такую статистику? Ибо руками, например, 12к снимков считать долго.

[trir]

SQL
http://msdn.microsoft.com/en-us/library/bb933832.aspx
http://msdn.microsoft.com/en-us/library/bb933923.aspx

[Boris]

программы безусловно есть. все зависит от того, как заданы полигоны и какой они сложности. если это просто прямоугольники в некой системе координат, ориентированные север-восток, то может хватить и стандартной электронной таблицы. Если что то более сложное, то нужна ГИС.
В любом случае придется использовать либо программирование, либо скриптовый язык ГИС. Для Mapinfo в ее MapBasic-е это около 10 строк:
1) найти пересечения
2) упорядочить их по возрастанию индекса, что бы не посчитать пересечения два раза
3) найти площадь пересечений.
Пункт № 1 - по времени квадратичен к количеству полигонов. Сортировка в зависимости от мощности алгоритма. Поиск площади - линейная по времени процедура, зависит от среднего количества пересечений.

[geo-aleksey]

Полигоны сложные, не прямоугольники. Ну т.е. у них по идее по 4 узла (надо проверить), но они в пространстве искажены после проективных преобразований. СК - WGS 84 / UTM 37 N. Насчет программирования задумывался, но сходу возможностей нет. Правда добрые люди уже помогли с конкретным примером, но способ не сказали. )))

[Boris]

Полигоны сложные, не прямоугольники. Ну т.е. у них по идее по 4 узла (надо проверить), но они в пространстве искажены после проективных преобразований. СК - WGS 84 / UTM 37 N. Насчет программирования задумывался, но сходу возможностей нет. Правда добрые люди уже помогли с конкретным примером, но способ не сказали. )))

Я в Mapinfo могу расписать как делать 1-ый, 2-й пункты стандартным образом, наверно и 3-й можно с помощью шаманства. Если не лень, то приложите образец (или все) данные в каком-нибудь ГИС (полу-ГИС) формате. До стадии список

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

"полигон А";"полигон Б";"площадь пересечения=xx га".

Скорее всего довести не трудно. Вот последующий анализ и т.п. - это уже потребует дополнительной обработки.

[geo-aleksey]

Полигоны можно скачать тут

[thegeo]

Результаты по вашим данным.
Area – в кв. м.
ListID – список ID листов пересекаемых текущим.
ListCross – список значений перекрытий в процентах по отношению к площади текущего листа (округлено до целого). Соответствует списку ID.

[geo-aleksey]

Спасибо большое

[geo-aleksey]

Результаты по вашим данным.

Шикарный результат. Можно узнать способ, это программирование или готовая ГИС?

[thegeo]

Это программа на MapBasic для MapInfo (для v.11.5, хотя нет никаких проблем откомпилировать ее и для предыдущих версий). Программу прилагаю.
Замечания по использованию:
1. принято, что все объекты в обрабатываемой таблице это области;
2. таблица должна иметь обязательные поля: ID (Integer), ListID и ListCross типа Char необходимой длины;
3. перед обработкой в поле ID д.б. записан идентификатор (через rowid())
Перед запуском нужно выделить любой объект в обрабатываемой таблице. Далее будет предложено указать значение минимального перекрытия учитываемого при обработке.
В ваших данных фигурировали также имена листов при необходимости их также можно вытащить. Для примера прилагаю таблицу с этими именами и ограничением по перекрытию 30%.

[geo-aleksey]

Спасибо огромное)