From e17391ca9f7ec9643a68593abc38018051985a8e Mon Sep 17 00:00:00 2001 From: jjangddu Date: Wed, 9 Sep 2026 10:14:52 +0900 Subject: [PATCH] =?UTF-8?q?feat(clinical):=20=EB=B3=91=EC=9B=90=20?= =?UTF-8?q?=EC=9E=84=EC=83=81=EB=A7=8C=20=EC=B5=9C=EC=A2=85=20=EC=95=8C?= =?UTF-8?q?=EA=B3=A0=EB=A6=AC=EC=A6=98=EC=9C=BC=EB=A1=9C=20=C2=B7=20?= =?UTF-8?q?=EC=A0=84=EC=97=AD=EC=9D=80=20=EA=B1=B4=EB=93=9C=EB=A6=AC?= =?UTF-8?q?=EC=A7=80=20=EC=95=8A=EB=8A=94=EB=8B=A4?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit 병원 임상 화면은 최종 보고본(레퍼런스)이 기본이어야 하는데, AlgoMode.reference 는 프로세스 전역이라 그냥 켜면 **일반 측정 화면의 BV 까지** 바뀐다. 그쪽은 현행 임상에서 쓰이고 있어 임의로 건드리면 안 된다. AlgoMode.withReference(ref) { } 범위 실행을 두고 ClinicalBv 가 그 안에서만 돈다. 계산이 끝나면 전역값은 원래대로다 — 임상 화면을 다녀왔다는 이유로 다른 화면의 값이 달라지지 않는다. 되돌리는지도 테스트로 고정했다. 임상 측정은 코루틴에서 돌아 두 계산이 겹칠 수 있으므로 동기화한다: 겹치면 한쪽의 복원이 다른 쪽의 설정을 지워 엉뚱한 경로로 계산된 값이 나오는데, 그 한 건이 검증 기록에 섞이면 나중에 찾아낼 방법이 없다. ClinicalBv.compute(reference = true) 가 기본. BV 측정 화면의 스위치는 이제 전역을 쓰지 않고 그 화면의 계산에만 적용된다(두 경로 비교용). **앞 커밋(21b0fca)의 숫자를 정정한다.** "align_cm1 기존 125.52 vs 레퍼런스 122.25 mL" 는 잘못된 측정이었다 — PiezoHW.activePreset 을 안 잡고 돌려 dps·delay 가 다른 프리셋 값으로 계산됐다(실측 자원은 VBT26050202=V1). 프리셋을 고정하면 이 데이터에서는 두 경로가 **같은 값**을 낸다(cm0 410.65 / cm1 490.63, 양쪽 동일). 두 경로가 같다는 뜻은 아니다. AlgoModeSwitchTest 실측으로 **cycle 44개 중 38개**의 검출·BV 가 갈린다. 여러 cycle 을 평균한 mean-scan 에서 차이가 묻힌 것이다. 즉 입력에 따라 같기도 다르기도 하므로, 값에 경로 표시를 붙이는 이유는 그대로 유효하다. 같은 함정을 테스트에도 반영했다 — ClinicalBvTest·AlgoPathReportTest 가 프리셋을 V1 로 고정한다. 다른 테스트가 남긴 전역 프리셋에 끌려가면 값이 통째로 흔들린다. Co-Authored-By: Claude Opus 5 --- .../medithings/vesiscan/managers/AlgoMode.kt | 29 +++++- .../vesiscan/managers/ClinicalBv.kt | 30 +++++- .../ui/views/clinical/BvMeasureSection.kt | 20 ++-- .../vesiscan/managers/AlgoPathReportTest.kt | 97 +++++++++++++------ .../vesiscan/managers/ClinicalBvTest.kt | 63 ++++++++---- 5 files changed, 177 insertions(+), 62 deletions(-) diff --git a/app/src/main/java/com/medithings/vesiscan/managers/AlgoMode.kt b/app/src/main/java/com/medithings/vesiscan/managers/AlgoMode.kt index 51cf772..a7a043a 100644 --- a/app/src/main/java/com/medithings/vesiscan/managers/AlgoMode.kt +++ b/app/src/main/java/com/medithings/vesiscan/managers/AlgoMode.kt @@ -47,5 +47,32 @@ object AlgoMode { var reference: Boolean = false /** 화면에 표시할 이름. */ - val label: String get() = if (reference) "레퍼런스 (Python 1:1)" else "기존 (현행 임상)" + val label: String get() = labelOf(reference) + + /** 전역 상태가 아니라 **그 계산에 실제로 쓰인** 값의 이름. */ + fun labelOf(ref: Boolean): String = + if (ref) "레퍼런스 (Python 1:1)" else "기존 (현행 임상)" + + private val lock = Any() + + /** + * [ref] 경로로 [block] 을 실행하고 **원래 값으로 되돌린다.** + * + * ## 왜 전역 토글을 그냥 켜지 않나 + * [reference] 는 프로세스 전역이라, 병원 임상 화면이 켜면 **일반 측정 화면의 BV 까지** + * 바뀐다. 그쪽은 현행 임상에서 쓰이고 있어 임의로 바꾸면 안 된다. + * + * 그래서 임상 경로([ClinicalBv])만 이 범위 안에서 돌린다. 계산이 끝나면 전역값은 + * 원래대로다 — 임상 화면을 다녀왔다는 이유로 다른 화면의 값이 달라지지 않는다. + * + * ## 동기화하는 이유 + * 임상 화면의 측정은 코루틴에서 돈다. 두 계산이 겹치면 한쪽의 복원이 다른 쪽의 + * 설정을 지워 **엉뚱한 경로로 계산된 값**이 나온다. 겹칠 일이 드물어도, 그때 나온 + * 한 건이 검증 기록에 섞이면 나중에 찾아낼 방법이 없다. + */ + fun withReference(ref: Boolean, block: () -> T): T = synchronized(lock) { + val prev = reference + reference = ref + try { block() } finally { reference = prev } + } } diff --git a/app/src/main/java/com/medithings/vesiscan/managers/ClinicalBv.kt b/app/src/main/java/com/medithings/vesiscan/managers/ClinicalBv.kt index 06bd557..832e395 100644 --- a/app/src/main/java/com/medithings/vesiscan/managers/ClinicalBv.kt +++ b/app/src/main/java/com/medithings/vesiscan/managers/ClinicalBv.kt @@ -56,9 +56,13 @@ object ClinicalBv { /** * 이 값을 낸 알고리즘 경로. **숫자와 함께 반드시 표시한다.** * - * 두 경로가 같은 입력에 다른 값을 낸다(실측 align_cm1: 기존 125.52 mL vs - * 레퍼런스 122.25 mL). 어느 쪽에서 나온 값인지 모르면 "정렬 위치의 BV 가 - * 타당한가"를 판단할 수 없다 — 검증 자체가 성립하지 않는다. + * 두 경로는 같은 신호에 다른 값을 낼 수 있다 — `AlgoModeSwitchTest` 실측에서 + * **cycle 44개 중 38개**의 검출·BV 가 갈렸다. 다만 여러 cycle 을 평균한 + * mean-scan 에서는 차이가 묻히기도 한다(align_cm0·cm1 은 두 경로가 동일). + * + * 즉 **입력에 따라 같기도 다르기도 하다.** 그래서 값만 남기면 나중에 해석할 수 + * 없다 — 어느 경로에서 나온 BV 인지 모르면 "정렬 위치가 타당한가"의 근거가 + * 되지 못한다. */ val algoLabel: String, /** 파형 마커용 (전벽, 후벽) 인덱스. 미검출 채널은 null. */ @@ -76,8 +80,26 @@ object ClinicalBv { signals: List, params: MethodDParams = MethodDParams.DEFAULT, supine: Boolean = false, + /** + * 병원 임상은 **최종 보고본(레퍼런스)** 이 기본이다. + * + * 전역 [AlgoMode.reference] 를 켜는 대신 [AlgoMode.withReference] 범위 안에서만 + * 돌린다 — 전역을 켜면 일반 측정 화면의 BV 까지 바뀌는데, 그쪽은 현행 임상에서 + * 쓰이고 있어 임의로 건드리면 안 된다. + * + * false 는 비교용이다(두 경로가 얼마나 다른지 그 자리에서 보기 위함). + */ + reference: Boolean = true, + ): Outcome = AlgoMode.withReference(reference) { + computeInScope(signals, params, supine, AlgoMode.labelOf(reference)) + } + + private fun computeInScope( + signals: List, + params: MethodDParams, + supine: Boolean, + algo: String, ): Outcome { - val algo = AlgoMode.label if (signals.size < PiezoHW.centerCh.size) { return Outcome(null, emptyList(), 0, "채널 부족 (${signals.size}/6)", algo, emptyList(), null) } diff --git a/app/src/main/java/com/medithings/vesiscan/ui/views/clinical/BvMeasureSection.kt b/app/src/main/java/com/medithings/vesiscan/ui/views/clinical/BvMeasureSection.kt index 4cab2cb..4feb9ba 100644 --- a/app/src/main/java/com/medithings/vesiscan/ui/views/clinical/BvMeasureSection.kt +++ b/app/src/main/java/com/medithings/vesiscan/ui/views/clinical/BvMeasureSection.kt @@ -16,7 +16,6 @@ import androidx.compose.ui.unit.dp import androidx.compose.ui.unit.sp import com.medithings.vesiscan.ble.BleManager import com.medithings.vesiscan.ble.PiezoChannelData -import com.medithings.vesiscan.managers.AlgoMode import com.medithings.vesiscan.managers.ClinicalBv import com.medithings.vesiscan.models.HospitalFixedParams import com.medithings.vesiscan.models.ProbeCycle @@ -79,13 +78,10 @@ fun BvMeasureSection( var error by remember { mutableStateOf(null) } var history by remember { mutableStateOf>(emptyList()) } - // 검증 목적이 **최종 보고된 알고리즘**이라 이 화면에서는 레퍼런스 경로를 기본으로 - // 켠다. AlgoMode 는 전역이라 일반 측정 화면에도 영향을 주므로, 켠 사실을 화면에 - // 드러내고 되돌릴 수 있게 둔다. - var useReference by remember { mutableStateOf(AlgoMode.reference) } - LaunchedEffect(Unit) { - if (!AlgoMode.reference) { AlgoMode.reference = true; useReference = true } - } + // 병원 임상은 **최종 보고본이 기본**이다. 전역 AlgoMode 를 켜지 않는다 — + // ClinicalBv 가 계산 동안만 그 경로로 돌리고 되돌리므로, 일반 측정 화면의 BV 는 + // 영향을 받지 않는다. 끄는 것은 두 경로를 그 자리에서 비교하기 위한 용도다. + var useReference by remember { mutableStateOf(true) } LaunchedEffect(mode) { val m = mode ?: return@LaunchedEffect @@ -142,7 +138,8 @@ fun BvMeasureSection( } if (!aborted && cycles.size == AVG_CYCLES) { error = null - val r = ClinicalBv.compute(ClinicalBv.meanScan(cycles), supine = supine) + val r = ClinicalBv.compute( + ClinicalBv.meanScan(cycles), supine = supine, reference = useReference) outcome = r r.volumeMl?.let { v -> history = (history + v).takeLast(HISTORY_MAX) @@ -173,7 +170,7 @@ fun BvMeasureSection( Row(verticalAlignment = Alignment.CenterVertically) { Switch( checked = useReference, - onCheckedChange = { useReference = it; AlgoMode.reference = it }, + onCheckedChange = { useReference = it }, enabled = !running, ) Spacer(Modifier.width(8.dp)) @@ -183,7 +180,8 @@ fun BvMeasureSection( fontSize = 13.sp, fontWeight = FontWeight.SemiBold, ) Text( - "끄면 기존 경로로 계산합니다. 두 경로는 같은 신호에 다른 값을 낼 수 있습니다.", + "이 화면에만 적용됩니다. 끄면 기존 경로로 계산합니다 — " + + "두 경로는 같은 신호에 다른 값을 낼 수 있습니다.", fontSize = 10.sp, color = MlSecondaryText, ) } diff --git a/app/src/test/java/com/medithings/vesiscan/managers/AlgoPathReportTest.kt b/app/src/test/java/com/medithings/vesiscan/managers/AlgoPathReportTest.kt index 6ab5e93..eef777c 100644 --- a/app/src/test/java/com/medithings/vesiscan/managers/AlgoPathReportTest.kt +++ b/app/src/test/java/com/medithings/vesiscan/managers/AlgoPathReportTest.kt @@ -1,6 +1,7 @@ package com.medithings.vesiscan.managers import org.junit.After +import org.junit.Before import org.junit.Test import java.io.File @@ -8,20 +9,31 @@ import java.io.File * 두 알고리즘 경로가 같은 입력에 어떤 값을 내는지 **숫자로 남긴다**. * * 검증 목적이 "정렬이 고른 위치의 BV 가 타당한가"인데, 그 BV 가 어느 경로에서 나온 - * 값인지 모르면 판단이 성립하지 않는다. [AlgoMode.reference] 는 기본이 **false**라 - * 앱은 기존 경로로 돈다 — 이 시험은 그 차이를 실측 데이터로 찍어 보고서에 쓴다. + * 값인지 모르면 판단이 성립하지 않는다. 병원 임상 경로([ClinicalBv])는 최종 보고본 + * (레퍼런스)이 기본이고, 일반 측정 화면은 기존 경로 그대로다 — 이 시험은 그 차이를 + * 실측 데이터로 찍어 보고서에 쓴다. * * 통과·실패를 가리는 시험이 아니다(단정문이 없다). 값을 출력하는 것이 목적이라 * 리포트에 그대로 옮길 수 있다. */ class AlgoPathReportTest { - @After fun restore() { AlgoMode.reference = false } + // 실측 자원은 VBT26050202(=V1) 로 받은 것이다. 프리셋이 다르면 dps·delay 가 달라져 + // 같은 sample index 가 다른 깊이로 매핑된다 — 프리셋을 안 잡으면 Python 과 어긋난 + // 값이 나오고, 그걸 "불일치"로 오해하게 된다(AnchorMeasureParityTest 와 같은 이유). + private val savedPreset = PiezoHW.activePreset - private fun scan(path: String): List? { + @Before fun setUp() { PiezoHW.activePreset = PiezoHW.DevicePreset.V1 } + + @After fun restore() { + AlgoMode.reference = false + PiezoHW.activePreset = savedPreset + } + + /** `cycle,ch,samples…` 를 cycle 목록으로. */ + private fun cyclesOf(path: String): List> { val f = File(path) - if (!f.exists()) return null - // cycle,ch,samples… — 전 cycle 을 채널별 평균(mean-scan)해 정렬과 같은 입력으로. + if (!f.exists()) return emptyList() val byCycle = linkedMapOf>() f.forEachLine { line -> val p = line.trim().split(',') @@ -31,32 +43,61 @@ class AlgoPathReportTest { byCycle.getOrPut(ci) { linkedMapOf() }[ch] = DoubleArray(p.size - 2) { p[it + 2].toDoubleOrNull() ?: 0.0 } } - val cycles = byCycle.values.filter { it.size >= 6 }.map { m -> (0..5).map { m.getValue(it) } } - return if (cycles.isEmpty()) null else ClinicalBv.meanScan(cycles) - } - - private fun report(label: String, sig: List) { - val out = StringBuilder("\n=== $label ===\n") - for (ref in listOf(false, true)) { - AlgoMode.reference = ref - val r = ClinicalBv.compute(sig) - out.append(" %-22s BV=%s 검출=%d/4 %s\n".format( - if (ref) "레퍼런스(Python 1:1)" else "기존(현행 임상)", - r.volumeMl?.let { "%.2f mL".format(it) } ?: "계산불가", - r.detectedCenter, - r.walls.filterNotNull().joinToString(" ") { "CH${it.channel}(${it.ant}~${it.post})" }, - )) - } - println(out) + return byCycle.values.filter { it.size >= 6 }.map { m -> (0..5).map { m.getValue(it) } } } @Test fun `두 경로의 BV 를 실측 데이터로 비교`() { - var any = false for (name in listOf("align_cm0.csv", "align_cm1.csv")) { - val sig = scan("src/test/resources/$name") ?: continue - any = true - report(name, sig) + val cycles = cyclesOf("src/test/resources/$name") + if (cycles.isEmpty()) continue + val sig = ClinicalBv.meanScan(cycles) + val sb = StringBuilder("\n=== BV · $name ===\n") + for (ref in listOf(false, true)) { + // ClinicalBv 가 범위 안에서 경로를 잡으므로 전역을 건드리지 않는다. + val r = ClinicalBv.compute(sig, reference = ref) + sb.append( + " %-22s BV=%-12s 검출=%d/4 %s\n".format( + if (ref) "레퍼런스(최종)" else "기존(현행)", + r.volumeMl?.let { "%.2f mL".format(it) } ?: "계산불가", + r.detectedCenter, + r.walls.filterNotNull() + .joinToString(" ") { "CH${it.channel}(${it.ant}~${it.post})" }, + ) + ) + } + println(sb) + } + // 범위 실행이 전역을 되돌렸는지 — 되돌리지 않으면 일반 측정 화면이 오염된다. + println("\n전역 AlgoMode.reference = ${AlgoMode.reference} (false 여야 한다)") + } + + @Test fun `정렬 지표를 두 경로로 비교`() { + val expect = File("src/test/resources/align_expect.tsv") + if (!expect.exists()) { println("\n(기대값 파일 없음)"); return } + val want = expect.readLines() + .filter { it.isNotBlank() && !it.startsWith("#") } + .associate { l -> l.split("\t").let { it[0].trim().toInt() to it } } + + for (ref in listOf(false, true)) { + AlgoMode.reference = ref + val sb = StringBuilder( + "\n=== AnchorGuide · ${if (ref) "레퍼런스(최종)" else "기존(현행)"} ===\n") + for (cm in 0..1) { + val cycles = cyclesOf("src/test/resources/align_cm$cm.csv") + if (cycles.isEmpty()) continue + val rec = AnchorGuide().measure(cycles, cm) + val w = want[cm] + sb.append( + " cm=%d nch=%d(py %s) ch3=%s(py %s) rate=%.10f(py %s) cap=%.10f(py %s)\n" + .format( + cm, rec.nch, w?.getOrNull(2), + if (rec.ch3Mean) "O" else "X", w?.getOrNull(3), + rec.ch3Rate, w?.getOrNull(4), + rec.capFrac, w?.getOrNull(5), + ) + ) + } + println(sb) } - if (!any) println("\n(실측 자원 없음 — 비교 생략)") } } diff --git a/app/src/test/java/com/medithings/vesiscan/managers/ClinicalBvTest.kt b/app/src/test/java/com/medithings/vesiscan/managers/ClinicalBvTest.kt index d6f53d7..25c055d 100644 --- a/app/src/test/java/com/medithings/vesiscan/managers/ClinicalBvTest.kt +++ b/app/src/test/java/com/medithings/vesiscan/managers/ClinicalBvTest.kt @@ -23,6 +23,15 @@ import java.io.File */ class ClinicalBvTest { + // 실측 자원은 VBT26050202(=V1) 로 받은 것이다. 프리셋이 다르면 dps·delay 가 달라져 + // 같은 sample index 가 다른 깊이로 매핑된다 — 다른 테스트가 남긴 전역 프리셋에 + // 끌려가면 여기 값이 통째로 흔들린다. + private val savedPreset = PiezoHW.activePreset + + @org.junit.Before fun pinPreset() { PiezoHW.activePreset = PiezoHW.DevicePreset.V1 } + + @org.junit.After fun restorePreset() { PiezoHW.activePreset = savedPreset } + /** * 정렬 검증용 실측 trace (`cycle,ch,samples…`). 파일이 없으면 **실패시킨다** — * 조용히 건너뛰면 이 클래스에서 제일 중요한 parity 검사가 통과한 것처럼 보인다. @@ -46,34 +55,52 @@ class ClinicalBvTest { private fun flat(n: Int = 120) = List(6) { DoubleArray(n) { 1000.0 } } + /** + * 두 경로 모두에서 대조한다. + * + * [ClinicalBv.compute] 는 병원 임상용이라 **레퍼런스가 기본**이다. 기대값을 기존 + * 경로로 만들어 비교하면 갈린다 — 실제로 갈린다(AlgoModeSwitchTest: cycle 44개 중 + * 38개). 그건 버그가 아니라 두 경로가 다르다는 사실이므로, 시험은 **같은 경로끼리** + * 맞춰 비교해야 의미가 있다. + */ @Test fun `AnchorGuide 의 cap_frac 과 같은 경로로 계산한다`() { val scan = realScan() - - // AnchorGuide.measure 가 하는 것과 **같은 호출**을 손으로 재현. - val outcome = MethodDRunner.detectMultichannelWithEvidence( - scan, MethodDParams.DEFAULT, applyTgc = true, applyCross = true) - val spans = outcome.walls.map { r -> - r?.let { WallWithSpan(it.antRefined, it.postRefined, it.lowStart, it.lowEnd) } + for (ref in listOf(false, true)) { + // AnchorGuide.measure 가 하는 것과 **같은 호출**을 손으로 재현. + val expected = AlgoMode.withReference(ref) { + val outcome = MethodDRunner.detectMultichannelWithEvidence( + scan, MethodDParams.DEFAULT, applyTgc = true, applyCross = true) + val spans = outcome.walls.map { r -> + r?.let { WallWithSpan(it.antRefined, it.postRefined, it.lowStart, it.lowEnd) } + } + estimateBv(spans, missingOutside = outcome.missingOutside) + } + val got = ClinicalBv.compute(scan, reference = ref) + if (expected == null) { + assertNull("ref=$ref 레퍼런스 경로가 null 인데 값을 냈다", got.volumeMl) + } else { + assertNotNull("ref=$ref 레퍼런스 경로는 값을 냈는데 null", got.volumeMl) + assertEquals("ref=$ref volumeMl", expected.volumeMl, got.volumeMl!!, 1e-9) + } } - val expected = estimateBv(spans, missingOutside = outcome.missingOutside) + } - val got = ClinicalBv.compute(scan) - - if (expected == null) { - assertNull("레퍼런스 경로가 null 인데 ClinicalBv 가 값을 냈다", got.volumeMl) - } else { - assertNotNull("레퍼런스 경로는 값을 냈는데 ClinicalBv 가 null", got.volumeMl) - assertEquals(expected.volumeMl, got.volumeMl!!, 1e-9) - } + @Test fun `범위 실행이 전역 상태를 되돌린다`() { + // 되돌리지 않으면 임상 화면을 다녀왔다는 이유로 **일반 측정 화면의 BV** 가 바뀐다. + val before = AlgoMode.reference + ClinicalBv.compute(realScan(), reference = true) + assertEquals("compute 뒤 전역값이 바뀌었다", before, AlgoMode.reference) } @Test fun `마커는 검출 인덱스를 그대로 준다`() { // 파형 위 세로선이 실제 검출 위치와 어긋나면, 눈으로 보는 것과 계산이 달라져 // 원인 분석이 성립하지 않는다. val scan = realScan() - val outcome = MethodDRunner.detectMultichannelWithEvidence( - scan, MethodDParams.DEFAULT, applyTgc = true, applyCross = true) - val got = ClinicalBv.compute(scan) + val outcome = AlgoMode.withReference(true) { + MethodDRunner.detectMultichannelWithEvidence( + scan, MethodDParams.DEFAULT, applyTgc = true, applyCross = true) + } + val got = ClinicalBv.compute(scan) // 기본 = 레퍼런스 assertEquals(outcome.walls.size, got.markers.size) outcome.walls.forEachIndexed { i, r ->