ogr2ogr

ogr2ogr - nám pomůže s konverzemi mezi formáty.

Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update]
               [-select field_list] [-where restricted_where|@filename]
               [-progress] [-sql <sql statement>|@filename] [-dialect dialect]
               [-preserve_fid] [-fid FID] [-limit nb_features]
               [-spat xmin ymin xmax ymax] [-spat_srs srs_def] [-geomfield field]
               [-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def] [-ct string]
               [-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]
               dst_datasource_name src_datasource_name
               [-lco NAME=VALUE] [-nln name]
               [-nlt type|PROMOTE_TO_MULTI|CONVERT_TO_LINEAR|CONVERT_TO_CURVE]
               [-dim XY|XYZ|XYM|XYZM|layer_dim] [layer [layer ...]]

Advanced options :
               [-gt n] [-ds_transaction]
               [[-oo NAME=VALUE] ...] [[-doo NAME=VALUE] ...]
               [-clipsrc [xmin ymin xmax ymax]|WKT|datasource|spat_extent]
               [-clipsrcsql sql_statement] [-clipsrclayer layer]
               [-clipsrcwhere expression]
               [-clipdst [xmin ymin xmax ymax]|WKT|datasource]
               [-clipdstsql sql_statement] [-clipdstlayer layer]
               [-clipdstwhere expression]
               [-wrapdateline][-datelineoffset val]
               [[-simplify tolerance] | [-segmentize max_dist]]
               [-makevalid]
               [-addfields] [-unsetFid] [-emptyStrAsNull]
               [-relaxedFieldNameMatch] [-forceNullable] [-unsetDefault]
               [-fieldTypeToString All|(type1[,type2]*)] [-unsetFieldWidth]
               [-mapFieldType srctype|All=dsttype[,srctype2=dsttype2]*]
               [-fieldmap identity | index1[,index2]*]
               [-splitlistfields] [-maxsubfields val]
               [-resolveDomains]
               [-explodecollections] [-zfield field_name]
               [-gcp ungeoref_x ungeoref_y georef_x georef_y [elevation]]* [-order n | -tps]
               [[-s_coord_epoch epoch] | [-t_coord_epoch epoch] | [-a_coord_epoch epoch]]
               [-nomd] [-mo "META-TAG=VALUE"]* [-noNativeData]

Note: ogr2ogr --long-usage for full help.

Příklad použití:

ogr2ogr _data/katastralni_uzemi.gpkg _data/659673/KATASTRALNI_UZEMI_L.shp

Poznámka

Všimněte si: pořadí souborů je obráceně oproti běžně užívanému. Nejprve výstupní a nakonec vstupní soubor.

Základní použití

Základní použití je převod z jednoho formátu na druhý. Formát může být specifikovaný koncovkou souboru, předponou (např. PG:) nebo pomocí parametru -f (output format). Jako první argument se udává výstup a jako poslední argument vstup.

ogr2ogr -f GPKG _data/katastralni_uzemi.gpkg _data/659673/KATASTRALNI_UZEMI_L.shp

Seznam podporovaných formátů samozřejmě získáme pomocí --formats

ogr2ogr --formats

Výstupní (a vstupní) speciality formátů se dočteme v dokumentaci k jednotlivým formátům, ať už na webových stránkách https://gdal.org/drivers/vector/index.html nebo pomocí přepínače --format

ogr2ogr --format "ESRI Shapefile"

Pokud chceme vidět postup zpracování, přidáme přepínač -progress

ogr2ogr -progress _data/katastralni_uzemi.gpkg _data/659673/KATASTRALNI_UZEMI_L.shp

Přidání prvků vs. jejich přepsání

Pomocí přepínačů -apend nebo -overwrite řídíme, zda se v cílovém datasetu bude aktualizovat celá vrstva nebo se pouze přidají další prvky na konec tabulky.

Celý soubor přepíšeme „natvrdo“ díky přepínači -overwrite

ogr2ogr _data/jihlava_parcely.gpkg _data/ruian_jihlava.gpkg parcely
ogr2ogr -append _data/jihlava_parcely.gpkg _data/ruian_polna.gpkg parcely
ogrinfo _data/jihlava_parcely.gpkg

Porovnejte s

ogr2ogr _data/jihlava_parcely.gpkg _data/ruian_jihlava.gpkg parcely -overwrite
ogr2ogr -update _data/jihlava_parcely.gpkg _data/ruian_polna.gpkg parcely
ogrinfo _data/jihlava_parcely.gpkg 57123

Výběr atributů ze zdrojového souboru

Výběr atributů můžeme ovlivnit parametrem -select

ogr2ogr -select Id,KmenoveCislo,PododdeleniCisla,VymeraParcely,ZpusobyVyuzitiPozemku,KatastralniUzemiKod _data/jihlava_parcely.gpkg _data/ruian_jihlava.gpkg parcely -overwrite
ogrinfo _data/jihlava_parcely.gpkg parcely -so

Výběr konkrétních prvků

Můžeme vybrat konkrétní provek na základě jeho ID -fid (a uložení do existujícího souboru, do nové vrstvy)

ogr2ogr -fid 23214 _data/jihlava_parcely.gpkg -nln prior _data/ruian_jihlava.gpkg parcely
ogrinfo _data/jihlava_parcely.gpkg prior

Můžeme omezit maximální počet prvků -limit

ogr2ogr -limit 10 _data/jihlava_parcely.gpkg -nln parcely_10 _data/ruian_jihlava.gpkg parcely
ogrinfo _data/jihlava_parcely.gpkg parcely_10

Můžeme uplatnit SQL podmínku -where

ogr2ogr -where "DruhPozemkuKod=14" _data/jihlava_parcely.gpkg -nln zahrady _data/ruian_jihlava.gpkg parcely
ogrinfo _data/jihlava_parcely.gpkg zahrady

Plná kontrola nad výběrem prvků pomocí SQL

S parametrem -sql můžeme využít veškeré možnosti, které nám tento jazyk nabízí

ogr2ogr -sql "SELECT * FROM parcely AS pa INNER JOIN katastralniuzemi AS ku ON pa.KatastralniUzemiKod = ku.Kod WHERE ku.nazev = 'Kosov u Jihlavy' AND ST_Area(pa.geom) > 10000" _data/jihlava_parcely.gpkg -nln kosov_1ha _data/ruian_jihlava.gpkg

Čitelněji zapsaný SQL výraz:

SELECT * FROM
        parcely AS pa
INNER JOIN
        katastralniuzemi AS ku
ON
        pa.KatastralniUzemiKod = ku.Kod
WHERE
        ku.nazev = 'Kosov u Jihlavy'
  AND
        ST_Area(pa.geom) > 10000"

Pomohli jsme si funkcí ST_Area, která vypočítá z geometri pa.geom plochu.

Plošné výřezy

ogr2ogr má značně univerzální možnosti jak plošně limitovat výběr prvků.

Nejjednodušší je možnost klasického boundingboxu (obálky). Můžeme dokonce specifikovat souřadnicový systém, ve kterém je filtr nastavený.

ogr2ogr -spat 15.5784 49.3919 15.6013 49.4050 -spat_srs EPSG:4326 _data/jihlava_parcely.gpkg -nln centrum _data/ruian_jihlava.gpkg parcely

Můžeme ale také místo výběru geometrie natvrdo oříznout pomocí -clip

To může být seznam hodnot minx miny maxx, maxy nebo přímo jiný datový zdroj, v našem případě natvrdo napsaný WKT (well known text) trojúhelník.

ogr2ogr -clipsrc "POLYGON ((-670036 -1130583, -668421 -1130579, -669195 -1129067, -670036 -1130583))" _data/jihlava_parcely.gpkg -nln orez _data/ruian_jihlava.gpkg parcely -nlt MULTIPOLYGON

Poznámka

Musíme ještě natvrdo přetypovat typ výstupní geometrie

Generalizace a segmentace

Pomocí parametrů -simplify a segmentize můžeme měnit charakter liniových (a plošných) geometrií

ogr2ogr -simplify 200 _data/jihlava_parcely.gpkg -nln ku_simple _data/ruian_jihlava.gpkg katastralniuzemi -nlt MULTIPOLYGON

Geometrii lze „zvalidnit“ pomocí -makevalid, ale na polygony se musí vzít jiný postup.

Změna typu geometrie

Volbu -nlt už jsme používali. Na tomto místě jenom zmíníme, že kromě přímo datového typu lze ovlivnit, bude-li výstup jako MULTI-geometrie, nebo například zápis pomocí křivek (u formátů, kde toto lze.

Změna souřadnicového systému

-t_srs - cílový souřadnicový systém, -s_srs zdrojový souřadnicový systém - pokud je potřeba (např. Shapefile bez .prj souboru).

ogr2ogr -s_srs EPSG:4326 _data/katastralni_uzemi-wgs84.gpkg _data/659673/KATASTRALNI_UZEMI_L.shp

nebo komplikovaněji

ogr2ogr -t_srs EPSG:3857 -s_srs "+proj=krovak +lat_0=49.5 +lon_0=24.8333333333333 +alpha=30.2881397527778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +towgs84=572.213,85.334,461.94,-4.9732,-1.529,-5.2484,3.5378 +units=m +no_defs +type=crs" _data/katastralni_uzemi-krovak.gpkg _data/659673/KATASTRALNI_UZEMI_L.shp

Specifické formáty

Většina formátů má specifické volby pro čtení a zápis, více lze nalézt v dokumentaci.

OGC GeoPackage

• …

ESRI Shapefile

• Max 2GB
• Kódování
• Délka názvů
• …

ogr2ogr -oo ENCODING=Windows-1250 -lco ENCODING=UTF-8 -lco 2GB_LIMIT=NO -lco SPATIAL_INDEX=YES _data/katastralni_uzemi-utf.shp _data/659673/KATASTRALNI_UZEMI_L.shp

CSV

• …