petekio 0.2.6

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

//! Scattered-XYZ CSV reader — wraps the `csv` crate behind petekIO types.
//!
//! Reads a headered CSV, taking X/Y/Z from the named columns. Every *other*
//! column becomes a `f64` attribute column **iff** all its cells parse as
//! numbers; a column with any non-numeric cell is silently dropped (it is text
//! metadata, not data). Rows with a non-numeric X/Y/Z are an error (readers
//! validate on load). Imports only from `foundation`.

use crate::foundation::{GeoError, Result};
use indexmap::IndexMap;
use std::path::Path;

/// Parsed XYZ coordinates plus named `f64` attribute columns (aligned 1:1).
pub type LoadedPoints = (Vec<[f64; 3]>, IndexMap<String, Vec<f64>>);

/// Read XYZ + attribute columns. Returns `(coords, attrs)` where each attribute
/// column is aligned 1:1 with `coords`.
pub fn load(path: &Path, x: &str, y: &str, z: &str) -> Result<LoadedPoints> {
    let mut rdr =
        csv::Reader::from_path(path).map_err(|e| GeoError::Parse(format!("points CSV: {e}")))?;
    let headers: Vec<String> = rdr
        .headers()
        .map_err(|e| GeoError::Parse(format!("points CSV: bad header: {e}")))?
        .iter()
        .map(str::to_string)
        .collect();
    let find = |col: &str| {
        headers
            .iter()
            .position(|h| h == col)
            .ok_or_else(|| GeoError::NotFound(format!("points CSV: column '{col}'")))
    };
    let (xi, yi, zi) = (find(x)?, find(y)?, find(z)?);

    // Candidate attribute columns: every column that is not X/Y/Z.
    let attr_cols: Vec<usize> = (0..headers.len())
        .filter(|i| *i != xi && *i != yi && *i != zi)
        .collect();

    let mut coords: Vec<[f64; 3]> = Vec::new();
    // Parsed cells per attribute column; `None` once a non-numeric cell appears.
    let mut raw: Vec<Option<Vec<f64>>> = vec![Some(Vec::new()); attr_cols.len()];

    let parse = |rec: &csv::StringRecord, idx: usize, what: &str, row: usize| -> Result<f64> {
        let cell = rec
            .get(idx)
            .ok_or_else(|| GeoError::Parse(format!("points CSV: row {row}: missing {what}")))?;
        cell.trim().parse::<f64>().map_err(|e| {
            GeoError::Parse(format!("points CSV: row {row}: bad {what} '{cell}': {e}"))
        })
    };

    for (row, rec) in rdr.records().enumerate() {
        let rec = rec.map_err(|e| GeoError::Parse(format!("points CSV: row {row}: {e}")))?;
        let xv = parse(&rec, xi, x, row)?;
        let yv = parse(&rec, yi, y, row)?;
        let zv = parse(&rec, zi, z, row)?;
        coords.push([xv, yv, zv]);
        for (slot, &col) in raw.iter_mut().zip(&attr_cols) {
            let Some(buf) = slot else { continue };
            match rec.get(col).map(str::trim).map(str::parse::<f64>) {
                Some(Ok(v)) => buf.push(v),
                _ => *slot = None, // non-numeric (or missing) → drop this column
            }
        }
    }

    let mut attrs = IndexMap::new();
    for (slot, &col) in raw.into_iter().zip(&attr_cols) {
        if let Some(buf) = slot {
            attrs.insert(headers[col].clone(), buf);
        }
    }
    Ok((coords, attrs))
}