wbraster 0.1.4

A pure-Rust library for reading and writing raster GIS formats
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
58
59
60
61
62
63

//! Raster reprojection example.
//!
//! Usage: cargo run --example reproject_io

use std::path::PathBuf;

use wbraster::{
    CrsInfo, DataType, DestinationFootprint, GridSizePolicy, NodataPolicy, Raster, RasterConfig,
    ReprojectOptions, ResampleMethod, Result,
};

fn data_dir() -> PathBuf {
    let dir = PathBuf::from(env!("CARGO_MANIFEST_DIR"))
        .join("data")
        .join("examples_reproject");
    std::fs::create_dir_all(&dir).unwrap();
    dir
}

fn main() -> Result<()> {
    let mut src = Raster::new(RasterConfig {
        cols: 120,
        rows: 80,
        bands: 1,
        x_min: -80.0,
        y_min: 35.0,
        cell_size: 0.01,
        nodata: -9999.0,
        data_type: DataType::F32,
        crs: CrsInfo::from_epsg(4326),
        ..Default::default()
    });

    for row in 0..src.rows {
        for col in 0..src.cols {
            let lon = src.col_center_x(col as isize);
            let lat = src.row_center_y(row as isize);
            let v = (lat.to_radians().sin() * lon.to_radians().cos()) * 1000.0;
            src.set(0, row as isize, col as isize, v)?;
        }
    }

    let opts = ReprojectOptions::new(3857, ResampleMethod::Bilinear)
        .with_grid_size_policy(GridSizePolicy::Expand)
        .with_destination_footprint(DestinationFootprint::SourceBoundary)
        .with_nodata_policy(NodataPolicy::PartialKernel)
        .with_resolution(1000.0, 1000.0);

    let projected = src.reproject_with_options(&opts)?;

    let out = data_dir().join("projected_3857.tif");
    projected.write(out.to_str().unwrap(), wbraster::RasterFormat::GeoTiff)?;

    println!("Source: {} x {}, EPSG {:?}", src.cols, src.rows, src.crs.epsg);
    println!(
        "Projected: {} x {}, EPSG {:?}",
        projected.cols, projected.rows, projected.crs.epsg
    );
    println!("Wrote {}", out.display());
    println!("reproject_io example OK");

    Ok(())
}