RasterLite: Porovnání verzí

Z GeoWikiCZ
Řádek 20: Řádek 20:


<source lang="sql">
<source lang="sql">
SELECT tile_id, AsText(geometry) FROM dem
SELECT tile_id, AsText(geometry) FROM landsat_b10_metadata
  WHERE NOT IsEmpty(geometry) LIMIT 3;
  WHERE NOT IsEmpty(geometry) LIMIT 3;
</source>
</source>


<pre>
<pre>
        tile_id = 0
tile_id|AsText(geometry)
astext(geometry) = POLYGON((-180 90, -177.904167 90, -177.904167 87.904167, -180 87.904167, -180 90))
0|POLYGON((403080 5679180, 433140 5679180, 433140 5652120, 403080 5652120, 403080 5679180))
 
1|POLYGON((433140 5679180, 463200 5679180, 463200 5652120, 433140 5652120, 433140 5679180))
        tile_id = 1
2|POLYGON((463200 5679180, 493260 5679180, 493260 5652120, 463200 5652120, 463200 5679180))
astext(geometry) = POLYGON((-177.904167 90, -175.808333 90, -175.808333 87.904167, -177.904167 87.904167, -177.904167 90))
 
        tile_id = 2
astext(geometry) = POLYGON((-175.808333 90, -173.7125 90, -173.7125 87.904167, -175.808333 87.904167, -175.808333 90))
</pre>
</pre>


V databázi se vytvoří dvě tabulky
V databázi se vytvoří dvě tabulky


* <tt>dem_rasters</tt> - binární data dlaždic
* <tt>landsat_b10_rasters</tt> - binární data dlaždic
* <tt>dem_metadata</tt> - metadata dlaždic
* <tt>landsat_b10_metadata</tt> - metadata dlaždic


<source lang=sql>
<source lang=sql>
SELECT * FROM dem_rasters LIMIT 3;
SELECT * FROM landsat_b10_rasters LIMIT 3;
</source>
</source>


Řádek 52: Řádek 48:


<source lang=sql>
<source lang=sql>
SELECT * FROM marble_metadata LIMIT 3;
SELECT * FROM landsat_b10_metadata LIMIT 3;
</source>
</source>


Řádek 59: Řádek 55:
----------  ---------------  ----------  ----------  ----------  ------------  ------------  ----------
----------  ---------------  ----------  ----------  ----------  ------------  ------------  ----------
0          raster metadata  0          0          0          0.0          0.0                     
0          raster metadata  0          0          0          0.0          0.0                     
1          /opt/vyuka/land 0          503         503         0.0041666666  0.0041666666           
1          L4191025_025198 0          501         451         60.0         60.0                  
2          /opt/vyuka/land 1          503         503         0.0041666666  0.0041666666         
2          L4191025_025198 1          501         451         60.0         60.0                  
</pre>
</pre>



Verze z 16. 4. 2012, 12:24

RasterLite je rozšíření SpatiaLite pro uložení, manipulaci a analýzu rastrových dat v prostředí SQLite.

Pod operačním systémem Debian GNU/Linux nainstalujeme SpatiaLite příkazem

sudo apt-get install librasterlite2 rasterlite-bin

Import dat

Nejprve stáhneme testovací data a vytvoříme SpatiaLite databázi.

wget http://geo102.fsv.cvut.cz/zodh/2011/tiff/pardubice-1980/L4191025_02519840414_B10.TIF

Data naimportujeme pomocí nástroje rasterlite_load.

rasterlite_load -d gisdb.sqlite -T landsat_b10 -f L4191025_02519840414_B10.TIF -i TIFF

Příklad jednoduchého dotazu.

SELECT tile_id, AsText(geometry) FROM landsat_b10_metadata
 WHERE NOT IsEmpty(geometry) LIMIT 3;
tile_id|AsText(geometry)
0|POLYGON((403080 5679180, 433140 5679180, 433140 5652120, 403080 5652120, 403080 5679180))
1|POLYGON((433140 5679180, 463200 5679180, 463200 5652120, 433140 5652120, 433140 5679180))
2|POLYGON((463200 5679180, 493260 5679180, 493260 5652120, 463200 5652120, 463200 5679180))

V databázi se vytvoří dvě tabulky

  • landsat_b10_rasters - binární data dlaždic
  • landsat_b10_metadata - metadata dlaždic
SELECT * FROM landsat_b10_rasters LIMIT 3;
id|raster
1 |II*
2 |II*
3 |II*
SELECT * FROM landsat_b10_metadata LIMIT 3;
id          source_name      tile_id     width       height      pixel_x_size  pixel_y_size  geometry  
----------  ---------------  ----------  ----------  ----------  ------------  ------------  ----------
0           raster metadata  0           0           0           0.0           0.0                     
1           L4191025_025198  0           501         451         60.0          60.0                    
2           L4191025_025198  1           501         451         60.0          60.0                    

Související články

Externí odkazy