GDAL/OGR

Schulungsunterlagen

Jakob Miksch

jakob@meggsimum.de
@jakobmiksch
@jakobmiksch

• Geoinformatiker
• Anwendungsentwickler
• OSGeo Foundation Charter Member

meggsimum

info@meggsimum.de
@meggsimum
@meggsimum

• meggsimum.de
• Dienstleistungen im Bereich GIS, Webmapping & GDI
• Maßgeschneiderte WebGIS-Lösungen
• Softwarekonzepte und Softwareentwicklung
• Geodaten
• Beratung und Schulung

## GDAL/OGR [www.gdal.org](https://gdal.org/) - C/C++ Bibliothek für Geo-Daten - Umwandlung - Verarbeitung - Analyse - alle gängigen Vektor- und Raster-Formate

## Entwicklung [github.com/OSGeo/gdal](https://github.com/OSGeo/gdal/) <img src="img/contributors.png" alt="" height="350px">

## Unterstützung [gdal.org/sponsors](https://gdal.org/sponsors/index.html) <img src="img/gdal-sponsors.png" alt="" height="400px">

## GDAL/OGR - integriert in [viele Programme](https://gdal.org/software_using_gdal.html#software-using-gdal): - QGIS, ArcGIS, Google Earth, ... - Anbindung zu: - Python, R, Java, ... - Kommandozeile

## Dokumentation <a href="https://gdal.org/" style="font-weight: bold; font-size:larger;">gdal.org</a> <img src="img/gdal-docs.png" alt="" height="400px">

## TLDR - "Spickzettel" - [gdalinfo](https://tldr.ostera.io/gdalinfo) - [gdal2tiles.py](https://tldr.ostera.io/gdal2tiles.py) - [gdalbuildvrt](https://tldr.ostera.io/gdalbuildvrt) - [ogr2ogr](https://tldr.ostera.io/ogr2ogr) - [ogrinfo](https://tldr.ostera.io/ogrinfo) - [ogrmerge.py](https://tldr.ostera.io/ogrmerge.py)

## Kommandozeile direkte Schnittstelle zum Computer <img src="img/terminal.svg" alt="" height="300px">

## Installation Windows [OSGeo4W](https://trac.osgeo.org/osgeo4w/) shell 

## Installation Linux über Paketmanager ```shell $ sudo apt install gdal-bin ```

## Installation Docker - offizielles [Docker image](https://hub.docker.com/r/osgeo/gdal) - Details auf [Blogpost](https://jakobmiksch.eu/post/gdal_ogr/#docker) ```shell $ docker run \ -t \ -i \ -v $(pwd):/data \ osgeo/gdal:alpine-small-latest \ sh ``` 

## QGIS integriert GDAL/OGR <img src="img/qgis_processing.png" alt="" height="400px">

Raster

## gdalinfo Geländemodell von Tirol <img src="img/tirol_dem.png" alt="" height="300px"> von [Land Tirol - data.tirol.gv.at](https://www.data.gv.at/katalog/dataset/land-tirol_tirolgelnde) - [CC BY 4.0](https://creativecommons.org/licenses/by/4.0/deed.de)

## gdalinfo ```shell $ gdalinfo tirol.tif ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text Driver: GTiff/GeoTIFF Files: tirol.tif tirol.tif.aux.xml Size is 43880, 24265 Coordinate System is: PROJCRS["MGI / Austria GK West", BASEGEOGCRS["MGI", DATUM["Militar-Geographische Institut", ELLIPSOID["Bessel 1841",6377397.155,299.1528128, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4312]], CONVERSION["Austria Gauss-Kruger West", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",0, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",10.3333333333333, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",1, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",0, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",-5000000, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["northing (X)",north, ORDER[1], LENGTHUNIT["metre",1]], AXIS["easting (Y)",east, ORDER[2], LENGTHUNIT["metre",1]], USAGE[ SCOPE["Engineering survey, topographic mapping."], AREA["Austria west of 11°50'E of Greenwich (29°30'E of Ferro)."], BBOX[46.77,9.53,47.61,11.84]], ID["EPSG",31254]] Data axis to CRS axis mapping: 2,1 Origin = (-18187.500000000000000,291477.500000000000000) Pixel Size = (5.000000000000000,-5.000000000000000) Metadata: AREA_OR_POINT=Area DataType=Elevation Image Structure Metadata: COMPRESSION=DEFLATE INTERLEAVE=BAND Corner Coordinates: Upper Left ( -18187.500, 291477.500) ( 10d 5'26.52"E, 47d45'43.96"N) Lower Left ( -18187.500, 170152.500) ( 10d 5'44.26"E, 46d40'14.89"N) Upper Right ( 201212.500, 291477.500) ( 13d 0'58.12"E, 47d43'51.96"N) Lower Right ( 201212.500, 170152.500) ( 12d57'42.23"E, 46d38'27.06"N) Center ( 91512.500, 230815.000) ( 11d32'29.27"E, 47d12'37.49"N) Band 1 Block=256x256 Type=Float32, ColorInterp=Gray Min=460.496 Max=3796.002 Minimum=460.496, Maximum=3796.002, Mean=1743.038, StdDev=649.372 NoData Value=-9999 Metadata: RepresentationType=ATHEMATIC STATISTICS_COVARIANCES=421683.5509713619 STATISTICS_MAXIMUM=3796.001953125 STATISTICS_MEAN=1743.0380317268 STATISTICS_MINIMUM=460.49600219727 STATISTICS_SKIPFACTORX=1 STATISTICS_SKIPFACTORY=1 STATISTICS_STDDEV=649.3716585834 ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdalinfo mit JSON Ausgabe ```shell $ gdalinfo \ -json \ tirol.tif ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```json { "description":"tirol.tif", "driverShortName":"GTiff", "driverLongName":"GeoTIFF", "files":[ "tirol.tif", "tirol.tif.aux.xml" ], "size":[ 43880, 24265 ], "coordinateSystem":{ "wkt":"PROJCRS[\"MGI / Austria GK West\",\n BASEGEOGCRS[\"MGI\",\n DATUM[\"Militar-Geographische Institut\",\n ELLIPSOID[\"Bessel 1841\",6377397.155,299.1528128,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n ID[\"EPSG\",4312]],\n CONVERSION[\"Austria Gauss-Kruger West\",\n METHOD[\"Transverse Mercator\",\n ID[\"EPSG\",9807]],\n PARAMETER[\"Latitude of natural origin\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8801]],\n PARAMETER[\"Longitude of natural origin\",10.3333333333333,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8802]],\n PARAMETER[\"Scale factor at natural origin\",1,\n SCALEUNIT[\"unity\",1],\n ID[\"EPSG\",8805]],\n PARAMETER[\"False easting\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8806]],\n PARAMETER[\"False northing\",-5000000,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8807]]],\n CS[Cartesian,2],\n AXIS[\"northing (X)\",north,\n ORDER[1],\n LENGTHUNIT[\"metre\",1]],\n AXIS[\"easting (Y)\",east,\n ORDER[2],\n LENGTHUNIT[\"metre\",1]],\n USAGE[\n SCOPE[\"Engineering survey, topographic mapping.\"],\n AREA[\"Austria west of 11°50'E of Greenwich (29°30'E of Ferro).\"],\n BBOX[46.77,9.53,47.61,11.84]],\n ID[\"EPSG\",31254]]", "dataAxisToSRSAxisMapping":[ 2, 1 ] }, "geoTransform":[ -18187.5, 5.0, 0.0, 291477.5, 0.0, -5.0 ], "metadata":{ "":{ "AREA_OR_POINT":"Area", "DataType":"Elevation" }, "IMAGE_STRUCTURE":{ "COMPRESSION":"DEFLATE", "INTERLEAVE":"BAND" } }, "cornerCoordinates":{ "upperLeft":[ -18187.5, 291477.5 ], "lowerLeft":[ -18187.5, 170152.5 ], "lowerRight":[ 201212.5, 170152.5 ], "upperRight":[ 201212.5, 291477.5 ], "center":[ 91512.5, 230815.0 ] }, "wgs84Extent":{ "type":"Polygon", "coordinates":[ [ [ 10.0903894, 47.7615684 ], [ 10.0953809, 46.6702868 ], [ 12.9610971, 46.6404207 ], [ 13.0154327, 47.7305462 ], [ 10.0903894, 47.7615684 ] ] ] }, "bands":[ { "band":1, "block":[ 256, 256 ], "type":"Float32", "colorInterpretation":"Gray", "min":460.496, "max":3796.002, "minimum":460.496, "maximum":3796.002, "mean":1743.038, "stdDev":649.372, "noDataValue":-9999.0, "metadata":{ "":{ "RepresentationType":"ATHEMATIC", "STATISTICS_COVARIANCES":"421683.5509713619", "STATISTICS_MAXIMUM":"3796.001953125", "STATISTICS_MEAN":"1743.0380317268", "STATISTICS_MINIMUM":"460.49600219727", "STATISTICS_SKIPFACTORX":"1", "STATISTICS_SKIPFACTORY":"1", "STATISTICS_STDDEV":"649.3716585834" } } } ] } ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdalinfo Histogramm ```shell $ gdalinfo \ -hist \ -json \ tirol.tif ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```json { "description":"tirol.tif", "driverShortName":"GTiff", "driverLongName":"GeoTIFF", "files":[ "tirol.tif", "tirol.tif.aux.xml" ], "size":[ 43880, 24265 ], "coordinateSystem":{ "wkt":"PROJCRS[\"MGI / Austria GK West\",\n BASEGEOGCRS[\"MGI\",\n DATUM[\"Militar-Geographische Institut\",\n ELLIPSOID[\"Bessel 1841\",6377397.155,299.1528128,\n LENGTHUNIT[\"metre\",1]]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n ID[\"EPSG\",4312]],\n CONVERSION[\"Austria Gauss-Kruger West\",\n METHOD[\"Transverse Mercator\",\n ID[\"EPSG\",9807]],\n PARAMETER[\"Latitude of natural origin\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8801]],\n PARAMETER[\"Longitude of natural origin\",10.3333333333333,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8802]],\n PARAMETER[\"Scale factor at natural origin\",1,\n SCALEUNIT[\"unity\",1],\n ID[\"EPSG\",8805]],\n PARAMETER[\"False easting\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8806]],\n PARAMETER[\"False northing\",-5000000,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8807]]],\n CS[Cartesian,2],\n AXIS[\"northing (X)\",north,\n ORDER[1],\n LENGTHUNIT[\"metre\",1]],\n AXIS[\"easting (Y)\",east,\n ORDER[2],\n LENGTHUNIT[\"metre\",1]],\n USAGE[\n SCOPE[\"Engineering survey, topographic mapping.\"],\n AREA[\"Austria west of 11°50'E of Greenwich (29°30'E of Ferro).\"],\n BBOX[46.77,9.53,47.61,11.84]],\n ID[\"EPSG\",31254]]", "dataAxisToSRSAxisMapping":[ 2, 1 ] }, "geoTransform":[ -18187.5, 5.0, 0.0, 291477.5, 0.0, -5.0 ], "metadata":{ "":{ "AREA_OR_POINT":"Area", "DataType":"Elevation" }, "IMAGE_STRUCTURE":{ "COMPRESSION":"DEFLATE", "INTERLEAVE":"BAND" } }, "cornerCoordinates":{ "upperLeft":[ -18187.5, 291477.5 ], "lowerLeft":[ -18187.5, 170152.5 ], "lowerRight":[ 201212.5, 170152.5 ], "upperRight":[ 201212.5, 291477.5 ], "center":[ 91512.5, 230815.0 ] }, "wgs84Extent":{ "type":"Polygon", "coordinates":[ [ [ 10.0903894, 47.7615684 ], [ 10.0953809, 46.6702868 ], [ 12.9610971, 46.6404207 ], [ 13.0154327, 47.7305462 ], [ 10.0903894, 47.7615684 ] ] ] }, "bands":[ { "band":1, "block":[ 256, 256 ], "type":"Float32", "colorInterpretation":"Gray", "min":460.496, "max":3796.002, "minimum":460.496, "maximum":3796.002, "mean":1743.038, "stdDev":649.372, "histogram":{ "count":256, "min":453.95579445035293, "max":3802.5421608719171, "buckets":[ 122307, 444451, 692463, 739296, [...] 62, 13 ] }, "noDataValue":-9999.0, "metadata":{ "":{ "RepresentationType":"ATHEMATIC", "STATISTICS_COVARIANCES":"421683.5509713619", "STATISTICS_MAXIMUM":"3796.001953125", "STATISTICS_MEAN":"1743.0380317268", "STATISTICS_MINIMUM":"460.49600219727", "STATISTICS_SKIPFACTORX":"1", "STATISTICS_SKIPFACTORY":"1", "STATISTICS_STDDEV":"649.3716585834" } } } ] } ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdalinfo entferne Datenquellen ```shell $ gdalinfo \ WMS:https://services.meggsimum.de/geoserver/ows ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text Driver: WMS/OGC Web Map Service Files: none associated Size is 512, 512 Subdatasets: SUBDATASET_1_NAME=WMS:https://services.meggsimum.de/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=mgsm-world%3Ane_10m_admin_0_countries&SRS=EPSG:4326&BBOX=-180.0,-90.0,180.0,83.6341 SUBDATASET_1_DESC=World Countries SUBDATASET_2_NAME=WMS:https://services.meggsimum.de/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=mgsm-world%3Ane_10m_populated_places&SRS=EPSG:4326&BBOX=-179.59,-90.0,179.383,82.4833 SUBDATASET_2_DESC=World Places SUBDATASET_3_NAME=WMS:https://services.meggsimum.de/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=mgsm-world%3Ane_10m_rivers_lake_centerlines&SRS=EPSG:4326&BBOX=-164.903,-52.1577,177.52,75.7935 SUBDATASET_3_DESC=World Rivers SUBDATASET_4_NAME=WMS:https://services.meggsimum.de/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=mgsm-world%3Ane_world&SRS=EPSG:4326&BBOX=-189.00892604139275,-114.50858369100003,189.00892604139275,114.508583691 SUBDATASET_4_DESC=World Image SUBDATASET_5_NAME=WMS:https://services.meggsimum.de/geoserver/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=mgsm-ger%3Apostal_codes_germany&SRS=EPSG:4326&BBOX=5.86632,47.2701,15.0419,55.0587 SUBDATASET_5_DESC=Postal Codes Germany Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 512.0) Upper Right ( 512.0, 0.0) Lower Right ( 512.0, 512.0) Center ( 256.0, 256.0) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdalinfo Auswahl eines bestimmten Layers ```shell $ gdalinfo \ WMS:https://services.meggsimum.de/geoserver/ows \ -sd 4 ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```json Driver: WMS/OGC Web Map Service Files: none associated Size is 1073741824, 650512428 Coordinate System is: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 Origin = (-189.008926041392755,114.508583690999998) Pixel Size = (0.000000352056559,-0.000000352056560) Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( -189.009, 114.509) (189d 0'32.13"W,114d30'30.90"N) Lower Left ( -189.009, -114.509) (189d 0'32.13"W,114d30'30.90"S) Upper Right ( 189.009, 114.509) (189d 0'32.13"E,114d30'30.90"N) Lower Right ( 189.009, -114.509) (189d 0'32.13"E,114d30'30.90"S) Center ( 0.0000000, -0.0000000) ( 0d 0' 0.01"E, 0d 0' 0.00"S) Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red Overviews: 536870912x325256214, 268435456x162628107, 134217728x81314054, 67108864x40657027, 33554432x20328513, 16777216x10164257, 8388608x5082128, 4194304x2541064, 2097152x1270532, 1048576x635266, 524288x317633, 262144x158817, 131072x79408, 65536x39704, 32768x19852, 16384x9926, 8192x4963, 4096x2482, 2048x1241, 1024x620 Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green Overviews: 536870912x325256214, 268435456x162628107, 134217728x81314054, 67108864x40657027, 33554432x20328513, 16777216x10164257, 8388608x5082128, 4194304x2541064, 2097152x1270532, 1048576x635266, 524288x317633, 262144x158817, 131072x79408, 65536x39704, 32768x19852, 16384x9926, 8192x4963, 4096x2482, 2048x1241, 1024x620 Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue Overviews: 536870912x325256214, 268435456x162628107, 134217728x81314054, 67108864x40657027, 33554432x20328513, 16777216x10164257, 8388608x5082128, 4194304x2541064, 2097152x1270532, 1048576x635266, 524288x317633, 262144x158817, 131072x79408, 65536x39704, 32768x19852, 16384x9926, 8192x4963, 4096x2482, 2048x1241, 1024x620 ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdal_translate Ein Koordinatensystem zuweisen ```shell $ gdal_translate \ -a_srs EPSG:4326 input.tif \ output.tif ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdal_translate Verkleinern ```shell $ gdal_translate \ -outsize 40% 40% \ earth.tif \ earth40.tif ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdalwarp Umprojizieren ```shell $ gdalwarp \ -t_srs EPSG:4326 \ input.tif \ output.tif ```

## gdalwarp <img src="img/earth.png" alt="" height="350px"> Rasterbild der Erde <img src="img/natural-earth-logo.png" alt="" height="100px">

## gdalwarp Extraktion über Koordinaten ```shell $ gdalwarp \ -te -29.7 35.0 47.1 72.4 \ -te_srs EPSG:4326 \ earth.tif \ europe.tif ```

## gdalwarp Extraktion über Vektor-Layer ```shell $ gdalwarp \ -cutline countries.shp \ -cl countries \ -cwhere "NAME = 'Germany'" \ -crop_to_cutline \ earth.tif \ germany.tif ``` <img src="img/germany.png" alt="" height="200px">

## gdalwarp Extraktion über einen entfernten Vektor-Layer ```shell $ gdalwarp \ -cutline https://raw.githubusercontent.com/martynafford/natural-earth-geojson/master/10m/cultural/ne_10m_admin_0_countries.json \ -cl ne_10m_admin_0_countries \ -cwhere "NAME = 'Switzerland'" \ -crop_to_cutline \ earth.tif \ switzerland.tif ``` <img src="img/switzerland.png" alt="" height="100px">

## gdalwarp Extraktion übere mehrere Vektor-Features ```shell $ gdalwarp \ -cutline countries.shp \ -cwhere "NAME IN ('Germany', 'Spain', 'Italy')" \ -crop_to_cutline \ earth.tif \ many_countries.tif ``` <img src="img/many_countries.png" alt="" height="200px">

## gdalwarp vs, gdal_translate | gdal_translate | gdalwarp | | ----- | ---- | | CRS zuweisen | CRS umwandeln | | Extraktion nach BBOX | Extrakt BBOX/Layer | | Formate umwandeln | | | Bänder tauschen | |

## Virtuelles Raster ```shell $ gdalbuildvrt \ output.vrt \ germany.tif \ switzerland.tif ``` `output.vrt`: <!-- .element: class="fragment" data-fragment-index="1" --> ```xml <VRTDataset rasterXSize="274" rasterYSize="276"> <SRS dataAxisToSRSAxisMapping="2,1">GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]</SRS> <GeoTransform> 5.8666666666481149e+00, 3.3333333333329968e-02, 0.0000000000000000e+00, 5.5033333333336827e+01, 0.0000000000000000e+00, -3.3333333333329995e-02</GeoTransform> <VRTRasterBand dataType="Byte" band="1"> <ColorInterp>Red</ColorInterp> <SimpleSource> <SourceFilename relativeToVRT="1">switzerland.tif</SourceFilename> <SourceBand>1</SourceBand> <SourceProperties RasterXSize="134" RasterYSize="59" DataType="Byte" BlockXSize="134" BlockYSize="20" /> <SrcRect xOff="0" yOff="0" xSize="134" ySize="59" /> <DstRect xOff="3" yOff="217" xSize="134" ySize="59" /> </SimpleSource> <SimpleSource> <SourceFilename relativeToVRT="1">germany.tif</SourceFilename> <SourceBand>1</SourceBand> <SourceProperties RasterXSize="274" RasterYSize="232" DataType="Byte" BlockXSize="274" BlockYSize="9" /> <SrcRect xOff="0" yOff="0" xSize="274" ySize="232" /> <DstRect xOff="0" yOff="0" xSize="274" ySize="232" /> </SimpleSource> </VRTRasterBand> <VRTRasterBand dataType="Byte" band="2"> <ColorInterp>Green</ColorInterp> <SimpleSource> <SourceFilename relativeToVRT="1">switzerland.tif</SourceFilename> <SourceBand>2</SourceBand> <SourceProperties RasterXSize="134" RasterYSize="59" DataType="Byte" BlockXSize="134" BlockYSize="20" /> <SrcRect xOff="0" yOff="0" xSize="134" ySize="59" /> <DstRect xOff="3" yOff="217" xSize="134" ySize="59" /> </SimpleSource> <SimpleSource> <SourceFilename relativeToVRT="1">germany.tif</SourceFilename> <SourceBand>2</SourceBand> <SourceProperties RasterXSize="274" RasterYSize="232" DataType="Byte" BlockXSize="274" BlockYSize="9" /> <SrcRect xOff="0" yOff="0" xSize="274" ySize="232" /> <DstRect xOff="0" yOff="0" xSize="274" ySize="232" /> </SimpleSource> </VRTRasterBand> <VRTRasterBand dataType="Byte" band="3"> <ColorInterp>Blue</ColorInterp> <SimpleSource> <SourceFilename relativeToVRT="1">switzerland.tif</SourceFilename> <SourceBand>3</SourceBand> <SourceProperties RasterXSize="134" RasterYSize="59" DataType="Byte" BlockXSize="134" BlockYSize="20" /> <SrcRect xOff="0" yOff="0" xSize="134" ySize="59" /> <DstRect xOff="3" yOff="217" xSize="134" ySize="59" /> </SimpleSource> <SimpleSource> <SourceFilename relativeToVRT="1">germany.tif</SourceFilename> <SourceBand>3</SourceBand> <SourceProperties RasterXSize="274" RasterYSize="232" DataType="Byte" BlockXSize="274" BlockYSize="9" /> <SrcRect xOff="0" yOff="0" xSize="274" ySize="232" /> <DstRect xOff="0" yOff="0" xSize="274" ySize="232" /> </SimpleSource> </VRTRasterBand> </VRTDataset> ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## Pyramiden / Übersichten <img src="img/build_pyramids.png" alt="" height="250px"> ```shell $ gdaladdo \ -r average \ input.tif \ 2 4 8 16 ``` <!-- .element: class="fragment" data-fragment-index="1" -->

### gdal2tiles ```shell $ gdal2tiles.py \ --zoom=2-5 \ input.tif \ output_folder ``` Ordnerstruktur: <!-- .element: class="fragment" data-fragment-index="1" --> ```text output_folder ├── 2 │ ├── 0 │ │ ├── 0.png │ │ ├── 1.png │ │ ├── 2.png │ │ └── 3.png │ ├── 1 │ │ ├── 0.png │ │ ├── 1.png │ │ ├── 2.png │ │ └── 3.png │ ├── 2 │ │ ├── 0.png │ │ ├── 1.png │ │ ├── 2.png │ │ └── 3.png │ └── 3 │ ├── 0.png │ ├── 1.png │ ├── 2.png │ └── 3.png ├── googlemaps.html ├── leaflet.html ├── openlayers.html └── tilemapresource.xml ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## gdaldem hillshade Schummerung von Erhebungen ```shell $ gdaldem hillshade \ reutte.tif \ hillshade.tif ``` <img src="img/hillshade.png" alt="" height="300px">

## gdaldem slope Hangneigung ```shell $ gdaldem slope \ reutte.tif \ slope.tif ``` <img src="img/slope.png" alt="" height="300px">

## gdaldem aspect Exposition ```shell $ gdaldem aspect \ reutte.tif \ aspect.tif ``` <img src="img/aspect.png" alt="" height="300px">

## gdal_viewshed Sichtfeld ```shell $ gdal_viewshed \ -ox 28648.097 \ -oy 261212.391 \ reutte.tif \ viewshed.tif ``` <img src="img/viewshed.png" alt="" height="300px">

## contour Isolinien ```shell $ gdal_contour \ -a myElevation \ -i 100.0 \ reutte.tif \ contour_lines.gpkg ``` <img src="img/contour-lines.png" alt="" height="300px">

## contour Iso-Polygone ```shell $ gdal_contour \ -i 100.0 \ -p \ reutte.tif \ contour_polygons.gpkg ``` <img src="img/contour-polygons.png" alt="" height="300px">

Vektor

## ogrinfo Informationen zu Vektordaten ```shell $ ogrinfo \ places.shp ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `places.shp' using driver `ESRI Shapefile' successful. 1: places (Point) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## ogrinfo ```shell $ ogrinfo \ -al \ places.shp ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `places.shp' using driver `ESRI Shapefile' successful. Layer name: places Metadata: DBF_DATE_LAST_UPDATE=2021-05-14 Geometry: Point Feature Count: 7343 Extent: (-179.589979, -90.000000) - (179.383304, 82.483323) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["latitude",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["longitude",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 NAME: String (100.0) POP_MAX: Integer (9.0) OGRFeature(places):0 NAME (String) = Colonia del Sacramento POP_MAX (Integer) = 21714 POINT (-57.8400024734013 -34.4799990054175) OGRFeature(places):1 NAME (String) = Trinidad POP_MAX (Integer) = 21093 POINT (-56.9009965601587 -33.5439989373607) OGRFeature(places):2 NAME (String) = Fray Bentos POP_MAX (Integer) = 23279 POINT (-58.3039974719095 -33.1389990288435) OGRFeature(places):3 NAME (String) = Canelones POP_MAX (Integer) = 19698 POINT (-56.2840014932431 -34.5380040066755) OGRFeature(places):4 NAME (String) = Florida POP_MAX (Integer) = 32234 POINT (-56.2149984479942 -34.0990020052172) OGRFeature(places):5 NAME (String) = Bassar POP_MAX (Integer) = 61845 POINT (0.789003573778132 9.26100006831723) OGRFeature(places):6 NAME (String) = Sotouboua POP_MAX (Integer) = 21054 POINT (0.98499646233131 8.5570021330571) OGRFeature(places):7 NAME (String) = Medenine POP_MAX (Integer) = 61705 POINT (10.4166995575967 33.3999990024673) OGRFeature(places):8 NAME (String) = Kebili POP_MAX (Integer) = 19875 POINT (8.9710025379012 33.6899970349206) OGRFeature(places):9 NAME (String) = Tataouine POP_MAX (Integer) = 62577 POINT (10.4667035894155 33.0000031543208) OGRFeature(places):10 NAME (String) = L'Ariana POP_MAX (Integer) = 97687 POINT (10.1999975474793 36.8666731534449) OGRFeature(places):11 NAME (String) = Jendouba POP_MAX (Integer) = 51408 POINT (8.74999861457911 36.5000040594261) OGRFeature(places):12 NAME (String) = Kasserine POP_MAX (Integer) = 76243 POINT (8.71669850332046 35.2167031048242) OGRFeature(places):13 NAME (String) = Sdid Bouzid POP_MAX (Integer) = 42098 POINT (9.50000448225563 35.0166960791497) OGRFeature(places):14 NAME (String) = Siliana POP_MAX (Integer) = 26960 POINT (9.38330159965864 36.0833041252355) OGRFeature(places):15 NAME (String) = Mahdia POP_MAX (Integer) = 45977 POINT (11.040876615322 35.4839130445027) [...] ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## ogrinfo ```shell $ ogrinfo \ places.shp \ -so \ -al ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `places.shp' using driver `ESRI Shapefile' successful. Layer name: places Metadata: DBF_DATE_LAST_UPDATE=2021-05-14 Geometry: Point Feature Count: 7343 Extent: (-179.589979, -90.000000) - (179.383304, 82.483323) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["latitude",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["longitude",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 NAME: String (100.0) POP_MAX: Integer (9.0) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## Shapefile Ordner von Shapefiles ```text shp ├── countries.dbf ├── countries.prj ├── countries.shp ├── countries.shx ├── places.dbf ├── places.prj ├── places.shp └── places.shx ```

## Shapefile Ordner wird als Datenquelle interpretiert ```shell $ ogrinfo shp ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `shp' using driver `ESRI Shapefile' successful. 1: places (Point) 2: countries (Polygon) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## GPX Austauschformat für GPS/GNSS Track ```xml <?xml version='1.0' encoding='UTF-8' standalone='yes' ?> <gpx version="1.1" creator="OsmAndRouterV2" xmlns="http://www.topografix.com/GPX/1/1" xmlns:osmand="https://osmand.net" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.topografix.com/GPX/1/1 http://www.topografix.com/GPX/1/1/gpx.xsd"> <metadata /> <trk> <name>Name of GPX track</name> <trkseg> <trkpt lat="34.7985776" lon="72.5726559"> <ele>2433.1</ele> <time>2021-04-01T06:35:40Z</time> <hdop>4.3</hdop> </trkpt> <trkpt lat="34.7986" lon="72.5726554"> <ele>2429.3</ele> <time>2021-04-01T06:36:00Z</time> <hdop>3.6</hdop> </trkpt> <trkpt lat="34.7985941" lon="72.5726613"> <ele>2433.6</ele> <time>2021-04-01T06:36:21Z</time> <hdop>3.9</hdop> </trkpt> <trkpt lat="34.7986255" lon="72.5727529"> <ele>2429.7</ele> <time>2021-04-01T06:36:41Z</time> <hdop>7.7</hdop> </trkpt> </trkseg> </trk> </gpx> ```

## GPX ```shell $ ogrinfo \ track.gpx ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `track.gpx' using driver `GPX' successful. 1: waypoints (Point) 2: routes (Line String) 3: tracks (Multi Line String) 4: route_points (Point) 5: track_points (Point) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## GPX ```shell $ ogrinfo \ track.gpx \ track_points \ -so ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `track.gpx' using driver `GPX' successful. Layer name: track_points Geometry: Point Feature Count: 902 Extent: (72.544079, 34.776858) - (72.574542, 34.798658) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 track_fid: Integer (0.0) track_seg_id: Integer (0.0) track_seg_point_id: Integer (0.0) ele: Real (0.0) time: DateTime (0.0) magvar: Real (0.0) geoidheight: Real (0.0) name: String (0.0) cmt: String (0.0) desc: String (0.0) src: String (0.0) link1_href: String (0.0) link1_text: String (0.0) link1_type: String (0.0) link2_href: String (0.0) link2_text: String (0.0) link2_type: String (0.0) sym: String (0.0) type: String (0.0) fix: String (0.0) sat: Integer (0.0) hdop: Real (0.0) vdop: Real (0.0) pdop: Real (0.0) ageofdgpsdata: Real (0.0) dgpsid: Integer (0.0) osmand_speed: Real (0.0) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## GeoJSON Austauschformat im Web ```json { "type": "FeatureCollection", "name": "places", "crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:OGC:1.3:CRS84" } }, "features": [ { "type": "Feature", "properties": { "NAME": "Colonia del Sacramento", "POP_MAX": 21714 }, "geometry": { "type": "Point", "coordinates": [ -57.84, -34.47 ] } }, { "type": "Feature", "properties": { "NAME": "Trinidad", "POP_MAX": 21093 }, "geometry": { "type": "Point", "coordinates": [ -56.90, -33.54 ] } }, { "type": "Feature", "properties": { "NAME": "Fray Bentos", "POP_MAX": 23279 }, "geometry": { "type": "Point", "coordinates": [ -58.30, -33.13 ] } }, { "type": "Feature", "properties": { "NAME": "Canelones", "POP_MAX": 19698 }, "geometry": { "type": "Point", "coordinates": [ -56.28, -34.53 ] } } ] } ```

## GeoJSON ```shell $ ogrinfo \ places.geojson \ -al \ -so ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `places.geojson' using driver `GeoJSON' successful. Layer name: places Geometry: Point Feature Count: 7343 Extent: (-179.589979, -90.000000) - (179.383304, 82.483323) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 NAME: String (0.0) POP_MAX: Integer (0.0) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## GeoJSON (entfernt) ```shell ogrinfo \ -so \ -al \ https://raw.githubusercontent.com/martynafford/natural-earth-geojson/master/10m/cultural/ne_10m_populated_places.json ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `https://raw.githubusercontent.com/martynafford/natural-earth-geojson/master/10m/cultural/ne_10m_populated_places.json' using driver `GeoJSON' successful. Layer name: ne_10m_populated_places Geometry: Point Feature Count: 7343 Extent: (-179.589979, -90.000000) - (179.383304, 82.483323) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 SCALERANK: Integer (0.0) NATSCALE: Integer (0.0) LABELRANK: Integer (0.0) FEATURECLA: String (0.0) NAME: String (0.0) NAMEPAR: String (0.0) NAMEALT: String (0.0) DIFFASCII: Integer (0.0) NAMEASCII: String (0.0) ADM0CAP: Real (0.0) CAPIN: String (0.0) WORLDCITY: Real (0.0) MEGACITY: Integer (0.0) SOV0NAME: String (0.0) SOV_A3: String (0.0) ADM0NAME: String (0.0) ADM0_A3: String (0.0) ADM1NAME: String (0.0) ISO_A2: String (0.0) NOTE: String (0.0) LATITUDE: Real (0.0) LONGITUDE: Real (0.0) CHANGED: Real (0.0) NAMEDIFF: Integer (0.0) DIFFNOTE: String (0.0) POP_MAX: Integer (0.0) POP_MIN: Integer (0.0) POP_OTHER: Integer (0.0) RANK_MAX: Integer (0.0) RANK_MIN: Integer (0.0) GEONAMEID: Real (0.0) MEGANAME: String (0.0) LS_NAME: String (0.0) LS_MATCH: Integer (0.0) CHECKME: Integer (0.0) MAX_POP10: Real (0.0) MAX_POP20: Real (0.0) MAX_POP50: Real (0.0) MAX_POP300: Real (0.0) MAX_POP310: Real (0.0) MAX_NATSCA: Real (0.0) MIN_AREAKM: Real (0.0) MAX_AREAKM: Real (0.0) MIN_AREAMI: Real (0.0) MAX_AREAMI: Real (0.0) MIN_PERKM: Real (0.0) MAX_PERKM: Real (0.0) MIN_PERMI: Real (0.0) MAX_PERMI: Real (0.0) MIN_BBXMIN: Real (0.0) MAX_BBXMIN: Real (0.0) MIN_BBXMAX: Real (0.0) MAX_BBXMAX: Real (0.0) MIN_BBYMIN: Real (0.0) MAX_BBYMIN: Real (0.0) MIN_BBYMAX: Real (0.0) MAX_BBYMAX: Real (0.0) MEAN_BBXC: Real (0.0) MEAN_BBYC: Real (0.0) COMPARE: Integer (0.0) GN_ASCII: String (0.0) FEATURE_CL: String (0.0) FEATURE_CO: String (0.0) ADMIN1_COD: Real (0.0) GN_POP: Real (0.0) ELEVATION: Real (0.0) GTOPO30: Real (0.0) TIMEZONE: String (0.0) GEONAMESNO: String (0.0) UN_FID: Integer (0.0) UN_ADM0: String (0.0) UN_LAT: Real (0.0) UN_LONG: Real (0.0) POP1950: Real (0.0) POP1955: Real (0.0) POP1960: Real (0.0) POP1965: Real (0.0) POP1970: Real (0.0) POP1975: Real (0.0) POP1980: Real (0.0) POP1985: Real (0.0) POP1990: Real (0.0) POP1995: Real (0.0) POP2000: Real (0.0) POP2005: Real (0.0) POP2010: Real (0.0) POP2015: Real (0.0) POP2020: Real (0.0) POP2025: Real (0.0) POP2050: Real (0.0) CITYALT: String (0.0) min_zoom: Real (0.0) wikidataid: String (0.0) wof_id: Integer (0.0) CAPALT: Integer (0.0) NAME_EN: String (0.0) NAME_DE: String (0.0) NAME_ES: String (0.0) NAME_FR: String (0.0) NAME_PT: String (0.0) NAME_RU: String (0.0) NAME_ZH: String (0.0) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## WFS Web Feature Service ```shell ogrinfo \ WFS:https://services.meggsimum.de/geoserver/ows ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```shell INFO: Open of `WFS:https://services.meggsimum.de/geoserver/ows' using driver `WFS' successful. Metadata: PROVIDER_NAME=meggsimum – Büro für Geoinformatik 1: mgsm-ger:postal_codes_germany (title: Postal Codes Germany) (Multi Surface) 2: mgsm-world:ne_10m_admin_0_countries (title: World Countries) (Multi Surface) 3: mgsm-world:ne_10m_populated_places (title: World Places) (Multi Point) 4: mgsm-world:ne_10m_rivers_lake_centerlines (title: World Rivers) (Multi Curve) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## WFS Web Feature Service ```shell ogrinfo \ -so \ WFS:https://services.meggsimum.de/geoserver/ows \ ne_10m_admin_0_countries ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ``` INFO: Open of `WFS:https://services.meggsimum.de/geoserver/ows' using driver `WFS' successful. Metadata: PROVIDER_NAME=meggsimum – Büro für Geoinformatik Layer name: mgsm-world:ne_10m_admin_0_countries Metadata: KEYWORD_1=features KEYWORD_2=ne_10m_admin_0_countries TITLE=World Countries Geometry: Multi Surface Feature Count: 255 Extent: (-180.000000, -90.000000) - (180.000000, 83.634100) Layer SRS WKT: GEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], USAGE[ SCOPE["Horizontal component of 3D system."], AREA["World."], BBOX[-90,-180,90,180]], ID["EPSG",4326]] Data axis to CRS axis mapping: 2,1 Geometry Column = geom gml_id: String (0.0) NOT NULL featurecla: String (0.0) scalerank: Integer (0.0) [...] ``` <!-- .element: class="fragment" data-fragment-index="1" -->

### CSV Werte durch Komma getrennt |id |name |breitengrad |laengengrad| |---|----------|------------|-----------| |1 |Salzburg | 47.8 | 13.0 | |2 |Heidelberg| 49.4 | 8.7 | |3 |Trondheim | 63.4 | 10.4 |

### CSV ```shell $ ogrinfo \ cities.csv \ -al ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `cities.csv' using driver `CSV' successful. Layer name: cities Geometry: None Feature Count: 3 Layer SRS WKT: (unknown) id: String (0.0) name: String (0.0) breitengrad: String (0.0) laengengrad: String (0.0) OGRFeature(cities):1 id (String) = 1 name (String) = Salzburg breitengrad (String) = 47.8 laengengrad (String) = 13.0 OGRFeature(cities):2 id (String) = 2 name (String) = Heidelberg breitengrad (String) = 49.4 laengengrad (String) = 8.7 OGRFeature(cities):3 id (String) = 3 name (String) = Trondheim breitengrad (String) = 63.4 laengengrad (String) = 10.4 ``` <!-- .element: class="fragment" data-fragment-index="1" -->

### CSV Geometriespalten zuweisen ```shell $ ogrinfo \ -al \ cities.csv \ -oo X_POSSIBLE_NAMES=laengengrad \ -oo Y_POSSIBLE_NAMES=breitengrad \ -oo KEEP_GEOM_COLUMNS=NO ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `cities.csv' using driver `CSV' successful. Layer name: cities Geometry: Point Feature Count: 3 Extent: (8.700000, 47.800000) - (13.000000, 63.400000) Layer SRS WKT: (unknown) ID: String (0.0) Name: String (0.0) OGRFeature(cities):1 ID (String) = 1 Name (String) = Salzburg POINT (13.0 47.8) OGRFeature(cities):2 ID (String) = 2 Name (String) = Heidelberg POINT (8.7 49.4) OGRFeature(cities):3 ID (String) = 3 Name (String) = Trondheim POINT (10.4 63.4) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

### PostgreSQL / PostGIS ```shell $ ogrinfo PG:dbname=demo ``` Ausgabe: <!-- .element: class="fragment" data-fragment-index="1" --> ```text INFO: Open of `PG:dbname=demo' using driver `PostgreSQL' successful. 1: places (Point) 2: tracks (Multi Line String) 3: track_points (Point) 4: waypoints (Point) 5: routes (Line String) 6: route_points (Point) 7: ne_10m_admin_1_states_provinces_shp (Multi Polygon) 8: ne_10m_admin_0_boundary_lines_disputed_areas (Multi Line String) 9: ne_10m_admin_0_scale_rank_minor_islands (Multi Polygon) 10: ne_10m_admin_0_map_units (Multi Polygon) 11: ne_10m_parks_and_protected_lands_area (Multi Polygon) 12: ne_10m_admin_1_states_provinces_lines (Multi Line String) 13: ne_10m_admin_0_boundary_lines_maritime_indicator (Multi Line String) 14: ne_10m_admin_0_map_subunits (Multi Polygon) 15: ne_10m_admin_0_boundary_lines_land (Multi Line String) 16: ne_10m_roads (Multi Line String) 17: ne_10m_admin_1_states_provinces_lines_shp (Multi Line String) 18: ne_10m_populated_places (Multi Point) 19: ne_10m_admin_0_boundary_lines_map_units (Multi Line String) 20: ne_10m_urban_areas (Multi Polygon) 21: ne_10m_parks_and_protected_lands_line (Multi Line String) 22: ne_10m_admin_1_states_provinces_lakes (Multi Polygon) 23: ne_10m_parks_and_protected_lands_point (Multi Point) 24: ne_10m_admin_0_disputed_areas (Multi Polygon) 25: ne_10m_admin_1_states_provinces (Multi Polygon) 26: ne_10m_admin_0_countries (Multi Polygon) ``` <!-- .element: class="fragment" data-fragment-index="1" -->

### PostgreSQL / PostGIS Verbindung über "Connection-String" ```shell $ ogrinfo \ PG:" dbname=demo host=localhost port=5432 user=postgres password=postgres schemas=public tables=places" ```

## ogr2ogr Umwandlung und Verarbeitung

## ogr2ogr Format umwandeln ```shell $ ogr2ogr \ -f GeoJSON \ output.geojson \ countries.shp ``` neuere GDAL Versionen erkennen Format anhand der Dateiendung<!-- .element: class="fragment" data-fragment-index="1" --> ```shell $ ogr2ogr \ output.geojson \ countries.shp ``` <!-- .element: class="fragment" data-fragment-index="1" -->

## ogr2ogr Umprojizieren ```shell $ ogr2ogr \ -t_srs EPSG:3857 \ output.gpkg \ input.gpkg ```

## ogr2ogr Umprojizieren inklusive Angabe des ursprünglichen Koordinatensystems ```shell $ ogr2ogr \ -s_srs EPSG:4326 \ -t_srs EPSG:3857 \ output.gpkg \ input.gpkg ```

## ogr2ogr Koordinatensystem zuweisen ```shell ogr2ogr \ output.gpkg \ cities.csv \ -oo X_POSSIBLE_NAMES=laengengrad \ -oo Y_POSSIBLE_NAMES=breitengrad \ -a_srs EPSG:4326 ```

## ogr2ogr räumliche Filterung ```shell ogr2ogr \ -spat -13.931 34.886 46.23 74.12 \ output.gpkg \ natural_earth_vector.gpkg ```

## ogr2ogr Filterung mit Bedingungen ```shell ogr2ogr \ -where "\"POP_EST\" < 1000000" \ output.gpkg \ natural_earth_vector.gpkg \ ne_10m_admin_0_countries ``` <img src="img/countries_pop_less_one_million.png" alt="" height="300px">

## ogr2ogr Filterung mit SQL ```shell ogr2ogr \ output.gpkg \ cities.csv \ -sql "SELECT * FROM cities WHERE latitude < 50" \ -oo X_POSSIBLE_NAMES=laengengrad \ -oo Y_POSSIBLE_NAMES=breitengrad \ -a_srs EPSG:4326 ```

## Links - [GDAL Website](https://gdal.org/) offizielle Dokumenation - [Blogpost](https://meggsimum.de/gdal-ogr-geodaten-automatisiert-verarbeiten/) von Jakob Miksch auf meggsimum.de - [Podcast über GDAL](https://podcasts.apple.com/us/podcast/gdal-geospatial-data-abstraction-library/id1452297085?i=1000520451921) von Mapscaping mit Paul Ramsey (auf Englisch) - [Blogpost](https://medium.com/planet-stories/a-gentle-introduction-to-gdal-part-1-a3253eb96082) von Robert Simmon (auf Englisch)

## Schulungen [meggsimum.de/schulungen](https://meggsimum.de/schulungen/) 

Impressum

Autor

Jakob Miksch
meggsimum - Büro für Geoinformatik (Christian Mayer)
Schillerstraße 2a
67122 Mutterstadt
info@meggsimum.de

Lizenz

Diese Folien sind unter CC BY-SA veröffentlicht.