TW202400076A - Evaluation method of sleep quality and computing apparatus related to sleep quality - Google Patents
Evaluation method of sleep quality and computing apparatus related to sleep quality Download PDFInfo
- Publication number
- TW202400076A TW202400076A TW112101793A TW112101793A TW202400076A TW 202400076 A TW202400076 A TW 202400076A TW 112101793 A TW112101793 A TW 112101793A TW 112101793 A TW112101793 A TW 112101793A TW 202400076 A TW202400076 A TW 202400076A
- Authority
- TW
- Taiwan
- Prior art keywords
- sleep quality
- sleep
- index
- characteristic data
- data
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Operations Research (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Entrepreneurship & Innovation (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
本發明是有關於一種資料分析技術,且特別是有關於一種睡眠品質的評估方法及相關於睡眠品質的運算裝置。The present invention relates to a data analysis technology, and in particular to a sleep quality evaluation method and a computing device related to sleep quality.
睡眠呼吸中止(Sleep apnea)是指睡眠時呼吸不自覺變弱甚至停止的症狀。呼吸中止通常難以察覺,直到身體嚴重缺氧才可能因不舒服而清醒醒來。然而,缺氧恐傷害身體,患者甚至可能因心血管疾病猝死。睡眠呼吸中止症的患者通常不知有症狀。除非患者自行到醫院透過特殊設備進行偵測診斷,才能發現症狀。Sleep apnea refers to the symptom of involuntary weakening or even stopping of breathing during sleep. The cessation of breathing is usually unnoticeable until the body is severely deprived of oxygen and may wake up from discomfort. However, hypoxia may harm the body, and patients may even die suddenly from cardiovascular disease. People with sleep apnea are often unaware of the symptoms. Symptoms can only be discovered unless the patient goes to the hospital for detection and diagnosis using special equipment.
有鑑於此,本發明實施例提供一種睡眠品質的評估方法及相關於睡眠品質的運算裝置,可輕易地檢測睡眠品質。In view of this, embodiments of the present invention provide a method for evaluating sleep quality and a computing device related to sleep quality, which can easily detect sleep quality.
本發明實施例的睡眠品質的評估方法包括(但不僅限於)下列步驟:取得感測資料。這感測資料是基於雷達回波所產生。將感測資料轉換成特徵資料。這特徵資料包括雷達回波在波形上的多個特徵點的統計量。依據特徵資料決定睡眠品質資訊。這睡眠品質資訊相關於睡眠品質的優劣程度。The sleep quality evaluation method according to the embodiment of the present invention includes (but is not limited to) the following steps: obtaining sensing data. This sensing data is generated based on radar echoes. Convert sensing data into feature data. This feature data includes the statistics of multiple feature points on the waveform of the radar echo. Sleep quality information is determined based on characteristic data. This sleep quality information is related to the degree of sleep quality.
本發明實施例的相關於睡眠品質的運算裝置包括(但不僅限於)記憶體及處理器。記憶體用以儲存程式碼。處理器耦接記憶體。處理器載入程式碼以執行:取得感測資料,將感測資料轉換成特徵資料,並依據特徵資料決定睡眠品質資訊。這感測資料是基於雷達回波所產生。這特徵資料包括雷達回波在波形上的多個特徵點的統計量。這睡眠品質資訊相關於睡眠品質的優劣程度。Computing devices related to sleep quality in embodiments of the present invention include (but are not limited to) memories and processors. Memory is used to store program code. The processor is coupled to the memory. The processor loads the program code to execute: obtain sensing data, convert the sensing data into characteristic data, and determine sleep quality information based on the characteristic data. This sensing data is generated based on radar echoes. This feature data includes the statistics of multiple feature points on the waveform of the radar echo. This sleep quality information is related to the degree of sleep quality.
基於上述,依據本發明實施例的睡眠品質的評估方法及相關於睡眠品質的運算裝置,使用基於雷達的感測資料預測睡眠品質資訊。自感測資料所取得的特徵資料對應於多項生理睡眠檢查(Polysomnography,PSG),多項生理睡眠檢查可反映於呼吸事件,呼吸事件又相關於睡眠品質的優劣。藉此,可透過非接觸式感測來評估睡眠品質。Based on the above, according to the sleep quality evaluation method and sleep quality-related computing device according to the embodiment of the present invention, radar-based sensing data is used to predict sleep quality information. The characteristic data obtained from the self-sensing data corresponds to multiple physiological sleep examinations (Polysomnography, PSG). Multiple physiological sleep examinations can be reflected in respiratory events, and respiratory events are related to the quality of sleep. In this way, sleep quality can be assessed through non-contact sensing.
為讓本發明的上述特徵和優點能更明顯易懂,下文特舉實施例,並配合所附圖式作詳細說明如下。In order to make the above-mentioned features and advantages of the present invention more obvious and easy to understand, embodiments are given below and described in detail with reference to the accompanying drawings.
圖1是依據本發明一實施例的運算裝置10及雷達50的元件方塊圖。請參照圖1,運算裝置10包括(但不僅限於)記憶體11及處理器12。運算裝置10可以是手機、平板電腦、筆記型電腦、桌上型電腦、語音助理裝置、智能家電、穿戴式裝置、車載裝置或其他電子裝置。FIG. 1 is a block diagram of components of the
記憶體11可以是任何型態的固定或可移動隨機存取記憶體(Radom Access Memory,RAM)、唯讀記憶體(Read Only Memory,ROM)、快閃記憶體(flash memory)、傳統硬碟(Hard Disk Drive,HDD)、固態硬碟(Solid-State Drive,SSD)或類似元件。在一實施例中,記憶體11用以儲存程式碼、軟體模組、組態配置、資料或檔案(例如,資料、事件、資訊、模型、或特徵),並待後續實施例詳述。The
處理器12耦接記憶體11。處理器12可以是中央處理單元(Central Processing Unit,CPU)、圖形處理單元(Graphic Processing unit,GPU),或是其他可程式化之一般用途或特殊用途的微處理器(Microprocessor)、數位信號處理器(Digital Signal Processor,DSP)、可程式化控制器、現場可程式化邏輯閘陣列(Field Programmable Gate Array,FPGA)、特殊應用積體電路(Application-Specific Integrated Circuit,ASIC)、神經網路加速器或其他類似元件或上述元件的組合。在一實施例中,處理器12用以執行運算裝置10的所有或部份作業,且可載入並執行記憶體11所儲存的各程式碼、軟體模組、檔案及資料。在一些實施例中,本發明實施例的方法中的部分作業可能透過不同或相同處理器12實現。The
在一實施例中,處理器12連接雷達50。例如,雷達50透過USB、Thunderbolt、Wi-Fi、藍芽或其他有線或無線通訊技術連接處理器12。又例如,運算裝置10內建雷達50,並透過內部線路連接處理器12與雷達50。雷達50可以是調頻連續波(Frequency Modulated Continuous Wave,FMCW)雷達或脈衝雷達(Impulse Radio,IR)-超寬頻(Ultra-Wideband,UWB)雷達。在一實施例中,雷達50用以產生感測資料。感測資料是基於雷達回波所產生。雷達回波是指雷達50的傳送訊號經物體(例如,人體、或衣服)反射的回波訊號。感測資料是雷達50的感測結果。例如,同相(In-phase)及/或正交(quadrature)訊號。In one embodiment,
在一實施例中,雷達50的傳送訊號的頻率可以是24GHz或其他可反射人體(例如,胸部或腹部)的頻率。In one embodiment, the frequency of the transmitted signal of the
在一應用情境中,雷達50可置於床頭、床邊或床尾,且其傳送訊號朝向人體的胸部或腹部,並據以檢測胸部或腹部的起伏情況。然而,雷達50的設置位置及朝向仍可依據實際需求而變更,且本發明實施例不加以限制。In an application scenario, the
下文中,將搭配運算裝置10及雷達50中的各項裝置、元件及模組說明本發明實施例所述之方法。本方法的各個流程可依照實施情形而隨之調整,且不僅限於此。In the following, the method described in the embodiment of the present invention will be described with reference to various devices, components and modules in the
圖2是依據本發明一實施例的睡眠品質的評估方法的流程圖。請參照圖2,處理器12取得感測資料(步驟S210)。具體而言,感測資料是基於雷達回波所產生。例如,雷達50發射連續波訊號,且連續波訊號受胸部或腹部反射而形成雷達回波。雷達50即可接收雷達回波,並據以產生感測資料。感測資料以同相(I)及正交(Q)訊號為例,同相訊號I1為[0.164144 0.179153 0.194716 ... 1.600188 1.590467 1.586891],且正交訊號Q1為[2.295545 2.278471 2.270613 ... 1.031502 1.027573 1.049331]。FIG. 2 is a flow chart of a sleep quality evaluation method according to an embodiment of the present invention. Referring to FIG. 2 , the
在一實施例中,處理器12可累積一段時間的感測資料。這一段時間例如是1、5或8小時。In one embodiment, the
處理器12將感測資料轉換成特徵資料(步驟S220)。在一實施例中,特徵資料包括感測資料中的二通道之間或單一通道內的變異數。這兩通道可以是同相及正交訊號。變異數的數學表示式為:
…(1)
Cov為變異數,X、Y為同相及正交訊號中的任一者,
為X的平均值,
為Y的平均值。以前述同相訊號I1及正交訊號Q1為例,其變異數為-0.21645484961728612。
The
在一實施例中,特徵資料包括感測資料的熵(entropy)。在資訊理論中,熵是指接收的每條訊息中包含的資訊的平均量,為不確定性的量度,且熵因訊息來源越隨機而越大。以熵為基礎的特徵例如是相對熵(relative entropy)、條件熵(conditional entropy)、互斥資訊(Mutual Information)、資訊熵(Information entropy)、夏農熵(Shannon entropy)、區塊熵(block entropy。In one embodiment, the characteristic data includes entropy of the sensed data. In information theory, entropy refers to the average amount of information contained in each message received. It is a measure of uncertainty, and entropy is larger as the source of the message is more random. Features based on entropy include relative entropy, conditional entropy, mutual information, information entropy, Shannon entropy, and block entropy. entropy.
以Shannon entropy為例,隨機變量X(具有x={x 1, …, x n}的值域)的熵H定義為: …(2), P X為隨機變數X的機率質量函數(probability mass function),且b為對數所用的底。 Taking Shannon entropy as an example, the entropy H of a random variable X (with the value range of x={x 1 , …, x n }) is defined as: ...(2), P X is the probability mass function of the random variable X, and b is the base used for logarithms.
以前述同相訊號I1及正交訊號Q1為例,其條件熵為1.4112874013149717。Taking the aforementioned in-phase signal I1 and quadrature signal Q1 as an example, their conditional entropy is 1.4112874013149717.
在一實施例中,特徵資料包括雷達回波在波形上的多個特徵點的統計量。特徵點可以是波形中的峰(peak)值及/或谷(valley)值。例如,圖3是依據本發明一實施例的波形相關特徵資料的示意圖。請參照圖3,這波形包括峰值P1、P2及谷值V1。峰值P1、P2可以是一個或多個週期內的最大值。而谷值V1可以是一個或多個週期內的最小值。In one embodiment, the feature data includes statistics of multiple feature points on the waveform of the radar echo. Feature points may be peak values and/or valley values in the waveform. For example, FIG. 3 is a schematic diagram of waveform-related characteristic data according to an embodiment of the present invention. Please refer to Figure 3. This waveform includes peak values P1, P2 and valley value V1. The peak values P1 and P2 can be the maximum values in one or more periods. The valley value V1 can be the minimum value within one or more periods.
此外,統計量可以是兩特徵點之間的間隔(interval)、間隔的變化量及/或那些特徵點的總數量。以圖3為例,統計量是兩峰值P1之間的間隔I PP。然而,統計量也可能是谷值V1與另一個相鄰谷值(圖未示)之間的間隔、峰值P1/P2與谷值V1之間的間隔或波形中指定兩點之間的間隔。間隔的變化量例如是兩個或更多個間隔之間的差異。例如,間隔I PP與另一個間隔I PP(圖未示,例如是峰值P2與下一個峰值之間的間隔)之間的差異。特徵點的總數量例如是一段時間(例如,1000個取樣點或3個小時)內的峰值及/或谷值的總數量。 Furthermore, the statistic may be an interval between two feature points, a variation of the interval, and/or the total number of those feature points. Taking Figure 3 as an example, the statistic is the interval I PP between the two peak values P1. However, the statistic may also be the interval between the valley value V1 and another adjacent valley value (not shown), the interval between the peak value P1/P2 and the valley value V1, or the interval between two specified points in the waveform. The variation of intervals is, for example, the difference between two or more intervals. For example, the difference between the interval I PP and another interval I PP (not shown in the figure, for example, the interval between the peak P2 and the next peak). The total number of feature points is, for example, the total number of peaks and/or valleys within a period of time (eg, 1000 sampling points or 3 hours).
在一實施例中,處理器12可分別決定同相及正交訊號的波形的統計量,也可取兩訊號的統計量的平均值,以作為特徵資料。In one embodiment, the
在一實施例中,特徵資料包括波形的趨勢,且這趨勢為波形在去除規律特性下的強度變化。規律特性可以是波形的週期性變化。例如,弦波訊號由零增加至最大值,由最大值下降至最小值,再由最小值增加至零的反覆變化。去除規律特性之後,即留下波形的走勢。也就是,強度變化。處理器12將趨勢抽象化(即,排除訊號的絕對強度干擾),即可作為描述睡眠品質(例如,呼吸事件、或睡眠事件)的特徵資料。In one embodiment, the characteristic data includes a trend of the waveform, and the trend is the intensity change of the waveform without regular characteristics. Regular characteristics can be periodic changes in the waveform. For example, a sine wave signal increases repeatedly from zero to a maximum value, drops from a maximum value to a minimum value, and then increases from a minimum value to zero. After removing the regular characteristics, the trend of the waveform is left. that is, intensity changes. The
例如,圖4是依據本發明一實施例的趨勢相關特徵資料的示意圖。請參照圖4,雷達回波可拆分成趨勢與規律。趨勢為線性函數,且線性函數的斜率為正,因此這波形的趨勢是強度逐漸增大。而規律為弦波。於另一實施例中,趨勢亦可為曲線。此外,以程式語言及前述同相訊號I1為例,python的package: seasonal_decompose可將同相訊號I1趨勢裡最大值的除以最小值,並可得到: 1.207502488(作為走勢的代表值)。For example, FIG. 4 is a schematic diagram of trend-related feature data according to an embodiment of the present invention. Please refer to Figure 4. Radar echoes can be divided into trends and patterns. The trend is a linear function, and the slope of the linear function is positive, so the trend of this waveform is that the intensity gradually increases. The pattern is a string wave. In another embodiment, the trend can also be a curve. In addition, taking the programming language and the aforementioned in-phase signal I1 as an example, python's package: seasonal_decompose can divide the maximum value in the trend of the in-phase signal I1 by the minimum value, and obtain: 1.207502488 (as a representative value of the trend).
在一實施例中,處理器12可分別決定同相及正交訊號的波形的趨勢,也可取兩訊號的趨勢的平均值,以作為特徵資料。In one embodiment, the
在一實施例中,處理器12可選擇前述感測資料的統計量、變異數、熵及趨勢中的一者或更多者作為特徵資料。In one embodiment, the
請參照圖2,處理器12依據特徵資料決定睡眠品質資訊(步驟S230)。具體而言,睡眠品質資訊相關睡眠品質的優劣程度。在一實施例中,睡眠品質資訊包括呼吸事件。例如,正常呼吸、呼吸不足/低通氣(hypopnea)、淺慢呼吸(Flow Limitation)、阻塞呼吸、清醒或無呼吸(apnea)事件。無呼吸定義為睡眠時呼吸完全停止,呼吸暫停,氣流中斷至少10秒以上,稱為一件事件(即,一次睡眠呼吸暫停)。低通氣呼吸定義為不正常的呼吸型態,對成年人而言,一次也至少超過10秒時間,氣流量及胸腹部的呼吸運動較正常的情況降低到只有30%~50%的程度,同時血液中的氧氣飽和度至少降低了4%。淺慢呼吸定義為不正常的呼吸型態,因呼吸道部份地受阻,氣流的通過量較正常的流量為低。Referring to FIG. 2 , the
而處理器12可依據特徵資料預測呼吸事件。本發明實施例的特徵資料是經比對多項生理睡眠檢查(Polysomnography,PSG)(例如,呼吸氣流、胸腔動作、腹肌行為或腦波圖)所得出較能區別呼吸事件的特徵。然而,有別於基於PSG辨識呼吸事件,本發明實施例是基於雷達的特徵資料辨識呼吸事件。The
在一實施例中,處理器12可透過機器學習模型預測呼吸事件。機器學習模型受訓練而理解特徵資料與呼吸事件之間的關聯性。機器學習模型例如是基於深度神經決策樹(Deep Neural Decision Tree,DNDT)、深度學習神經網路(Temporal Convolutional Network,TCN)、決策樹(Decision Tree)、隨機森林(Random Forest)或其他機器學習演算法。深度學習神經網路例如是時間卷積網路(Temporal Convolutional Network,TCN)及卷積神經網路(Convolutional Neural Network,CNN)。DNDT是混合深度學習及決策樹的策略。機器學習演算法可分析訓練樣本以自中獲得規律,從而透過規律對未知資料預測。例如,機器學習模型依據已標記樣本(例如,已知呼吸低通氣事件的特徵資料、或已知正常呼吸事件的特徵資料)建立特徵資料(即,模型的輸入)與呼吸事件(即,模型的輸出)之間的隱藏層中各節點之間的關聯。而機器學習模型即是經學習後所建構出的模型,並可據以對待評估資料(例如,特徵資料)推論。In one embodiment, the
例如,CNN可學習影像相關特徵,且TCN可學習時序上的特徵。DNDT結合機器學習領域中的決策樹以及深度學習的概念。傳統決策樹無法最佳化每個樹節點,從而造成判斷上的限制。然而,DNDT可以讓決策樹的每個節點進行深度學習的學習機器決定權重。For example, CNN can learn image-related features, and TCN can learn temporal features. DNDT combines the concepts of decision trees and deep learning in the field of machine learning. Traditional decision trees cannot optimize each tree node, resulting in limitations in judgment. However, DNDT allows the learning machine to perform deep learning to determine the weight of each node of the decision tree.
例如,圖5是依據本發明一實施例的深度神經決策樹(DNDT)的示意圖。請參照圖5,在模型的訓練過程中,透過深度學習不斷改變決策樹中的各節點的權重(如圖中輸入至丟棄層(Binning Layer)之間的權重或克羅內克積(Kron Product)至輸出之間的權重)來進一步降低損失(loss)值,從而達成訓練權重的目的。丟棄層用以實現可微分的分群函數,且克羅內克積用於決定決策樹中的哪個葉值(leaf value)進行預測。因此,透過DNDT的架構學習特徵資料,可進一步用於判斷是否為正常睡眠或是有睡眠/呼吸事件的發生。For example, FIG. 5 is a schematic diagram of a deep neural decision tree (DNDT) according to an embodiment of the present invention. Please refer to Figure 5. During the training process of the model, deep learning is used to continuously change the weight of each node in the decision tree (as shown in the figure, the weight between the input to the discarding layer (Binning Layer) or the Kronecker Product (Kron Product) ) to the output) to further reduce the loss value, thereby achieving the purpose of training weights. The dropout layer is used to implement differentiable grouping functions, and the Kronecker product is used to determine which leaf value in the decision tree to make predictions. Therefore, the characteristic data learned through the DNDT architecture can be further used to determine whether it is normal sleep or whether a sleep/breathing event occurs.
以前述同相訊號I1為例,處理器12可將同相訊號I1直接轉換成二維的矩陣(例如,矩陣大小為40 x 50或 30 x 50)(作為模型的輸入)對機器學習模型進行學習。或者,特徵資料的形式可能依不同時間大小有所變化(例如,100、1500、2000或2500個取樣點但不以此為限)。在一些應用情境中,時間越長或取樣點數越多,則準確度較高,但不以此為限。又或者,處理器12可利用表格(table)類型的特徵資料(作為模型的輸入)對機器學習模型進行學習。例如,將這些時間序列的特徵資料轉換成數值後整理成表格如下:
表(1)
圖6是依據本發明一實施例的感測資料與事件的示意圖。請參照圖6,事件E1為呼吸低通氣事件,而感測資料及PSG可反映出明顯的快速起伏變化。例如,事件E1的末段有較高或較低數值。事件E2為正常睡眠事件,因此感測資料及PSG的波形大致規律變化。Figure 6 is a schematic diagram of sensing data and events according to an embodiment of the present invention. Please refer to Figure 6. Event E1 is a respiratory hypoventilation event, and the sensing data and PSG can reflect obvious rapid fluctuations. For example, the end of event E1 has a higher or lower value. Event E2 is a normal sleep event, so the waveforms of the sensing data and PSG change roughly regularly.
處理器12可統計一段時間(例如,2、5或8小時)內特定呼吸事件的次數及/或期間,以作為睡眠品質資訊。若正常睡眠事件的統計量越高,則睡眠品質較好(例如,優劣程度越高,且高代表優);若諸如呼吸低通氣及/或無呼吸的統計量越高,則睡眠品質較差(優劣程度越低,且低代表劣)。The
在一實施例中,睡眠品質資訊包括睡眠統計指標。睡眠統計指標為睡眠呼吸障礙指數(Respiratory Disturbance Index,RDI)或呼吸暫停低通氣指數(Apnea-Hypopnea Index,AHI)。RDI是睡眠時呼吸中斷的次數,也有人直接使用AHI。同樣情況的量測下,RDI的指數會稍微大於AHI的指數。依據美國睡眠協會的標準,AHI低於5為正常情況,AHI為5~14屬於輕度,AHI為15~29屬於中度,且AHI為30以上屬於重度睡眠呼吸障礙。也就是說,睡眠統計指標越低,則睡眠品質的優劣程度越高,且高代表優;睡眠統計指標越高,則睡眠品質的優劣程度越低,且低代表劣。In one embodiment, the sleep quality information includes sleep statistical indicators. The sleep statistical index is the Respiratory Disturbance Index (RDI) or the Apnea-Hypopnea Index (AHI). RDI is the number of interrupted breathing during sleep, and some people use AHI directly. Under the same measurement conditions, the RDI index will be slightly larger than the AHI index. According to the standards of the American Sleep Association, an AHI below 5 is normal, an AHI of 5 to 14 is mild, an AHI of 15 to 29 is moderate, and an AHI of 30 or above is severe sleep apnea. That is to say, the lower the sleep statistical index, the higher the degree of sleep quality, and high means excellent; the higher the sleep statistical index, the lower the degree of sleep quality, and low means poor.
在一實施例中,處理器12可依據預測的呼吸事件決定睡眠統計指標。處理器12可統計一段時間(例如,3小時、5小時或8小時)內先前呼吸事件的預測結果(例如,機器學習模型的輸出),並產生預測的睡眠統計指標。例如,特定呼吸事件的次數除以統計時間所得的值。In one embodiment, the
例如,表(2)是時間點(例如,每分鐘、每30分鐘或每小時)與預測結果的對應關係:
表(2)
為了驗證本發明實施例所產出的類RDI值(即,睡眠統計指標)是否能接近真實的RDI數值,實際蒐集臨床研究案中103位在睡眠中心進行睡眠檢測的資料,並用這些資料來進行驗證。圖7是依據本發明一實施例的指標驗證的示意圖。請參照圖7,本發明實施例所得出的類RDI(以Y軸的bRDI表示)與睡眠技師判斷的RDI(以X軸的RDI)進行趨勢及準確度分析(表(3))。由圖7可知,本發明實施例所得出的類RDI數值與真實的RDI數值呈現正相關且關聯度(例如,0.7481)有大於0.7。In order to verify whether the RDI-like value (i.e., sleep statistical index) produced by the embodiment of the present invention can be close to the real RDI value, the data of 103 people who underwent sleep testing in the sleep center in the clinical research case were actually collected, and these data were used to conduct Verify. Figure 7 is a schematic diagram of indicator verification according to an embodiment of the present invention. Please refer to Figure 7. The trend and accuracy analysis of the similar RDI (expressed as bRDI on the Y-axis) obtained by the embodiment of the present invention and the RDI judged by the sleep technician (expressed as the RDI on the X-axis) are performed (Table (3)). It can be seen from Figure 7 that the RDI-like value obtained by the embodiment of the present invention is positively correlated with the real RDI value, and the correlation degree (for example, 0.7481) is greater than 0.7.
另外,在臨床上會以RDI大於或等於15/小時(h)及大於或等於30/h的定義為有中度及重度以上的呼吸中止症狀,而表(3)可得出比對結果:
表(3)
在一實施例中,處理器12可依據特徵資料預測睡眠統計指標。例如,處理器12另外訓練另一個機器學習模型,並據以理解特徵資料與預測睡眠統計指標之間的關聯性。機器學習模型的介紹可參酌前述說明,於此不再贅述。例如,機器學習模型依據已標記樣本(例如,已知RDI的特徵資料、或已知AHI的特徵資料)建立特徵資料(即,模型的輸入)與睡眠統計指標(即,模型的輸出)之間的隱藏層中各節點之間的關聯。由於本發明實施例的特徵資料可用於區別呼吸事件且睡眠統計指標是基於呼吸事件(例如,特定的一個或更多個呼吸事件的次數除以統計時間)得出,因此可證明特徵資料可用於預測睡眠統計指標。In one embodiment, the
綜上所述,在本發明實施例的睡眠品質的評估方法及相關於睡眠品質的運算裝置中,依據雷達的感測資料所轉換的特徵資料(例如,相關於變異數、熵、波形及/或趨勢)判斷睡眠品質的好壞。藉此,可利用非接觸式感測的方式評估睡眠品質。To sum up, in the sleep quality evaluation method and sleep quality-related computing device according to the embodiment of the present invention, the characteristic data (for example, related to variation, entropy, waveform and/or related to the radar sensing data) is converted. or trend) to determine the quality of sleep. In this way, sleep quality can be assessed using non-contact sensing.
雖然本發明已以實施例揭露如上,然其並非用以限定本發明,任何所屬技術領域中具有通常知識者,在不脫離本發明的精神和範圍內,當可作些許的更動與潤飾,故本發明的保護範圍當視後附的申請專利範圍所界定者為準。Although the present invention has been disclosed above through embodiments, they are not intended to limit the present invention. Anyone with ordinary knowledge in the technical field may make some modifications and modifications without departing from the spirit and scope of the present invention. Therefore, The protection scope of the present invention shall be determined by the appended patent application scope.
10:運動裝置 11:記憶體 12:處理器 50:雷達 S210~S230:步驟 P1、P2:峰值 V1:谷值 I PP:間隔 E1、E2:事件 10: Motion device 11: Memory 12: Processor 50: Radar S210~S230: Steps P1, P2: Peak value V1: Valley value I PP : Interval E1, E2: Event
圖1是依據本發明一實施例的運算裝置及雷達的元件方塊圖。 圖2是依據本發明一實施例的睡眠品質的評估方法的流程圖。 圖3是依據本發明一實施例的波形相關特徵資料的示意圖。 圖4是依據本發明一實施例的趨勢相關特徵資料的示意圖。 圖5是依據本發明一實施例的深度神經決策樹(Deep Neural Decision Tree,DNDT)的示意圖。 圖6是依據本發明一實施例的感測資料與事件的示意圖。 圖7是依據本發明一實施例的指標驗證的示意圖。 FIG. 1 is a component block diagram of a computing device and a radar according to an embodiment of the present invention. FIG. 2 is a flow chart of a sleep quality evaluation method according to an embodiment of the present invention. FIG. 3 is a schematic diagram of waveform-related characteristic data according to an embodiment of the present invention. FIG. 4 is a schematic diagram of trend-related feature data according to an embodiment of the present invention. Figure 5 is a schematic diagram of a deep neural decision tree (Deep Neural Decision Tree, DNDT) according to an embodiment of the present invention. Figure 6 is a schematic diagram of sensing data and events according to an embodiment of the present invention. Figure 7 is a schematic diagram of indicator verification according to an embodiment of the present invention.
S210~S230:步驟 S210~S230: steps
Claims (20)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US18/190,125 US20230404474A1 (en) | 2022-06-16 | 2023-03-27 | Evaluation method of sleep quality and computing apparatus related to sleep quality |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US202263352644P | 2022-06-16 | 2022-06-16 | |
| US63/352,644 | 2022-06-16 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| TW202400076A true TW202400076A (en) | 2024-01-01 |
| TWI849690B TWI849690B (en) | 2024-07-21 |
Family
ID=89135642
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW111137595A TWI865935B (en) | 2022-06-16 | 2022-10-03 | Data processing method and apparatus |
| TW112101793A TWI849690B (en) | 2022-06-16 | 2023-01-16 | Evaluation method of sleep quality and computing apparatus related to sleep quality |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW111137595A TWI865935B (en) | 2022-06-16 | 2022-10-03 | Data processing method and apparatus |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN117252268A (en) |
| TW (2) | TWI865935B (en) |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TW201223505A (en) * | 2010-12-14 | 2012-06-16 | Univ Nat Cheng Kung | Medical analysis method and apparatus |
| TW201612845A (en) * | 2014-08-28 | 2016-04-01 | Resmed Ltd | Method, system and apparatus for diagnosis, monitoring and treatment of respiratory disorders |
| EP3356560A4 (en) * | 2015-09-30 | 2019-06-12 | Inform Genomics, Inc. | Systems and methods for predicting treatment-regimen-related outcomes |
| WO2018022986A1 (en) * | 2016-07-29 | 2018-02-01 | The Dun & Bradstreet Corporation | Diagnostic engine to enhance derandomized entity behavior identification and classification |
| CN108065916B (en) * | 2017-12-14 | 2021-04-09 | 中国人民解放军国防科技大学 | A non-contact sleep quality monitoring method based on bioradar |
| KR102658390B1 (en) * | 2017-12-22 | 2024-04-17 | 레스메드 센서 테크놀로지스 리미티드 | Devices, systems, and methods for health and medical sensing |
| EP3883468A2 (en) * | 2018-11-19 | 2021-09-29 | ResMed Sensor Technologies Limited | Methods and apparatus for detection of disordered breathing |
| CN109480783B (en) * | 2018-12-20 | 2022-02-18 | 深圳和而泰智能控制股份有限公司 | Apnea detection method and device and computing equipment |
| CN111142102B (en) * | 2019-12-26 | 2022-07-19 | 华为技术有限公司 | Respiratory data calculation method and related equipment |
| CN111210023B (en) * | 2020-01-13 | 2023-04-11 | 哈尔滨工业大学 | Automatic selection system and method for data set classification learning algorithm |
| CN111783093B (en) * | 2020-06-28 | 2024-10-22 | 南京航空航天大学 | A malware classification and detection method based on soft dependency |
-
2022
- 2022-10-03 TW TW111137595A patent/TWI865935B/en active
- 2022-10-28 CN CN202211333180.6A patent/CN117252268A/en active Pending
-
2023
- 2023-01-16 TW TW112101793A patent/TWI849690B/en active
Also Published As
| Publication number | Publication date |
|---|---|
| TWI849690B (en) | 2024-07-21 |
| TW202401306A (en) | 2024-01-01 |
| TWI865935B (en) | 2024-12-11 |
| CN117252268A (en) | 2023-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116602642B (en) | Heart rate monitoring method, device and equipment | |
| Zhuang et al. | Accurate contactless sleep apnea detection framework with signal processing and machine learning methods | |
| CN112043252B (en) | Emotion recognition system and method based on respiratory component in pulse signal | |
| TW202038141A (en) | Living body detection method and living body detection system | |
| Li et al. | Mo-sleep: Unobtrusive sleep and movement monitoring via Wi-Fi signal | |
| WO2022119811A1 (en) | Health monitoring system | |
| CN111387936A (en) | Sleep stage identification method, device and equipment | |
| He et al. | Optimum target range bin selection method for time-frequency analysis to detect falls using wideband radar and a lightweight network | |
| TWI849690B (en) | Evaluation method of sleep quality and computing apparatus related to sleep quality | |
| Chen et al. | Wi-Fi-based human activity recognition for continuous, whole-room monitoring of motor functions in Parkinson’s disease | |
| TWM593240U (en) | Detection device for apnea based on unipolar ecg | |
| US20230404474A1 (en) | Evaluation method of sleep quality and computing apparatus related to sleep quality | |
| US20220104724A1 (en) | Computer-based dynamic rating of ataxic breathing | |
| Liu et al. | The Multivariate Largest Lyapunov Exponent as an Age‐Related Metric of Quiet Standing Balance | |
| CN120323951A (en) | Non-contact sleeping posture detection and vital signs monitoring method based on millimeter wave radar | |
| KR20250012708A (en) | Server, method and computer program for determining disease using radar | |
| Yin et al. | Unobtrusive sleep health assessment using impulse radar: a pilot study in older people | |
| Pascual et al. | Analytical modelling of raw data for flow-guided in-body nanoscale localization | |
| KR102731278B1 (en) | Server, method and computer program for analysis of sleep breathing using radar | |
| US20220313156A1 (en) | Monitoring system and monitoring method for sleep apnea | |
| Zhao et al. | Rfid-based respiration monitoring under daily body motion disturbances | |
| Han | Respiratory patterns classification using UWB radar | |
| CN114287922A (en) | Method, device and equipment for detecting human body out-of-bed state and storage medium | |
| CN119864140B (en) | A health scanning method and system based on multidimensional slow waves and AI | |
| CN113658695B (en) | Fitness exercise data monitoring method, system and computer readable storage medium |