a `§»eMã@sXddlmZddlmZddlmZddlmZddlmZdd„Zdd „Z d d „Z d S) é)Úosr)Úgdal)Ú GA_Update)Úio)ÚlogcCsdS)Nz>+proj=utm +zone=30 +ellps=WGS84 +datum=WGS84 +units=m +no_defs©rrrú/code/opendm/pseudogeo.pyÚget_pseudogeo_utmsr cCsdS)Ngš™™™™™¹?rrrrrÚget_pseudogeo_scale sr c CsÈt |¡st d|¡dSzlt d|¡t |t¡}t  ¡}|  tj ¡|  t ƒ¡| | ¡¡| dtƒdddtƒ g¡d}Wn:tyÂ}z"t d|t|ƒf¡WYd}~n d}~00dS)Nz3Cannot add pseudo georeferencing, %s does not existzIAdding pseudo georeferencing (raster should show up at the equator) to %sgz8Cannot add pseudo georeferencing to %s (%s), skipping...)rÚ file_existsrÚ ODM_WARNINGÚODM_INFOrÚOpenrrÚSpatialReferenceÚSetAxisMappingStrategyÚOAMS_TRADITIONAL_GIS_ORDERÚImportFromProj4r Ú SetProjectionÚ ExportToWktÚSetGeoTransformr Ú ExceptionÚstr)ÚgeotiffÚdst_dsÚsrsÚerrrÚadd_pseudo_georeferencing s    rN) ÚosgeorrÚosgeo.gdalconstrÚopendmrrr r rrrrrÚs