Difference between revisions of "Croatia/Kratki vodici/digitalni-model-terena-iz-csv-datoteke-pomocu-gdala"

From OSGeo
Jump to navigation Jump to search
Line 1: Line 1:
 +
[[Category:Croatia/Kratki_vodici|D]]
 +
 
= Digitalni model terena iz CSV datoteke pomoću GDALa =
 
= Digitalni model terena iz CSV datoteke pomoću GDALa =
  
−
[[Category:Croatia/Kratki_vodici|D]]
+
Kako bi uspješno pratili kratki vodič morate zadovoljiti sljedeće:
 +
# instalirati GDAL programski paket za vaš operativni sustav
 +
# preuzeti i otpakirati ZIP datoteku koja sadrži podatke u CSV formatu
 +
 
 +
Datoteka ''podaci.csv'' sadrži sljedeće zapise:
 +
<pre>
 +
X,Y,Z
 +
5585884.50,5063582.05,63.45
 +
5585884.67,5063582.95,63.90
 +
5585884.75,5063583.81,63.97
 +
5585884.65,5063584.70,64.18
 +
5585884.27,5063585.74,64.30
 +
5585883.73,5063586.89,64.37
 +
...
 +
</pre>
 +
 
 +
Zaglavlje datoteke definira koje informacije sadrži koja kolona, u ovom slučaju '''X''' i '''Y''' koordinatu te visinu '''Z'''.
 +
 
 +
Tijek procesa:
 +
# napraviti raster pomoću ''gdal_grid''
 +
## definirati CSV kao virtualni format za OGR
 +
## odabranom metodom interpolacije napraviti raster
 +
# izraditi slojnice pomoću ''gdal_contour''
 +
 
 +
 +
== ''Gridding'' pomoću gdal_grid ==
 +
 
 +
=== VRT virtualni format ===
 +
 
 +
[http://www.gdal.org/gdal_vrttut.html VRT] je tzv. virtualni format, odnosno tekstualna datoteka s metapodacima koji opisuju neki set podataka.
 +
 
 +
<pre>
 +
<OGRVRTDataSource>
 +
  <OGRVRTLayer name="podaci">
 +
    <SrcDataSource>podaci.csv</SrcDataSource>
 +
    <GeometryType>wkbPoint</GeometryType>
 +
    <GeometryField encoding="PointFromColumns" x="X" y="Y" z="Z"/>
 +
  </OGRVRTLayer>
 +
</OGRVRTDataSource>
 +
</pre>
 +
 
 +
Zapis ''<OGRVRTLayer name="podaci">'' definira naziv sloja podataka, ''<GeometryType>'' tip podataka, koji je u ovom slučaju točka, i ''<GeometryField'' koji je definiran zapis, odnosno u kojoj koloni se nalazi koji podatak. Naziv sloja podataka bi trebao biti isti kao naziv datoteke bez ekstenzije, dakle ''podaci'' i ''podaci''.
 +
 
 +
Tekstualnu datoteku možemo nazvati kako želimo, ali za potrebe ovog vodiča neka bude ''podaci.vrt''.
 +
 
 +
Nakon spremanja datoteke možemo iskoristiti alat ''ogrinfo'' kako bi saznali osnovne informacije o ovo skupu podataka:
 +
<pre>
 +
ogrinfo podaci.vrt podaci -so
 +
INFO: Open of `podaci.vrt'
 +
      using driver `VRT' successful.
 +
 
 +
Layer name: podaci
 +
Geometry: Point
 +
Feature Count: 4448
 +
Layer SRS WKT:
 +
(unknown)
 +
X: String (0.0)
 +
Y: String (0.0)
 +
Z: String (0.0)
 +
</pre>
 +
 
 +
Primjerice znamo da skup podataka ima 4448 točaka i da nema definiran prostorni referentni sustav (SRS). SRS možemo definirati tako da dodamo:
 +
<pre>
 +
<LayerSRS>+proj=tmerc +ellps=bessel +k_o=0.9999 +x_0=5500000 +lon_0=15 +units=m</LayerSRS>
 +
</pre>
 +
čime bi definirali 5 zonu GaussKruger projekcije bez transformacijskih parametara prema wgs84 u PROJ.4 formatu.
 +
 
 +
=== ''Gridding'' i odabir metode interploacije ===
 +
 
 +
Već je spomenuto kako ''gdal_grid'' podržava tri metode interpolacije (link). Osim toga za svaku se metodu mogu kontrolirati [http://www.gdal.org/gdal_grid.html različiti parametri]. Najkorisnije je spomenuti da bez definicije parametara ''gdal_grid'' u obzir uzima sve točke zadanog područja, što vrlo vjerojatno nije željno ponašanje.
 +
 
 +
Od parametara su najzanimljiviji:

Revision as of 06:18, 7 December 2010


Digitalni model terena iz CSV datoteke pomoću GDALa

Kako bi uspješno pratili kratki vodič morate zadovoljiti sljedeće:

  1. instalirati GDAL programski paket za vaš operativni sustav
  2. preuzeti i otpakirati ZIP datoteku koja sadrži podatke u CSV formatu

Datoteka podaci.csv sadrži sljedeće zapise: 

X,Y,Z
5585884.50,5063582.05,63.45
5585884.67,5063582.95,63.90
5585884.75,5063583.81,63.97
5585884.65,5063584.70,64.18
5585884.27,5063585.74,64.30
5585883.73,5063586.89,64.37
...

Zaglavlje datoteke definira koje informacije sadrži koja kolona, u ovom slučaju X i Y koordinatu te visinu Z.

Tijek procesa:

  1. napraviti raster pomoću gdal_grid
    1. definirati CSV kao virtualni format za OGR
    2. odabranom metodom interpolacije napraviti raster
  2. izraditi slojnice pomoću gdal_contour


Gridding pomoću gdal_grid

VRT virtualni format

VRT je tzv. virtualni format, odnosno tekstualna datoteka s metapodacima koji opisuju neki set podataka. 

<OGRVRTDataSource>
  <OGRVRTLayer name="podaci">
     <SrcDataSource>podaci.csv</SrcDataSource> 
     <GeometryType>wkbPoint</GeometryType>
     <GeometryField encoding="PointFromColumns" x="X" y="Y" z="Z"/> 
  </OGRVRTLayer>
</OGRVRTDataSource>

Zapis <OGRVRTLayer name="podaci"> definira naziv sloja podataka, <GeometryType> tip podataka, koji je u ovom slučaju točka, i <GeometryField koji je definiran zapis, odnosno u kojoj koloni se nalazi koji podatak. Naziv sloja podataka bi trebao biti isti kao naziv datoteke bez ekstenzije, dakle podaci i podaci.

Tekstualnu datoteku možemo nazvati kako želimo, ali za potrebe ovog vodiča neka bude podaci.vrt.

Nakon spremanja datoteke možemo iskoristiti alat ogrinfo kako bi saznali osnovne informacije o ovo skupu podataka: 

ogrinfo podaci.vrt podaci -so
INFO: Open of `podaci.vrt'
      using driver `VRT' successful.

Layer name: podaci
Geometry: Point
Feature Count: 4448
Layer SRS WKT:
(unknown)
X: String (0.0)
Y: String (0.0)
Z: String (0.0)

Primjerice znamo da skup podataka ima 4448 točaka i da nema definiran prostorni referentni sustav (SRS). SRS možemo definirati tako da dodamo: 

<LayerSRS>+proj=tmerc +ellps=bessel +k_o=0.9999 +x_0=5500000 +lon_0=15 +units=m</LayerSRS>

čime bi definirali 5 zonu GaussKruger projekcije bez transformacijskih parametara prema wgs84 u PROJ.4 formatu.

Gridding i odabir metode interploacije

Već je spomenuto kako gdal_grid podržava tri metode interpolacije (link). Osim toga za svaku se metodu mogu kontrolirati različiti parametri. Najkorisnije je spomenuti da bez definicije parametara gdal_grid u obzir uzima sve točke zadanog područja, što vrlo vjerojatno nije željno ponašanje.

Od parametara su najzanimljiviji:

Retrieved from "http://wiki.osgeo.org/w/index.php?title=Croatia/Kratki_vodici/digitalni-model-terena-iz-csv-datoteke-pomocu-gdala&oldid=51872"