lasrs 0.1.3

A library for parsing geophysical well log (.las) files
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
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

use lasrs::Las;

#[test]
fn version_test() {
    let las = Las::new("./sample/example.las");
    assert_eq!(las.version(), 2.0);
}
#[test]
fn wrap_test() {
    let las = Las::new("./sample/example.las");
    assert_eq!(las.wrap(), false);
}
#[test]
fn headers_test() {
    let las = Las::new("./sample/example.las");
    assert_eq!(
        vec!["DEPT", "DT", "RHOB", "NPHI", "SFLU", "SFLA", "ILM", "ILD"],
        las.headers()
    );
    let las = Las::new("./sample/A10.las");
    assert_eq!(
        vec![
            "DEPT",
            "Perm",
            "Gamma",
            "Porosity",
            "Fluvialfacies",
            "NetGross"
        ],
        las.headers()
    );
}

#[test]
fn data_test() {
    let las = Las::new("./sample/example.las");
    let expected: Vec<Vec<f64>> = vec![
        vec![1670.0, 123.45, 2550.0, 0.45, 123.45, 123.45, 110.2, 105.6],
        vec![1669.875, 123.45, 2550.0, 0.45, 123.45, 123.45, 110.2, 105.6],
        vec![1669.75, 123.45, 2550.0, 0.45, 123.45, 123.45, 110.2, 105.6],
        vec![
            1669.745, 123.45, 2550.0, -999.25, 123.45, 123.45, 110.2, 105.6,
        ],
    ];
    assert_eq!(expected, las.data());
    let las = Las::new("./sample/A10.las");
    let expected: Vec<Vec<f64>> = vec![
        vec![1499.879, -999.25, -999.25, -999.25, -999.25, 0.0],
        vec![1500.129, -999.25, -999.25, -999.25, -999.25, 0.0],
        vec![1500.629, -999.25, -999.25, -999.25, -999.25, 0.0],
        vec![1501.129, -999.25, -999.25, 0.270646, 0.0, 0.0],
        vec![1501.629, 124.5799, 78.869453, 0.267428, 0.0, 0.0],
    ];
    assert_eq!(expected, &las.data()[0..5]);
}

#[test]
fn test_column() {
    let las = Las::new("./sample/example.las");
    assert_eq!(
        vec![1670.0, 1669.875, 1669.75, 1669.745],
        las.column("DEPT")
    );
}

#[test]
fn test_counts() {
    let las = Las::new("./sample/example.las");
    assert_eq!(las.column_count(), 8);
    assert_eq!(las.row_count(), 4);
}
#[test]
fn test_header_and_desc() {
    let las = Las::new("./sample/example.las");
    let mut expected = vec![
        ("DEPT", "DEPTH"),
        ("DT", "SONIC TRANSIT TIME"),
        ("RHOB", "BULK DENSITY"),
        ("NPHI", "NEUTRON POROSITY"),
        ("SFLU", "SHALLOW RESISTIVITY"),
        ("SFLA", "SHALLOW RESISTIVITY"),
        ("ILM", "MEDIUM RESISTIVITY"),
        ("ILD", "DEEP RESISTIVITY"),
    ]
    .into_iter()
    .map(|a| (a.0.to_string(), a.1.to_string()))
    .collect::<Vec<_>>();
    assert_eq!(
        expected.sort_by(|a, b| a.0.partial_cmp(&b.0).unwrap()),
        las.headers_and_desc()
            .sort_by(|a, b| a.0.partial_cmp(&b.0).unwrap())
    );
}

#[test]
fn other_test() {
    let las = Las::new("./sample/example.las");
    let expected = [
        "Note: The logging tools became stuck at 625 metres causing the data",
        "between 625 metres and 615 metres to be invalid.",
    ];
    assert_eq!(las.other(), expected.join("\n").to_string());

    let las = Las::new("./sample/A10.las");
    let expected = "".to_string();
    assert_eq!(las.other(), expected);
}

#[test]
fn csv_test() {
    use std::fs;
    let las = Las::new("./sample/example.las");
    las.to_csv("example");
    let result = fs::read_to_string("example.csv").expect("File was not created");
    let expected = [
        "DEPT,DT,RHOB,NPHI,SFLU,SFLA,ILM,ILD",
        "1670,123.45,2550,0.45,123.45,123.45,110.2,105.6",
        "1669.875,123.45,2550,0.45,123.45,123.45,110.2,105.6",
        "1669.75,123.45,2550,0.45,123.45,123.45,110.2,105.6",
        "1669.745,123.45,2550,-999.25,123.45,123.45,110.2,105.6",
    ];
    assert_eq!(expected.join("\n").to_string(), result);
    fs::remove_file("example.csv").expect("Could not clean up after test");
}