Cómo generar archivos de mundo para la trama de X?

Question

Dada una lista de referencias geográficas rásteres, alrededor de un centenar, lo que es el método más sencillo para generar archivos de mundo de cada uno de ellos?

En este caso, el origen de los rásteres son geotiff con la inserción de la georreferenciación o .BIL con el acompañamiento .PRJ, pero la solución ideal no importa lo que el formato es. Tengo la esperanza de evitar la duplicación de datos.

Tengo a mi disposición Arcgis 10 (Arcinfo Y Spatial Analyst), GDAL Utilidades (gdal v1.8.0), y de Quantum GIS, 1.8.

Answer 1

Si estás usando una versión reciente de QGIS, puede que ya tenga acceso al Extracto herramienta de Proyección mencionado aquí.

Del artículo: "Lo que esto le permite hacer es averiguar lo que la proyección de su rásteres. Luego de las exportaciones de este a un .WLD archivo. También le permite crear .Archivos PRJ. Si usted tiene GDAL instaladas las Herramientas (y activa), que encontrarás en esta nueva herramienta en Raster > Proyecciones > Extracto de la Proyección."

Incluye un procesamiento por lotes opción, lo que es ideal para sus propósitos.

Answer 2

Se desencadenó por Mike Toews puntero a GetGeoTransform, me las arreglé para crear un pequeño gdal secuencia de comandos de python que construye el mundo de los archivos para cualquier ráster georreferenciado (creo). El código completo está aquí: gdal-makeworld.py. Lo esencial bits son:

geotransform = dataset.GetGeoTransform()

if geotransform is not None:
    x, x_size, x_rot, y, y_rot, y_size = geotransform
    world_file.write('%s\n' % x_size)
    world_file.write('%s\n' % x_rot)
    world_file.write('%s\n' % y_rot)
    world_file.write('%s\n' % y_size)
    world_file.write('%s\n' % x)
    world_file.write('%s\n' % y)
    world_file.close()

''' geotransform tuple key:
    [0] /* top left x */
    [1] /* w-e pixel resolution */
    [2] /* rotation, 0 if image is "north up" */
    [3] /* top left y */
    [4] /* rotation, 0 if image is "north up" */
    [5] /* n-s pixel resolution */
'''

Adicional gracias a Schuyler Erle para escribir gdalcopyproj.py que he utilizado como punto de partida.