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
[email protected]
http://lists.osgeo.org/mailman/listinfo/portugal
