<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Boa tarde.<br>
<br>
Habitualmente utilizo dados vectoriais, pelo que não tenho muita
experiência com grandes volumes de dados em formato raster.<br>
<br>
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.<br>
<br>
Fiz alguma pesquisa e tentei 3 caminhos:<br>
<ol>
<li>"retile"</li>
<ol>
<li>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...
<pre>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 </pre>
</li>
</ol>
<li>Mosaico + pirâmide</li>
<ol>
<li>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....<br>
</li>
</ol>
<li>Pirâmide + mosaico</li>
<ol>
<li>Criei primeiro as pirâmides para todos os ficheiros com um
script (cerca de 1 hora a correr):
<pre>#!/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</pre>
</li>
<li>Criei depois o mosaico:
<pre>gdalbuildvrt mosaic.vrt *.jpg</pre>
</li>
<li>Carreguei para o QGis (deixei a abrir durante cerca de 30
min)</li>
</ol>
</ol>
<p><br>
Neste ultimo caso até consegui que abrisse, mas fica praticamente
inutilizavel porque demora demasiado tempo a arrastar ou fazer
zoom.<br>
</p>
<p>Existe mais alguma coisa que possa fazer ou outro formato que
possa utilizar para as imagens (envio informação gdalinfo abaixo)?
<br>
</p>
<p>Ou o melhor mesmo é usar menos imagens de cada vez?<br>
</p>
<p>Peço desculpa por alguma imprecisão ao nível dos termos usados.<br>
</p>
<p>Obrigado desde já.<br>
</p>
<p>Cumprimentos,<br>
Rui Henriques.<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>O resultado do gdalinfo (depois de criar a pirâmide) é o
seguinte:<br>
<font face="monospace">Driver: JPEG/JPEG JFIF<br>
Files: ortofoto_1.jpg<br>
ortofoto_1.jpg.ovr<br>
ortofoto_1.jgw<br>
Size is 10200, 10200<br>
Coordinate System is `'<br>
Origin = (XXXX.000000000000000,-XXXXX.000000000000000)<br>
Pixel Size = (0.100000000000000,-0.100000000000000)<br>
Image Structure Metadata:<br>
SOURCE_COLOR_SPACE=YCbCr<br>
INTERLEAVE=PIXEL<br>
COMPRESSION=JPEG<br>
Corner Coordinates:<br>
Upper Left ( XXXX.000, -XXXXX.000) <br>
Lower Left ( XXXX.000, -XXXXX.000) <br>
Upper Right ( XXXX.000, -XXXXX.000) <br>
Lower Right ( XXXX.000, -XXXXX.000) <br>
Center ( XXXX.000, -XXXXX.000) <br>
Band 1 Block=10200x1 Type=Byte, ColorInterp=Red<br>
Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
160x160, 80x80, 40x40, 20x20, 10x10<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 2 Block=10200x1 Type=Byte, ColorInterp=Green<br>
Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
160x160, 80x80, 40x40, 20x20, 10x10<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 3 Block=10200x1 Type=Byte, ColorInterp=Blue<br>
Overviews: 5100x5100, 2550x2550, 1275x1275, 638x638, 319x319,
160x160, 80x80, 40x40, 20x20, 10x10<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG</font><br>
</p>
</body>
</html>