[Portugal] Carregamento de dados raster

[Giovanni Manghi]

Olá Rui,


eu normalmente uso as indicações descritas aqui

http://mapserver.org/input/raster.html#preprocessing-rasters

e sempre tive bons resultados: mapas de varios GB rapidos a abrir no
QGIS e a serem visualizados em aplicações webgis.

Abraços

-- Giovanni --


On Fri, 2012-03-09 at 16:34 +0000, Rui Pedro Henriques wrote:
> Boa tarde.
> 
> Habitualmente utilizo dados vectoriais, pelo que não tenho muita
> experiência com grandes volumes de dados em formato raster.
> 
> Recentemente tive necessidade de utilizar uma quantidade grande de
> ortofotomapas (283 imagens, 5,3 GB) que, naturalmente, tive grande
> dificuldade em carregar para onde fosse.
> 
> Fiz alguma pesquisa e tentei 3 caminhos:
>      1. "retile"
>              1. Utilizei o seguinte comando mas desisti quando
>                 verifiquei (ao fim de 2 horas) que já tinha 20 GB de
>                 dados novos e ainda só ia no segundo ortofoto... 
>                 ls *.jpg > files.txt
>                 mkdir tiles
>                 gdal_retile.py -levels 3 -ps 2048 2048 -co "TILED=YES" -co "BLOCKXSIZE=256" -co "BLOCKYSIZE=256" -s_srs EPSG:27493 -tileIndex tiles_index -targetDir tiles --optfile files.txt 
>      2. Mosaico + pirâmide
>              1. A tentativa era criar primeiro um mosaico com todos os
>                 ortofotos e depois criar as pirâmides. Desisti depois
>                 de 2 horas e ainda a 70% com um ficheiro de "overview"
>                 de 23GB....
>      3. Pirâmide + mosaico
>              1. Criei primeiro as pirâmides para todos os ficheiros
>                 com um script (cerca de 1 hora a correr): 
>                 #!/bin/bash
>                 
>                 for FILE in *.jpg
>                 do
>                   gdaladdo -ro --config INTERLEAVE_OVERVIEW PIXEL --config COMPRESS_OVERVIEW JPEG $FILE 2 4 8 16 32 64 128 256 512 1024
>                 done
>              2. Criei depois o mosaico: 
>                 gdalbuildvrt mosaic.vrt *.jpg
>              3. Carreguei para o QGis (deixei a abrir durante cerca de
>                 30 min)
> 
> 
> Neste ultimo caso até consegui que abrisse, mas fica praticamente
> inutilizavel porque demora demasiado tempo a arrastar ou fazer zoom.
> 
> 
> Existe mais alguma coisa que possa fazer ou outro formato que possa
> utilizar para as imagens (envio informação gdalinfo abaixo)? 
> 
> 
> Ou o melhor mesmo é usar menos imagens de cada vez?
> 
> 
> Peço desculpa por alguma imprecisão ao nível dos termos usados.
> 
> 
> Obrigado desde já.
> 
> 
> Cumprimentos,
> Rui Henriques.
> 
> 
> 
> 
> O resultado do gdalinfo (depois de criar a pirâmide) é o seguinte:
> Driver: JPEG/JPEG JFIF
> Files: ortofoto_1.jpg
>        ortofoto_1.jpg.ovr
>        ortofoto_1.jgw
> Size is 10200, 10200
> Coordinate System is `'
> Origin = (XXXX.000000000000000,-XXXXX.000000000000000)
> Pixel Size = (0.100000000000000,-0.100000000000000)
> Image Structure Metadata:
>   SOURCE_COLOR_SPACE=YCbCr
>   INTERLEAVE=PIXEL
>   COMPRESSION=JPEG
> Corner Coordinates:
> Upper Left  (    XXXX.000,  -XXXXX.000) 
> Lower Left  (    XXXX.000,  -XXXXX.000) 
> Upper Right (    XXXX.000,  -XXXXX.000) 
> Lower Right (    XXXX.000,  -XXXXX.000) 
> Center      (    XXXX.000,  -XXXXX.000) 
> Band 1 Block=10200x1 Type=Byte, ColorInterp=Red
>   Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
> 160x160, 80x80, 40x40, 20x20, 10x10
>   Image Structure Metadata:
>     COMPRESSION=JPEG
> Band 2 Block=10200x1 Type=Byte, ColorInterp=Green
>   Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
> 160x160, 80x80, 40x40, 20x20, 10x10
>   Image Structure Metadata:
>     COMPRESSION=JPEG
> Band 3 Block=10200x1 Type=Byte, ColorInterp=Blue
>   Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
> 160x160, 80x80, 40x40, 20x20, 10x10
>   Image Structure Metadata:
>     COMPRESSION=JPEG
> 
> 
> _______________________________________________
> Portugal mailing list
> Portugal at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/portugal