petekio 0.3.15

Subsurface data ingestion + structure layer: surfaces, wells, points, polygons with loading, interpolation, and statistics.
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
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

//! Human-readable export writers for scattered points and polygon rings —
//! GeoJSON + CSV. The write side of `io::vector` / `io::csv_points`, so an
//! element round-trips through its interchange format. `NaN` → `null` (GeoJSON)
//! or the literal `NaN` (CSV), both of which parse back to `f64::NAN`.

use crate::foundation::{GeoError, Result};
use indexmap::IndexMap;
use serde_json::{json, Map, Value};
use std::io::Write;
use std::path::Path;

fn num(v: f64) -> Value {
    if v.is_finite() {
        json!(v)
    } else {
        Value::Null // NaN/Inf aren't representable in JSON
    }
}

/// Write scattered points (`[x, y, z]` + attribute columns) as a GeoJSON
/// `FeatureCollection` of `Point`s (attributes become feature properties).
pub fn write_points_geojson(
    path: &Path,
    coords: &[[f64; 3]],
    attrs: &IndexMap<String, Vec<f64>>,
) -> Result<()> {
    let features: Vec<Value> = coords
        .iter()
        .enumerate()
        .map(|(i, c)| {
            let mut props = Map::new();
            for (k, col) in attrs {
                props.insert(k.clone(), num(col.get(i).copied().unwrap_or(f64::NAN)));
            }
            json!({
                "type": "Feature",
                "geometry": {"type": "Point", "coordinates": [c[0], c[1], c[2]]},
                "properties": Value::Object(props),
            })
        })
        .collect();
    let fc = json!({"type": "FeatureCollection", "features": features});
    write_json(path, &fc)
}

/// Write polygon `rings` (each `[x, y, z]`) as a GeoJSON `FeatureCollection` of
/// `Polygon`s (2-D exteriors; z is dropped, matching the areal model).
pub fn write_polygons_geojson(path: &Path, rings: &[Vec<[f64; 3]>]) -> Result<()> {
    let features: Vec<Value> = rings
        .iter()
        .map(|ring| {
            let ext: Vec<Value> = ring.iter().map(|p| json!([p[0], p[1]])).collect();
            json!({
                "type": "Feature",
                "geometry": {"type": "Polygon", "coordinates": [ext]},
                "properties": {},
            })
        })
        .collect();
    let fc = json!({"type": "FeatureCollection", "features": features});
    write_json(path, &fc)
}

/// Write scattered points as CSV with `x,y,z` + one column per attribute.
pub fn write_points_csv(
    path: &Path,
    coords: &[[f64; 3]],
    attrs: &IndexMap<String, Vec<f64>>,
) -> Result<()> {
    let mut out = String::from("x,y,z");
    for k in attrs.keys() {
        out.push(',');
        out.push_str(k);
    }
    out.push('\n');
    for (i, c) in coords.iter().enumerate() {
        out.push_str(&format!("{},{},{}", c[0], c[1], c[2]));
        for col in attrs.values() {
            out.push(',');
            out.push_str(&col.get(i).copied().unwrap_or(f64::NAN).to_string());
        }
        out.push('\n');
    }
    let mut f = std::fs::File::create(path)?;
    f.write_all(out.as_bytes())?;
    Ok(())
}

fn write_json(path: &Path, value: &Value) -> Result<()> {
    let s = serde_json::to_string(value).map_err(|e| GeoError::Parse(e.to_string()))?;
    std::fs::File::create(path)?.write_all(s.as_bytes())?;
    Ok(())
}