Обработка по «A Practical Introduction to Wavelet Analysis in Electroretinography» (Shwetar et al., 2025)
Авторский код: CWT amor, линейный detrend, нормированный модуль |W| (0–1), DWT Haar с энергией |D|² (0–1), общий максимум для пары сигналов и палитра jet. Здесь воспроизведена эта последовательность для CWT/DWT. FFT обрабатывается по прежней методике отдельно.
Адаптация браузера: Morlet ω₀ = 6, L1-нормировка, симметричное продолжение границ; частоты 5–100 Гц, 48 строк для Roland PERG. Сетка, продолжение и нормировка не верифицированы на идентичность MATLAB cwt: численного совпадения с MATLAB не заявляем. Краевая область отмечена штрихами. Полосы DWT dyadic, подписи — условные центральные частоты; они не являются физиологическими генераторами ответа.
Статья рассматривает full-field ERG. Перенос на PERG — исследовательская адаптация, без диагностических порогов. Фильтр 0,3–300 Гц относится к регистрации данных авторов: повторно к Roland не применяется. Сохранена возможность выбрать прежний режим CWT.
Авторский MATLAB-код · Опубликованная версия статьи
Загрузить исследование Roland
Выберите один или несколько экспортов Roland RETIport RAW или mfERG FOK 61. XLSX, XLS и CSV; до 20 МБ на файл. Файл обрабатывается в вашем браузере.
Пациент 001 · демонстрационное исследование
Распознавание методики Roland
Поддерживаются полные экспорты RETIport RAW (PERG, VEP, ERG и другие именованные протоколы) и mfERG FOK 61. Отчёты без исходных кривых и иные форматы требуют отдельного парсера. Неопределённую методику уточните вручную.
Независимый SNRНедоступен
Подготовка графиков…
OD правый глаз
FFT · амплитудный спектр
CWT Morlet
Мощность выбранной частоты во времени
DWT · Haar
OS левый глаз
FFT · амплитудный спектр
CWT Morlet
Мощность выбранной частоты во времени
DWT · Haar
Методика и ограничения
FFT: удаление среднего всей эпохи, симметричное окно Hann, односторонний спектр; амплитуда d|FFT|/Σw, d = 2 для внутренних бинов и 1 для DC/Nyquist. Полный спектр и PSD доступны в JSON; на графике 0–100 Гц. Шаг сетки 5,54 Гц не равен полной разрешающей способности окна. Для transient спектр описывает форму эпохи; точность гармоник steady-state ограничена длительностью записи.
Прежний режим CWT: комплексный Morlet, ω₀ = 6, 48 логарифмических частот 5–100 Гц. Прямая дискретная свёртка, нормировка 1/√s, удаление среднего; нулевое продолжение вне эпохи. В прежнем режиме цвет: 10 log₁₀(|W|² / максимум), диапазон −40…0 dB. Это относительная мощность коэффициентов CWT, не PSD. Штриховка обозначает влияние края: расстояние до края < √2s. Низкие частоты в короткой эпохе особенно ненадёжны.
Для демонстрационного файла фильтры приобретения 1–50 Гц. В новых файлах параметры показаны после загрузки; пунктир 50 Гц — ориентир. Область выше 50 Гц показана только для контроля, с пунктирной границей. Дополнительная цифровая фильтрация не применяется. Transient и steady-state анализируются отдельно; частота StimFreq не интерпретируется автоматически как частота реверсий. Название исходного файла «норма» не заменяет совместимые лабораторные референсы. Заключение о норме или патологии здесь не формируется.
Morlet и cone of influence: PyCWTYousif Shwetar1, David Lalush1, Jason McAnany2, Brett Jeffrey3, Melissa Haendel4
PMID: 40778145
PMCID: PMC12330432
DOI: 10.1101/2025.07.25.25331915
Организации авторов
- Joint Department of Biomedical Engineering, University of North Carolina and North Carolina State University, Chapel Hill, NC, United States.
- Department of Ophthalmology and Visual Sciences, University of Illinois Chicago, Chicago, Illinois, United States.
- Ophthalmic Genetics and Visual Function Branch, National Eye Institute, NIH, Bethesda, MD, United States.
- Department of Genetics, University of North Carolina, Chapel Hill, NC, United States.
Открыть статью в PubMed