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

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
> >> [email protected]
> >> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz
> >>
> >
> >
> >
> > _______________________________________________
> > FreeGeoCZ mailing list
> > [email protected]
> > http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz
>
>
>
>
> _______________________________________________
> FreeGeoCZ mailing list
> [email protected]
> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz
>

-- 
Jachym Cepicky
e-mail: [email protected]
twitter: @jachymc

_______________________________________________
FreeGeoCZ mailing list
[email protected]
http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz