[Portugal] Carregamento de dados raster

[Rui Pedro Henriques]

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

-------------- próxima parte ----------
Um anexo em HTML foi limpo...
URL: http://lists.osgeo.org/pipermail/portugal/attachments/20120309/576198fa/attachment.html