GDAL/OGR извлечение самого большого полигона

Ответить
taksadev
Новоприбывший
Сообщения: 8
Зарегистрирован: 18 апр 2018, 18:23
Репутация: 0
Откуда: Земля

GDAL/OGR извлечение самого большого полигона

Сообщение taksadev » 15 июл 2018, 12:17

Здравствуйте!

Есть файл Shp в котором находятся полигоны, один большой и много полигонов значительно меньшей площади. Мне нужно используя команды GDAL/OGR выделить самый большой полигон и сохранить его в отдельный файл.

Подскажите, пожалуйста, как это сделать?
Последний раз редактировалось taksadev 20 июл 2018, 05:02, всего редактировалось 1 раз.

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: GDAL/OGR извлечение самого большого полигона

Сообщение rhot » 15 июл 2018, 13:08

Подсказываю. Используйте OGR SQL.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: GDAL/OGR извлечение самого большого полигона

Сообщение Александр Мурый » 15 июл 2018, 13:21

Пример:

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

ogr2ogr -f "ESRI Shapefile" -dialect SQLite -sql 'SELECT *,ST_Area(geometry) AS area FROM test ORDER BY area DESC LIMIT 1' test1.shp test.shp
Редактор материалов, модератор форума

Аватара пользователя
Denis Rykov
Гуру
Сообщения: 3376
Зарегистрирован: 11 апр 2008, 21:09
Репутация: 529
Ваше звание: Author
Контактная информация:

Re: GDAL/OGR извлечение самого большого полигона

Сообщение Denis Rykov » 15 июл 2018, 18:57

С помощью командных утилит fio и jq:

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

fio cat input.shp | fio calc area 'shape(f.geometry).area' | jq -s 'max_by(.properties.area)'
Spatial is now, more than ever, just another column- The Geometry Column.

taksadev
Новоприбывший
Сообщения: 8
Зарегистрирован: 18 апр 2018, 18:23
Репутация: 0
Откуда: Земля

Re: GDAL/OGR извлечение самого большого полигона

Сообщение taksadev » 20 июл 2018, 05:02

Большое спасибо!

Ответить

Вернуться в «GDAL/OGR»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 3 гостя