#gis #geography #geospatial

inverse_distance_weight

Inverse Distance Weighting (IDW) algorithm for spatial interpolation

2 releases

0.1.1 Sep 9, 2023
0.1.0 Sep 9, 2023

#219 in Geospatial

MIT license

17KB
268 lines

inverse_distance_weight

Latest version Documentation MIT

An implementation of the Inverse Distance Weighting (IDW) algorithm.

The crate supports points of 1 to 3 dimension to perform the interpolation.

The weighted function used in the algorithm is weightᵢ = 1 / distance(pointᵢ, position).

Examples

use inverse_distance_weight::IDW;

// 1 dimension
let points = vec![0.0, 1.0];
let values = vec![0.0, 1.0];
let idw = IDW::new(points, values);

let result = idw.evaluate(0.5);

// 2 dimension
let points = vec![(0.0, 0.0), (1.0, 1.0)];
let values = vec![0.0, 1.0];
let idw = IDW::new(points, values);

let result = idw.evaluate((0.5, 0.5));

// 3 dimension
let points = vec![(0.0, 0.0, 0.0), (1.0, 1.0, 1.0)];
let values = vec![0.0, 1.0];
let idw = IDW::new(points, values);

let result = idw.evaluate((0.5, 0.5, 0.5));

// Customize
let points = vec![0.0, 1.0];
let values = vec![0.0, 1.0];
let idw = IDW::new(points, values)
    // Sets a power parameter. Default is 2.
    .power(0.5)
    // Sets a transform function for weights.
    .weighted_function(|weight| (1.0 + (4.0 * std::f64::consts::PI * weight).sin()) * 0.5);

Dependencies

~150KB