def arbitraryMapEquals(self, firstMap, secondMap): import math def floatEquals(float1, float2): threshold = 1e-6 return math.fabs(float1 - float2) <= threshold if((firstMap.dataType() == _pcraster.Scalar) | \ (firstMap.dataType() == _pcraster.Directional)): result = _pcraster._closeAtTolerance(firstMap, secondMap) else: result = pcraster.pcreq(firstMap, secondMap) minMap = pcraster.mapminimum(pcraster.spatial(pcraster.scalar(result))) value, isValid = pcraster.cellvalue(minMap, 1) assert isValid return floatEquals(value, 1.0)
def mapEqualsValidated(self, result, filename): import math def floatEquals(float1, float2): threshold = 1e-6 return math.fabs(float1 - float2) <= threshold validatedMap = self.readValidatedField(filename) if (result.dataType() == _pcraster.Scalar) | (result.dataType() == _pcraster.Directional): result = _pcraster._closeAtTolerance(result, validatedMap) else: result = pcraster.pcreq(result, validatedMap) min = pcraster.mapminimum(pcraster.scalar(result)) value, isValid = pcraster.cellvalue(min, 1) assert isValid return floatEquals(value, 1.0)
def mapEqualsValidated(self, result, filename): import math def floatEquals(float1, float2): threshold = 1e-6 return math.fabs(float1 - float2) <= threshold validatedMap = self.readValidatedField(filename) if((result.dataType() == _pcraster.Scalar) | \ (result.dataType() == _pcraster.Directional)): result = _pcraster._closeAtTolerance(result, validatedMap) else: result = pcraster.pcreq(result, validatedMap) min = pcraster.mapminimum(pcraster.scalar(result)) value, isValid = pcraster.cellvalue(min, 1) assert isValid return floatEquals(value, 1.0)