surtgis 0.7.0

High-performance geospatial analysis CLI
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147

//! CLI handlers for vector geoprocessing commands.

use std::time::Instant;

use anyhow::{Context, Result};

use surtgis_algorithms::vector::overlay;
use surtgis_algorithms::vector::{buffer_points, BufferParams};
use surtgis_core::vector::{read_vector, FeatureCollection};

use crate::commands::VectorCommands;
use crate::helpers::spinner;

pub fn handle(action: VectorCommands) -> Result<()> {
    match action {
        VectorCommands::Intersection { input_a, input_b, output } => {
            let start = Instant::now();
            let a = read_vector(&input_a).context("Failed to read layer A")?;
            let b = read_vector(&input_b).context("Failed to read layer B")?;
            println!("A: {} features, B: {} features", a.features.len(), b.features.len());

            let pb = spinner("Computing intersection...");
            let result = overlay::intersection(&a, &b);
            pb.finish_and_clear();

            println!("Result: {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }

        VectorCommands::Union { input_a, input_b, output } => {
            let start = Instant::now();
            let a = read_vector(&input_a)?;
            let b = read_vector(&input_b)?;

            let pb = spinner("Computing union...");
            let result = overlay::union(&a, &b);
            pb.finish_and_clear();

            println!("Result: {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }

        VectorCommands::Difference { input_a, input_b, output } => {
            let start = Instant::now();
            let a = read_vector(&input_a)?;
            let b = read_vector(&input_b)?;

            let pb = spinner("Computing difference...");
            let result = overlay::difference(&a, &b);
            pb.finish_and_clear();

            println!("Result: {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }

        VectorCommands::SymDifference { input_a, input_b, output } => {
            let start = Instant::now();
            let a = read_vector(&input_a)?;
            let b = read_vector(&input_b)?;

            let pb = spinner("Computing symmetric difference...");
            let result = overlay::symmetric_difference(&a, &b);
            pb.finish_and_clear();

            println!("Result: {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }

        VectorCommands::Dissolve { input, output } => {
            let start = Instant::now();
            let fc = read_vector(&input)?;
            println!("Input: {} features", fc.features.len());

            let pb = spinner("Dissolving...");
            let result = overlay::dissolve(&fc);
            pb.finish_and_clear();

            println!("Dissolved to {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }

        VectorCommands::Buffer { input, distance, segments, output } => {
            let start = Instant::now();
            let fc = read_vector(&input)?;
            println!("Input: {} features, distance: {}", fc.features.len(), distance);

            let pb = spinner("Buffering...");
            let params = BufferParams {
                distance,
                segments,
            };

            let mut result = FeatureCollection::new();
            for feature in &fc.features {
                if let Some(geom) = &feature.geometry {
                    if let geo::Geometry::Point(pt) = geom {
                        let buffered = buffer_points(pt, &params);
                        result.push(surtgis_core::vector::Feature {
                            geometry: Some(geo::Geometry::Polygon(buffered)),
                            properties: feature.properties.clone(),
                            id: feature.id.clone(),
                        });
                    }
                }
            }
            pb.finish_and_clear();

            println!("Buffered {} features in {:.1?}", result.features.len(), start.elapsed());
            write_vector(&result, &output)?;
        }
    }
    Ok(())
}

fn write_vector(fc: &FeatureCollection, path: &std::path::Path) -> Result<()> {
    use geo::algorithm::coords_iter::CoordsIter;

    let features: Vec<serde_json::Value> = fc.features.iter().filter_map(|f| {
        let geom = f.geometry.as_ref()?;
        let coords: Vec<Vec<Vec<f64>>> = match geom {
            geo::Geometry::Polygon(p) => {
                let mut rings = vec![];
                rings.push(p.exterior().0.iter().map(|c| vec![c.x, c.y]).collect());
                for interior in p.interiors() {
                    rings.push(interior.0.iter().map(|c| vec![c.x, c.y]).collect());
                }
                rings
            }
            _ => return None,
        };
        Some(serde_json::json!({
            "type": "Feature",
            "geometry": { "type": "Polygon", "coordinates": coords },
            "properties": {}
        }))
    }).collect();

    let geojson = serde_json::json!({
        "type": "FeatureCollection",
        "features": features
    });

    std::fs::write(path, serde_json::to_string_pretty(&geojson)?)
        .context("Failed to write output")?;
    println!("Written: {}", path.display());
    Ok(())
}