RasterLite: Porovnání verzí
mBez shrnutí editace |
|||
Řádek 35: | Řádek 35: | ||
tile_id = 2 | tile_id = 2 | ||
astext(geometry) = POLYGON((-175.808333 90, -173.7125 90, -173.7125 87.904167, -175.808333 87.904167, -175.808333 90)) | astext(geometry) = POLYGON((-175.808333 90, -173.7125 90, -173.7125 87.904167, -175.808333 87.904167, -175.808333 90)) | ||
</pre> | |||
V databázi se vytvoří dvě tabulky | |||
* <tt>marble_rasters</tt> - binární data dlaždic | |||
* <tt>marble_metadata</tt> - metadata dlaždic | |||
<source lang=sql> | |||
select * from marble_rasters limit 3; | |||
</source> | |||
<pre> | |||
id|raster | |||
1 |II* | |||
2 |II* | |||
3 |II* | |||
</pre> | |||
<source lang=sql> | |||
select * from marble_metadata limit 3; | |||
</source> | |||
<pre> | |||
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 /opt/vyuka/land 0 503 503 0.0041666666 0.0041666666 | |||
2 /opt/vyuka/land 1 503 503 0.0041666666 0.0041666666 | |||
</pre> | </pre> | ||
Verze z 13. 4. 2011, 18:27
RasterLite je rozšíření SpatiaLite pro uložení, manipulaci a analýzu rastrových dat.
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://ueod-globe.net/globe/TrueMarble_GeoTIFF/TrueMarble.500m.21600x21600.A1.tif.gz gzip -d TrueMarble.500m.21600x21600.A1.tif.gz spatialite raster.sqlite < init_spatialite-2.3.sql
Data naimportujeme pomocí nástroje rasterlite_load.
rasterlite_load -d raster.sqlite -T marble -f TrueMarble.500m.21600x21600.A1.tif -i TIFF -e 4326
Příklad jednoduchého dotazu.
SELECT tile_id, AsText(geometry) FROM marble_metadata
WHERE NOT IsEmpty(geometry) LIMIT 3;
tile_id = 0 astext(geometry) = POLYGON((-180 90, -177.904167 90, -177.904167 87.904167, -180 87.904167, -180 90)) tile_id = 1 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))
V databázi se vytvoří dvě tabulky
- marble_rasters - binární data dlaždic
- marble_metadata - metadata dlaždic
select * from marble_rasters limit 3;
id|raster 1 |II* 2 |II* 3 |II*
select * from marble_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 /opt/vyuka/land 0 503 503 0.0041666666 0.0041666666 2 /opt/vyuka/land 1 503 503 0.0041666666 0.0041666666