Re: [FreeGeoCZ] PostGIS raster - export vice tiles do jednoho TIFFu
Antonin Orlik <[email protected]> Tue, 15 Nov 2016 22:38:50 +0100
| Newsgroups | gmane.comp.gis.grass.cz |
|---|---|
| Message-ID | <[email protected]> |
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