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