Non-contact estimation of respiratory pattern (RP) and respiration rate (RR) has multiple applications. Existing methods for RP and RR measurement fall into one of the three categories - (i) estimation through nasal air flow measurement, (ii) estimation from video-based remote photoplethysmography, and (iii) estimation by measurement of motion induced by respiration using motion detectors. These methods, however, require specialized sensors, are computationally expensive and/or critically depend on selection of a region of interest (ROI) for processing. In this paper a general framework is described for estimating a periodic signal driving noisy LTI channels connected in parallel with unknown dynamics. The method is then applied to derive a computationally inexpensive method for estimating RP using 2D cameras that does not critically depend on ROI. Specifically, RP is estimated by imaging the changes in the reflected light caused by respiration-induced motion. Each spatial location in the field of view of the camera is modeled as a noise-corrupted linear time-invariant (LTI) measurement channel with unknown system dynamics, driven by a single generating respiratory signal. Estimation of RP is cast as a blind deconvolution problem and is solved through a method comprising subspace projection and statistical aggregation. Experiments are carried out on 31 healthy human subjects by generating multiple RPs and comparing the proposed estimates with simultaneously acquired ground truth from an impedance pneumograph device. The proposed estimator agrees well with the ground truth device in terms of correlation measures, despite variability in clothing pattern, angle of view and ROI.
翻译:对呼吸模式(RP)和呼吸率(RR)的非接触估计具有多种应用。现有RP和RR测量方法可分为三种类别之一:(一) 通过鼻气流测量估算,(二) 通过视频远程光谱采集估算,(三) 通过使用运动探测器测量呼吸变化产生的运动量度估算,但这些方法需要专门的传感器,在计算上成本昂贵和(或)关键取决于对处理感兴趣的区域(ROI)的选择。本文描述了一个总体框架,用于估计一个定期信号,该信号驱动着与未知的动态同时连接的噪音LTI频道。然后,该方法用于利用不严重依赖ROI的2D相机估算RP的计算,(二)基于视频的远程光谱谱谱采集,(三)通过对反映的光的变化进行成像来估算。这些方法在摄影场的每个空间位置都建模成一个噪音-线性时间变压(ROI)测量频道,由生成的单一呼吸空间信号驱动的系统变化模式驱动,然后用来计算出一种计算结果变压(三),在模拟中,在模拟中,在模拟变压中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟变压中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟中,在模拟变相变相变压中,在模拟中,将,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,在模拟,进行,进行,在模拟,在模拟,在模拟,在模拟,进行,进行,进行,进行,进行中,进行中,进行着,进行着,在模拟,在模拟,进行,进行中,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行着,进行,