Inverse Distance Weighting interpolation algorithm with elevation detrending/reprojection.
The closest n stations to each pixel are used to compute the local lapse rate, allowing to project the contributions of these n stations to the local pixel with an inverse distance weight. It therefore takes the following arguments:
- NEIGHBORS: how many neighbouring stations should be used;
- MAX_DISTANCE: maximum allowed distance (in meters) between the stations and grid points to interpolate to;
- SCALE: this is a scaling parameter to smooth the IDW distribution. In effect, this is added to the distance in order to move into the tail of the 1/d distribution (default: 1000m);
- ALPHA: this is an exponent to the 1/d distribution (default: 1);
Either NEIGHBORS or MAX_DISTANCE needs to be specified. When both are specified, both restrictions are used. In such cases, only stations less than MAX_DISTANCE away from the grid point to interpolate to are used, up to a maximum of n stations. This can lead to less than n stations to be included in the interpolation.
- Note
- Beware, this method sometimes produces very sharp transitions as it spatially moves from one station's area of influence to another one!
TA::algorithms = LIDW_LAPSE
TA::lidw_lapse::neighbors = 6
|
const std::string | algo |
|
std::vector< double > | getData (const Date &i_date, const std::string &i_param) |
|
size_t | getData (const Date &i_date, const std::string &i_param, std::vector< double > &o_vecData, std::vector< StationData > &o_vecMeta) |
|
TimeSeriesManager & | tsmanager |
|
Date | date |
|
std::vector< MeteoData > | vecMeteo |
|
std::vector< double > | vecData |
| store the measurement for the given parameter
|
|
std::vector< StationData > | vecMeta |
| store the station data for the given parameter
|
|
std::ostringstream | info |
| to store some extra information about the interplation process
|
|
const std::string | param |
| the parameter that we will interpolate
|
|
size_t | nrOfMeasurments |
| Number of stations that have been used, so this can be reported to the user.
|
|