mappers_warp 0.1.0

Very simplistic tool for reprojecting maps, based on GdalWarp
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57

#![cfg(feature = "multithreading")]

use anyhow::Result;
use mappers::{
    Ellipsoid,
    projections::{LambertConformalConic, LongitudeLatitude},
};
use mappers_warp::{CubicBSpline, RasterBoundsDefinition, Warper};
use ndarray::Array2;

mod utils;
use utils::*;

#[test]
fn parallel_equivalence() -> Result<()> {
    let src_proj = LongitudeLatitude;
    let tgt_proj = LambertConformalConic::builder()
        .ref_lonlat(80., 24.)
        .standard_parallels(12.472955, 35.1728044444444)
        .ellipsoid(Ellipsoid::WGS84)
        .initialize_projection()?;

    let source_bounds =
        RasterBoundsDefinition::new((60.00, 68.25), (31.75, 40.0), 0.25, 0.25, src_proj)?;
    let target_bounds = RasterBoundsDefinition::new(
        (2_320_000. - 4_000_000., 2_740_000. - 4_000_000.),
        (5_090_000. - 4_000_000., 5_640_000. - 4_000_000.),
        10_000.,
        10_000.,
        tgt_proj,
    )?;

    let source_raster: Array2<f64> = open_nc_data("./tests/data/waves_34.nc")?;

    let serial_warper = Warper::initialize::<CubicBSpline, _, _>(&source_bounds, &target_bounds)?;
    let prllel_warper =
        Warper::initialize_parallel::<CubicBSpline, _, _>(&source_bounds, &target_bounds)?;
    assert_eq!(serial_warper, prllel_warper);

    let serial_raster = serial_warper.warp_unchecked(&source_raster);
    let prllel_raster = prllel_warper.warp_unchecked_parallel(&source_raster);
    assert_eq!(serial_raster, prllel_raster); // we don't want any nans in result so assert_eq is fine

    let serial_raster = serial_warper.warp_reject_nodata(&source_raster)?;
    let prllel_raster = prllel_warper.warp_reject_nodata_parallel(&source_raster)?;
    assert_eq!(serial_raster, prllel_raster); // we don't want any nans in result so assert_eq is fine

    let serial_raster = serial_warper.warp_ignore_nodata(&source_raster)?;
    let prllel_raster = prllel_warper.warp_ignore_nodata_parallel(&source_raster)?;
    assert_eq!(serial_raster, prllel_raster); // we don't want any nans in result so assert_eq is fine

    let serial_raster = serial_warper.warp_discard_nodata(&source_raster)?;
    let prllel_raster = prllel_warper.warp_discard_nodata_parallel(&source_raster)?;
    assert_eq!(serial_raster, prllel_raster); // we don't want any nans in result so assert_eq is fine

    Ok(())
}