GDAL/OGR

Kätevin tapa kuvata haluttua koordinaattijärjestelmää on EPSG koodi. Suomalaisten koordinaattijärjestelmien EPSG koodit ovat:

Lisää Suomen EPSG koodeja löytyy MML:n sivuilta. Kaikien EPSG koodien tietokanta.

Joskus proj4 kirjasto ei ymmärrä alkuperäisen tiedoston koordinaattijärjestelmää tai tiedosto ei sisällä tätä tietoa. Siinä tapauksessa kannattaa komennossa määritellä myös alkuperäisen aineiston koordinaattijärjestelmä -s_srs muuttujalla.

gdalwarp -s_srs EPSG:3386 -t_srs EPSG:2393 input_kkj0.tif output_kkj3.tif

Jos tiedosto ei sisällä tietoa koordinaattijärjestelmästä, sen voi määritellä a_srs muuttujalla, esimerkkejä löytyy täältä. Huomaa, että näiden komentojon tuloksena tiedostolle vain annetaan koordinaattijärjestelmän tiedot, mutta aineistoa ei muuteta koordinaattijärjestelmään.

Formaattimuunnokset

Rasteriaineistot

Rasteriaineistojen tiedostomuotoja voi vaihtaa gdal_translate komennolla. Luotavan tiedoston formaatin määrittelee -of muuttuja. Sen oletusarvo on GeoTiff, joten sitä ei ole pakko erikseen määritellä. Toisaalta oletusarvoisesti GeoTiff tiedosto luodaan pakkaamattomana, jos haluaa käyttää esimerkiksi LCW pakkausta, se tulisi lisätä -co muuttujan taakse.
gdal_translate input.ecw output.tif
gdal_translate -co "compress=LCW" input.ecw output.tif
gdal_translate -of ECW input.tif output.ecw

Vektoriaineistot

Vektoriaineistojen formaatin vaihtamiseen soveltuu FWToolsissa ogr2ogr komento. Tiedostomuunnosta varten tulee määritellä -f muuttuja. Oletuksena muunnokselle on ESRI shape, joten sitä ei ole pakko erikseen määritellä. MapInfo mif tai tab formaatti valitaan tiedoston nimen mukaan.
ogr2ogr output.shp input.mif
ogr2ogr -f "MapInfo File" output.tab input.shp
ogr2ogr -f "MapInfo File" output.mif input.shp

Jos samassa .tab tiedostossa on sekä pisteitä, viivoja että alueita, ne tallennetaan useampaan .shp tiedostoon. Ks. ohje.

Oletuksena Shape tiedoston String kenttien koko on 80 merkkiä, jota tarvittaessa pitennetään. Isompien paljon teksti-kenttiä sisältävien aineistojen osalta on suositeltava on lisätä komennon loppuun vielä: –lco RESIZE=YES. Silloin String ja Integer kenttien koot säätetään tarkaan sisällön mukaan.

ogr2ogr output.shp input.mif –lco RESIZE=YES

Jo olemassa olevan Shape tiedoston .dbf tiedostoa voi pienentää ogrinfo työkalulla.

ogrinfo file1.dbf -sql "RESIZE file1"

Shape tiedoston merksitö määritellään "-lco ENCODING=XX" asetuksella.

ogr2ogr output.shp input.gml -lco ENCODING=UTF-8

OGR työkaluilla voi tallentaa myös WFS rajapinnalta saatavila olevia aineistoja.

ogr2ogr -f 'ESRI Shapefile' output.shp WFS:"http://geohub.jrc.ec.europa.eu/effis/ows" layername –lco RESIZE=YES

OGR pystyy lukemaan myös OpenStreetMap projektissa käytettäviä formaatteja. Esim. halutessaan ladata omalle koneelle koko Suomen OSM aineiston, yksi helppo tapa on ladata ensin Geofabrik sivuilta Suomen aineisto pbf muodossa. Ja sitten tallentaa se ogr2ogr työkalulla esim. Shape formaattiin (tai SpatialLite tietokantaan ym).

ogr2ogr OSM_aineisto finland-latest.osm.pbf

Jos halutaan kuitenkin vain tiettyjä kohteita OSM-kannassa, sekin on mahdollista. Rajauksen tekemiseksi on hyödyllinen tutustua aineiston rakenteeseen ogrinfo työkalulla ja OpenStreetMap Map Features sivuun. Seuraavat esimerkkikomennot tallentavat Shape tiedostoon kaikki järvet ja kuntarajat.

ogr2ogr -where "natural like 'water'" fin_lakes.shp finland-latest.osm.pbf multipolygons
ogr2ogr -where "admin_level='8'" fin_municipalities.shp finland-latest.osm.pbf multipolygons

OGR työkaluja voidaan käyttää myös CSV tiedostossa olevien koordinaattien siirtämiseksi paikkatietoformaattiin. CSV tiedoston ensimmäisellä rivillä pitäisi olla kenttien nimet. Esimerkiksi CSV tiedosto (vain UTF8), jossa on koordinaatit CoordX ja CoordY nimisissä sarakkeissa voidaan muuttaa Shape formaattiin näin (GDAL 2.0 ja uudempi):

ogr2ogr output.shp input.csv -oo X_POSSIBLE_NAMES=CoordX -oo Y_POSSIBLE_NAMES=CoordY -a_srs EPSG:3067 –lco RESIZE=YES

3D shape > 2D shape

FOR /R %i IN (*.tif) DO gdal_translate %i F:\compress95%~pni.jp2 -of JP2ECW -co target=95

for i in *.shp; do echo ogrinfo $i -sql "CREATE SPATIAL INDEX ON ${i/.shp/}"; done

#!/bin/bash
for i in $(find -name '*.tif')
do
    echo $i
    gdal_translate $i ../infrared_lossless/${i/tif/jp2} -of JP2ECW -co target=0
done