Gigantesca immagine JP2 da 2,5 GB: come affiancarla?


11

Ho una singola immagine JP2 da 2,5 GB di una città che devo dividere in tessere, che poi elaborerò per creare scenari di simulazione di volo.

Google come potrei, non riesco a trovare la soluzione.

A quanto pare JP2 ha i dati XML GIS incorporati, il che non mi riguarda in alcun modo, poiché non ne ho bisogno. Tutto quello che sto cercando di fare è suddividere l'immagine in numero x di tessere preferibilmente 2048 x 2048.

Un'utilità per farlo in OSX o Win7 sarebbe ottima.

Idee?

Grazie!


Sei interessato ad un approccio QGIS o ArcGIS?
Aaron

5
prova gdal2tiles.py
mdsumner

Risposte:


14

È possibile scrivere un semplice script di shell e utilizzare l'utilità gdal_translate per tagliare i riquadri. Su Mac, ottieni la build GDAL dalla sua e. Su Windows, il modo migliore per ottenere GDAL è tramite il programma di installazione OSGeo4W . Nel seguente script, è necessario sostituire le seguenti cose

  • XDIM e YDIM con le dimensioni X e Y jp2 di origine
  • gdal_translate con il percorso effettivo sul tuo sistema
  • input.jp2 con percorso al file di input.

Salvare lo script in un file di testo chiamato script.sh. Apri un Terminale ed esegui "bash script.sh". Ciò creerà i riquadri in una directory denominata riquadri /. Lo script viene copiato da qui con una piccola modifica per il tuo scopo.

#!/bin/bash
mkdir tiles

XDIM=24000
YDIM=24000
BLOCKSIZE=2048
XPOS=0
YPOS=0
BLOCKNO=0
while [ $YPOS -le $YDIM ]
do
while [ $XPOS -le $XDIM ]
  do
    echo "$XPOS $YPOS : ${BLOCKNO}.tif"
    gdal_translate -of GTiff -srcwin $XPOS $YPOS $BLOCKSIZE $BLOCKSIZE input.jp2 \
      tiles/${BLOCKNO}.tif
    BLOCKNO=`echo "$BLOCKNO + 1" | bc`
    XPOS=`echo "$XPOS + $BLOCKSIZE" | bc`
  done
YPOS=`echo "$YPOS + $BLOCKSIZE" | bc`
XPOS=0
done

2
Ah, nessun problema. Il Finder non è a conoscenza del formato JP2, quindi non aiuta. Hai bisogno della libreria GDAL per leggere JP2. Se hai installato GDAL dal link nella mia risposta, puoi usare il comando 'gdalinfo' per ottenere quelle informazioni. Sul mio sistema, GDAL è su /Library/Frameworks/GDAL.framework/Versions/1.9/Programs/gdalinfo. esegui gdalinfo <percorso / to / tuo / jp2> e riporterà le dimensioni. Allo stesso modo aggiungi il percorso completo a gdal_translate nello script affinché funzioni.
pensieri spaziali
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.