nmr-processing 8.1.0 → 8.2.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/lib/assignment/utils/getAssignment/checkIDs.d.ts +1 -1
- package/lib/assignment/utils/getAssignment/checkIDs.js.map +1 -1
- package/lib/index.d.ts +2 -0
- package/lib/index.js +2 -0
- package/lib/index.js.map +1 -1
- package/lib/peaks/peaksToRanges.d.ts +16 -0
- package/lib/peaks/peaksToRanges.js +2 -1
- package/lib/peaks/peaksToRanges.js.map +1 -1
- package/lib/peaks/solventSuppression.d.ts +1 -1
- package/lib/peaks/solventSuppression.js +34 -26
- package/lib/peaks/solventSuppression.js.map +1 -1
- package/lib/peaks/util/jAnalyzer.js +3 -2
- package/lib/peaks/util/jAnalyzer.js.map +1 -1
- package/lib/ranges/markSolventSignal.d.ts +3 -0
- package/lib/ranges/markSolventSignal.js +107 -0
- package/lib/ranges/markSolventSignal.js.map +1 -0
- package/lib/xyz/util/formatZone.d.ts +3 -0
- package/lib/xyz/util/formatZone.js +38 -0
- package/lib/xyz/util/formatZone.js.map +1 -0
- package/lib/xyz/xyzAutoZonesPicking.d.ts +1 -2
- package/lib/xyz/xyzAutoZonesPicking.js +2 -34
- package/lib/xyz/xyzAutoZonesPicking.js.map +1 -1
- package/lib/xyz/xyzJResAnalyzer.d.ts +6 -1
- package/lib/xyz/xyzJResAnalyzer.js +9 -8
- package/lib/xyz/xyzJResAnalyzer.js.map +1 -1
- package/lib-esm/assignment/utils/getAssignment/checkIDs.js.map +1 -1
- package/lib-esm/index.js +2 -0
- package/lib-esm/index.js.map +1 -1
- package/lib-esm/peaks/peaksToRanges.js +1 -1
- package/lib-esm/peaks/peaksToRanges.js.map +1 -1
- package/lib-esm/peaks/solventSuppression.js +32 -21
- package/lib-esm/peaks/solventSuppression.js.map +1 -1
- package/lib-esm/peaks/util/jAnalyzer.js +3 -2
- package/lib-esm/peaks/util/jAnalyzer.js.map +1 -1
- package/lib-esm/ranges/markSolventSignal.js +100 -0
- package/lib-esm/ranges/markSolventSignal.js.map +1 -0
- package/lib-esm/xyz/util/formatZone.js +34 -0
- package/lib-esm/xyz/util/formatZone.js.map +1 -0
- package/lib-esm/xyz/xyzAutoZonesPicking.js +1 -33
- package/lib-esm/xyz/xyzAutoZonesPicking.js.map +1 -1
- package/lib-esm/xyz/xyzJResAnalyzer.js +9 -8
- package/lib-esm/xyz/xyzJResAnalyzer.js.map +1 -1
- package/package.json +2 -2
- package/src/assignment/utils/getAssignment/checkIDs.ts +1 -1
- package/src/index.ts +3 -0
- package/src/peaks/peaksToRanges.ts +1 -1
- package/src/peaks/solventSuppression.ts +57 -34
- package/src/peaks/util/jAnalyzer.ts +4 -2
- package/src/ranges/markSolventSignal.ts +121 -0
- package/src/xyz/util/formatZone.ts +36 -0
- package/src/xyz/xyzAutoZonesPicking.ts +1 -35
- package/src/xyz/xyzJResAnalyzer.ts +14 -7
|
@@ -0,0 +1,100 @@
|
|
|
1
|
+
import generateID from '../assignment/utils/generateID';
|
|
2
|
+
import { addIDs } from '../assignment/utils/getAssignment/checkIDs';
|
|
3
|
+
import { detectSignals } from '../peaks/peaksToRanges';
|
|
4
|
+
import { solventSuppression } from '../peaks/solventSuppression';
|
|
5
|
+
export function markSolventPeaks(input, solventSignals, options = {}) {
|
|
6
|
+
const { frequency = 400 } = options;
|
|
7
|
+
const ranges = addIDs([...input]);
|
|
8
|
+
let peakList = [];
|
|
9
|
+
for (const range of ranges) {
|
|
10
|
+
const rangeID = range.id;
|
|
11
|
+
for (const signal of range.signals || []) {
|
|
12
|
+
const signalID = signal.id;
|
|
13
|
+
for (let peak of signal.peaks || []) {
|
|
14
|
+
peakList.push({
|
|
15
|
+
...peak,
|
|
16
|
+
rangeID,
|
|
17
|
+
signalID,
|
|
18
|
+
});
|
|
19
|
+
}
|
|
20
|
+
}
|
|
21
|
+
}
|
|
22
|
+
const peaks = solventSuppression(peakList, solventSignals, {
|
|
23
|
+
markSolventPeaks: true,
|
|
24
|
+
});
|
|
25
|
+
const signalsWithSolvent = getSignalIDsWithSolvent(peaks);
|
|
26
|
+
for (let range of ranges) {
|
|
27
|
+
if (!range.signals)
|
|
28
|
+
continue;
|
|
29
|
+
for (let signal of range.signals) {
|
|
30
|
+
if (!signalsWithSolvent.includes(signal.id))
|
|
31
|
+
continue;
|
|
32
|
+
const signalPeaks = [];
|
|
33
|
+
const solventSignalPeaks = [];
|
|
34
|
+
for (let peak of peaks) {
|
|
35
|
+
if (peak.signalID === signal.id) {
|
|
36
|
+
const { signalID, rangeID, ...newPeak } = peak;
|
|
37
|
+
if (peak.kind === 'solvent') {
|
|
38
|
+
solventSignalPeaks.push(newPeak);
|
|
39
|
+
}
|
|
40
|
+
else {
|
|
41
|
+
signalPeaks.push(newPeak);
|
|
42
|
+
}
|
|
43
|
+
}
|
|
44
|
+
}
|
|
45
|
+
const newSignals = signalPeaks.length > 0
|
|
46
|
+
? detectSignals({ x: [], y: [] }, signalPeaks, {
|
|
47
|
+
integralType: 'peak',
|
|
48
|
+
frequency,
|
|
49
|
+
})
|
|
50
|
+
: [];
|
|
51
|
+
solventSignalPeaks.sort((a, b) => a.x - b.x);
|
|
52
|
+
const lowPPMValue = solventSignalPeaks[0].x;
|
|
53
|
+
const highPPMValue = solventSignalPeaks[solventSignalPeaks.length - 1].x;
|
|
54
|
+
const newSolventSignals = detectSignals({ x: [], y: [] }, solventSignalPeaks, {
|
|
55
|
+
integralType: 'peak',
|
|
56
|
+
frequencyCluster: (highPPMValue - lowPPMValue) * frequency,
|
|
57
|
+
frequency,
|
|
58
|
+
});
|
|
59
|
+
range.signals = [
|
|
60
|
+
...adaptSignals(newSignals),
|
|
61
|
+
...adaptSignals(newSolventSignals),
|
|
62
|
+
...range.signals.filter((currSignal) => currSignal.id !== signal.id),
|
|
63
|
+
];
|
|
64
|
+
}
|
|
65
|
+
}
|
|
66
|
+
return ranges;
|
|
67
|
+
}
|
|
68
|
+
function getSignalIDsWithSolvent(solventPeaks) {
|
|
69
|
+
const uniqueSignalID = new Set();
|
|
70
|
+
for (let peak of solventPeaks) {
|
|
71
|
+
if (peak.kind !== 'solvent')
|
|
72
|
+
continue;
|
|
73
|
+
uniqueSignalID.add(peak.signalID);
|
|
74
|
+
}
|
|
75
|
+
return Array.from(uniqueSignalID);
|
|
76
|
+
}
|
|
77
|
+
function adaptSignals(signals) {
|
|
78
|
+
const newSignals = [];
|
|
79
|
+
for (let signal of signals) {
|
|
80
|
+
const signalResult = {
|
|
81
|
+
id: generateID(),
|
|
82
|
+
delta: signal.delta,
|
|
83
|
+
kind: signal.kind || 'signal',
|
|
84
|
+
multiplicity: signal.multiplicity,
|
|
85
|
+
integration: signal.integralData.value,
|
|
86
|
+
};
|
|
87
|
+
signalResult.peaks = signal.peaks.map((peak) => {
|
|
88
|
+
const newResult = {
|
|
89
|
+
y: peak.intensity,
|
|
90
|
+
...peak,
|
|
91
|
+
};
|
|
92
|
+
delete newResult.intensity;
|
|
93
|
+
return newResult;
|
|
94
|
+
});
|
|
95
|
+
signalResult.js = signal.nmrJs || [];
|
|
96
|
+
newSignals.push(signalResult);
|
|
97
|
+
}
|
|
98
|
+
return newSignals;
|
|
99
|
+
}
|
|
100
|
+
//# sourceMappingURL=markSolventSignal.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"markSolventSignal.js","sourceRoot":"","sources":["../../src/ranges/markSolventSignal.ts"],"names":[],"mappings":"AACA,OAAO,UAAU,MAAM,gCAAgC,CAAC;AACxD,OAAO,EAAE,MAAM,EAAE,MAAM,4CAA4C,CAAC;AAEpE,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AACvD,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AAUjE,MAAM,UAAU,gBAAgB,CAC9B,KAAU,EACV,cAA6B,EAC7B,UAAe,EAAE;IAEjB,MAAM,EAAE,SAAS,GAAG,GAAG,EAAE,GAAG,OAAO,CAAC;IACpC,MAAM,MAAM,GAAG,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,CAAsB,CAAC;IAEvD,IAAI,QAAQ,GAA4B,EAAE,CAAC;IAC3C,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE;QAC1B,MAAM,OAAO,GAAG,KAAK,CAAC,EAAE,CAAC;QACzB,KAAK,MAAM,MAAM,IAAI,KAAK,CAAC,OAAO,IAAI,EAAE,EAAE;YACxC,MAAM,QAAQ,GAAG,MAAM,CAAC,EAAE,CAAC;YAC3B,KAAK,IAAI,IAAI,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,EAAE;gBACnC,QAAQ,CAAC,IAAI,CAAC;oBACZ,GAAG,IAAI;oBACP,OAAO;oBACP,QAAQ;iBACT,CAAC,CAAC;aACJ;SACF;KACF;IAED,MAAM,KAAK,GAAG,kBAAkB,CAAC,QAAQ,EAAE,cAAc,EAAE;QACzD,gBAAgB,EAAE,IAAI;KACvB,CAAC,CAAC;IAEH,MAAM,kBAAkB,GAAG,uBAAuB,CAAC,KAAK,CAAC,CAAC;IAC1D,KAAK,IAAI,KAAK,IAAI,MAAM,EAAE;QACxB,IAAI,CAAC,KAAK,CAAC,OAAO;YAAE,SAAS;QAC7B,KAAK,IAAI,MAAM,IAAI,KAAK,CAAC,OAAO,EAAE;YAChC,IAAI,CAAC,kBAAkB,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;gBAAE,SAAS;YACtD,MAAM,WAAW,GAAG,EAAE,CAAC;YACvB,MAAM,kBAAkB,GAAG,EAAE,CAAC;YAC9B,KAAK,IAAI,IAAI,IAAI,KAAK,EAAE;gBACtB,IAAI,IAAI,CAAC,QAAQ,KAAK,MAAM,CAAC,EAAE,EAAE;oBAC/B,MAAM,EAAE,QAAQ,EAAE,OAAO,EAAE,GAAG,OAAO,EAAE,GAAG,IAAI,CAAC;oBAC/C,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS,EAAE;wBAC3B,kBAAkB,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;qBAClC;yBAAM;wBACL,WAAW,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;qBAC3B;iBACF;aACF;YACD,MAAM,UAAU,GACd,WAAW,CAAC,MAAM,GAAG,CAAC;gBACpB,CAAC,CAAC,aAAa,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,WAAW,EAAE;oBAC3C,YAAY,EAAE,MAAM;oBACpB,SAAS;iBACV,CAAC;gBACJ,CAAC,CAAC,EAAE,CAAC;YACT,kBAAkB,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7C,MAAM,WAAW,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5C,MAAM,YAAY,GAAG,kBAAkB,CAAC,kBAAkB,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACzE,MAAM,iBAAiB,GAAG,aAAa,CACrC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAChB,kBAAkB,EAClB;gBACE,YAAY,EAAE,MAAM;gBACpB,gBAAgB,EAAE,CAAC,YAAY,GAAG,WAAW,CAAC,GAAG,SAAS;gBAC1D,SAAS;aACV,CACF,CAAC;YAEF,KAAK,CAAC,OAAO,GAAG;gBACd,GAAG,YAAY,CAAC,UAAU,CAAC;gBAC3B,GAAG,YAAY,CAAC,iBAAiB,CAAC;gBAClC,GAAG,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,UAAU,EAAE,EAAE,CAAC,UAAU,CAAC,EAAE,KAAK,MAAM,CAAC,EAAE,CAAC;aACrE,CAAC;SACH;KACF;IACD,OAAO,MAAa,CAAC;AACvB,CAAC;AAED,SAAS,uBAAuB,CAAC,YAAqC;IACpE,MAAM,cAAc,GAAG,IAAI,GAAG,EAAE,CAAC;IACjC,KAAK,IAAI,IAAI,IAAI,YAAY,EAAE;QAC7B,IAAI,IAAI,CAAC,IAAI,KAAK,SAAS;YAAE,SAAS;QACtC,cAAc,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;KACnC;IACD,OAAO,KAAK,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;AACpC,CAAC;AAED,SAAS,YAAY,CAAC,OAAuB;IAC3C,MAAM,UAAU,GAAG,EAAE,CAAC;IACtB,KAAK,IAAI,MAAM,IAAI,OAAO,EAAE;QAC1B,MAAM,YAAY,GAAsB;YACtC,EAAE,EAAE,UAAU,EAAE;YAChB,KAAK,EAAE,MAAM,CAAC,KAAK;YACnB,IAAI,EAAE,MAAM,CAAC,IAAI,IAAI,QAAQ;YAC7B,YAAY,EAAE,MAAM,CAAC,YAAY;YACjC,WAAW,EAAE,MAAM,CAAC,YAAY,CAAC,KAAK;SACvC,CAAC;QACF,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE;YAC7C,MAAM,SAAS,GAAQ;gBACrB,CAAC,EAAE,IAAI,CAAC,SAAS;gBACjB,GAAG,IAAI;aACR,CAAC;YACF,OAAO,SAAS,CAAC,SAAS,CAAC;YAC3B,OAAO,SAAsB,CAAC;QAChC,CAAC,CAAC,CAAC;QACH,YAAY,CAAC,EAAE,GAAG,MAAM,CAAC,KAAK,IAAI,EAAE,CAAC;QACrC,UAAU,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;KAC/B;IACD,OAAO,UAAU,CAAC;AACpB,CAAC"}
|
|
@@ -0,0 +1,34 @@
|
|
|
1
|
+
export function formatZones(signals) {
|
|
2
|
+
let zones = [];
|
|
3
|
+
for (const signal of signals) {
|
|
4
|
+
let minMax1 = [Number.MAX_VALUE, 0];
|
|
5
|
+
let minMax2 = [Number.MAX_VALUE, 0];
|
|
6
|
+
for (const peak of signal.peaks || []) {
|
|
7
|
+
if (peak.minX < minMax1[0]) {
|
|
8
|
+
minMax1[0] = peak.minX;
|
|
9
|
+
}
|
|
10
|
+
if (peak.maxX > minMax1[1]) {
|
|
11
|
+
minMax1[1] = peak.maxX;
|
|
12
|
+
}
|
|
13
|
+
if (peak.minY < minMax2[0]) {
|
|
14
|
+
minMax2[0] = peak.minY;
|
|
15
|
+
}
|
|
16
|
+
if (peak.maxY > minMax2[1]) {
|
|
17
|
+
minMax2[1] = peak.maxY;
|
|
18
|
+
}
|
|
19
|
+
}
|
|
20
|
+
zones.push({
|
|
21
|
+
x: {
|
|
22
|
+
from: minMax1[0],
|
|
23
|
+
to: minMax1[1],
|
|
24
|
+
},
|
|
25
|
+
y: {
|
|
26
|
+
from: minMax2[0],
|
|
27
|
+
to: minMax2[1],
|
|
28
|
+
},
|
|
29
|
+
signals: [signal],
|
|
30
|
+
});
|
|
31
|
+
}
|
|
32
|
+
return zones;
|
|
33
|
+
}
|
|
34
|
+
//# sourceMappingURL=formatZone.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"formatZone.js","sourceRoot":"","sources":["../../../src/xyz/util/formatZone.ts"],"names":[],"mappings":"AAGA,MAAM,UAAU,WAAW,CAAC,OAAsB;IAChD,IAAI,KAAK,GAAc,EAAE,CAAC;IAC1B,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE;QAC5B,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACpC,KAAK,MAAM,IAAI,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,EAAE;YACrC,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;YACD,IAAI,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;gBAC1B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;aACxB;SACF;QACD,KAAK,CAAC,IAAI,CAAC;YACT,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,CAAC,EAAE;gBACD,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC;gBAChB,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC;aACf;YACD,OAAO,EAAE,CAAC,MAAM,CAAC;SAClB,CAAC,CAAC;KACJ;IACD,OAAO,KAAK,CAAC;AACf,CAAC"}
|
|
@@ -5,6 +5,7 @@ import simpleClustering from 'ml-simple-clustering';
|
|
|
5
5
|
import { determineRealTop } from '../peaks/util/determineRealTop';
|
|
6
6
|
import { getKernel } from '../peaks/util/getKernel';
|
|
7
7
|
import * as PeakOptimizer from '../peaks/util/peakOptimizer';
|
|
8
|
+
import { formatZones } from './util/formatZone';
|
|
8
9
|
const smallFilter = [
|
|
9
10
|
[0, 0, 1, 2, 2, 2, 1, 0, 0],
|
|
10
11
|
[0, 1, 4, 7, 7, 7, 4, 1, 0],
|
|
@@ -80,39 +81,6 @@ export function xyzAutoZonesPicking(spectraData, options) {
|
|
|
80
81
|
}
|
|
81
82
|
return formatZones(signals);
|
|
82
83
|
}
|
|
83
|
-
function formatZones(signals) {
|
|
84
|
-
let zones = [];
|
|
85
|
-
for (const signal of signals) {
|
|
86
|
-
let minMax1 = [Number.MAX_VALUE, 0];
|
|
87
|
-
let minMax2 = [Number.MAX_VALUE, 0];
|
|
88
|
-
for (const peak of signal.peaks || []) {
|
|
89
|
-
if (peak.minX < minMax1[0]) {
|
|
90
|
-
minMax1[0] = peak.minX;
|
|
91
|
-
}
|
|
92
|
-
if (peak.maxX > minMax1[1]) {
|
|
93
|
-
minMax1[1] = peak.maxX;
|
|
94
|
-
}
|
|
95
|
-
if (peak.minY < minMax2[0]) {
|
|
96
|
-
minMax2[0] = peak.minY;
|
|
97
|
-
}
|
|
98
|
-
if (peak.maxY > minMax2[1]) {
|
|
99
|
-
minMax2[1] = peak.maxY;
|
|
100
|
-
}
|
|
101
|
-
}
|
|
102
|
-
zones.push({
|
|
103
|
-
x: {
|
|
104
|
-
from: minMax1[0],
|
|
105
|
-
to: minMax1[1],
|
|
106
|
-
},
|
|
107
|
-
y: {
|
|
108
|
-
from: minMax2[0],
|
|
109
|
-
to: minMax2[1],
|
|
110
|
-
},
|
|
111
|
-
signals: [signal],
|
|
112
|
-
});
|
|
113
|
-
}
|
|
114
|
-
return zones;
|
|
115
|
-
}
|
|
116
84
|
const createSignals2D = (peaks, options) => {
|
|
117
85
|
let { nCols, nRows, absoluteData, originalData, observedFrequencies, tolerances, nuclei, realTopDetection, minY, maxY, minX, maxX, } = options;
|
|
118
86
|
let [nucleusX, nucleusY] = nuclei;
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyzAutoZonesPicking.js","sourceRoot":"","sources":["../../src/xyz/xyzAutoZonesPicking.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AACnC,OAAO,KAAK,WAAW,MAAM,uBAAuB,CAAC;AACrD,OAAO,KAAK,iBAAiB,MAAM,wBAAwB,CAAC;AAE5D,OAAO,gBAAgB,MAAM,sBAAsB,CAAC;
|
|
1
|
+
{"version":3,"file":"xyzAutoZonesPicking.js","sourceRoot":"","sources":["../../src/xyz/xyzAutoZonesPicking.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AACnC,OAAO,KAAK,WAAW,MAAM,uBAAuB,CAAC;AACrD,OAAO,KAAK,iBAAiB,MAAM,wBAAwB,CAAC;AAE5D,OAAO,gBAAgB,MAAM,sBAAsB,CAAC;AAEpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAEpD,OAAO,KAAK,aAAa,MAAM,6BAA6B,CAAC;AAG7D,OAAO,EAAE,WAAW,EAAE,MAAM,mBAAmB,CAAC;AAEhD,MAAM,WAAW,GAAG;IAClB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CAC5B,CAAC;AAkEF,MAAM,UAAU,mBAAmB,CACjC,WAAmB,EACnB,OAAmC;IAEnC,IAAI,EACF,SAAS,GAAG,EAAE,EACd,gBAAgB,GAAG,IAAI,EACvB,eAAe,GAAG,GAAG,EACrB,MAAM,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,EACrB,mBAAmB,EACnB,eAAe,GAAG,KAAK,EACvB,KAAK,GAAG,IAAI,EACZ,gBAAgB,GAAG,IAAI,EACvB,UAAU,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,EACrB,gBAAgB,GAAG,IAAI,EACvB,MAAM,EAAE,aAAa,GACtB,GAAG,OAAO,CAAC;IAEZ,IACE,CAAC,KAAK,CAAC,OAAO,CAAC,mBAAmB,CAAC;QACnC,CAAC,WAAW,CAAC,MAAM,CAAC,mBAAmB,CAAC,EACxC;QACA,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;KACrD;IAED,eAAe,GAAG,eAAe,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC;IAExE,IAAI,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACvC,IAAI,YAAY,GAAG,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC;IAExC,IAAI,YAAY,GAAG,SAAS,EAAE;QAC5B,WAAW,GAAG,OAAO,CAAC,WAAW,EAAE,EAAE,KAAK,EAAE,SAAS,EAAE,CAAC,CAAC;QACzD,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QACnC,YAAY,GAAG,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC;KACrC;IAED,IAAI,YAAY,GAAG,IAAI,YAAY,CAAC,QAAQ,GAAG,YAAY,CAAC,CAAC;IAC7D,IAAI,YAAY,GAAG,IAAI,YAAY,CAAC,QAAQ,GAAG,YAAY,CAAC,CAAC;IAE7D,KAAK,IAAI,WAAW,GAAG,CAAC,EAAE,WAAW,GAAG,YAAY,EAAE,WAAW,EAAE,EAAE;QACnE,IAAI,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC;QAC1C,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,QAAQ,EAAE,IAAI,EAAE,EAAE;YAC1C,IAAI,KAAK,GAAG,WAAW,GAAG,QAAQ,GAAG,IAAI,CAAC;YAC1C,YAAY,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC;YAC/C,YAAY,CAAC,KAAK,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,8CAA8C;SACrF;KACF;IAED,MAAM,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,SAAS,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC;IAEtE,IAAI,kBAAkB,GAAG,gBAAgB;QACvC,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,YAAY,EAAE,MAAM,EAAE;YACpC,IAAI,EAAE,YAAY;YAClB,IAAI,EAAE,QAAQ;SACf,CAAC;QACJ,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,YAAY,EAAE,MAAM,EAAE;YACvC,IAAI,EAAE,YAAY;YAClB,IAAI,EAAE,QAAQ;SACf,CAAC,CAAC;IAEP,IAAI,QAAQ,GAAG,iBAAiB,CAAC,iBAAiB,CAAC,YAAY,EAAE;QAC/D,YAAY;QACZ,YAAY,EAAE,kBAAkB;QAChC,IAAI,EAAE,YAAY;QAClB,IAAI,EAAE,QAAQ;QACd,OAAO,EAAE,eAAe;KACzB,CAAC,CAAC;IAEH,IAAI,KAAK,EAAE;QACT,sEAAsE;QACtE,QAAQ,GAAG,aAAa,CAAC,KAAK,CAAC,QAAQ,EAAE,gBAAgB,CAAC,CAAC;KAC5D;IAED,IAAI,OAAO,GAAG,eAAe,CAAC,QAAQ,EAAE;QACtC,KAAK,EAAE,YAAY;QACnB,KAAK,EAAE,QAAQ;QACf,IAAI,EAAE,WAAW,CAAC,IAAI;QACtB,IAAI,EAAE,WAAW,CAAC,IAAI;QACtB,IAAI,EAAE,WAAW,CAAC,IAAI;QACtB,IAAI,EAAE,WAAW,CAAC,IAAI;QACtB,YAAY;QACZ,YAAY;QACZ,UAAU;QACV,MAAM;QACN,mBAAmB;QACnB,gBAAgB;KACjB,CAAC,CAAC;IAEH,IAAI,eAAe,EAAE;QACnB,OAAO,GAAG,aAAa,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;KAClD;IAED,OAAO,WAAW,CAAC,OAAO,CAAC,CAAC;AAC9B,CAAC;AAsBD,MAAM,eAAe,GAAG,CAAC,KAAe,EAAE,OAA+B,EAAE,EAAE;IAC3E,IAAI,EACF,KAAK,EACL,KAAK,EACL,YAAY,EACZ,YAAY,EACZ,mBAAmB,EACnB,UAAU,EACV,MAAM,EACN,gBAAgB,EAChB,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,GACL,GAAG,OAAO,CAAC;IAEZ,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,GAAG,MAAM,CAAC;IAClC,IAAI,CAAC,UAAU,EAAE,UAAU,CAAC,GAAG,UAAU,CAAC;IAC1C,IAAI,CAAC,iBAAiB,EAAE,iBAAiB,CAAC,GAAG,mBAAmB,CAAC;IAEjE,IAAI,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACrC,IAAI,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IAErC,IAAI,gBAAgB,EAAE;QACpB,KAAK,GAAG,gBAAgB,CAAC,KAAK,EAAE;YAC9B,KAAK;YACL,YAAY;YACZ,YAAY;YACZ,IAAI;YACJ,IAAI;YACJ,IAAI;YACJ,IAAI;SACL,CAAC,CAAC;KACJ;IAED,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC1C,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAC1C,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAC1C,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAC1C,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAE1C,oGAAoG;KACrG;IACD,0GAA0G;IAC1G,kBAAkB;IAClB,IAAI,YAAY,GAAG,EAAE,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACrC,IACE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,iBAAiB,GAAG,UAAU;gBAClE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,iBAAiB,GAAG,UAAU,EAClE;gBACA,yEAAyE;gBACzE,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aACtB;iBAAM;gBACL,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aACtB;SACF;KACF;IACD,IAAI,QAAQ,GAAG,gBAAgB,CAAC,YAAY,CAAC,CAAC;IAE9C,IAAI,OAAO,GAAkB,EAAE,CAAC;IAChC,IAAI,QAAQ,EAAE;QACZ,KAAK,MAAM,OAAO,IAAI,QAAQ,EAAE;YAC9B,IAAI,MAAM,GAAQ;gBAChB,CAAC,EAAE;oBACD,KAAK,EAAE,CAAC;oBACR,OAAO,EAAE,QAAQ;oBACjB,UAAU,EAAE,EAAE;iBACf;gBACD,CAAC,EAAE;oBACD,KAAK,EAAE,CAAC;oBACR,OAAO,EAAE,QAAQ;oBACjB,UAAU,EAAE,EAAE;iBACf;aACF,CAAC;YACF,IAAI,OAAO,GAAa,EAAE,CAAC;YAC3B,IAAI,IAAI,GAAG,CAAC,CAAC;YAEb,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE;gBACnD,IAAI,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE;oBACxB,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;oBAC3B,MAAM,CAAC,CAAC,CAAC,KAAK,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAClD,MAAM,CAAC,CAAC,CAAC,KAAK,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAClD,IAAI,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;iBACxB;aACF;YAED,MAAM,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC;YACvB,MAAM,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC;YACvB,MAAM,CAAC,KAAK,GAAG,OAAO,CAAC;YACvB,OAAO,CAAC,IAAI,CAAC,MAAqB,CAAC,CAAC;SACrC;KACF;IACD,OAAO,OAAO,CAAC;AACjB,CAAC,CAAC;AAEF,MAAM,OAAO,GAAG,CAAC,WAAmB,EAAE,OAA0B,EAAE,EAAE;IAClE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,WAAW,CAAC;IAE7C,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IAE5B,IAAI,QAAQ,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IACvC,IAAI,YAAY,GAAG,WAAW,CAAC,CAAC,CAAC,MAAM,CAAC;IAExC,IAAI,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;IACnD,IAAI,SAAS,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;IAE/C,IAAI,KAAK,GAAG,KAAK,GAAG,YAAY,CAAC;IACjC,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,QAAQ,EAAE,CAAC,CAAC,CAAC;IAC1C,IAAI,KAAK,GAAG,CAAC;QAAE,KAAK,EAAE,CAAC;IACvB,IAAI,KAAK,GAAG,CAAC;QAAE,KAAK,EAAE,CAAC;IAEvB,IAAI,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IACxB,IAAI,OAAO,GAAG,KAAK,GAAG,CAAC,CAAC;IACxB,IAAI,SAAS,GAAG,MAAM,CAAC,KAAK,CAAC,YAAY,GAAG,KAAK,EAAE,QAAQ,GAAG,KAAK,CAAC,CAAC;IACrE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,CAAC,EAAE,EAAE;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE;YACjC,SAAS,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,EAAE,CAAC,GAAG,OAAO,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC9D;KACF;IAED,OAAO;QACL,CAAC,EAAE,SAAS,CAAC,SAAS,EAAE;QACxB,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;QAChC,IAAI,EAAE,IAAI,GAAG,OAAO,GAAG,SAAS;KACjC,CAAC;AACJ,CAAC,CAAC"}
|
|
@@ -1,24 +1,25 @@
|
|
|
1
1
|
import simpleClustering from 'ml-simple-clustering';
|
|
2
2
|
import jAnalyzer from '../peaks/util/jAnalyzer';
|
|
3
|
+
import { formatZones } from './util/formatZone';
|
|
3
4
|
export function xyzJResAnalyzer(signals, options = {}) {
|
|
4
|
-
let { reference = 0, referenceMaxShiftError = 0.08, tolerances = [10, 100], nuclei = ['1H', '1H'], observedFrequencies = [400, 400], jAnalyzer = {
|
|
5
|
+
let { reference = 0, getZones = false, referenceMaxShiftError = 0.08, tolerances = [10, 100], nuclei = ['1H', '1H'], observedFrequencies = [400, 400], jAnalyzer = {
|
|
5
6
|
jAxisKey: { jAxis: 'y', intensity: 'z' },
|
|
6
7
|
}, } = options;
|
|
7
|
-
let temporalSignals = compilePattern(signals, {
|
|
8
|
+
let temporalSignals = compilePattern([...signals], {
|
|
8
9
|
observedFrequencies,
|
|
9
10
|
tolerances,
|
|
10
11
|
nuclei,
|
|
11
12
|
jAnalyzer,
|
|
12
13
|
});
|
|
13
14
|
//check if the signal are symmetric around the reference
|
|
14
|
-
let
|
|
15
|
+
let signals2D = [];
|
|
15
16
|
for (const tempSignal of temporalSignals) {
|
|
16
17
|
let delta = tempSignal.y.delta;
|
|
17
18
|
if (Math.abs(delta - reference) > referenceMaxShiftError)
|
|
18
19
|
continue;
|
|
19
|
-
|
|
20
|
+
signals2D.push(tempSignal);
|
|
20
21
|
}
|
|
21
|
-
return
|
|
22
|
+
return getZones ? formatZones(signals2D) : signals2D;
|
|
22
23
|
}
|
|
23
24
|
function compilePattern(signals, options) {
|
|
24
25
|
let { observedFrequencies, tolerances, nuclei, jAnalyzer: jAnalyzerOptions, } = options;
|
|
@@ -67,9 +68,9 @@ function compilePattern(signals, options) {
|
|
|
67
68
|
}
|
|
68
69
|
if (peaksO.length > 0) {
|
|
69
70
|
peaksO.reverse();
|
|
70
|
-
let
|
|
71
|
-
for (const
|
|
72
|
-
newSignals.push(
|
|
71
|
+
let signals2D = createSignals2D(peaksO, signalOptions);
|
|
72
|
+
for (const signal of signals2D) {
|
|
73
|
+
newSignals.push(signal);
|
|
73
74
|
}
|
|
74
75
|
}
|
|
75
76
|
}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"xyzJResAnalyzer.js","sourceRoot":"","sources":["../../src/xyz/xyzJResAnalyzer.ts"],"names":[],"mappings":"AACA,OAAO,gBAAgB,MAAM,sBAAsB,CAAC;AAGpD,OAAO,SAAS,MAAM,yBAAyB,CAAC;
|
|
1
|
+
{"version":3,"file":"xyzJResAnalyzer.js","sourceRoot":"","sources":["../../src/xyz/xyzJResAnalyzer.ts"],"names":[],"mappings":"AACA,OAAO,gBAAgB,MAAM,sBAAsB,CAAC;AAGpD,OAAO,SAAS,MAAM,yBAAyB,CAAC;AAIhD,OAAO,EAAE,WAAW,EAAE,MAAM,mBAAmB,CAAC;AAgChD,MAAM,UAAU,eAAe,CAC7B,OAAsB,EACtB,UAAkC,EAAE;IAEpC,IAAI,EACF,SAAS,GAAG,CAAC,EACb,QAAQ,GAAG,KAAK,EAChB,sBAAsB,GAAG,IAAI,EAC7B,UAAU,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EACtB,MAAM,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,EACrB,mBAAmB,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,EAChC,SAAS,GAAG;QACV,QAAQ,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,SAAS,EAAE,GAAG,EAAE;KACzC,GACF,GAAG,OAAO,CAAC;IACZ,IAAI,eAAe,GAAG,cAAc,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE;QACjD,mBAAmB;QACnB,UAAU;QACV,MAAM;QACN,SAAS;KACV,CAAC,CAAC;IACH,wDAAwD;IACxD,IAAI,SAAS,GAAG,EAAE,CAAC;IACnB,KAAK,MAAM,UAAU,IAAI,eAAe,EAAE;QACxC,IAAI,KAAK,GAAG,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;QAC/B,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,SAAS,CAAC,GAAG,sBAAsB;YAAE,SAAS;QACnE,SAAS,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;KAC5B;IACD,OAAO,QAAQ,CAAC,CAAC,CAAC,WAAW,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;AACvD,CAAC;AAED,SAAS,cAAc,CACrB,OAAsB,EACtB,OAAuC;IAEvC,IAAI,EACF,mBAAmB,EACnB,UAAU,EACV,MAAM,EACN,SAAS,EAAE,gBAAgB,GAC5B,GAAG,OAAO,CAAC;IAEZ,IAAI,aAAa,GAAG;QAClB,mBAAmB;QACnB,UAAU;QACV,MAAM;QACN,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,UAAU;QAC3B,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,UAAU;KAC5B,CAAC;IAEF,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;IACvD,uBAAuB;IACvB,KAAK,MAAM,MAAM,IAAI,UAAU,EAAE;QAC/B,IAAI,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC;QACzB,MAAM,CAAC,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC;QACrC,MAAM,CAAC,YAAY,GAAG,EAAE,CAAC;QACzB,MAAM,CAAC,OAAO,GAAG,EAAE,CAAC;QACpB,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;QAC9B,MAAM,CAAC,OAAO,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,CAAC,YAAY,GAAG;YACpB,IAAI,EAAE,MAAM,CAAC,gBAAgB;YAC7B,EAAE,EAAE,MAAM,CAAC,gBAAgB;SAC5B,CAAC;QACF,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE;YACxB,IAAI,CAAC,IAAI,CAAC,KAAK;gBAAE,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;SACpC;QACD,KAAK,CAAC,IAAI,CAAC,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;KACjD;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACvC,SAAS,CAAC,cAAc,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC;QAE1D,IACE,UAAU,CAAC,CAAC,CAAC,CAAC,WAAW;YACzB,UAAU,CAAC,CAAC,CAAC,CAAC,YAAY,KAAK,GAAG;YAClC,UAAU,CAAC,CAAC,CAAC,CAAC,YAAY,KAAK,EAAE,EACjC;YACA,6CAA6C;YAC7C,IAAI,MAAM,GAAG,EAAE,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC9D,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,KAAK,KAAK,EAAE;oBAC1C,IAAI,KAAK,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAChD,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;oBACnB,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBAChC,UAAU,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBACjC,UAAU,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBACvC,UAAU,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC;iBACzB;aACF;YACD,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE;gBACrB,MAAM,CAAC,OAAO,EAAE,CAAC;gBACjB,IAAI,SAAS,GAAG,eAAe,CAAC,MAAM,EAAE,aAAa,CAAC,CAAC;gBAEvD,KAAK,MAAM,MAAM,IAAI,SAAS,EAAE;oBAC9B,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;iBACzB;aACF;SACF;KACF;IAED,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACpB,OAAO,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,OAAO,OAAO,CAAC;AACjB,CAAC;AAcD,SAAS,eAAe,CAAC,UAAoB,EAAE,OAAY;IACzD,IAAI,EAAE,mBAAmB,EAAE,UAAU,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,OAAO,CAAC;IAElE,MAAM,KAAK,GAAiB,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,CAAC,CAAC;IAEnE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,GAAG,MAAM,CAAC;IAClC,IAAI,CAAC,UAAU,EAAE,UAAU,CAAC,GAAG,UAAU,CAAC;IAC1C,IAAI,CAAC,iBAAiB,EAAE,iBAAiB,CAAC,GAAG,mBAAmB,CAAC;IAEjE,0GAA0G;IAC1G,kBAAkB;IAClB,IAAI,YAAY,GAAG,EAAE,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACrC,IACE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,iBAAiB,GAAG,UAAU;gBAClE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,iBAAiB,GAAG,UAAU,EAClE;gBACA,yEAAyE;gBACzE,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aACtB;iBAAM;gBACL,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aACtB;SACF;KACF;IACD,IAAI,QAAQ,GAAG,gBAAgB,CAAC,YAAY,CAAC,CAAC;IAE9C,IAAI,OAAO,GAAG,EAAE,CAAC;IACjB,IAAI,QAAQ,EAAE;QACZ,KAAK,IAAI,OAAO,IAAI,QAAQ,EAAE;YAC5B,IAAI,MAAM,GAAQ;gBAChB,QAAQ;gBACR,QAAQ;gBACR,YAAY,EAAE;oBACZ,IAAI,EAAE,MAAM,CAAC,gBAAgB;oBAC7B,EAAE,EAAE,MAAM,CAAC,gBAAgB;iBAC5B;gBACD,OAAO,EAAE,CAAC;gBACV,YAAY,EAAE,EAAE;gBAChB,OAAO,EAAE,EAAE;gBACX,OAAO,EAAE,iBAAiB;gBAC1B,WAAW,EAAE,EAAE;gBACf,WAAW,EAAE,EAAE;gBACf,MAAM,EAAE,CAAC;gBACT,MAAM,EAAE,CAAC;aACV,CAAC;YACF,IAAI,OAAO,GAAG,EAAE,CAAC;YACjB,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;YACpC,IAAI,OAAO,GAAG,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;YACpC,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,KAAK,IAAI,KAAK,GAAG,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,KAAK,IAAI,CAAC,EAAE,KAAK,EAAE,EAAE;gBACxD,IAAI,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE;oBACxB,MAAM,CAAC,OAAO,EAAE,CAAC;oBACjB,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,KAAK;wBAAE,KAAK,CAAC,KAAK,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC;oBACnD,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;oBAC3B,MAAM,CAAC,MAAM,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,MAAM,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBACjD,IAAI,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBACvB,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;wBAClC,OAAO,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;qBAChC;oBACD,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;wBAClC,OAAO,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;qBAChC;oBACD,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;wBAClC,OAAO,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;qBAChC;oBACD,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,EAAE;wBAClC,OAAO,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC;qBAChC;iBACF;aACF;YAED,MAAM,CAAC,MAAM,GAAG;gBACd,EAAE,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE;gBACpC,EAAE,IAAI,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE;aACrC,CAAC;YACF,MAAM,CAAC,MAAM,IAAI,IAAI,CAAC;YACtB,MAAM,CAAC,MAAM,IAAI,IAAI,CAAC;YACtB,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;YAC9B,MAAM,CAAC,KAAK,GAAG,OAAO,CAAC;YACvB,OAAO,CAAC,IAAI,CAAC,MAAsB,CAAC,CAAC;SACtC;KACF;IACD,OAAO,OAAO,CAAC;AACjB,CAAC"}
|
package/package.json
CHANGED
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
{
|
|
2
2
|
"name": "nmr-processing",
|
|
3
|
-
"version": "8.
|
|
3
|
+
"version": "8.2.0",
|
|
4
4
|
"description": "Pure functions allowing to process NMR spectra.",
|
|
5
5
|
"main": "./lib/index.js",
|
|
6
6
|
"module": "./lib-esm/index.js",
|
|
@@ -60,7 +60,7 @@
|
|
|
60
60
|
"form-data": "^4.0.0",
|
|
61
61
|
"gyromagnetic-ratio": "^1.0.0",
|
|
62
62
|
"is-any-array": "^2.0.0",
|
|
63
|
-
"linear-sum-assignment": "^1.0.
|
|
63
|
+
"linear-sum-assignment": "^1.0.4",
|
|
64
64
|
"ml-array-mean": "^1.1.6",
|
|
65
65
|
"ml-array-rescale": "^1.3.7",
|
|
66
66
|
"ml-array-sequential-fill": "^1.1.8",
|
|
@@ -42,7 +42,7 @@ export function hasIDs(
|
|
|
42
42
|
}
|
|
43
43
|
}
|
|
44
44
|
|
|
45
|
-
export function addIDs
|
|
45
|
+
export function addIDs<T extends NMRRange[] | NMRZone[]>(data: T) {
|
|
46
46
|
for (const element of data) {
|
|
47
47
|
if (!element.id) element.id = generateID();
|
|
48
48
|
for (let signal of element.signals || []) {
|
package/src/index.ts
CHANGED
|
@@ -37,6 +37,9 @@ export * from './databases/getDatabase';
|
|
|
37
37
|
export * from './databases/carbonImpurities';
|
|
38
38
|
export * from './databases/protonImpurities';
|
|
39
39
|
|
|
40
|
+
export * from './peaks/solventSuppression';
|
|
41
|
+
export * from './ranges/markSolventSignal';
|
|
42
|
+
|
|
40
43
|
export type { NMRSignal1D } from './signals/NMRSignal1D';
|
|
41
44
|
export type { NMRSignal2D, Signal2DProjection } from './xyz/NMRSignal2D';
|
|
42
45
|
export type { NMRRange } from './xy/NMRRange';
|
|
@@ -1,5 +1,5 @@
|
|
|
1
1
|
import { PointXY } from 'cheminfo-types';
|
|
2
|
-
import linearSumAssignment from 'linear-sum-assignment';
|
|
2
|
+
import { linearSumAssignment } from 'linear-sum-assignment';
|
|
3
3
|
import { Matrix } from 'ml-matrix';
|
|
4
4
|
import { gaussianFct } from 'ml-peak-shape-generator';
|
|
5
5
|
import { xFindClosestIndex, xMaxValue } from 'ml-spectra-processing';
|
|
@@ -12,45 +12,47 @@ import { splitSpinSystem } from '../signals/simulation/splitSpinSystem';
|
|
|
12
12
|
|
|
13
13
|
import { NMRPeak1D } from './NMRPeak1D';
|
|
14
14
|
|
|
15
|
-
export function solventSuppression(
|
|
16
|
-
peakList:
|
|
15
|
+
export function solventSuppression<T extends NMRPeak1D>(
|
|
16
|
+
peakList: T[],
|
|
17
17
|
solvent: NMRSignal1D[],
|
|
18
18
|
options: any = {},
|
|
19
19
|
) {
|
|
20
|
-
const peaks = [...peakList];
|
|
21
|
-
sortAscending(peaks);
|
|
22
|
-
const xValues = peaks.map((peak) => peak.x);
|
|
20
|
+
const peaks = [...peakList].sort((a, b) => a.x - b.x);
|
|
23
21
|
|
|
24
|
-
const
|
|
22
|
+
const xValues = peaks.map((peak) => peak.x);
|
|
25
23
|
|
|
26
|
-
|
|
27
|
-
// solventSignals = solventSignals.slice(0, solventSignals.findIndex((signal) => signal.delta > 400));
|
|
24
|
+
const { markSolventPeaks = false, solventZoneExtension = 1.2 } = options;
|
|
28
25
|
|
|
29
|
-
for (const solventSignal of
|
|
30
|
-
|
|
26
|
+
for (const solventSignal of solvent) {
|
|
27
|
+
let solventXYPeaks = solventSignal.peaks
|
|
31
28
|
? solventSignal.peaks
|
|
32
29
|
: getSolventPeaks(solventSignal);
|
|
33
|
-
|
|
30
|
+
solventXYPeaks.sort((a, b) => a.x - b.x);
|
|
31
|
+
|
|
34
32
|
let upIndex = xFindClosestIndex(
|
|
35
33
|
xValues,
|
|
36
|
-
solventXYPeaks[solventXYPeaks.length - 1].x +
|
|
34
|
+
solventXYPeaks[solventXYPeaks.length - 1].x + solventZoneExtension,
|
|
37
35
|
);
|
|
38
|
-
let lowIndex = xFindClosestIndex(
|
|
36
|
+
let lowIndex = xFindClosestIndex(
|
|
37
|
+
xValues,
|
|
38
|
+
solventXYPeaks[0].x - solventZoneExtension,
|
|
39
|
+
);
|
|
40
|
+
|
|
39
41
|
if (upIndex === lowIndex) continue;
|
|
40
42
|
const nearPeaks = peaks.slice(lowIndex, upIndex + 1);
|
|
43
|
+
|
|
41
44
|
const amplitudeResiduals = [];
|
|
42
45
|
const deltaResiduals = [];
|
|
43
46
|
const positionResiduals = [];
|
|
44
47
|
for (let peak of nearPeaks) {
|
|
45
|
-
const shiftedSolventPeaks =
|
|
46
|
-
peak,
|
|
47
|
-
|
|
48
|
-
solventXYPeaks,
|
|
49
|
-
);
|
|
48
|
+
const { peaks: shiftedSolventPeaks, delta: currentDelta } =
|
|
49
|
+
getShiftedSolventPeaks(peak, solventSignal, solventXYPeaks);
|
|
50
|
+
|
|
50
51
|
const closestPeaks = getClosestPeaks(shiftedSolventPeaks, nearPeaks);
|
|
51
52
|
let deltaResidual = 0;
|
|
52
53
|
let amplitudeResidual = 0;
|
|
53
54
|
let positionResidual = 0;
|
|
55
|
+
|
|
54
56
|
for (let i = 0; i < closestPeaks.length; i++) {
|
|
55
57
|
amplitudeResidual += Math.abs(
|
|
56
58
|
shiftedSolventPeaks[i].y - closestPeaks[i].y,
|
|
@@ -63,8 +65,12 @@ export function solventSuppression(
|
|
|
63
65
|
amplitudeResidual = Number.MAX_SAFE_INTEGER;
|
|
64
66
|
positionResidual = Number.MAX_SAFE_INTEGER;
|
|
65
67
|
} else {
|
|
66
|
-
positionResidual = gaussianFct(
|
|
68
|
+
positionResidual = gaussianFct(
|
|
69
|
+
Math.abs(solventSignal.delta - currentDelta),
|
|
70
|
+
0.5,
|
|
71
|
+
);
|
|
67
72
|
}
|
|
73
|
+
|
|
68
74
|
amplitudeResiduals.push(amplitudeResidual);
|
|
69
75
|
deltaResiduals.push(deltaResidual);
|
|
70
76
|
positionResiduals.push(positionResidual);
|
|
@@ -76,18 +82,28 @@ export function solventSuppression(
|
|
|
76
82
|
|
|
77
83
|
let minIndex = -1;
|
|
78
84
|
let minScore = Number.MAX_SAFE_INTEGER;
|
|
85
|
+
let score = [];
|
|
79
86
|
for (let i = 0; i < deltaResiduals.length; i++) {
|
|
80
87
|
const value =
|
|
81
88
|
(amplitudeResiduals[i] / maxAmplitude +
|
|
82
89
|
deltaResiduals[i] / maxDelta +
|
|
83
|
-
|
|
90
|
+
1 -
|
|
91
|
+
positionResiduals[i] / maxPosition) /
|
|
84
92
|
3;
|
|
93
|
+
|
|
94
|
+
score.push(value);
|
|
85
95
|
if (minScore > value) {
|
|
86
96
|
minIndex = i;
|
|
87
97
|
minScore = value;
|
|
88
98
|
}
|
|
89
99
|
}
|
|
90
|
-
|
|
100
|
+
|
|
101
|
+
if (minScore < 0) {
|
|
102
|
+
new Error('There is not a correct match with the pattern');
|
|
103
|
+
return peaks;
|
|
104
|
+
}
|
|
105
|
+
|
|
106
|
+
const { peaks: shiftedSolventPeaks } = getShiftedSolventPeaks(
|
|
91
107
|
nearPeaks[minIndex],
|
|
92
108
|
solventSignal,
|
|
93
109
|
solventXYPeaks,
|
|
@@ -96,7 +112,11 @@ export function solventSuppression(
|
|
|
96
112
|
const { rowAssignments, gain } = linearSumAssignment(diff, {
|
|
97
113
|
maximaze: false,
|
|
98
114
|
});
|
|
99
|
-
|
|
115
|
+
|
|
116
|
+
if (gain < 0) {
|
|
117
|
+
new Error('The gain is below to zero');
|
|
118
|
+
return peaks;
|
|
119
|
+
}
|
|
100
120
|
|
|
101
121
|
if (markSolventPeaks) {
|
|
102
122
|
for (let index of rowAssignments) {
|
|
@@ -131,27 +151,34 @@ function getSolventPeaks(
|
|
|
131
151
|
return peaks.filter((peak) => peak.x < 1000);
|
|
132
152
|
}
|
|
133
153
|
|
|
134
|
-
function getShiftedSolventPeaks(
|
|
135
|
-
peak:
|
|
154
|
+
function getShiftedSolventPeaks<T extends PointXY>(
|
|
155
|
+
peak: T,
|
|
136
156
|
solventSignal: any,
|
|
137
|
-
solventXYPeaks:
|
|
157
|
+
solventXYPeaks: T[],
|
|
138
158
|
) {
|
|
139
159
|
const shiftedSolventPeaks: PointXY[] = JSON.parse(
|
|
140
160
|
JSON.stringify(solventXYPeaks),
|
|
141
161
|
);
|
|
142
162
|
// shift x values of solventPeaks to center it to the current peak.
|
|
143
|
-
|
|
163
|
+
let deltaPPM = peak.x - solventSignal.delta;
|
|
164
|
+
if (solventXYPeaks.length > 1 && solventXYPeaks.length % 2 === 0) {
|
|
165
|
+
deltaPPM += solventXYPeaks[0].x;
|
|
166
|
+
}
|
|
144
167
|
const maxIntensity = shiftedSolventPeaks.reduce(
|
|
145
168
|
(max, current) => (current.y > max ? current.y : max),
|
|
146
169
|
shiftedSolventPeaks[0].y,
|
|
147
170
|
);
|
|
171
|
+
let currentDelta = 0;
|
|
148
172
|
for (let shiftedSolventPeak of shiftedSolventPeaks) {
|
|
149
173
|
shiftedSolventPeak.x += deltaPPM;
|
|
150
174
|
shiftedSolventPeak.y /= maxIntensity;
|
|
175
|
+
currentDelta += shiftedSolventPeak.x;
|
|
151
176
|
}
|
|
152
|
-
return
|
|
177
|
+
return {
|
|
178
|
+
peaks: shiftedSolventPeaks,
|
|
179
|
+
delta: currentDelta / shiftedSolventPeaks.length,
|
|
180
|
+
};
|
|
153
181
|
}
|
|
154
|
-
|
|
155
182
|
function getClosestPeaks(shiftedSolventPeaks: PointXY[], nearPeaks: PointXY[]) {
|
|
156
183
|
const diff = getDiffMatrix(shiftedSolventPeaks, nearPeaks);
|
|
157
184
|
const { rowAssignments, gain } = linearSumAssignment(diff, {
|
|
@@ -169,11 +196,7 @@ function getClosestPeaks(shiftedSolventPeaks: PointXY[], nearPeaks: PointXY[]) {
|
|
|
169
196
|
return assignmentPeaks;
|
|
170
197
|
}
|
|
171
198
|
|
|
172
|
-
function
|
|
173
|
-
return peaks[0].x > peaks[1].x ? peaks.sort((a, b) => a.x - b.x) : peaks;
|
|
174
|
-
}
|
|
175
|
-
|
|
176
|
-
function getDiffMatrix(rows: PointXY[], columns: PointXY[]) {
|
|
199
|
+
function getDiffMatrix<T extends PointXY>(rows: T[], columns: T[]) {
|
|
177
200
|
const nbColumns = columns.length;
|
|
178
201
|
const nbRows = rows.length;
|
|
179
202
|
const diff = new Matrix(nbRows, nbColumns);
|
|
@@ -62,10 +62,12 @@ type SignalInternMandatory = MakeMandatory<
|
|
|
62
62
|
function checkSignalForCompilePattern(
|
|
63
63
|
signal: SignalIntern,
|
|
64
64
|
): asserts signal is SignalInternMandatory {
|
|
65
|
-
if (!signal
|
|
65
|
+
if (!('symRank' in signal)) {
|
|
66
66
|
throw new Error('Internal error, symRank was not calculated');
|
|
67
67
|
}
|
|
68
|
-
if (!
|
|
68
|
+
if (!('mask' in signal)) {
|
|
69
|
+
throw new Error('Internal Error, mask was not added');
|
|
70
|
+
}
|
|
69
71
|
}
|
|
70
72
|
|
|
71
73
|
export interface OpitonsCompilePatter {
|