def __runFastGeoreferencing__(self, georefObject, messtischblatt, tmpDir, destPath): try: gcps = createGCPs(georefObject.clipparameter, messtischblatt.BoundingBoxObj.getCornerPointsAsList(), messtischblatt.zoomify_height) return georeference(messtischblatt.original_path, destPath, tmpDir, gcps, DATABASE_SRID, DATABASE_SRID, 'polynom', self.logger) except: self.logger.error('Something went wrong while trying to process a fast georefercing process') raise
def __runStableGeoreferencing__(self, georefObject, messtischblatt, tmpDir, destPath): try: # create a shapefile which represents the boundingbox of the messtischblatt and is latery used for clipping shpPath = messtischblatt.BoundingBoxObj.asShapefile( os.path.join(tmpDir, "shape")) gcps = createGCPs( georefObject.clipparameter, messtischblatt.BoundingBoxObj.getCornerPointsAsList(), messtischblatt.zoomify_height) georeference_result_file = georeference( messtischblatt.original_path, destPath, tmpDir, gcps, DATABASE_SRID, DATABASE_SRID, 'polynom', self.logger, shpPath) # add overviews if self.__executeCommands__( "gdaladdo --config GDAL_CACHEMAX 500 -r average %s %s" % (georeference_result_file, GEOREFERENCE_OVERVIEW_LEVELS)): return destPath else: raise GeoreferenceProcessRunningError( 'Something went wrong while trying to process a georefercing process' ) except: self.logger.error( 'Something went wrong while trying to process a georefercing process' ) raise
def testCreateGCPs(self): print 'Test clip raster with shapfile function ...' clipparameters = '985.27322387696:8628.600769043,7056.9926300049:8619.2794799805,7043.318145752:2149.7222900391,969.47880554199:2167.4909973145' boundingbox = BoundingBox(20.9999980926514, 55.7999992370605, 21.1666679382324, 55.9000015258789, 4314) response = createGCPs(clipparameters, boundingbox.getCornerPointsAsList(), 9498) for gcp in response: print gcp self.assertTrue(isinstance(response, list), 'Response is not like expected.') print '=============='
def __runStableGeoreferencing__(self, georefObject, messtischblatt, tmpDir, destPath): try: # create a shapefile which represents the boundingbox of the messtischblatt and is latery used for clipping shpPath = messtischblatt.BoundingBoxObj.asShapefile(os.path.join(tmpDir,"shape")) gcps = createGCPs(georefObject.clipparameter, messtischblatt.BoundingBoxObj.getCornerPointsAsList(), messtischblatt.zoomify_height) georeference_result_file = georeference(messtischblatt.original_path, destPath, tmpDir, gcps, DATABASE_SRID, DATABASE_SRID, 'polynom', self.logger, shpPath) # add overviews if self.__executeCommands__("gdaladdo --config GDAL_CACHEMAX 500 -r average %s %s"%(georeference_result_file,GEOREFERENCE_OVERVIEW_LEVELS)): return destPath else: raise GeoreferenceProcessRunningError('Something went wrong while trying to process a georefercing process') except: self.logger.error('Something went wrong while trying to process a georefercing process') raise
def __runFastGeoreferencing__(self, georefObject, messtischblatt, tmpDir, destPath): try: gcps = createGCPs( georefObject.clipparameter, messtischblatt.BoundingBoxObj.getCornerPointsAsList(), messtischblatt.zoomify_height) return georeference(messtischblatt.original_path, destPath, tmpDir, gcps, DATABASE_SRID, DATABASE_SRID, 'polynom', self.logger) except: self.logger.error( 'Something went wrong while trying to process a fast georefercing process' ) raise