¿Puedo desenfocar una imagen con GDAL?

Question

Tengo una imagen GeoTiff mymap.tif https://www.dropbox.com/s/51chrzv9vp9olg8/mymap.tif

$ gdalinfo mymap.tif Driver: GTiff/GeoTIFF Files: mymap.tif Size is 1085, 281 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (-122.509475852140255,37.769141730303630) Pixel Size = (0.000010000000000,-0.000010000000000) Metadata: AREA_OR_POINT=Area Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=PIXEL Corner Coordinates: Upper Left (-122.5094759, 37.7691417) (122d30'34.11"W, 37d46' 8.91"N) Lower Left (-122.5094759, 37.7663317) (122d30'34.11"W, 37d45'58.79"N) Upper Right (-122.4986259, 37.7691417) (122d29'55.05"W, 37d46' 8.91"N) Lower Right (-122.4986259, 37.7663317) (122d29'55.05"W, 37d45'58.79"N) Center (-122.5040509, 37.7677367) (122d30'14.58"W, 37d46' 3.85"N) Band 1 Block=1085x32 Type=Byte, ColorInterp=Red Mask Flags: PER_DATASET ALPHA Band 2 Block=1085x32 Type=Byte, ColorInterp=Green Mask Flags: PER_DATASET ALPHA Band 3 Block=1085x32 Type=Byte, ColorInterp=Blue Mask Flags: PER_DATASET ALPHA Band 4 Block=1085x32 Type=Byte, ColorInterp=Alpha

Si obtengo los límites de esta imagen desde Google Map, obtengo esta imagen de la verdad sobre el terreno https://www.dropbox.com/s/mrb3d1ho33mw9bw/truth.tif . Las dos imágenes no están del todo alineadas.

¿Es posible desenfocar mi imagen para que al superponerla a la de Google Maps se alineen correctamente?

Answer 1

Si cargo su mymap.tif en QGIS con el plugin Openlayers, se alinea bien con la imagen aérea de Google y Bing. El fondo de Openstreetmap está un poco desplazado, pero no sé qué precisión tiene.

Tu truth.tif no tiene información de límites, por eso no se alinea. Ambos tienen el mismo tamaño de píxel, pero no se puede saber la resolución sin los límites en grados o cualquier otra unidad.

Después de georreferenciarla manualmente, ambas imágenes aparecen en el lugar correcto, pero con tamaños diferentes:

Convirtiendo truth.tif a WGS84, obtengo esta extensión reportada por gdalinfo:

Upper Left  (-122.5099500,  37.7696580) (122d30'35.82"W, 37d46'10.77"N)
Lower Left  (-122.5099500,  37.7652832) (122d30'35.82"W, 37d45'55.02"N)
Upper Right (-122.4975252,  37.7696580) (122d29'51.09"W, 37d46'10.77"N)
Lower Right (-122.4975252,  37.7652832) (122d29'51.09"W, 37d45'55.02"N)
Center      (-122.5037376,  37.7674706) (122d30'13.46"W, 37d46' 2.89"N)

No sé cómo has conseguido coger la foto de Google Map, pero evidentemente no se han conservado los límites. Una razón podría ser que no comparten la misma proyección.