def test_clipRaster(self): """Raster layers can be clipped """ # Create a raster layer myName = 'shake' myRasterLayer = QgsRasterLayer(RASTERPATH, myName) myMessage = 'Did not find layer "%s" in path "%s"' % \ (myName, RASTERPATH) assert myRasterLayer is not None, myMessage # Create a bounding box myRect = [97, -3, 104, 1] # Clip the vector to the bbox myResult = clipLayer(myRasterLayer, myRect) # Check the output is valid assert os.path.exists(myResult) # Clip and give a desired resolution for the output mySize = 0.05 myResult = clipLayer(myRasterLayer, myRect, mySize) myNewRasterLayer = QgsRasterLayer(myResult, myName) assert myNewRasterLayer.isValid(), 'Resampled raster is not valid' myMessage = ('Resampled raster has incorrect pixel size.' 'Expected: %f, Actual: %f' % (mySize, myNewRasterLayer.rasterUnitsPerPixel())) assert myNewRasterLayer.rasterUnitsPerPixel() == mySize, myMessage
def test_clipRaster(self): """Raster layers can be clipped """ # Create a raster layer myName = 'shake' myRasterLayer = QgsRasterLayer(RASTERPATH, myName) msg = 'Did not find layer "%s" in path "%s"' % (myName, RASTERPATH) assert myRasterLayer is not None, msg # Create a bounding box myRect = [97, -3, 104, 1] # Clip the vector to the bbox myResult = clipLayer(myRasterLayer, myRect) # Check the output is valid assert os.path.exists(myResult) # Clip and give a desired resolution for the output mySize = 0.05 myResult = clipLayer(myRasterLayer, myRect, mySize) myNewRasterLayer = QgsRasterLayer(myResult, myName) assert myNewRasterLayer.isValid(), 'Resampled raster is not valid' msg = ('Resampled raster has incorrect pixel size.' 'Expected: %f, Actual: %f' % (mySize, myNewRasterLayer.rasterUnitsPerPixel())) assert myNewRasterLayer.rasterUnitsPerPixel() == mySize, msg