w-geo 1.1.21 → 1.1.23

package/src/dtIntrpPsdContentIso.mjs

@@ -0,0 +1,109 @@
+import get from 'lodash-es/get.js'
+import find from 'lodash-es/find.js'
+import round from 'lodash-es/round.js'
+import size from 'lodash-es/size.js'
+import isnum from 'wsemi/src/isnum.mjs'
+import isestr from 'wsemi/src/isestr.mjs'
+import cdbl from 'wsemi/src/cdbl.mjs'
+import dtIntrpPsdBySize from './dtIntrpPsdBySize.mjs'
+
+
+function dtIntrpPsdContentIso(dt, opt = {}) {
+
+    //keySize
+    let keySize = get(opt, 'keySize', '')
+    if (!isestr(keySize)) {
+        keySize = 'size'
+    }
+
+    //keyFraction
+    let keyFraction = get(opt, 'keyFraction', '')
+    if (!isestr(keyFraction)) {
+        keyFraction = 'fraction'
+    }
+
+    //dtIntrpPsdBySize
+    let psizes = [
+        630,
+        200,
+        63,
+        20,
+        6.3,
+        2,
+        0.63,
+        0.2,
+        0.063,
+        0.02,
+        0.0063,
+        0.002,
+    ]
+    let psds = dtIntrpPsdBySize(dt, psizes, { keySize, keyFraction })
+    // console.log('psds', psds)
+
+    //check
+    if (size(psds) === 0) {
+        console.log('dt', dt)
+        console.log('psizes', psizes)
+        console.log('keySize', keySize)
+        console.log('keyFraction', keyFraction)
+        throw new Error('invalid PSD')
+    }
+
+    //1e20->2: 礫石
+    //2->0.063: 砂
+    //0.063->0.02: 粉土
+    //0.02->0: 黏土
+    let p2 = find(psds, { size: 2 })
+    let p0063 = find(psds, { size: 0.063 })
+    // let p0005 = find(psds, { size: 0.005 })
+    let p0002 = find(psds, { size: 0.002 })
+    let f2 = get(p2, 'fraction')
+    let f0063 = get(p0063, 'fraction')
+    // let f0005 = get(p0005, 'fraction')
+    let f0002 = get(p0002, 'fraction')
+    // console.log('p2', p2, 'p0063', p0063, 'p0002', p0002)
+    // console.log('f2', f2, 'f0063', f0063, 'f0002', f0002)
+
+    //check
+    if (!isnum(f2)) {
+        throw new Error('invalid PSD in 2mm')
+    }
+    if (!isnum(f0063)) {
+        throw new Error('invalid PSD in 0.063mm')
+    }
+    if (!isnum(f0002)) {
+        throw new Error('invalid PSD in 0.002mm')
+    }
+    f2 = cdbl(f2)
+    f0063 = cdbl(f0063)
+    f0002 = cdbl(f0002)
+    // console.log('f2', f2, 'f0063', f0063, 'f0002', f0002)
+
+    //ctGravelISO, ctSandISO, ctFineISO, ctClayISO, 四捨五入可能會導致4項總和非100%
+    let ctGravelISO = 100 - f2
+    ctGravelISO = round(ctGravelISO)
+    let ctSandISO = f2 - f0063
+    ctSandISO = round(ctSandISO)
+    let ctSiltISO = f0063 - f0002
+    ctSiltISO = round(ctSiltISO)
+    let ctClayISO = f0002
+    ctClayISO = round(ctClayISO)
+    // console.log('ctGravelISO', ctGravelISO, 'ctSandISO', ctSandISO, 'ctSiltISO', ctSiltISO, 'ctClayISO', ctClayISO)
+
+    //add prop.
+    dt.ctCoarseISO = ctGravelISO + ctSandISO
+    dt.ctGravelISO = ctGravelISO
+    dt.ctSandISO = ctSandISO
+    dt.ctFineISO = ctSiltISO + ctClayISO
+    dt.ctSiltISO = ctSiltISO
+    dt.ctClayISO = ctClayISO
+    // console.log('ctGravel Ori', dt.ctGravel, 'New', ctGravelISO)
+    // console.log('ctSand Ori', dt.ctSand, 'New', ctSandISO)
+    // console.log('ctFine Ori', (cdbl(dt.ctSilt) + cdbl(dt.ctClay)), 'New', ctFineISO)
+    // console.log('ctClay Ori', dt.ctClay, 'New', ctClayISO)
+
+    return dt
+}
+
+
+export default dtIntrpPsdContentIso

package/test/buildInterpFun.test.mjs

@@ -0,0 +1,65 @@
+import assert from 'assert'
+import _ from 'lodash-es'
+import w from 'wsemi'
+import buildInterpFun from '../src/buildInterpFun.mjs'
+
+
+describe(`buildInterpFun`, function() {
+
+    let rows = [
+        {
+            depth: 0,
+            value: 1,
+        },
+        {
+            depth: 1,
+            value: 2,
+        },
+        {
+            depth: 3,
+            value: 3,
+        },
+        {
+            depth: 10,
+            value: 4,
+        },
+    ]
+
+    let fun1 = buildInterpFun(rows, 'depth', 'value')
+
+    let x1 = 0.5
+    // let y1 = fun1(x1)
+    // console.log(y1)
+    let y1 = 1
+
+    let xs1 = [0.5, 1, 1.9, 2, 2.1]
+    // let ys1 = fun1(xs1)
+    // console.log(ys1)
+    let ys1 = [1, 2, 2, 2, 3]
+
+    let fun2 = buildInterpFun(rows, 'depth', 'value', { mode: 'linear' })
+
+    let xs2 = [0.5, 1, 1.9, 2, 2.1]
+    // let ys2 = fun2(xs2)
+    // console.log(ys2)
+    let ys2 = [1.5, 2, 2.45, 2.5, 2.55]
+
+    it(`should return ${y1} when fun1(${x1}) by buildInterpFun(rows, 'depth', 'value')`, function() {
+        let r = fun1(x1)
+        let rr = y1
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+    it(`should return ${ys1} when fun1(${xs1}) by buildInterpFun(rows, 'depth', 'value')`, function() {
+        let r = fun1(xs1)
+        let rr = ys1
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+    it(`should return ${ys2} when fun1(${xs2}) by buildInterpFun(rows, 'depth', 'value', { mode: 'linear' })`, function() {
+        let r = fun2(xs2)
+        let rr = ys2
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+})

package/test/calcDepthMaxMin.test.mjs

@@ -0,0 +1,29 @@
+import assert from 'assert'
+import { calcDepthMaxMin } from '../src/calcDepthMaxMin.mjs'
+import rowsIn1 from './calcDepthMaxMin-rowsIn1.json' assert { type: "json" }
+
+
+describe(`calcDepthMaxMin`, function() {
+
+    let ropt = {
+        1: { approximate: false },
+        2: { approximate: true },
+    }
+    let rout = {
+        1: { depthMax: 99.98, depthMin: 0 },
+        2: { depthMax: 100, depthMin: 0 },
+    }
+
+    it(`should return ${JSON.stringify(rout[1])} when calcDepthMaxMin(rowsIn1, ${JSON.stringify(ropt[1])})`, function() {
+        let r = calcDepthMaxMin(rowsIn1, ropt[1])
+        let rr = rout[1]
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+    it(`should return ${JSON.stringify(rout[2])} when calcDepthMaxMin(rowsIn1, ${JSON.stringify(ropt[2])})`, function() {
+        let r = calcDepthMaxMin(rowsIn1, ropt[2])
+        let rr = rout[2]
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+})

package/test/dtIntrpPsdContentIso.test.mjs

@@ -0,0 +1,30 @@
+import assert from 'assert'
+import _ from 'lodash-es'
+import w from 'wsemi'
+import dtIntrpPsdContentIso from '../src/dtIntrpPsdContentIso.mjs'
+
+
+describe(`dtIntrpPsdContentIso`, function() {
+
+    let dt = {
+        GSD: `[101.6, 76.2, 50.8, 25.4, 19.0, 9.5, 4.75, 2.0, 0.85, 0.425, 0.25, 0.15, 0.106, 0.075, 0.027, 0.018, 0.011, 0.008, 0.006, 0.003, 0.001, 0.0, 0.0, 0.0]`,
+        GSP: `[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, -1.0, 100.0, 99.0, 96.0, 86.0, 74.0, 61.0, 51.0, 32.0, 16.0, 0.0, 0.0, 0.0]`,
+    }
+
+    let out1 = {
+        GSD: '[101.6, 76.2, 50.8, 25.4, 19.0, 9.5, 4.75, 2.0, 0.85, 0.425, 0.25, 0.15, 0.106, 0.075, 0.027, 0.018, 0.011, 0.008, 0.006, 0.003, 0.001, 0.0, 0.0, 0.0]',
+        GSP: '[100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, -1.0, 100.0, 99.0, 96.0, 86.0, 74.0, 61.0, 51.0, 32.0, 16.0, 0.0, 0.0, 0.0]',
+        ctCoarseISO: 2,
+        ctGravelISO: 0,
+        ctSandISO: 2,
+        ctFineISO: 98,
+        ctSiltISO: 72,
+        ctClayISO: 26
+    }
+    it(`should return ${out1} when dtIntrpPsdContentIso(dt)`, function() {
+        let r = dtIntrpPsdContentIso(dt)
+        let rr = out1
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+})

package/test/getDepthMaxMin.test.mjs

@@ -1,29 +0,0 @@
-import assert from 'assert'
-import { getDepthMaxMin } from '../src/getDepthMaxMin.mjs'
-import rowsIn1 from './getDepthMaxMin-rowsIn1.json' assert { type: "json" }
-
-
-describe(`getDepthMaxMin`, function() {
-
-    let ropt = {
-        1: { approximate: false },
-        2: { approximate: true },
-    }
-    let rout = {
-        1: { depthMax: 99.98, depthMin: 0 },
-        2: { depthMax: 100, depthMin: 0 },
-    }
-
-    it(`should return ${JSON.stringify(rout[1])} when getDepthMaxMin(rowsIn1, ${JSON.stringify(ropt[1])})`, function() {
-        let r = getDepthMaxMin(rowsIn1, ropt[1])
-        let rr = rout[1]
-        assert.strict.deepStrictEqual(r, rr)
-    })
-
-    it(`should return ${JSON.stringify(rout[2])} when getDepthMaxMin(rowsIn1, ${JSON.stringify(ropt[2])})`, function() {
-        let r = getDepthMaxMin(rowsIn1, ropt[2])
-        let rr = rout[2]
-        assert.strict.deepStrictEqual(r, rr)
-    })
-
-})