w-geo 1.1.20 → 1.1.21

package/docs/-_class.mjs.html

@@ -62,7 +62,7 @@
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcDepthByDepthStartEnd.mjs.html

@@ -233,7 +233,7 @@ export default calcDepthByDepthStartEnd
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcDepthStartEndByConnect.mjs.html

@@ -322,7 +322,7 @@ export default calcDepthStartEndByConnect
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcDepthStartEndByDepth.mjs.html

@@ -347,7 +347,7 @@ export default calcDepthStartEndByDepth
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcDepthStartEndByGroup.mjs.html

@@ -566,7 +566,7 @@ export default calcDepthStartEndByGroup
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcLiquefaction.mjs.html

@@ -7801,7 +7801,7 @@ export default calcLiquefaction
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/calcVerticalStress.mjs.html

@@ -733,7 +733,7 @@ export default { //整合輸出預設得要有default
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/checkDepth.mjs.html

@@ -214,7 +214,7 @@ export default checkDepth
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/checkDepthStartEnd.mjs.html

@@ -273,7 +273,7 @@ export default checkDepthStartEnd
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/complementDepthData.mjs.html

@@ -506,7 +506,7 @@ export default complementDepthData
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/groupByDepthStartEnd.mjs.html

@@ -492,7 +492,7 @@ export default groupByDepthStartEnd
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/index.html

@@ -71,7 +71,7 @@
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/mergeByDepthStartEnd.mjs.html

@@ -907,7 +907,7 @@ export default mergeByDepthStartEnd
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/relaPlasticityParams.mjs.html

@@ -334,7 +334,7 @@ export default relaPlasticityParams
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/relaPorousParams.mjs.html

@@ -524,7 +524,7 @@ export default relaPorousParams
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/sepDepthStartEndByDepth.mjs.html

@@ -324,7 +324,7 @@ export default sepDepthStartEndByDepth
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/smoothDepthByKey.mjs.html

@@ -280,7 +280,7 @@ export default smoothDepthByKey
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/docs/w-geo.html

@@ -9658,7 +9658,7 @@
 <br class="clear">
 
 <footer>
-    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Sat May 04 2024 00:06:28 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
+    Documentation generated by <a href="https://github.com/jsdoc3/jsdoc">JSDoc 4.0.2</a> on Wed May 08 2024 22:17:07 GMT+0800 (台北標準時間) using the <a href="https://github.com/clenemt/docdash">docdash</a> theme.
 </footer>
 
 <script>prettyPrint();</script>

package/g_7_1-intrpPsdBySize.mjs

@@ -0,0 +1,52 @@
+import fs from 'fs'
+import _ from 'lodash-es'
+import w from 'wsemi'
+import intrpPsdBySize from './src/intrpPsdBySize.mjs'
+
+
+let psd = [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]
+let psf = [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 pfss = _.map(psd, (_v, k) => {
+    let s = _.get(psd, k)
+    let f = _.get(psf, k)
+    return {
+        size: s,
+        fraction: f,
+    }
+})
+console.log('pfss', pfss)
+// => pfss [
+//   { size: 101.6, fraction: 100 },
+//   { size: 76.2, fraction: 100 },
+//   { size: 50.8, fraction: 100 },
+//   { size: 25.4, fraction: 100 },
+//   { size: 19, fraction: 100 },
+//   { size: 9.5, fraction: 100 },
+//   { size: 4.75, fraction: 100 },
+//   { size: 2, fraction: 100 },
+//   { size: 0.85, fraction: 100 },
+//   { size: 0.425, fraction: 100 },
+//   { size: 0.25, fraction: 100 },
+//   { size: 0.15, fraction: -1 },
+//   { size: 0.106, fraction: 100 },
+//   { size: 0.075, fraction: 99 },
+//   { size: 0.027, fraction: 96 },
+//   { size: 0.018, fraction: 86 },
+//   { size: 0.011, fraction: 74 },
+//   { size: 0.008, fraction: 61 },
+//   { size: 0.006, fraction: 51 },
+//   { size: 0.003, fraction: 32 },
+//   { size: 0.001, fraction: 16 },
+//   { size: 0, fraction: 0 },
+//   { size: 0, fraction: 0 },
+//   { size: 0, fraction: 0 }
+// ]
+
+let psizes = 0.063
+
+let r = intrpPsdBySize(pfss, psizes)
+console.log(r)
+// => 98.4880247416152
+
+
+//node --experimental-modules g_7_1-intrpPsdBySize.mjs

package/{g_7_9-calcPropInterfaceFrictionAngle.mjs → g_7_a1-calcPropInterfaceFrictionAngle.mjs}

@@ -105,4 +105,4 @@ await w.pmSeries(ms, async (m, i) => {
 
 })
 
-//node --experimental-modules g_7_9-calcPropInterfaceFrictionAngle.mjs
+//node --experimental-modules g_7_a1-calcPropInterfaceFrictionAngle.mjs

package/g_8_1-dtIntrpPsdBySize.mjs

@@ -0,0 +1,50 @@
+import fs from 'fs'
+import _ from 'lodash-es'
+import w from 'wsemi'
+import dtIntrpPsdBySize from './src/dtIntrpPsdBySize.mjs'
+
+
+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]`,
+}
+console.log('dt', dt)
+// => 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 psizes = [
+    630,
+    200,
+    63,
+    20,
+    6.3,
+    2,
+    0.63,
+    0.2,
+    0.063,
+    0.02,
+    0.0063,
+    0.002,
+]
+
+let fss = dtIntrpPsdBySize(dt, psizes)
+console.log('fss', fss)
+// => [
+//   { size: 630, fraction: 100 },
+//   { size: 200, fraction: 100 },
+//   { size: 63, fraction: 99.99999999999999 },
+//   { size: 20, fraction: 99.99999999999999 },
+//   { size: 6.3, fraction: 100.00000000000001 },
+//   { size: 2, fraction: 100 },
+//   { size: 0.63, fraction: 100.00000000000001 },
+//   { size: 0.2, fraction: 100 },
+//   { size: 0.063, fraction: 98.4880247416152 },
+//   { size: 0.02, fraction: 88.59851004564663 },
+//   { size: 0.0063, fraction: 52.69597513510718 },
+//   { size: 0.002, fraction: 26.094876057143324 }
+// ]
+
+
+//node --experimental-modules g_8_1-dtIntrpPsdBySize.mjs

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "w-geo",
-  "version": "1.1.20",
+  "version": "1.1.21",
   "main": "dist/w-geo.umd.js",
   "dependencies": {
     "decimal.js": "^10.4.3",

package/src/WGeo.mjs

@@ -25,10 +25,12 @@ export { default as checkDepth } from './checkDepth.mjs'
 export { default as checkDepthStartEnd } from './checkDepthStartEnd.mjs'
 export { default as cnst } from './cnst.mjs'
 export { default as complementDepthData } from './complementDepthData.mjs'
+export { default as dtIntrpPsdBySize } from './dtIntrpPsdBySize.mjs'
 export { default as getCptDepthHcOpt } from './getCptDepthHcOpt.mjs'
 export { default as getDepthMaxMin } from './getDepthMaxMin.mjs'
 export { default as groupByDepthStartEnd } from './groupByDepthStartEnd.mjs'
 export { default as intrpDefParam } from './intrpDefParam.mjs'
+export { default as intrpPsdBySize } from './intrpPsdBySize.mjs'
 export { default as judge } from './judge.mjs'
 export { default as mergeByDepthStartEnd } from './mergeByDepthStartEnd.mjs'
 export { default as relaPlasticityParams } from './relaPlasticityParams.mjs'

package/src/dtIntrpPsdBySize.mjs

@@ -0,0 +1,137 @@
+import each from 'lodash-es/each.js'
+import map from 'lodash-es/map.js'
+import get from 'lodash-es/get.js'
+import j2o from 'wsemi/src/j2o.mjs'
+import isearr from 'wsemi/src/isearr.mjs'
+import isestr from 'wsemi/src/isestr.mjs'
+import isnum from 'wsemi/src/isnum.mjs'
+import cdbl from 'wsemi/src/cdbl.mjs'
+import intrpPsdBySize from './intrpPsdBySize.mjs'
+
+
+//parseStrArray
+let parseStrArray = (c) => {
+    let r = j2o(c)
+    if (!isearr(r)) {
+        console.log('c', c)
+        throw new Error('無法解析JSON字串型陣列')
+    }
+    return r
+}
+
+
+function dtIntrpPsdBySize(dt, psizes, opt = {}) {
+
+    //keySize
+    let keySize = get(opt, 'keySize', '')
+    if (!isestr(keySize)) {
+        keySize = 'size'
+    }
+
+    //keyFraction
+    let keyFraction = get(opt, 'keyFraction', '')
+    if (!isestr(keyFraction)) {
+        keyFraction = 'fraction'
+    }
+
+    //GSD
+    let GSD = get(dt, 'GSD', '')
+    if (!isestr(GSD) && !isearr(GSD)) {
+        throw new Error('invalid GSD')
+    }
+    if (isestr(GSD)) {
+        GSD = parseStrArray(GSD)
+    }
+    // console.log('GSD', GSD)
+
+    //GSP
+    let GSP = get(dt, 'GSP', '')
+    if (!isestr(GSP) && !isearr(GSP)) {
+        throw new Error('invalid GSP')
+    }
+    if (isestr(GSP)) {
+        GSP = parseStrArray(GSP)
+    }
+    // console.log('GSP', GSP)
+
+    //psizes
+    if (!isearr(psizes)) {
+        //ASTM
+        psizes = [
+            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.034,
+            0.022,
+            0.013,
+            0.009,
+            0.007,
+            0.003,
+            0.001,
+        ]
+        // //ISO
+        // psizes = [
+        //     630,
+        //     200,
+        //     63,
+        //     20,
+        //     6.3,
+        //     2,
+        //     0.63,
+        //     0.2,
+        //     0.063,
+        //     0.02,
+        //     0.0063,
+        //     0.002,
+        // ]
+    }
+
+    //fss
+    let fss = []
+    each(GSD, (_v, k) => {
+        let size = get(GSD, k, '')
+        if (!isnum(size)) {
+            console.log(k, 'size', size)
+            throw new Error(`size[${size}] is not an effective number`)
+        }
+        size = cdbl(size)
+        let fraction = get(GSP, k, '')
+        if (!isnum(fraction)) {
+            console.log(k, 'fraction', fraction)
+            throw new Error(`fraction[${fraction}] is not an effective number`)
+        }
+        fraction = cdbl(fraction)
+        fss.push({
+            [keySize]: size,
+            [keyFraction]: fraction,
+        })
+    })
+    // console.log('fss', fss)
+
+    //重算各size處之通過百分比
+    let rs = intrpPsdBySize(fss, psizes, { keySize, keyFraction })
+    rs = map(rs, (fraction, k) => {
+        let size = get(psizes, k)
+        return {
+            [keySize]: size,
+            [keyFraction]: fraction,
+        }
+    })
+    // console.log('rs', rs)
+
+    return rs
+}
+
+
+export default dtIntrpPsdBySize

package/src/intrpPsdBySize.mjs

@@ -0,0 +1,125 @@
+import get from 'lodash-es/get.js'
+import size from 'lodash-es/size.js'
+import filter from 'lodash-es/filter.js'
+import map from 'lodash-es/map.js'
+import maxBy from 'lodash-es/maxBy.js'
+import minBy from 'lodash-es/minBy.js'
+import isnum from 'wsemi/src/isnum.mjs'
+import isearr from 'wsemi/src/isearr.mjs'
+import isestr from 'wsemi/src/isestr.mjs'
+import cdbl from 'wsemi/src/cdbl.mjs'
+import interp1 from 'wsemi/src/interp1.mjs'
+
+
+function intrpPsdBySize(pfss, psizes, opt = {}) {
+
+    //keySize
+    let keySize = get(opt, 'keySize', '')
+    if (!isestr(keySize)) {
+        keySize = 'size'
+    }
+
+    //keyFraction
+    let keyFraction = get(opt, 'keyFraction', '')
+    if (!isestr(keyFraction)) {
+        keyFraction = 'fraction'
+    }
+
+    //check
+    let pfs0 = get(pfss, 0)
+    let size0 = get(pfs0, keySize, '')
+    let fraction0 = get(pfs0, keyFraction, '')
+    if (!isnum(size0)) {
+        throw new Error(`pfssp[0][${keySize}] is not an effective number`)
+    }
+    if (!isnum(fraction0)) {
+        throw new Error(`fraction0[0][${keyFraction}] is not an effective number`)
+    }
+
+    //check
+    if (!isnum(psizes) && !isearr(psizes)) {
+        console.log('psizes', psizes)
+        throw new Error('psizes is not a number or an effective array')
+    }
+
+    //one
+    let one = false
+    if (isnum(psizes)) {
+        one = true
+        psizes = [cdbl(psizes)]
+    }
+
+    //取粒徑size與fraction為數字, 且size大於0(因要轉log內插), fraction要大於等於0, 且無施篩孔時值為-1亦須剔除
+    let rs = filter(pfss, (v) => {
+        let bs = cdbl(v[keySize]) > 0
+        let bf = cdbl(v[keyFraction]) >= 0
+        // console.log('bs', bs, v[keySize])
+        // console.log('bf', bf, v[keyFraction])
+        return bs && bf
+    })
+    // console.log('rs', rs)
+
+    //check
+    if (size(rs) === 0) {
+        console.log('pfss', pfss)
+        throw new Error(`invalid rs`)
+    }
+
+    //lrs
+    let lrs = map(rs, (v) => {
+        v[keySize] = Math.log(v[keySize])
+        return v
+    })
+
+    //sizeLogMax, sizeLogMin
+    let sizeLogMax = maxBy(lrs, keySize)
+    let sizeLogMin = minBy(lrs, keySize)
+
+    //ifs
+    let ifs = map(psizes, (psize) => {
+
+        //sizeLog
+        let sizeLog = Math.log(psize)
+        // console.log('psize', psize)
+        // console.log('sizeLog', sizeLog)
+
+        //check
+        if (sizeLog >= sizeLogMax[keySize]) {
+            return sizeLogMax[keyFraction]
+        }
+        if (sizeLog <= sizeLogMin[keySize]) {
+            return sizeLogMin[keyFraction]
+        }
+
+        //interp1
+        let optP = {
+            keyX: keySize,
+            keyY: keyFraction,
+        }
+        let pfraction = interp1(lrs, sizeLog, optP)
+        // console.log('pfraction', pfraction)
+
+        //check
+        if (get(pfraction, 'err')) {
+            console.log('lrs', lrs)
+            console.log('sizeLog', sizeLog)
+            console.log('optP', optP)
+            console.log('pfraction', pfraction)
+            return null
+        }
+
+        return pfraction
+    })
+    // console.log('ifs', ifs)
+
+    //r
+    let r = ifs
+    if (one) {
+        r = ifs[0]
+    }
+
+    return r
+}
+
+
+export default intrpPsdBySize

package/test/dtIntrpPsdBySize.test.mjs

@@ -0,0 +1,50 @@
+import assert from 'assert'
+import _ from 'lodash-es'
+import w from 'wsemi'
+import dtIntrpPsdBySize from '../src/dtIntrpPsdBySize.mjs'
+
+
+describe(`dtIntrpPsdBySize`, 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 psizes = [
+        630,
+        200,
+        63,
+        20,
+        6.3,
+        2,
+        0.63,
+        0.2,
+        0.063,
+        0.02,
+        0.0063,
+        0.002,
+    ]
+    // let fss = dtIntrpPsdBySize(dt, psizes)
+
+    let out1 = [
+        { size: 630, fraction: 100 },
+        { size: 200, fraction: 100 },
+        { size: 63, fraction: 99.99999999999999 },
+        { size: 20, fraction: 99.99999999999999 },
+        { size: 6.3, fraction: 100.00000000000001 },
+        { size: 2, fraction: 100 },
+        { size: 0.63, fraction: 100.00000000000001 },
+        { size: 0.2, fraction: 100 },
+        { size: 0.063, fraction: 98.4880247416152 },
+        { size: 0.02, fraction: 88.59851004564663 },
+        { size: 0.0063, fraction: 52.69597513510718 },
+        { size: 0.002, fraction: 26.094876057143324 }
+    ]
+    it(`should return ${out1} when dtIntrpPsdBySize(dt, ${psizes})`, function() {
+        let r = dtIntrpPsdBySize(dt, psizes)
+        let rr = out1
+        assert.strict.deepStrictEqual(r, rr)
+    })
+
+})