Compare commits
1 Commits
fix/36
...
fix/131-vo
| Author | SHA1 | Date | |
|---|---|---|---|
| dd38f362d6 |
@@ -1,195 +1 @@
|
||||
"""Crisis synthesizer — learn from anonymized crisis interactions.
|
||||
|
||||
This is deliberately simple and privacy-preserving. It does not train a model or
|
||||
modify detection rules automatically. It only logs metadata, summarizes patterns,
|
||||
and suggests human-reviewed keyword weight adjustments.
|
||||
"""
|
||||
|
||||
from __future__ import annotations
|
||||
|
||||
import argparse
|
||||
import json
|
||||
import time
|
||||
from collections import Counter, defaultdict
|
||||
from pathlib import Path
|
||||
from typing import Iterable
|
||||
|
||||
DEFAULT_LOG_PATH = Path.home() / ".the-door" / "crisis-interactions.jsonl"
|
||||
LEVELS = ("NONE", "LOW", "MEDIUM", "HIGH", "CRITICAL")
|
||||
|
||||
|
||||
def build_interaction_event(
|
||||
level: str,
|
||||
indicators: list[str],
|
||||
response_given: str,
|
||||
continued_conversation: bool,
|
||||
false_positive: bool,
|
||||
*,
|
||||
now: float | None = None,
|
||||
) -> dict:
|
||||
return {
|
||||
"timestamp": float(time.time() if now is None else now),
|
||||
"level": level,
|
||||
"indicators": list(indicators),
|
||||
"indicator_count": len(indicators),
|
||||
"response_given": response_given,
|
||||
"continued_conversation": bool(continued_conversation),
|
||||
"false_positive": bool(false_positive),
|
||||
}
|
||||
|
||||
|
||||
def append_interaction_event(
|
||||
log_path: str | Path,
|
||||
*,
|
||||
level: str,
|
||||
indicators: list[str],
|
||||
response_given: str,
|
||||
continued_conversation: bool,
|
||||
false_positive: bool,
|
||||
now: float | None = None,
|
||||
) -> dict:
|
||||
event = build_interaction_event(
|
||||
level,
|
||||
indicators,
|
||||
response_given,
|
||||
continued_conversation,
|
||||
false_positive,
|
||||
now=now,
|
||||
)
|
||||
path = Path(log_path)
|
||||
path.parent.mkdir(parents=True, exist_ok=True)
|
||||
with path.open("a", encoding="utf-8") as handle:
|
||||
handle.write(json.dumps(event) + "\n")
|
||||
return event
|
||||
|
||||
|
||||
def load_interaction_events(log_path: str | Path) -> list[dict]:
|
||||
path = Path(log_path)
|
||||
if not path.exists():
|
||||
return []
|
||||
events = []
|
||||
for line in path.read_text(encoding="utf-8").splitlines():
|
||||
if not line.strip():
|
||||
continue
|
||||
events.append(json.loads(line))
|
||||
return events
|
||||
|
||||
|
||||
def summarize_keywords(events: Iterable[dict]) -> list[dict]:
|
||||
counts: Counter[str] = Counter()
|
||||
for event in events:
|
||||
counts.update(event.get("indicators", []))
|
||||
return [{"keyword": keyword, "count": count} for keyword, count in counts.most_common(10)]
|
||||
|
||||
|
||||
def suggest_keyword_adjustments(events: Iterable[dict], *, min_observations: int = 5) -> list[dict]:
|
||||
stats: dict[str, dict[str, int]] = defaultdict(lambda: {
|
||||
"observations": 0,
|
||||
"true_positive_count": 0,
|
||||
"false_positive_count": 0,
|
||||
"continued_conversation_count": 0,
|
||||
})
|
||||
|
||||
for event in events:
|
||||
for keyword in event.get("indicators", []):
|
||||
bucket = stats[keyword]
|
||||
bucket["observations"] += 1
|
||||
if event.get("false_positive"):
|
||||
bucket["false_positive_count"] += 1
|
||||
else:
|
||||
bucket["true_positive_count"] += 1
|
||||
if event.get("continued_conversation"):
|
||||
bucket["continued_conversation_count"] += 1
|
||||
|
||||
suggestions = []
|
||||
for keyword, bucket in sorted(stats.items()):
|
||||
if bucket["observations"] < min_observations:
|
||||
continue
|
||||
fp = bucket["false_positive_count"]
|
||||
tp = bucket["true_positive_count"]
|
||||
if fp >= min_observations and tp == 0:
|
||||
adjustment = "lower_weight"
|
||||
rationale = "Observed only false positives across the sample window."
|
||||
elif tp >= min_observations and fp == 0:
|
||||
adjustment = "raise_weight"
|
||||
rationale = "Observed repeated genuine crises with no false positives."
|
||||
else:
|
||||
adjustment = "observe"
|
||||
rationale = "Mixed evidence; keep monitoring before changing weights."
|
||||
suggestions.append(
|
||||
{
|
||||
"keyword": keyword,
|
||||
**bucket,
|
||||
"suggested_adjustment": adjustment,
|
||||
"rationale": rationale,
|
||||
}
|
||||
)
|
||||
return suggestions
|
||||
|
||||
|
||||
def build_weekly_report(
|
||||
events: Iterable[dict],
|
||||
*,
|
||||
now: float | None = None,
|
||||
window_days: int = 7,
|
||||
min_observations: int = 3,
|
||||
) -> dict:
|
||||
current_time = float(time.time() if now is None else now)
|
||||
cutoff = current_time - (window_days * 86400)
|
||||
filtered = [event for event in events if float(event.get("timestamp", 0)) >= cutoff]
|
||||
|
||||
detections_per_level = {level: 0 for level in LEVELS}
|
||||
detected_events = []
|
||||
continued_after_intervention = 0
|
||||
for event in filtered:
|
||||
level = event.get("level", "NONE")
|
||||
detections_per_level[level] = detections_per_level.get(level, 0) + 1
|
||||
if level != "NONE":
|
||||
detected_events.append(event)
|
||||
if event.get("continued_conversation"):
|
||||
continued_after_intervention += 1
|
||||
|
||||
false_positive_count = sum(1 for event in detected_events if event.get("false_positive"))
|
||||
false_positive_estimate = false_positive_count / len(detected_events) if detected_events else 0.0
|
||||
|
||||
return {
|
||||
"window_days": window_days,
|
||||
"total_events": len(filtered),
|
||||
"detections_per_level": detections_per_level,
|
||||
"most_common_keywords": summarize_keywords(filtered),
|
||||
"false_positive_estimate": false_positive_estimate,
|
||||
"continued_after_intervention": continued_after_intervention,
|
||||
"keyword_weight_suggestions": suggest_keyword_adjustments(filtered, min_observations=min_observations),
|
||||
}
|
||||
|
||||
|
||||
def render_weekly_report(summary: dict) -> str:
|
||||
return json.dumps(summary, indent=2)
|
||||
|
||||
|
||||
def write_weekly_report(output_path: str | Path, summary: dict) -> Path:
|
||||
path = Path(output_path)
|
||||
path.parent.mkdir(parents=True, exist_ok=True)
|
||||
path.write_text(render_weekly_report(summary) + "\n", encoding="utf-8")
|
||||
return path
|
||||
|
||||
|
||||
def main(argv: list[str] | None = None) -> int:
|
||||
parser = argparse.ArgumentParser(description="Summarize anonymized crisis interactions")
|
||||
parser.add_argument("--log-path", default=str(DEFAULT_LOG_PATH), help="JSONL crisis interaction log")
|
||||
parser.add_argument("--days", type=int, default=7, help="Lookback window in days")
|
||||
parser.add_argument("--min-observations", type=int, default=3, help="Minimum observations before suggesting keyword adjustments")
|
||||
parser.add_argument("--output", help="Optional file to write the weekly report JSON")
|
||||
args = parser.parse_args(argv)
|
||||
|
||||
events = load_interaction_events(args.log_path)
|
||||
summary = build_weekly_report(events, window_days=args.days, min_observations=args.min_observations)
|
||||
rendered = render_weekly_report(summary)
|
||||
print(rendered)
|
||||
if args.output:
|
||||
write_weekly_report(args.output, summary)
|
||||
return 0
|
||||
|
||||
|
||||
if __name__ == "__main__":
|
||||
raise SystemExit(main())
|
||||
...
|
||||
@@ -680,7 +680,7 @@ html, body {
|
||||
|
||||
<!-- Footer -->
|
||||
<footer id="footer">
|
||||
<a href="/about.html" aria-label="About The Door">about</a>
|
||||
<a href="/about" aria-label="About The Door">about</a>
|
||||
<button id="safety-plan-btn" aria-label="Open My Safety Plan">my safety plan</button>
|
||||
<button id="clear-chat-btn" aria-label="Clear chat history">clear chat</button>
|
||||
</footer>
|
||||
|
||||
@@ -1,111 +0,0 @@
|
||||
"""Tests for evolution/crisis_synthesizer.py (issue #36)."""
|
||||
|
||||
from __future__ import annotations
|
||||
|
||||
import importlib.util
|
||||
import json
|
||||
import pathlib
|
||||
import sys
|
||||
import tempfile
|
||||
import unittest
|
||||
|
||||
ROOT = pathlib.Path(__file__).resolve().parents[1]
|
||||
SCRIPT = ROOT / 'evolution' / 'crisis_synthesizer.py'
|
||||
|
||||
spec = importlib.util.spec_from_file_location('crisis_synthesizer', str(SCRIPT))
|
||||
mod = importlib.util.module_from_spec(spec)
|
||||
sys.modules['crisis_synthesizer'] = mod
|
||||
spec.loader.exec_module(mod)
|
||||
|
||||
|
||||
class TestCrisisSynthesizerEvent(unittest.TestCase):
|
||||
def test_build_interaction_event_is_privacy_preserving(self):
|
||||
event = mod.build_interaction_event(
|
||||
level='CRITICAL',
|
||||
indicators=['want_to_die', 'no_way_out'],
|
||||
response_given='guardian',
|
||||
continued_conversation=True,
|
||||
false_positive=False,
|
||||
now=1700000000,
|
||||
)
|
||||
self.assertEqual(event['timestamp'], 1700000000)
|
||||
self.assertEqual(event['level'], 'CRITICAL')
|
||||
self.assertEqual(event['response_given'], 'guardian')
|
||||
self.assertTrue(event['continued_conversation'])
|
||||
self.assertFalse(event['false_positive'])
|
||||
self.assertEqual(event['indicators'], ['want_to_die', 'no_way_out'])
|
||||
for forbidden in ['text', 'message', 'content', 'ip', 'session_id', 'user_id']:
|
||||
self.assertNotIn(forbidden, event)
|
||||
|
||||
|
||||
class TestCrisisSynthesizerStorage(unittest.TestCase):
|
||||
def test_append_and_load_events_round_trip(self):
|
||||
with tempfile.TemporaryDirectory() as tmp:
|
||||
log_path = pathlib.Path(tmp) / 'crisis-events.jsonl'
|
||||
mod.append_interaction_event(
|
||||
log_path,
|
||||
level='HIGH',
|
||||
indicators=['hopeless'],
|
||||
response_given='companion',
|
||||
continued_conversation=False,
|
||||
false_positive=True,
|
||||
now=1700000100,
|
||||
)
|
||||
events = mod.load_interaction_events(log_path)
|
||||
self.assertEqual(len(events), 1)
|
||||
self.assertEqual(events[0]['level'], 'HIGH')
|
||||
self.assertEqual(events[0]['indicators'], ['hopeless'])
|
||||
|
||||
|
||||
class TestCrisisSynthesizerSummary(unittest.TestCase):
|
||||
def test_weekly_report_contains_required_metrics(self):
|
||||
events = [
|
||||
mod.build_interaction_event('CRITICAL', ['want_to_die'], 'guardian', True, False, now=1700000000),
|
||||
mod.build_interaction_event('HIGH', ['hopeless'], 'companion', False, True, now=1700000100),
|
||||
mod.build_interaction_event('LOW', ['rough_day'], 'friend', False, False, now=1700000200),
|
||||
mod.build_interaction_event('CRITICAL', ['want_to_die'], 'guardian', False, False, now=1700000300),
|
||||
mod.build_interaction_event('NONE', [], 'friend', False, False, now=1700000400),
|
||||
]
|
||||
summary = mod.build_weekly_report(events, now=1700000500, window_days=7)
|
||||
self.assertEqual(summary['detections_per_level']['CRITICAL'], 2)
|
||||
self.assertEqual(summary['detections_per_level']['HIGH'], 1)
|
||||
self.assertEqual(summary['detections_per_level']['LOW'], 1)
|
||||
self.assertEqual(summary['detections_per_level']['NONE'], 1)
|
||||
self.assertEqual(summary['continued_after_intervention'], 1)
|
||||
self.assertAlmostEqual(summary['false_positive_estimate'], 0.25)
|
||||
self.assertEqual(summary['most_common_keywords'][0]['keyword'], 'want_to_die')
|
||||
self.assertEqual(summary['most_common_keywords'][0]['count'], 2)
|
||||
|
||||
|
||||
class TestCrisisSynthesizerSuggestions(unittest.TestCase):
|
||||
def test_suggests_weight_adjustments_from_interactions(self):
|
||||
events = []
|
||||
for ts in range(3):
|
||||
events.append(mod.build_interaction_event('CRITICAL', ['want_to_die'], 'guardian', True, False, now=1700000000 + ts))
|
||||
for ts in range(3):
|
||||
events.append(mod.build_interaction_event('LOW', ['rough_day'], 'friend', False, True, now=1700000100 + ts))
|
||||
suggestions = mod.suggest_keyword_adjustments(events, min_observations=3)
|
||||
by_keyword = {s['keyword']: s for s in suggestions}
|
||||
self.assertEqual(by_keyword['want_to_die']['suggested_adjustment'], 'raise_weight')
|
||||
self.assertEqual(by_keyword['rough_day']['suggested_adjustment'], 'lower_weight')
|
||||
|
||||
|
||||
class TestCrisisSynthesizerRendering(unittest.TestCase):
|
||||
def test_render_weekly_report_outputs_json(self):
|
||||
summary = {
|
||||
'detections_per_level': {'NONE': 0, 'LOW': 1, 'MEDIUM': 0, 'HIGH': 0, 'CRITICAL': 0},
|
||||
'most_common_keywords': [{'keyword': 'rough_day', 'count': 1}],
|
||||
'false_positive_estimate': 0.0,
|
||||
'continued_after_intervention': 0,
|
||||
'keyword_weight_suggestions': [],
|
||||
'window_days': 7,
|
||||
'total_events': 1,
|
||||
}
|
||||
rendered = mod.render_weekly_report(summary)
|
||||
parsed = json.loads(rendered)
|
||||
self.assertEqual(parsed['window_days'], 7)
|
||||
self.assertEqual(parsed['most_common_keywords'][0]['keyword'], 'rough_day')
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
unittest.main()
|
||||
350
voice_analysis.py
Normal file
350
voice_analysis.py
Normal file
@@ -0,0 +1,350 @@
|
||||
"""
|
||||
voice_analysis.py — Voice message distress analysis via paralinguistic features.
|
||||
|
||||
Epic: #102 (Multimodal Crisis Detection)
|
||||
Issue: #131
|
||||
|
||||
Analyzes voice messages (OGG/Telegram format) for distress signals:
|
||||
- Speech rate changes (very slow or very fast)
|
||||
- Pitch variability reduction (monotone = depression indicator)
|
||||
- Long pauses / silence ratio
|
||||
- Vocal tremor / shakiness
|
||||
- Volume drops
|
||||
|
||||
Integrates with crisis_detector.py text-based detection for multimodal coverage.
|
||||
"""
|
||||
|
||||
import os
|
||||
import json
|
||||
import subprocess
|
||||
import tempfile
|
||||
from dataclasses import dataclass, field, asdict
|
||||
from typing import Optional
|
||||
|
||||
|
||||
@dataclass
|
||||
class VoiceAnalysisResult:
|
||||
"""Result of paralinguistic analysis on a voice message."""
|
||||
transcript: str = ""
|
||||
speech_rate: float = 0.0 # words per minute
|
||||
pitch_mean: float = 0.0 # Hz, average fundamental frequency
|
||||
pitch_variability: float = 0.0 # std dev of pitch (low = monotone)
|
||||
silence_ratio: float = 0.0 # 0-1, fraction of audio that is silence
|
||||
tremor_score: float = 0.0 # 0-1, vocal shakiness estimate
|
||||
volume_drop_score: float = 0.0 # 0-1, sudden volume decreases
|
||||
distress_score: float = 0.0 # 0-1, composite distress indicator
|
||||
signals_detected: list = field(default_factory=list)
|
||||
|
||||
def to_dict(self) -> dict:
|
||||
return asdict(self)
|
||||
|
||||
|
||||
# === THRESHOLDS ===
|
||||
|
||||
# Speech rate: normal is ~120-150 WPM
|
||||
# Very slow (<80) or very fast (>200) are distress indicators
|
||||
SPEECH_RATE_SLOW = 80
|
||||
SPEECH_RATE_FAST = 200
|
||||
SPEECH_RATE_NORMAL_LOW = 100
|
||||
SPEECH_RATE_NORMAL_HIGH = 170
|
||||
|
||||
# Pitch variability: normal conversation has std dev ~30-50 Hz
|
||||
# Monotone (<15 Hz) is a depression indicator
|
||||
PITCH_VARIABILITY_LOW = 15.0 # Hz — monotone threshold
|
||||
PITCH_VARIABILITY_NORMAL = 30.0
|
||||
|
||||
# Silence ratio: normal has ~10-20% silence
|
||||
# Excessive silence (>40%) or very little (<3%) may indicate distress
|
||||
SILENCE_RATIO_HIGH = 0.4
|
||||
SILENCE_RATIO_LOW = 0.03
|
||||
|
||||
# Composite thresholds
|
||||
DISTRESS_LOW = 0.3
|
||||
DISTRESS_MEDIUM = 0.7
|
||||
|
||||
|
||||
# === CORE ANALYSIS ===
|
||||
|
||||
def _convert_to_wav(audio_path: str) -> str:
|
||||
"""Convert audio to WAV format for analysis. Returns path to temp WAV file."""
|
||||
wav_path = tempfile.mktemp(suffix='.wav')
|
||||
try:
|
||||
subprocess.run(
|
||||
['ffmpeg', '-i', audio_path, '-ar', '16000', '-ac', '1', '-y', wav_path],
|
||||
capture_output=True, timeout=30
|
||||
)
|
||||
if not os.path.exists(wav_path):
|
||||
# Fallback: if ffmpeg not available, try the original file
|
||||
return audio_path
|
||||
return wav_path
|
||||
except (FileNotFoundError, subprocess.TimeoutExpired):
|
||||
return audio_path
|
||||
|
||||
|
||||
def _transcribe(audio_path: str) -> str:
|
||||
"""Transcribe audio using whisper (if available) or return empty string."""
|
||||
try:
|
||||
import whisper
|
||||
model = whisper.load_model("base")
|
||||
result = model.transcribe(audio_path)
|
||||
return result.get("text", "").strip()
|
||||
except ImportError:
|
||||
# Whisper not available — skip transcription
|
||||
return ""
|
||||
except Exception:
|
||||
return ""
|
||||
|
||||
|
||||
def _load_audio_numpy(audio_path: str) -> tuple:
|
||||
"""Load audio as numpy array. Returns (samples, sample_rate) or (None, None)."""
|
||||
try:
|
||||
import librosa
|
||||
samples, sr = librosa.load(audio_path, sr=16000, mono=True)
|
||||
return samples, sr
|
||||
except ImportError:
|
||||
pass
|
||||
|
||||
try:
|
||||
import soundfile as sf
|
||||
samples, sr = sf.read(audio_path)
|
||||
if len(samples.shape) > 1:
|
||||
samples = samples.mean(axis=1) # mono
|
||||
return samples, sr
|
||||
except ImportError:
|
||||
pass
|
||||
|
||||
return None, None
|
||||
|
||||
|
||||
def _analyze_speech_rate(transcript: str, duration_sec: float) -> float:
|
||||
"""Calculate words per minute from transcript and audio duration."""
|
||||
if not transcript or duration_sec <= 0:
|
||||
return 0.0
|
||||
words = len(transcript.split())
|
||||
minutes = duration_sec / 60.0
|
||||
return words / minutes if minutes > 0 else 0.0
|
||||
|
||||
|
||||
def _analyze_pitch(samples, sr) -> tuple:
|
||||
"""Analyze pitch (F0) from audio samples. Returns (mean_hz, variability_hz)."""
|
||||
try:
|
||||
import librosa
|
||||
f0, voiced_flag, _ = librosa.pyin(
|
||||
samples, fmin=librosa.note_to_hz('C2'),
|
||||
fmax=librosa.note_to_hz('C7'), sr=sr
|
||||
)
|
||||
import numpy as np
|
||||
f0_clean = f0[~np.isnan(f0)]
|
||||
if len(f0_clean) == 0:
|
||||
return 0.0, 0.0
|
||||
return float(np.mean(f0_clean)), float(np.std(f0_clean))
|
||||
except (ImportError, Exception):
|
||||
return 0.0, 0.0
|
||||
|
||||
|
||||
def _analyze_silence(samples, sr, threshold_db: float = -40.0) -> float:
|
||||
"""Calculate ratio of silence in audio (0-1)."""
|
||||
try:
|
||||
import librosa
|
||||
import numpy as np
|
||||
rms = librosa.feature.rms(y=samples)[0]
|
||||
rms_db = librosa.amplitude_to_db(rms, ref=np.max)
|
||||
silence_frames = np.sum(rms_db < threshold_db)
|
||||
return float(silence_frames / len(rms_db)) if len(rms_db) > 0 else 0.0
|
||||
except (ImportError, Exception):
|
||||
return 0.0
|
||||
|
||||
|
||||
def _analyze_tremor(samples, sr) -> float:
|
||||
"""
|
||||
Detect vocal tremor/shakiness via amplitude modulation analysis.
|
||||
|
||||
Tremor manifests as periodic amplitude fluctuations (3-12 Hz range).
|
||||
Returns 0-1 score where 1 = strong tremor detected.
|
||||
"""
|
||||
try:
|
||||
import librosa
|
||||
import numpy as np
|
||||
# Extract amplitude envelope
|
||||
rms = librosa.feature.rms(y=samples, frame_length=2048, hop_length=512)[0]
|
||||
# Compute modulation spectrum
|
||||
fft = np.abs(np.fft.rfft(rms))
|
||||
freqs = np.fft.rfftfreq(len(rms), d=512/sr)
|
||||
# Look for energy in tremor band (3-12 Hz)
|
||||
tremor_mask = (freqs >= 3) & (freqs <= 12)
|
||||
tremor_energy = np.sum(fft[tremor_mask])
|
||||
total_energy = np.sum(fft[1:]) # skip DC
|
||||
if total_energy == 0:
|
||||
return 0.0
|
||||
ratio = tremor_energy / total_energy
|
||||
return float(min(1.0, ratio * 5)) # normalize — typical tremor is 0.1-0.3 of total
|
||||
except (ImportError, Exception):
|
||||
return 0.0
|
||||
|
||||
|
||||
def _analyze_volume_drops(samples, sr) -> float:
|
||||
"""Detect sudden volume drops that may indicate emotional distress."""
|
||||
try:
|
||||
import librosa
|
||||
import numpy as np
|
||||
rms = librosa.feature.rms(y=samples, frame_length=2048, hop_length=512)[0]
|
||||
if len(rms) < 2:
|
||||
return 0.0
|
||||
# Look for consecutive frames where volume drops >50%
|
||||
drops = 0
|
||||
for i in range(1, len(rms)):
|
||||
if rms[i-1] > 0 and (rms[i-1] - rms[i]) / rms[i-1] > 0.5:
|
||||
drops += 1
|
||||
return float(min(1.0, drops / (len(rms) * 0.1)))
|
||||
except (ImportError, Exception):
|
||||
return 0.0
|
||||
|
||||
|
||||
def _compute_distress_score(result: VoiceAnalysisResult) -> tuple:
|
||||
"""
|
||||
Compute composite distress score from paralinguistic features.
|
||||
|
||||
Returns (score, signals_detected).
|
||||
"""
|
||||
signals = []
|
||||
score = 0.0
|
||||
weights = 0
|
||||
|
||||
# Speech rate (0.2 weight)
|
||||
if result.speech_rate > 0:
|
||||
if result.speech_rate < SPEECH_RATE_SLOW:
|
||||
signals.append(f"very_slow_speech ({result.speech_rate:.0f} WPM)")
|
||||
score += 0.8 * 0.2
|
||||
elif result.speech_rate > SPEECH_RATE_FAST:
|
||||
signals.append(f"very_fast_speech ({result.speech_rate:.0f} WPM)")
|
||||
score += 0.6 * 0.2
|
||||
elif result.speech_rate < SPEECH_RATE_NORMAL_LOW:
|
||||
score += 0.3 * 0.2
|
||||
weights += 0.2
|
||||
|
||||
# Pitch variability (0.25 weight — monotone is strong depression indicator)
|
||||
if result.pitch_variability > 0:
|
||||
if result.pitch_variability < PITCH_VARIABILITY_LOW:
|
||||
signals.append(f"monotone_voice (variability={result.pitch_variability:.1f} Hz)")
|
||||
score += 0.9 * 0.25
|
||||
elif result.pitch_variability < PITCH_VARIABILITY_NORMAL:
|
||||
signals.append(f"reduced_pitch_variability ({result.pitch_variability:.1f} Hz)")
|
||||
score += 0.5 * 0.25
|
||||
weights += 0.25
|
||||
|
||||
# Silence ratio (0.2 weight)
|
||||
if result.silence_ratio > 0:
|
||||
if result.silence_ratio > SILENCE_RATIO_HIGH:
|
||||
signals.append(f"excessive_silence ({result.silence_ratio:.0%})")
|
||||
score += 0.7 * 0.2
|
||||
elif result.silence_ratio < SILENCE_RATIO_LOW:
|
||||
signals.append(f"minimal_pauses ({result.silence_ratio:.0%})")
|
||||
score += 0.3 * 0.2
|
||||
weights += 0.2
|
||||
|
||||
# Tremor (0.2 weight)
|
||||
if result.tremor_score > 0:
|
||||
if result.tremor_score > 0.5:
|
||||
signals.append(f"vocal_tremor (score={result.tremor_score:.2f})")
|
||||
score += result.tremor_score * 0.2
|
||||
weights += 0.2
|
||||
|
||||
# Volume drops (0.15 weight)
|
||||
if result.volume_drop_score > 0:
|
||||
if result.volume_drop_score > 0.4:
|
||||
signals.append(f"volume_drops (score={result.volume_drop_score:.2f})")
|
||||
score += result.volume_drop_score * 0.15
|
||||
weights += 0.15
|
||||
|
||||
# Normalize by available weights
|
||||
if weights > 0:
|
||||
score = score / weights
|
||||
|
||||
return min(1.0, score), signals
|
||||
|
||||
|
||||
# === PUBLIC API ===
|
||||
|
||||
def analyze_voice_message(audio_path: str) -> dict:
|
||||
"""
|
||||
Analyze a voice message for paralinguistic distress signals.
|
||||
|
||||
Args:
|
||||
audio_path: Path to audio file (OGG, WAV, MP3, etc.)
|
||||
|
||||
Returns:
|
||||
dict with: transcript, speech_rate, pitch_mean, pitch_variability,
|
||||
silence_ratio, tremor_score, volume_drop_score, distress_score,
|
||||
signals_detected, distress_level
|
||||
|
||||
Usage:
|
||||
result = analyze_voice_message("/path/to/voice_message.ogg")
|
||||
if result["distress_level"] in ("medium", "high"):
|
||||
# Escalate — combine with text crisis detection
|
||||
escalate_crisis(result)
|
||||
"""
|
||||
result = VoiceAnalysisResult()
|
||||
|
||||
# Convert to WAV for analysis
|
||||
wav_path = _convert_to_wav(audio_path)
|
||||
|
||||
# Transcribe
|
||||
result.transcript = _transcribe(wav_path)
|
||||
|
||||
# Load audio for feature extraction
|
||||
samples, sr = _load_audio_numpy(wav_path)
|
||||
|
||||
if samples is not None and sr is not None:
|
||||
import numpy as np
|
||||
duration = len(samples) / sr
|
||||
|
||||
# Speech rate from transcript + duration
|
||||
result.speech_rate = _analyze_speech_rate(result.transcript, duration)
|
||||
|
||||
# Pitch analysis
|
||||
result.pitch_mean, result.pitch_variability = _analyze_pitch(samples, sr)
|
||||
|
||||
# Silence ratio
|
||||
result.silence_ratio = _analyze_silence(samples, sr)
|
||||
|
||||
# Tremor detection
|
||||
result.tremor_score = _analyze_tremor(samples, sr)
|
||||
|
||||
# Volume drops
|
||||
result.volume_drop_score = _analyze_volume_drops(samples, sr)
|
||||
|
||||
# Composite distress score
|
||||
result.distress_score, result.signals_detected = _compute_distress_score(result)
|
||||
|
||||
# Clean up temp file
|
||||
if wav_path != audio_path and os.path.exists(wav_path):
|
||||
os.unlink(wav_path)
|
||||
|
||||
# Classify distress level
|
||||
if result.distress_score >= DISTRESS_MEDIUM:
|
||||
distress_level = "high"
|
||||
elif result.distress_score >= DISTRESS_LOW:
|
||||
distress_level = "medium"
|
||||
elif result.distress_score > 0:
|
||||
distress_level = "low"
|
||||
else:
|
||||
distress_level = "none"
|
||||
|
||||
output = result.to_dict()
|
||||
output["distress_level"] = distress_level
|
||||
return output
|
||||
|
||||
|
||||
def get_audio_duration(audio_path: str) -> float:
|
||||
"""Get audio duration in seconds."""
|
||||
try:
|
||||
import librosa
|
||||
duration = librosa.get_duration(path=audio_path)
|
||||
return float(duration)
|
||||
except (ImportError, Exception):
|
||||
try:
|
||||
import soundfile as sf
|
||||
info = sf.info(audio_path)
|
||||
return float(info.duration)
|
||||
except (ImportError, Exception):
|
||||
return 0.0
|
||||
Reference in New Issue
Block a user