Re: [FreeGeoCZ] PostGIS raster - export vice tiles do jednoho TIFFu

[Antonin Orlik <antonin.orlik AT gmail.com>]

Nakonec jsem na to přiąel - přehlédl jsem parametr "mode" pro PG driver GDALu - ten musí být v tomto případě nastaven na hodnotu 2 - viz: https://trac.osgeo.org/gdal/wiki/frmts_wtkraster.html

Správně tedy:
gdal_translate "PG:dbname='test' table=tile where='rid IN (18,19,20,21)' mode=2" /tmp/test.tif

Jinak co do rychlosti - před časem jsem prováděl benchmark - porovnával jsem různé operace mezi PostGIS Raster, GRASS, RASDAMAN a souborovým přístupem přes GDAL API. Hodně záleľí, co chceą s těmi rastry dělat - PostGIS Raster sice nedosahuje rychlosti jako souborový přístup přes GDAL API, ale v mnoha ohledech zjednoduąuje ľivot a zavádí v datech pořádek. V mém případě provádím nad velkým mnoľstvím časových řad časté prostorové dotazy, různé agregace, statistiky a exportuji kompozice - zde je PostGIS Raster dobrou volbou...
Jinak souborový přístup přes GDAL API + numpy + scipy je bezkonkurenčně nejrychlejąí, pokud provádím jednorázové zpracování dat.

I tak díky!

Zdraví
Tonda

On 15.11.2016 21:12, Jáchym Čepický wrote:

Ahoj,

nemám moc odpověď, jenom dám do placu co jsem slyąel od autorů: postgis raster sice funguje, ale je to příąerně neefektivní a pomalý a moc by se to pouľívat nemělo.

No ... nějaký vacuum analyze nenapoví? Moc toho o databázích nevím, jenom tak malinko, asi to máą uľ načtený..

J

Dne 15.11.2016 v 10:59 Antonin Orlik napsal(a):

Zdravim,

nemuzu najit efektivni zpusob, jak z PostGIS Raster DB ulozit do
jedineho souboru (treba GTiff) vice tiles (rastr byl do DB importovan s
parametrem -t).
Pokud exportuji jeden tile z cele mozaiky, pak funguje napr.:
gdal_translate "PG:dbname='test' table=tile where='rid=18'" /tmp/test.tif

ale pokud chci exportovat cely rastr nebo jen nekolik tiles z cele
mozaiky, ktera je slozena treba z 500 tiles, pak toto nefunguje:
gdal_translate "PG:host=localhost port=5432 dbname='test' user='user'
password='password' schema='public' table=tile where='rid IN
(18,19,20,21)'" /tmp/test.tif
(Input file contains subdatasets. Please, select one of them for reading.)

Zkousel jsem to i pres St_AsTIFF, kde jsem tiles spojil pomoci ST_Union,
ale pokud je tiles hodne, trva to zatracene dlouho.
Zkousel jsem i gdal_merge.py, ale neumi s tim pracovat.

Mate s tim nekdo zkusenosti? Ve skriptech od Martina L. jsem videl, ze
uvadi utilitu pgsql2raster, ale zadnou takovou jsem nenasel (jen
samozrejme raster2pgsql) - pravdepodobne omyl.

Predem diky!
Tonda




_______________________________________________
FreeGeoCZ mailing list
FreeGeoCZ AT fsv.cvut.cz
http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz

_______________________________________________
FreeGeoCZ mailing list
FreeGeoCZ AT fsv.cvut.cz
http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz