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
