Страница 1 из 1

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

Добавлено: 15 июл 2018, 12:17
taksadev
Здравствуйте!

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

Подскажите, пожалуйста, как это сделать?

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

Добавлено: 15 июл 2018, 13:08
rhot
Подсказываю. Используйте OGR SQL.

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

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

Добавлено: 15 июл 2018, 18:57
Denis Rykov
С помощью командных утилит fio и jq:

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

fio cat input.shp | fio calc area 'shape(f.geometry).area' | jq -s 'max_by(.properties.area)'

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

Добавлено: 20 июл 2018, 05:02
taksadev
Большое спасибо!