To mě zajímá,

co ten RASDAMAN? Nějaké praktické postřehy?

Dík

J

Dne 15.11.2016 v 22:38 Antonin Orlik napsal(a):
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@fsv.cvut.cz
>> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz
>>
>
>
>
> _______________________________________________
> FreeGeoCZ mailing list
> FreeGeoCZ@fsv.cvut.cz
> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz




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


--
Jachym Cepicky
e-mail: jachym.cepi...@gmail.com
twitter: @jachymc

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

Reply via email to