TWM574469U - QRS wave instant detection device - Google Patents

QRS wave instant detection device Download PDF

Info

Publication number
TWM574469U
TWM574469U TW107216587U TW107216587U TWM574469U TW M574469 U TWM574469 U TW M574469U TW 107216587 U TW107216587 U TW 107216587U TW 107216587 U TW107216587 U TW 107216587U TW M574469 U TWM574469 U TW M574469U
Authority
TW
Taiwan
Prior art keywords
value
ecg signal
signal
unit
detection device
Prior art date
Application number
TW107216587U
Other languages
Chinese (zh)
Inventor
林俊成
黃彥樺
Original Assignee
國立勤益科技大學
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 國立勤益科技大學 filed Critical 國立勤益科技大學
Priority to TW107216587U priority Critical patent/TWM574469U/en
Publication of TWM574469U publication Critical patent/TWM574469U/en

Links

Abstract

本新型QRS波即時檢測裝置,其藉由接收單元接收待檢測的心電訊號,並透過擷取單元藉由擷取窗格擷取該心電訊號,再配合處理單元消除基線飄移,以及轉換得到所有尺度的高頻訊號及低頻訊號,接著計算單元選用分析尺度,以及計算出正、負門檻值,並推算該心電訊號的對應座標,最後紀錄單元判斷該座標與前一座標間之距離並予以紀錄或刪除,藉以有效精準判斷QRS波的位置,不僅降低雜訊所造成的誤判,同時大幅提升判斷準確度、降低計算量,以及提升分析速度,進而有效減少人力資源的消耗,更加符合經濟效益。The new QRS wave real-time detection device receives a ECG signal to be detected through a receiving unit, and captures the ECG signal through an acquisition unit through an acquisition unit, and cooperates with a processing unit to eliminate baseline drift and conversion. The high-frequency signal and low-frequency signal of all scales, then the calculation unit selects the analysis scale, calculates the positive and negative threshold values, and calculates the corresponding coordinates of the ECG signal. Finally, the recording unit judges the distance between the coordinates and the previous one. It can be recorded or deleted to effectively and accurately determine the position of the QRS wave, which not only reduces the misjudgment caused by noise, but also greatly improves the accuracy of judgment, reduces the calculation amount, and improves the analysis speed, thereby effectively reducing the consumption of human resources and more economically. benefit.

Description

QRS波即時檢測裝置QRS wave instant detection device

本新型是有關於一種檢測裝置,特別是一種QRS波即時檢測裝置。The present invention relates to a detection device, particularly a QRS wave instant detection device.

隨著生活水準與品質的提升,重大疾病的預防及檢測已經逐漸受到現代人所重視,同時疾病的型態也由急性轉為慢性居多,例如癌症、心臟疾病、腦性血管疾病等,受助於科技的發展,藉由量測相關的生理訊號並進行分析,使得重大疾病的早期徵兆得以被察覺,以利患者及早接受治療,同時有助於控制病情並提升治癒率。With the improvement of living standards and quality, the prevention and detection of major diseases have gradually been valued by modern people. At the same time, the type of disease has also changed from acute to chronic, such as cancer, heart disease, and cerebrovascular disease. With the development of science and technology, by measuring and analyzing relevant physiological signals, early signs of major diseases can be detected, so that patients can be treated early, and it can help control the disease and improve the cure rate.

查,心臟疾病普遍利用心電圖進行病情的追蹤、分析及診斷,其採用非侵入性的量測方式,並提供良好的訊息,使得醫療人員或研究人員得以藉由心電圖訊號所呈現的心臟整體的電位變化,進而清楚地瞭解患者的心臟狀況,以及檢測出有潛在危險的患者,亦或是幫助監測心臟有缺陷的患者,以便降低患者心臟疾病發作的風險,而在進行心電圖訊號的判斷及分類之前,醫療人員或研究人員必須要在大量的心電圖訊號中找出QRS波型,進而藉由P、Q、R、S、T波的位置、大小及形狀做為診斷心臟疾病的依據。It is commonly used to track, analyze, and diagnose heart disease in cardiac diseases. It uses non-invasive measurement methods and provides good information, so that medical personnel or researchers can use the electrocardiogram signal to show the overall potential of the heart. Changes to further understand the patient ’s heart condition and detect potentially dangerous patients, or to help monitor patients with heart defects in order to reduce the risk of heart attacks in patients, before judging and classifying ECG signals Medical personnel or researchers must find out the QRS wave shape in a large number of ECG signals, and then use the position, size, and shape of P, Q, R, S, and T waves as the basis for diagnosing heart disease.

惟,為求提升心電圖訊號分析判斷的精準性,其通常需要花費相當長的時間進行量測,藉以收集足夠的訊號進行分析,同時提高紀錄到病症發作的成功率,所以大量的心電圖訊號造成分析判讀作業需要耗費很長的時數,且其需經由專業的醫療人員或是研究人員來進行判讀,故在人力資源的消耗是相當龐大的,因此,如何快速且精準地判斷心電圖訊號中的QRS波,藉以讓專業的醫療人員或是研究人員快速評估及診斷病人的病症,進而有效減少人力資源的消耗,以及提升整體分析的速度,仍是本領域技術人員努力的目標之一。However, in order to improve the accuracy of ECG signal analysis and judgment, it usually takes a considerable time to measure, so as to collect enough signals for analysis, and at the same time improve the success rate of records to the onset of the disease, so a large number of ECG signals cause analysis The interpretation operation takes a long time, and it needs to be interpreted by professional medical personnel or researchers, so the consumption of human resources is quite huge. Therefore, how to quickly and accurately determine the QRS in the ECG signal Wave, so that professional medical personnel or researchers can quickly evaluate and diagnose the patient's disease, thereby effectively reducing the consumption of human resources and improving the speed of overall analysis, is still one of the goals of those skilled in the art.

因此,本新型之目的,是在提供一種QRS波即時檢測裝置,其可大幅提升判斷準確度、降低計算量及提升分析速度,以達到病症的快速評估及診斷,進而有效減少人力資源的消耗。Therefore, the purpose of the present invention is to provide a QRS wave real-time detection device, which can greatly improve the accuracy of judgment, reduce the amount of calculation, and increase the speed of analysis, so as to achieve rapid assessment and diagnosis of diseases, and thereby effectively reduce the consumption of human resources.

於是,本新型QRS波即時檢測裝置,其包含有接收單元、擷取單元、處理單元、計算單元,以及紀錄單元;其中,該處理單元包含有濾波器與轉換器,而該紀錄單元包含有判斷組與紀錄組,藉由該接收單元接收待檢測之心跳訊號,接著該擷取單元使二擷取窗格重疊並透過該擷取窗格擷取該心電訊號,再配合該濾波器消除該心電訊號之基線飄移訊號,並藉由該轉換器轉換該心電訊號得到所有尺度的高頻訊號與低頻訊號,並將波型轉換為正、負特徵極值點,然後該計算單元選用其一尺度計算出正、負門檻值,並依據該正、負門檻值標示出特徵極值對,再推算出該心電訊號在時間領域對應的座標,最後配合該判斷組進一步依據該等特徵極值對間的距離判斷,並將符合的特徵極值對傳送至該紀錄組進行紀錄,如此不僅提升判斷準確度、降低計算量及提升分析速度,以達到病症的快速評估及診斷,同時有效減少人力資源的消耗,更加符合經濟效益。Therefore, the new QRS wave real-time detection device includes a receiving unit, a capturing unit, a processing unit, a calculation unit, and a recording unit; wherein the processing unit includes a filter and a converter, and the recording unit includes a judgment Group and record group, the heartbeat signal to be detected is received by the receiving unit, then the acquisition unit overlaps two acquisition panes and acquires the ECG signal through the acquisition pane, and cooperates with the filter to eliminate the The baseline drift signal of the ECG signal, and the ECG signal is converted by the converter to obtain high-frequency signals and low-frequency signals of all scales, and the waveform is converted into positive and negative characteristic extreme points, and then the calculation unit selects it. Calculate the positive and negative threshold values on a scale, mark the characteristic extreme value pairs according to the positive and negative threshold values, and then calculate the coordinates of the ECG signal in the time domain. Finally, cooperate with the judgment group to further rely on the characteristic poles. The distance between value pairs is judged, and the matching characteristic extreme value pairs are transmitted to the record group for recording, which not only improves the accuracy of judgment, reduces the calculation amount, and improves the analysis. Degree in order to achieve rapid assessment and diagnosis of disease, while effectively reducing the consumption of human resources, more cost-effective.

有關本新型之前述及其他技術內容、特點與功效,在以下配合參考圖式之較佳實施例的詳細說明中,將可清楚的明白。The foregoing and other technical contents, features, and effects of the present invention will be clearly understood in the following detailed description of the preferred embodiments with reference to the drawings.

參閱圖1,本新型之第一較佳實施例,該QRS波即時檢測裝置3包含有一與外部心電訊號量測裝置4或心電訊號資料庫5連接之接收單元31,一連接該接收單元31之擷取單元32,一連接該擷取單元32之處理單元33,一連接該處理單元33之計算單元34,以及一連接該計算單元34之紀錄單元35;其中,該處理單元33包含有一濾波器331,以及一連接該濾波器331之轉換器332,另外,該紀錄單元35包含有一判斷組351,以及一與該判斷組351連接之紀錄組352。Referring to FIG. 1, a first preferred embodiment of the present invention. The QRS wave real-time detection device 3 includes a receiving unit 31 connected to an external ECG signal measuring device 4 or an ECG signal database 5, and a receiving unit connected to the receiving unit 31. An extraction unit 32 of 31, a processing unit 33 connected to the extraction unit 32, a calculation unit 34 connected to the processing unit 33, and a recording unit 35 connected to the calculation unit 34; wherein the processing unit 33 includes a The filter 331 and a converter 332 connected to the filter 331. In addition, the recording unit 35 includes a judgment group 351 and a recording group 352 connected to the judgment group 351.

參閱圖1及圖2,進行QRS波檢測作業時,首先該接收單元31由該心電訊號量測裝置4或該心電訊號資料庫5接收一待檢測的心電訊號S,而該心電訊號量測裝置4可直接穿戴於患者身上,藉以直接接收患者的心電訊號S並進行後續分析處理,接著該擷取單元32透過小波轉換公式如下列:Referring to FIG. 1 and FIG. 2, when performing a QRS wave detection operation, first, the receiving unit 31 receives an ECG signal S to be detected by the ECG signal measuring device 4 or the ECG signal database 5, and the ECG signal The measurement device 4 can be directly worn on the patient, so as to directly receive the patient's ECG signal S and perform subsequent analysis and processing. Then, the acquisition unit 32 uses the wavelet conversion formula as follows:

利用母小波Ψ(t-b/a)中的係數a可以讓母小波伸縮,而係數b可以讓母小波平移,該擷取單元32進而透過該係數a、b調整讓母小波按照輸入的參數伸縮與平移即可得到類似可以自由調整的擷取窗格,參閱圖3,並使任二擷取窗格間重疊複數個取樣點,不同數量的取樣點重疊亦會影響後續偵測正確性,而本實施例中係以任二擷取窗格間重疊150個取樣點且每個擷取窗格為4096個取樣點為例加以說明,該擷取單元32接著藉由該擷取窗格擷取該待檢測之心電訊號S。The coefficient a in the mother wavelet Ψ (tb / a) can be used to expand and contract the mother wavelet, and the coefficient b can be used to translate the mother wavelet. The acquisition unit 32 further adjusts the mother wavelet to expand and contract according to the input parameters through the coefficients a and b. Panning can obtain similar acquisition panes that can be adjusted freely, see Figure 3, and make multiple sampling points overlap between any two acquisition panes. The overlapping of different numbers of sampling points will also affect the subsequent detection accuracy. In the embodiment, an example is described in which 150 sampling points overlap between any two acquisition panes and each acquisition pane has 4096 sampling points. The acquisition unit 32 then acquires the acquisition panes through the acquisition panes. ECG signal S to be detected.

配合參閱圖4,接著該處理單元33藉由該濾波器331針對該擷取之心電訊號S進行濾波,而本實施例中係使用200ms中值濾波器331與600ms中值濾波器331進行濾波處理,該200ms中值濾波器331首先針對該擷取之心電訊號S進行初步濾波處理,藉以濾除QRS波及P波,接著該600ms中值濾波器331再針對經過該200ms中值濾波器331濾波之心電訊號S進行濾波處理,藉以濾除T波,並得到不包含QRS波、P波及T波的基線飄移訊號,而該處理單元33再將該心電訊號S減去該基線飄移訊號,即得到修正之心電訊號S,即如圖5-2所示,而本實施例中,該濾波器331不包含消除高頻雜訊之動作,藉以避免該心電訊號S於消除高頻雜訊的過程中,連同頻率較高的QRS波型一起被影響而降低特徵值,然後該轉換器332進行該修正之心電訊號S的轉換作業,其採用的雙正交樣條小波的母小波為下列公式:With reference to FIG. 4, the processing unit 33 then filters the captured ECG signal S through the filter 331. In this embodiment, the 200 ms median filter 331 and the 600 ms median filter 331 are used for filtering. Processing, the 200ms median filter 331 first performs preliminary filtering processing on the captured ECG signal S, thereby filtering out QRS waves and P waves, then the 600ms median filter 331 and then the 200ms median filter 331 The filtered ECG signal S is filtered to remove the T wave and obtain a baseline drift signal that does not include QRS, P and T waves, and the processing unit 33 subtracts the baseline drift signal from the ECG signal S. That is, the corrected ECG signal S is obtained, as shown in FIG. 5-2. In this embodiment, the filter 331 does not include the action of eliminating high frequency noise, so as to avoid the ECG signal S from eliminating high frequency. In the process of noise, the high frequency QRS waveform is affected to reduce the eigenvalue, and then the converter 332 performs the conversion operation of the modified ECG signal S. The mother of the bi-orthogonal spline wavelet is used. The wavelet is the following formula:

該轉換器332藉由上述公式轉換該修正之心電訊號S,並得到所有尺度的高頻訊號D與低頻訊號A,而該轉換器332轉換尺度的數量可依實際需求與情況進行調整,而本實施例中係以該轉換器332透過小波轉換得到四個尺度的高頻訊號D1、D2、D3、D4及低頻訊號A1、A2、A3、A4,即如圖6至圖9所示,同時該轉換器332將該等高頻訊號D1、D2、D3、D4中之波型皆轉換為正、負特徵極值點,即如圖10所示,同時本實施例中,該轉換器332係將該等高頻訊號D1、D2、D3、D4之擷取窗格前50個取樣點及後50個取樣點予以去除,不納入檢測運算。The converter 332 converts the modified ECG signal S by the above formula, and obtains the high-frequency signal D and low-frequency signal A of all scales, and the number of scales converted by the converter 332 can be adjusted according to actual needs and conditions, and In this embodiment, the converter 332 is used to obtain high-frequency signals D1, D2, D3, and D4 and low-frequency signals A1, A2, A3, and A4 through wavelet transformation, as shown in FIGS. 6 to 9. The converter 332 converts the waveforms of the high-frequency signals D1, D2, D3, and D4 into positive and negative characteristic extreme points, as shown in FIG. 10. In the embodiment, the converter 332 is The first 50 sampling points and the last 50 sampling points of the acquisition pane of these high-frequency signals D1, D2, D3, and D4 are removed, and they are not included in the detection calculation.

仍續前述,該計算單元34選用其一尺度進行計算,而本實施例中,該計算單元34係選用第三尺度進行計算為例加以說明,該計算單元34先將該第三尺度之高頻訊號D3分成四段,並在該四段訊號內分別取正的最大值與負的最小值,然後將四個最大值取平均數得到平均最大值,以及將四個最小值取平均數得到平均最小值,接著取平均最大值的四分之一的數值,即得到正門檻值,另取平均最小值的四分之一的數值,即得到負門檻值,參閱圖11及圖12,該計算單元34進一步標示出一超過該正門檻值之正特徵極值點,以及標示出一與該正特徵極值點極性相反且超過該負門檻值之負特徵極值點,以形成特徵極值對,且該正、負特徵極值點間之距離不超過預設值,而本實施例中,該預設值為45個取樣點,當距離超過45個取樣點則表示波型較圓滑,不太符合QRS波型高聳尖銳的特徵,接著該計算單元34將正特徵極值點座標(X2, Y2)、負特徵極值點座標 (X1, Y1),以及將該心電訊號S做小波轉換後對時間領域超前或落後的量值帶入一推算公式,而本實施例中,該值依不同使用尺度帶入,其於該第二尺度為5,該第三尺度為10,該第四尺度為25,而該推算公式如下列:Continuing from the foregoing, the calculation unit 34 selects one of the scales for calculation. In this embodiment, the calculation unit 34 selects the third scale for calculation. For example, the calculation unit 34 first selects the high frequency of the third scale. The signal D3 is divided into four segments, and the positive maximum value and the negative minimum value are respectively taken in the four segment signals, and then the four maximum values are averaged to obtain the average maximum value, and the four minimum values are averaged to obtain the average value. The minimum value, then a quarter of the average maximum value, to obtain the positive threshold value, and another quarter value of the average minimum value, to obtain the negative threshold value, see Figure 11 and Figure 12, the calculation Unit 34 further marks a positive characteristic extreme point that exceeds the positive threshold value, and a negative characteristic extreme point that has the opposite polarity to the positive characteristic extreme point and exceeds the negative threshold value to form a characteristic extreme value pair. And the distance between the positive and negative characteristic extreme points does not exceed a preset value. In this embodiment, the preset value is 45 sampling points. When the distance exceeds 45 sampling points, it means that the waveform is smoother. Too high QRS wave shape Then, the calculation unit 34 sets the positive characteristic extreme point coordinates (X2, Y2), the negative characteristic extreme point coordinates (X1, Y1), and performs wavelet transformation of the ECG signal S to lead or lag the time domain. The magnitude of is brought into an estimation formula, and in this embodiment, the value is brought in according to different usage scales, which is 5 in the second scale, 10 in the third scale, 25 in the fourth scale, and the calculation The formula is as follows:

該計算單元34進而藉由上述推算公式推算出該心電訊號S在時間領域對應的座標P(n)。The calculation unit 34 further calculates the coordinates P (n) of the ECG signal S in the time domain by using the above-mentioned estimation formula.

仍續前述,該紀錄單元35之判斷組351將前述座標P(n)與前一特徵極值對所推算出之座標P(n-1)距離做比較,因年紀、性別、生活習慣等不同因素的影響,該距離的判斷標準值亦會有所不同,不應以此為限,當該距離小於第一數值,而本實施例中該第一數值為100個取樣點,該判斷組351取該座標P(n)與P(n-1)波型較尖銳者,亦即該判斷組351將該P(n)與P(n-1)個別的正特徵極值點座標(X2, Y2)及負特徵極值點座標(X1, Y1)帶入一判斷公式,該判斷公式如下列:Continuing the foregoing, the judging group 351 of the recording unit 35 compares the distance between the aforementioned coordinate P (n) and the calculated coordinate P (n-1) of the previous characteristic extreme value pair, due to differences in age, gender, living habits, etc. The judgment standard value of the distance will also be different due to factors, and should not be limited to this. When the distance is smaller than the first value, and the first value is 100 sampling points in this embodiment, the judgment group 351 Take the coordinates P (n) and P (n-1) with sharper waveforms, that is, the judgment group 351 sets the individual positive characteristic extreme point coordinates of the P (n) and P (n-1) coordinates (X2, Y2) and negative characteristic extreme point coordinates (X1, Y1) are brought into a judgment formula, which is as follows:

該判斷組351進而透過該判斷公式分別計算該P(n)及P(n-1)後取結果值(α)較大者(即為R波波峰),並將結果值較大的特徵極值對傳送至該紀錄組352進行記錄,配合參閱圖13,框起處出現兩對特徵極值對,該判斷組351即透過該判斷公式進行計算並取結果值較大者,亦即R波波峰,再傳送至該紀錄組352進行記錄;當該距離小於第二數值且大於該第一數值,代表受測者心率較快,亦或該心電訊號S中含有較陡峭的T波或P波,必須進一步判斷,而本實施例中該第二數值為130個取樣點,是以,該判斷組351加入另一尺度之高頻訊號D協助判斷,而本實施例中,該判斷組351加入第二尺度之高頻訊號D2進行判斷,當該第三尺度之該特徵極值對亦出現在該第二尺度之高頻訊號D2內且位置相同,則該判斷組351將該P(n)座標位置傳送至該紀錄組352進行記錄,亦即當座標P(n-1)與P(n)皆出現在該第三尺度與該第二尺度的高頻訊號D3、D2當中,則代表該二座標P(n-1)、P(n)皆為R波波峰,因此該二座標P(n-1)、 P(n)皆須傳送至該紀錄組352紀錄,反之,未同時出現在該第三尺度與該第二尺度高頻訊號D3、D2中的座標代表不是R波波峰,則不傳送至該紀錄組352紀錄,參閱圖14,該第三尺度高頻訊號D3框起處之特徵極值對亦出現在該第二尺度高頻訊號D2中,則該判斷組351將其傳送至該紀錄組352加入紀錄,另外參閱圖15,圖15所示之框起處在該第三尺度的高頻訊號D3出現特徵極值對,然而在該第二尺度的高頻訊號D2並沒有該特徵極值對,則該判斷組351不將其傳送至該紀錄組352,亦即不予以紀錄;當該距離大於該第二數值,亦即大於130個取樣點,則該判斷組351直接將其傳送至該紀錄組352進行紀錄,如此即可得到所有QRS波的位置,不僅減少檢測時的干擾,有效增加檢測的準確度,進而快速且準確做到QRS波的定位,同時大大提升判斷準確度、降低整體計算量,以及有效提升分析速度,更能減少人力資源的消耗,以及符合經濟效益。The judging group 351 further calculates the P (n) and P (n-1) by using the judging formula, and then takes the one with the larger result value (α) (that is, the R wave peak), and sets the characteristic pole with the larger result value. The value pair is transmitted to the record group 352 for recording. With reference to FIG. 13, two characteristic extreme value pairs appear at the frame. The judgment group 351 calculates through the judgment formula and takes the larger value, that is, the R wave. The wave peak is transmitted to the record group 352 for recording; when the distance is less than the second value and greater than the first value, it means that the subject's heart rate is faster, or the ECG signal S contains a steeper T wave or P It is necessary to further judge the wave. In this embodiment, the second value is 130 sampling points. Therefore, the judging group 351 adds a high-frequency signal D of another scale to assist in judging. In this embodiment, the judging group 351 The high-frequency signal D2 of the second scale is added for judgment. When the characteristic extreme value pair of the third scale also appears in the high-frequency signal D2 of the second scale and has the same position, the judgment group 351 sets the P (n ) The coordinate position is transmitted to the record group 352 for recording, that is, when the coordinates P (n-1) and P (n) both appear Among the high-frequency signals D3 and D2 of the third scale and the second scale, it means that the two coordinates P (n-1) and P (n) are both R wave peaks, so the two coordinates P (n-1 ), P (n) must be transmitted to the record group 352 record. Conversely, if the coordinates that do not appear in the third scale and the second scale high frequency signals D3 and D2 are not R wave peaks, they are not transmitted to The record group 352 records, refer to FIG. 14. The characteristic extreme pair at the beginning of the frame of the third-scale high-frequency signal D3 also appears in the second-scale high-frequency signal D2, and then the judgment group 351 transmits it to the record. Group 352 is added to the record, and referring to FIG. 15, the frame shown in FIG. 15 has characteristic extreme value pairs at the high-frequency signal D3 at the third scale, but the high-frequency signal D2 at the second scale does not have the characteristic pole. Value pair, the judgment group 351 does not transmit it to the recording group 352, that is, it does not record it; when the distance is greater than the second value, that is, greater than 130 sampling points, the judgment group 351 directly transmits it Go to the record group 352 for recording. In this way, the positions of all QRS waves can be obtained, which not only reduces interference during detection, but also effectively increases Measurement accuracy, and then quickly and accurately locate the QRS complex to do, while greatly enhance the accuracy of judgment, reducing the overall amount of calculation, and effectively enhance the speed of analysis, and help reduce the consumption of human resources, as well as cost-effective.

本新型採用MIT-BIH心律不整資料庫的48筆心電圖紀錄,總共109,488個QRS波,進一步測試本案所提出的QRS波即時檢測裝置3,並採用下列算式計算本案所提之QRS波即時檢測裝置3的靈敏度(Se)、陽性預測值(+P)及錯誤率(Der)等數值:The new model uses 48 ECG records of the MIT-BIH arrhythmia database, a total of 109,488 QRS waves, to further test the QRS wave real-time detection device 3 proposed in this case, and uses the following formula to calculate the QRS wave real-time detection device 3 proposed in this case Sensitivity (Se), positive predictive value (+ P), and error rate (Der):

TP:正確的QRS波總數TP: Total number of correct QRS waves

FP:非QRS波卻被偵測到的總數FP: Total number of non-QRS waves detected

FN:是QRS波卻沒被偵測到的總數FN: The total number of QRS waves that were not detected

FP+FN:總錯誤數量FP + FN: Total number of errors

參下列表1,實驗結果顯示QRS波的偵測正確性可達到99.84%的靈敏度和99.91%的陽性預測值,錯誤率僅為0.2471%,因此,本案所提出之QRS波即時檢測裝置3具有相當的準確度且也有實際應用的價值。See Table 1 below. The experimental results show that the QRS wave detection accuracy can reach 99.84% sensitivity and 99.91% positive predictive value, and the error rate is only 0.2471%. Therefore, the QRS wave real-time detection device 3 proposed in this case has comparable The accuracy is also of practical value.

表1 編號 TP FP FN FP+FN 編號 TP FP FN FP+FN 100 2,273 0 0 0 201 1,963 6 2 8 101 1,865 0 0 0 202 2,136 0 3 3 102 2,187 0 0 0 203 2,980 13 33 46 103 2,084 0 0 0 205 2,656 0 6 6 104 2,229 1 1 2 207 1,860 1 18 19 105 2,572 27 4 31 208 2,955 2 16 18 106 2,027 0 5 5 209 3,005 0 2 2 107 2,137 0 1 1 210 2,650 4 8 12 108 1,763 10 13 23 212 2,748 0 0 0 109 2,532 0 0 0 213 3,251 0 0 0 111 2,124 0 1 1 214 2,262 1 3 4 112 2,539 0 0 0 215 3,363 0 1 1 113 1,789 6 0 6 217 2,208 3 7 10 114 1,879 1 8 9 219 2,154 0 0 0 115 1,953 0 0 0 220 2,048 0 0 0 116 2,412 0 20 20 221 2,427 0 2 2 117 1,535 0 0 0 222 2,483 1 2 3 118 2,278 0 0 0 223 2,605 0 3 3 119 1,987 0 0 0 228 2,053 8 13 21 121 1,863 0 1 1 230 2,256 1 0 1 122 2,476 0 0 0 231 1,571 0 0 0 123 1,518 0 0 0 232 1,780 4 0 4 124 1,619 0 0 0 233 3,079 0 2 2 200 2,601 2 5 7 234 2,753 0 0 0 Table 1 Numbering TP FP FN FP + FN Numbering TP FP FN FP + FN 100 2,273 0 0 0 201 1,963 6 2 8 101 1,865 0 0 0 202 2,136 0 3 3 102 2,187 0 0 0 203 2,980 13 33 46 103 2,084 0 0 0 205 2,656 0 6 6 104 2,229 1 1 2 207 1,860 1 18 19 105 2,572 27 4 31 208 2,955 2 16 18 106 2,027 0 5 5 209 3,005 0 2 2 107 2,137 0 1 1 210 2,650 4 8 12 108 1,763 10 13 twenty three 212 2,748 0 0 0 109 2,532 0 0 0 213 3,251 0 0 0 111 2,124 0 1 1 214 2,262 1 3 4 112 2,539 0 0 0 215 3,363 0 1 1 113 1,789 6 0 6 217 2,208 3 7 10 114 1,879 1 8 9 219 2,154 0 0 0 115 1,953 0 0 0 220 2,048 0 0 0 116 2,412 0 20 20 221 2,427 0 2 2 117 1,535 0 0 0 222 2,483 1 2 3 118 2,278 0 0 0 223 2,605 0 3 3 119 1,987 0 0 0 228 2,053 8 13 twenty one 121 1,863 0 1 1 230 2,256 1 0 1 122 2,476 0 0 0 231 1,571 0 0 0 123 1,518 0 0 0 232 1,780 4 0 4 124 1,619 0 0 0 233 3,079 0 2 2 200 2,601 2 5 7 234 2,753 0 0 0

參閱圖16及圖17所示,本新型之第二較佳實施例,其仍包含有前一實施例所述之構件,且其目的功效均與前一實施例相同,在此不再贅述,而本實施例特別在於:該紀錄單元35另連接有一排除單元36,藉由該排除單元36透過前述判斷公式計算該每一特徵極值對之結果值(α),當該結果值小於所有結果值之平均值的五分之一,該排除單元36即將該點由紀錄中刪除,藉以有效提升偵測的正確性。Referring to FIG. 16 and FIG. 17, the second preferred embodiment of the present invention still includes the components described in the previous embodiment, and its purpose and effect are the same as those in the previous embodiment, and will not be repeated here. This embodiment is particularly that the recording unit 35 is further connected with an exclusion unit 36, and the result unit (α) of each characteristic extreme value pair is calculated by the exclusion unit 36 through the foregoing judgment formula. When the result value is less than all results One-fifth of the average value of the value, the exclusion unit 36 deletes the point from the record, thereby effectively improving the accuracy of detection.

歸納前述,本新型QRS波即時檢測裝置,其藉由該接收單元接收待檢測的心電訊號,進一步由該擷取單元調整擷取窗格並透過該擷取窗格擷取該心電訊號後,再配合該濾波器消除基線飄移,以及該轉換器轉換得到所有尺度的高頻訊號及低頻訊號,接著該計算單元選用分析尺度,以及計算出正、負門檻值,並推算該心電訊號的對應座標,最後該判斷組依據該座標與前一點座標間之距離判斷是否傳送至該紀錄組進行紀錄,如此即完成QRS波的檢測,不僅降低雜訊以及大振幅的P波與T波所造成的誤判,同時大幅提升判斷準確度、降低計算量,以及提升分析速度,進而有效減少人力資源的消耗,更加符合經濟效益。In summary, the new QRS wave real-time detection device receives the ECG signal to be detected by the receiving unit, further adjusts the acquisition pane by the acquisition unit, and acquires the ECG signal through the acquisition pane. , Together with the filter to eliminate baseline drift, and the converter converts high-frequency signals and low-frequency signals of all scales, then the calculation unit selects the analysis scale, calculates the positive and negative thresholds, and calculates the ECG signal. Corresponds to the coordinates. Finally, the judgment group judges whether it is transmitted to the record group for recording according to the distance between the coordinate and the previous point. In this way, the QRS wave detection is completed, which not only reduces the noise and the large amplitude of P waves and T waves. At the same time, it greatly improves the accuracy of judgment, reduces the amount of calculation, and increases the speed of analysis, which effectively reduces the consumption of human resources and is more in line with economic benefits.

惟以上所述者,僅為說明本新型之較佳實施例而已,當不能以此限定本新型實施之範圍,即大凡依本新型申請專利範圍及新型說明書內容所作之簡單的等效變化與修飾,皆應仍屬本新型專利涵蓋之範圍內。However, the above are only for explaining the preferred embodiments of the present invention. When the scope of the implementation of the present invention cannot be limited by this, that is, the simple equivalent changes and modifications made according to the scope of the patent application of the new application and the content of the new specification , All should still fall within the scope of this new patent.

(本創作)(This creation)

3‧‧‧QRS波即時檢測裝置 3‧‧‧QRS wave real-time detection device

31‧‧‧接收單元 31‧‧‧Receiving unit

32‧‧‧擷取單元 32‧‧‧Capture unit

33‧‧‧處理單元 33‧‧‧Processing unit

34‧‧‧計算單元 34‧‧‧ Computing Unit

35‧‧‧紀錄單元 35‧‧‧Record Unit

36‧‧‧排除單元 36‧‧‧Excluded units

331‧‧‧濾波器 331‧‧‧Filter

332‧‧‧轉換器 332‧‧‧ converter

351‧‧‧判斷組 351‧‧‧Judgment Team

352‧‧‧紀錄組 352‧‧‧Record Team

4‧‧‧心電訊號量測裝置 4‧‧‧ ECG signal measuring device

5‧‧‧心電訊號資料庫 5‧‧‧ ECG Signal Database

S‧‧‧心電訊號 S‧‧‧heart signal

A1‧‧‧第一尺度之低頻訊號 A1‧‧‧First-level low-frequency signal

A2‧‧‧第二尺度之低頻訊號 A2‧‧‧The second low frequency signal

A3‧‧‧第三尺度之低頻訊號 A3‧‧‧ Third-frequency low-frequency signal

A4‧‧‧第四尺度之低頻訊號 A4‧‧‧Fourth scale low frequency signal

D1‧‧‧第一尺度之高頻訊號 D1‧‧‧The first high-frequency signal

D2‧‧‧第二尺度之高頻訊號 D2‧‧‧Second-scale high-frequency signal

D3‧‧‧第三尺度之高頻訊號 D3‧‧‧ Third-frequency high-frequency signal

D4‧‧‧第四尺度之高頻訊號 D4‧‧‧Four-scale high-frequency signal

圖1是本新型第一較佳實施例之示意圖。 圖2是該第一較佳實施例之檢測流程圖。 圖3是調整擷取窗格之示意圖 圖4是該第一較佳實施例之局部流程圖。 圖5-1是基線飄移之心電訊號。 圖5-2是修正基線飄移之心電訊號。 圖6是小波多重分解示意圖。 圖7是小波一階轉換波型圖。 圖8是心電訊號經四階小波轉換後的低頻訊號。 圖9是心電訊號經四階小波轉換後的高頻訊號。 圖10是高頻訊號呈現正、負特徵極值點。 圖11是特徵極值對示意圖。 圖12是該第一較佳實施例之局部流程圖。 圖13是出現兩對特徵極值對則取波型較尖銳者。 圖14-15是加入第二尺度的高頻成分判斷。 圖16是本新型第二較佳實施例之示意圖。 圖17是該第二較佳實施例之檢測流程圖。FIG. 1 is a schematic diagram of a first preferred embodiment of the present invention. FIG. 2 is a detection flowchart of the first preferred embodiment. FIG. 3 is a schematic diagram of adjusting the capture pane. FIG. 4 is a partial flowchart of the first preferred embodiment. Figure 5-1 is the baseline ECG signal. Figure 5-2 is the ECG signal corrected for baseline drift. FIG. 6 is a schematic diagram of wavelet multiple decomposition. Fig. 7 is a wavelet first-order conversion pattern. Figure 8 shows the low-frequency signal of the ECG signal after fourth-order wavelet conversion. Figure 9 shows the high-frequency signal of the ECG signal after fourth-order wavelet conversion. Figure 10 shows the extreme points of high-frequency signals showing positive and negative characteristics. FIG. 11 is a schematic diagram of characteristic extreme value pairs. FIG. 12 is a partial flowchart of the first preferred embodiment. Figure 13 shows the sharper waveforms when two pairs of characteristic extreme pairs appear. Figure 14-15 shows the high-frequency component judgment added to the second scale. FIG. 16 is a schematic diagram of a second preferred embodiment of the present invention. FIG. 17 is a detection flowchart of the second preferred embodiment.

Claims (9)

一種QRS波即時檢測裝置,其包含有: 一接收單元,其與外部之心電訊號資料庫或心電訊號量測裝置連接,藉以接收一待檢測的心電訊號; 一擷取單元,其與該接收單元連接,而該擷取單元可透過小波轉換公式調整擷取窗格的位置,以使任二擷取窗格間重疊數個取樣點,並透過該等擷取窗格擷取該心電訊號; 一處理單元,其與該擷取單元連接,而該處理單元包含有一針對擷取之心電訊號進行濾波處理之濾波器,以及一連接該濾波器並針對經濾波之心電訊號進行轉換之轉換器;其中,該濾波器針對擷取之心電訊號進行濾波並得到基線飄移訊號,再將該心電訊號減去該基線飄移訊號,以得到修正之心電訊號,而該轉換器將修正之心電訊號轉換得到所有尺度的高頻訊號與低頻訊號,同時將該心電訊號中的波型轉換為正、負特徵極值點; 一計算單元,其與該處理單元連接,而該計算單元可選用前述其一尺度進行計算,以得到正門檻值及負門檻值,並依該正門檻值及該負門檻值為基礎,標示出超過該正門檻值之正特徵極值點,以及超過該負門檻值之負特徵極值點,且該正、負特徵極值點間之距離需小於預設值,以形成一特徵極值對,並推算出該心電訊號在時間領域對應的座標;以及 一紀錄單元,其包含有一判斷組,以及一與該判斷組連接之紀錄組;其中,該判斷組比較二特徵極值對所推算出之座標間的距離,當該距離小於第一數值時,該判斷組分別計算該等座標並將結果值較大者傳送至該紀錄組進行記錄;當該距離小於第二數值且大於該第一數值,該判斷組加入另一尺度之高頻訊號做判斷,若該特徵極值對亦出現在該另一尺度之高頻訊號內且位置相同,則將該特徵極值對傳送至該紀錄組進行記錄;反之,當該距離大於該第二數值,該判斷組將該特徵極值對傳送至該紀錄組進行紀錄。A QRS wave real-time detection device includes: a receiving unit, which is connected to an external ECG signal database or an ECG signal measurement device, so as to receive an ECG signal to be detected; an acquisition unit, and The receiving unit is connected, and the capturing unit can adjust the position of the capturing pane through a wavelet transformation formula, so that there are several sampling points overlapping between any two capturing panes, and the heart is captured through the capturing panes. Electrical signal; a processing unit connected to the acquisition unit, and the processing unit includes a filter for filtering the captured ECG signal, and a filter connected to the filter for ECG signal filtering A conversion converter; wherein, the filter filters the captured ECG signal and obtains a baseline drift signal, and then subtracts the baseline drift signal from the ECG signal to obtain a modified ECG signal, and the converter The modified ECG signal is converted into high-frequency signals and low-frequency signals of all scales, and the waveforms in the ECG signal are converted into positive and negative characteristic extreme points; a calculation unit, which The processing unit is connected, and the calculation unit may use one of the aforementioned scales to perform calculations to obtain a positive threshold value and a negative threshold value, and based on the positive threshold value and the negative threshold value, mark those that exceed the positive threshold value. Positive characteristic extreme points and negative characteristic extreme points that exceed the negative threshold, and the distance between the positive and negative characteristic extreme points must be less than a preset value to form a characteristic extreme value pair, and calculate the heart The coordinates corresponding to the time signal in the telecommunication signal; and a recording unit, which includes a judgment group and a recording group connected to the judgment group; wherein the judgment group compares the distances between the coordinates calculated by the two characteristic extreme values When the distance is less than the first value, the judgment group calculates the coordinates separately and transmits the larger value to the record group for recording; when the distance is less than the second value and greater than the first value, the judgment group Add a high-frequency signal of another scale for judgment. If the characteristic extreme value pair also appears in the high-frequency signal of the other scale and has the same position, then transmit the characteristic extreme value pair to the record group for recording; Conversely, when the distance is greater than the second value, the judgment group transmits the characteristic extreme value pair to the record group for recording. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該二擷取窗格間重疊150個取樣點。According to the QRS wave real-time detection device described in item 1 of the scope of the patent application, wherein the two acquisition panes overlap 150 sampling points. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該濾波器為中值濾波器。According to the QRS wave real-time detection device described in item 1 of the patent application scope, wherein the filter is a median filter. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該轉換器將該心電訊號轉換為四個尺度的高頻訊號與低頻訊號,且轉換完成的高頻訊號前、後各去除50個取樣點。According to the QRS wave real-time detection device described in item 1 of the patent application scope, wherein the converter converts the ECG signal into high-frequency signals and low-frequency signals of four scales, and the converted high-frequency signals are respectively front and rear. Remove 50 sampling points. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該計算單元選用第三尺度之高頻訊號進行分析計算。According to the QRS wave real-time detection device described in item 1 of the scope of the patent application, wherein the calculation unit selects a third-scale high-frequency signal for analysis and calculation. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該第一數值為100個取樣點,而該第二數值為130個取樣點。According to the QRS wave real-time detection device described in item 1 of the scope of patent application, the first value is 100 sampling points, and the second value is 130 sampling points. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該預設值為45個取樣點。According to the QRS wave real-time detection device described in item 1 of the scope of patent application, the preset value is 45 sampling points. 根據申請專利範圍第1項或第6項所述之QRS波即時檢測裝置,其中,該判斷組加入第二尺度之高頻訊號做判斷。According to the QRS wave real-time detection device described in item 1 or item 6 of the scope of patent application, wherein the judging group adds a high-frequency signal of the second scale for judgment. 根據申請專利範圍第1項所述之QRS波即時檢測裝置,其中,該紀錄單元另連接有一排除單元,該排除單元計算該每一特徵極值對之結果值,並依據該結果值進一步排除不符合的特徵極值對。According to the QRS wave real-time detection device described in item 1 of the scope of the patent application, wherein the recording unit is further connected with an exclusion unit, the exclusion unit calculates the result value of each characteristic extreme value pair, and further excludes the inaccuracies based on the result value. Consistent characteristic extreme pairs.
TW107216587U 2018-12-06 2018-12-06 QRS wave instant detection device TWM574469U (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW107216587U TWM574469U (en) 2018-12-06 2018-12-06 QRS wave instant detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW107216587U TWM574469U (en) 2018-12-06 2018-12-06 QRS wave instant detection device

Publications (1)

Publication Number Publication Date
TWM574469U true TWM574469U (en) 2019-02-21

Family

ID=66215121

Family Applications (1)

Application Number Title Priority Date Filing Date
TW107216587U TWM574469U (en) 2018-12-06 2018-12-06 QRS wave instant detection device

Country Status (1)

Country Link
TW (1) TWM574469U (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112022122A (en) * 2020-09-29 2020-12-04 深圳职业技术学院 Sleep monitoring earphone

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112022122A (en) * 2020-09-29 2020-12-04 深圳职业技术学院 Sleep monitoring earphone

Similar Documents

Publication Publication Date Title
CN110327036B (en) Method for extracting respiratory signal and respiratory frequency from wearable electrocardiogram
JPH1170089A (en) Cardiac waveform characteristics displaying
Saadi et al. Automatic real-time embedded QRS complex detection for a novel patch-type electrocardiogram recorder
US10912479B2 (en) Method for accurately extracting abnormal potential within QRS
WO2019100563A1 (en) Method for assessing electrocardiogram signal quality
WO2012092766A1 (en) Method and system for automated detection and analysis in pediatric electrocardiography
CN110123304B (en) Dynamic electrocardio noise filtering method based on multi-template matching and correlation coefficient matrix
CN110236508A (en) A kind of non-invasive blood pressure continuous monitoring method
US20200100697A1 (en) Methods and system for processing an emg signal
Bsoul et al. Detection of P, QRS, and T components of ECG using wavelet transformation
Liu et al. Systematic methods for fetal electrocardiographic analysis: Determining the fetal heart rate, RR interval and QT interval
TWI527560B (en) Method and apparatus for monitoring heartbeat signal based on empirical mode decomposition
Di Maria et al. An algorithm for the analysis of fetal ECGs from 4-channel non-invasive abdominal recordings
Tan et al. EMD-based electrocardiogram delineation for a wearable low-power ECG monitoring device
CN106236041A (en) A kind of measure in real time and accurately heart rate and the algorithm of breathing rate and system
CN111419219A (en) PPG heart beat signal preprocessing method and device and atrial fibrillation detection equipment
CN112001862B (en) Non-contact type apparent heart rate detection method for eliminating motion noise of video heart impact signal
Choudhary et al. Delineation and analysis of seismocardiographic systole and diastole profiles
TWM574469U (en) QRS wave instant detection device
CN110491504B (en) Method for acquiring medical index data of heart sound signal
TWI696191B (en) Algorithm of qrs detection capable of reducing noise effects
CN110584638A (en) Non-contact heart rate measurement method based on CMOR wavelet
Golpaygani et al. Detection and identification of S1 and S2 heart sounds using wavelet decomposition method
TWI672127B (en) Algorithm of qrs detection based on wavelet transformation
CN203898296U (en) Novel electrocardiogram detection device

Legal Events

Date Code Title Description
MM4K Annulment or lapse of a utility model due to non-payment of fees