...
The estimation formula used:
where:
Pest(t) is the estimated value at the test station at time t
Pi(t) is the measured value at neighbour station i at time t
Di is the distance between the test station and the neighbour station i
N is the number of neighbour stations taken into account
b is the power of distance D, (default b=2)
Test criteria
The test criterion with an absolute threshold is exceeded when the following condition fails:
The test criterion with relative threshold is exceeded when the following condition fails:
with:
b the admissable absolute difference
Xrel the multiplier of the standard deviation
SPest(t)
the standard deviation of neighbouring values
Configuration
- id: identifier of the check.
- variableDefinition: embedded variable definition (see above).
- inputVariableId: One or more identifiers for variables of which the flags have to be used.
- outputVariableId: One or more identifiers for variables for which the flags have to be modified.
- searchRadius: The maximum radius for selecting reference locations in meters.
...
A configuration example for the spatialHomogeneityCheck is given below:
Code Block |
---|
<?xml version="1.0" encoding="UTF-8"?>
<secondaryValidation xmlns="http://www.wldelft.nl/fews" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.wldelft.nl/fews http://fews.wldelft.nl/schemas/version1.0/secondaryValidation.xsd">
<variableDefinition>
<variableId>spatialHomogeneityCheck_location1</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location1</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location2</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location2</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location3</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location3</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location4</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location4</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location5</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location5</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location6</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location6</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location7</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location7</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location8</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location8</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location9</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location9</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location10</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location10</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<variableDefinition>
<variableId>spatialHomogeneityCheck_location11</variableId>
<timeSeriesSet>
<moduleInstanceId>SpatialHomogeneityCheck</moduleInstanceId>
<valueType>scalar</valueType>
<parameterId>H.obs</parameterId>
<locationId>location11</locationId>
<timeSeriesType>external historical</timeSeriesType>
<timeStep unit="minute" multiplier="15"/>
<readWriteMode>read complete forecast</readWriteMode>
</timeSeriesSet>
</variableDefinition>
<spatialHomogeneityCheck id="spatialHomogeneityCheck">
<inputVariableId>spatialHomogeneityCheck_location1</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location2</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location3</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location4</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location5</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location6</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location7</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location8</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location9</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location10</inputVariableId>
<inputVariableId>spatialHomogeneityCheck_location11</inputVariableId>
<outputVariableId>spatialHomogeneityCheck_location1</outputVariableId>
<outputVariableId>spatialHomogeneityCheck_location3</outputVariableId>
<outputVariableId>spatialHomogeneityCheck_location5</outputVariableId>
<outputVariableId>spatialHomogeneityCheck_location7</outputVariableId>
<outputVariableId>spatialHomogeneityCheck_location9</outputVariableId>
<outputVariableId>spatialHomogeneityCheck_location11</outputVariableId>
<searchRadius>100000</searchRadius>
<numberOfPoints>4</numberOfPoints>
<distancePower>2</distancePower>
<threshold>
<absolute>3</absolute>
<outputFlag>unreliable</outputFlag>
<logLevel>WARN</logLevel>
<logEventCode>SecondaryValidation.spatialHomogeneityCheck</logEventCode>
<logMessage>%AMOUNT_CHANGED_FLAGS% flags set to %OUTPUT_FLAG% by %CHECK_ID%, header=%HEADER%, location(s)=%LOCATION_NAME%</logMessage>
</threshold>
<threshold>
<relative>2</relative>
<outputFlag>doubtful</outputFlag>
<logLevel>INFO</logLevel>
<logEventCode>SecondaryValidation.spatialHomogeneityCheck</logEventCode>
<logMessage>%AMOUNT_CHANGED_FLAGS% flags set to %OUTPUT_FLAG% by %CHECK_ID%, header=%HEADER%, location(s)=%LOCATION_NAME%</logMessage>
</threshold>
</spatialHomogeneityCheck>
</secondaryValidation>
|