CN109741827A - Prediction method of respiratory signal cycle in thoracic-abdominal surface area combined with double-cycle judgment - Google Patents

Prediction method of respiratory signal cycle in thoracic-abdominal surface area combined with double-cycle judgment Download PDF

Info

Publication number
CN109741827A
CN109741827A CN201910018163.5A CN201910018163A CN109741827A CN 109741827 A CN109741827 A CN 109741827A CN 201910018163 A CN201910018163 A CN 201910018163A CN 109741827 A CN109741827 A CN 109741827A
Authority
CN
China
Prior art keywords
signal
period
breathing
cycle
pause
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
CN201910018163.5A
Other languages
Chinese (zh)
Other versions
CN109741827B (en
Inventor
于晓洋
王明
韩岫君
樊琪
胡亚欣
徐雅婷
赵烟桥
张琴
王妍
王淼
朱子桐
李宇潇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201910018163.5A priority Critical patent/CN109741827B/en
Publication of CN109741827A publication Critical patent/CN109741827A/en
Application granted granted Critical
Publication of CN109741827B publication Critical patent/CN109741827B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明结合双周期判断的胸腹表面区域呼吸信号周期预测方法属于肿瘤医学、精密仪器、工程技术和数学等技术领域;该方法首先构建理想单周期呼吸信号,然后进行周期延拓,得到理想多周期呼吸信号,最后对理想多周期呼吸信号进行周期提取,整个提取过程在数据量相同的两个假想周期间进行,并在两个周期内即可完成周期提取;该方法方法能够准确提取理想多周期呼吸信号的周期,为准确预测呼吸运动奠定仿真实验基础。

The method for predicting the breathing signal period of the chest and abdomen surface area combined with the double-period judgment of the invention belongs to the technical fields of oncology medicine, precision instruments, engineering technology, mathematics, etc. Periodic breathing signal, and finally the ideal multi-period breathing signal is periodically extracted. The entire extraction process is performed between two imaginary cycles with the same amount of data, and the periodic extraction can be completed within two cycles; this method can accurately extract the ideal multi-cycle breathing signal. The period of the periodic breathing signal lays a simulation experiment foundation for accurately predicting the breathing movement.

Description

The chest and abdomen surface region breath signal period forecasting method judged in conjunction with binary cycle
Technical field
The chest and abdomen surface region breath signal period forecasting method of present invention combination binary cycle judgement belongs to medical oncology, essence The technical fields such as close instrument, engineering technology and mathematics.
Background technique
Radiotherapy (being hereinafter radiotherapy) is one of most important means for the treatment of cancer.About 70% cancer patient Radiotherapy can be used during treating cancer, there are about 40% cancers can cure by radiotherapy.The effect of radiation therapy with put The precision of radiation exposure to tumor region is related with dosage.It is intraperitoneal that alimentary system malignant tumour is located at thoracic cavity, by respiratory movement It influences, the position of tumour and volume can all change over time, and cause radioactive ray irradiation precision to reduce, to influence radiotherapy effect Fruit.
In order to solve the problems, such as that respiratory movement reduces radiotherapeutic effect, many scholars have carried out target area flared end, have held one's breath, breathe door Control, passive pressurization, the research work of the four-dimensional techniques or methods such as radiotherapy and real-time tracking.These work effectively reduce breathing fortune The dynamic influence to radiotherapeutic effect, however, since these work are all the lag compensation work carried out after respiratory movement generation, because This can not fundamentally overcome the problems, such as that respiratory movement reduces radiotherapeutic effect.
To solve the above-mentioned problems, some scholars begin one's study lead compensation method, and this method is transported firstly the need of to breathing It is dynamic to be predicted.Currently, respiratory movement predicting method can be divided into following two major classes: the first, model prediction method, it passes through monitoring Respiratory movement early period data establish respiratory movement model, then following respiratory movement is speculated based on model;The second, model-free Prediction technique, it passes through the respiratory movement mode observed and is trained, and didactic learning algorithm is formed, by learning come pre- Survey following respiratory movement.
Although respiratory movement constantly repeats air-breathing and expiration movement, respiratory movement with air-breathing and is not exhaled as one Period is simply repeated, and during each air-breathing and expiration, by itself and external influence, period and amplitude all can There is change, model prediction method can not look after this variation, therefore precision of prediction will receive restriction, so, more scholars open Beginning is attempted in model-free prediction technique, for example, predicting respiratory movement with Gaussian process recurrence.
It is a kind of pervasive prediction algorithm that Gaussian process, which returns, predicts Future Data by study given data, is predicted As a result it is provided in the form of mean value and variance.Since Gaussian process regression algorithm itself is not to predict for respiratory movement merely And it proposes, therefore do not take into account that respirometric some constraint conditions.If respirometric constraint condition taken into account, Estimation range will be reduced, prediction result also will be more accurate.
For respiratory movement, frequency can be limited within the scope of some, therefore can be by obtaining breath signal Period or frequency carry out constrained forecast as a result, improve precision of prediction.Be conducive to improve prediction as it can be seen that obtaining the respirometric period Precision.
However, since respiratory movement is not monocyclic simple repetition, be difficult with two same characteristic features values away from From judging the period.So, if obtaining the respirometric period has become the key technical problem for improving precision of prediction.
Summary of the invention
In order to predict the period of breath signal, the present invention starts with from ideal breath signal, proposes a kind of breath signal Periodicity extraction method, this method not only include according to the work of breathing rule building identical breaths signal, but also including to ideal The work that the breath signal period extracts.
The object of the present invention is achieved like this:
In conjunction with the chest and abdomen surface region breath signal period forecasting method that binary cycle judges, comprise the steps of:
Step a, ideal monocycle breath signal is constructed;
Step b, periodic extension is carried out to the ideal monocycle breath signal that step a is obtained, obtains ideal multicycle breathing Signal;
Step c, periodicity extraction is carried out to the ideal multicycle breath signal that step b is obtained, comprising:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
The chest and abdomen surface region breath signal period forecasting method of above-mentioned combination binary cycle judgement, building described in step a Ideal monocycle breath signal, comprising the following steps:
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal.
In conjunction with the chest and abdomen surface region breath signal period forecasting method that binary cycle judges, comprise the steps of:
Step a, ideal monocycle breath signal is constructed;
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal;
Step b, periodic extension is carried out to the ideal monocycle breath signal that step a is obtained, obtains ideal multicycle breathing Signal;
Step c, periodicity extraction is carried out to the ideal multicycle breath signal that step b is obtained, comprising:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
A kind of ideal monocycle breath signal construction method, comprising the following steps:
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal.
A kind of ideal multicycle breath signal periodicity extraction method, comprising the following steps:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
The utility model has the advantages that
Although first, the present invention is under the jurisdiction of medical oncology technical field, the present invention be pair from the technical point of view Signal characteristic quantity extracts, and singly cannot achieve and diagnose to disease from the point of view of technical purpose and effect of the invention And treatment, therefore it is not belonging to the diagnostic and therapeutic method of disease described in Patent Law Article 25, object is not present in the present invention Problem.
The second, the present invention discloses abundant enough, as long as the source code that specification part is recorded is transported in MATLAB software Row, can be obtained operation result, either those skilled in the art are also non-those skilled in the art, can realize this Shen Please.
Third, the present invention also provides a kind of ideal monocycle breath signal construction method, this method only needs to use three Section SIN function can fit the breath signal in a cycle, method is simple, function is simple, and fitting result with really exhale It is corresponding to inhale movement.
4th, the method for the present invention carries out during data volume identical two imaginations week, and with this research team Shen on the same day The patent of invention " the human body chest and abdomen surface region breath signal period forecasting method calculated based on variance " of report is compared, as long as one Periodicity extraction can be completed in a true period, therefore calculation amount is smaller, from simulation result it can also be seen that, runing time from 0.40s is increased to 0.14s.
5th, the method for the present invention carries out during data volume identical two imaginations week, declares on the same day with this research team Patent of invention " in conjunction with three periods judge chest and abdomen surface region breath signal period forecasting method ", it is contemplated that respiratory movement is exhaled Gas and breathing process are asymmetric, therefore data can be realized within two imaginary periods and be mutually authenticated, and pseudoperiodicity is being avoided to occur, Prediction result more accurately meanwhile, it is capable to runing time is further increased, from simulation result it can also be seen that, runing time 0.14s is increased to from 0.32s.
Detailed description of the invention
Fig. 1 is ideal multicycle breath signal figure.
Fig. 2 is program runnable interface.
Specific embodiment
The specific embodiment of the invention is described in further detail with reference to the accompanying drawing.
Specific embodiment one
The present embodiment is the chest and abdomen surface region breath signal period forecasting embodiment of the method in conjunction with binary cycle judgement.
The chest and abdomen surface region breath signal period forecasting method of the combination binary cycle judgement of the present embodiment, by following steps Composition:
Step a, ideal monocycle breath signal is constructed;
Step b, periodic extension is carried out to the ideal monocycle breath signal that step a is obtained, obtains ideal multicycle breathing Signal;
Step c, periodicity extraction is carried out to the ideal multicycle breath signal that step b is obtained, comprising:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
Specific embodiment two
The present embodiment is the chest and abdomen surface region breath signal period forecasting embodiment of the method in conjunction with binary cycle judgement.
The chest and abdomen surface region breath signal period forecasting method of the combination binary cycle judgement of the present embodiment, is being embodied On the basis of example one, the ideal monocycle breath signal of building described in step a is further limited, comprising the following steps:
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal.
Specific embodiment three
The present embodiment is the chest and abdomen surface region breath signal period forecasting embodiment of the method in conjunction with binary cycle judgement.
The chest and abdomen surface region breath signal period forecasting method of the combination binary cycle judgement of the present embodiment, by following steps Composition:
Step a, ideal monocycle breath signal is constructed;
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal;
Step b, periodic extension is carried out to the ideal monocycle breath signal that step a is obtained, obtains ideal multicycle breathing Signal;
Step c, periodicity extraction is carried out to the ideal multicycle breath signal that step b is obtained, comprising:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
Specific embodiment four
The present embodiment is a kind of ideal monocycle breath signal construction method embodiment.
A kind of ideal monocycle breath signal construction method of the present embodiment, comprising the following steps:
Step a1, according to the process of mankind's eupnea, respiratory movement is divided into breathing process, exhalation process and class and is suspended Process three phases;
Step a2, the duration T of breathing process is determined respectively1, exhalation process duration T2With class pause process Duration T3
It step a3, is 4T with the period1Phase is the sine function simulation breathing process of [0, pi/2], obtains simulation air-breathing letter Number;It is 4T with the period2Phase is the sine function simulation exhalation process of [pi/2, π], obtains simulation breath signal;It is 4T with the period3 Phase is that the sine function simulation class of [π, 3 pi/2s] suspends process, obtains simulation class pause signal;
Step a4, the relative magnitude for simulating inspiration signal is adjusted between [a, b], the opposite width of breath signal will be simulated Value is adjusted between [c, b], and the relative magnitude for simulating class pause signal is adjusted between [a, c];
Step a5, obtained simulation inspiration signal, simulation breath signal and simulation class pause signal is end to end, it obtains Ideal monocycle breath signal.
Specific embodiment five
The present embodiment is a kind of ideal multicycle breath signal periodicity extraction embodiment of the method.
A kind of ideal multicycle breath signal periodicity extraction method of the present embodiment, comprising the following steps:
Step c1, cycle data variable is set as 1, enters step c2;
Step c2, according to cycle data variable, the sum of period 1 all data and second round all data are calculated The sum of, enter step c3;
Step c3, judge whether the sum of period 1 all data are equal to the sum of second round all data, if:
It is to enter step c4;
It is no, cycle data variable+1, return step c2;
Step c4, cycle data variable is the period of multicycle breath signal.
Specific embodiment seven
The present embodiment is the chest and abdomen surface region breath signal period forecasting embodiment of the method in conjunction with binary cycle judgement.
In order to verify the method for the present invention, this method is run on MATLAB R2014a software.Applied computer is matched It sets as follows:
Allocation of computer table
Operating system 7 Ultimate X86 of Windows
Processor Intel(R)Core(TM)i5-8250OU@1.60GHz 1.80GHz
Memory is installed 8.00GB
The MATLAB program write is as follows:
Program operation result difference is as depicted in figs. 1 and 2.Wherein, Fig. 1 is the signal graph of ideal multicycle breath signal, Fig. 2 is program runnable interface.
Program operation result shows that the method for the present invention has accurately extracted the signal period from respiratory movement signal.
Finally, thank state natural sciences fund general project " in radiotherapy human body chest and abdomen surface dynamic 3 D measurement and Region breath Motion Estimation and space-time integrated are predicted " (project number: 61571168) and Students' Innovation project is " swollen towards lung The chest and abdomen surface respiratory movement forecasting research of tumor radiotherapy " (project number: 201810214267) to the fund branch of this patent It holds.

Claims (5)

1.结合双周期判断的胸腹表面区域呼吸信号周期预测方法,其特征在于,由以下步骤组成:1. the method for predicting the breathing signal period of the chest and abdomen surface area judged in conjunction with the double period, is characterized in that, is made up of the following steps: 步骤a、构建理想单周期呼吸信号;Step a, constructing an ideal single-cycle breathing signal; 步骤b、对步骤a得到的理想单周期呼吸信号进行周期延拓,得到理想多周期呼吸信号;Step b, performing periodic extension on the ideal single-cycle breathing signal obtained in step a, to obtain an ideal multi-cycle breathing signal; 步骤c、对步骤b得到的理想多周期呼吸信号进行周期提取,包括:Step c. Periodically extract the ideal multi-period breathing signal obtained in step b, including: 步骤c1、设定周期数据变量为1,进入步骤c2;Step c1, set the period data variable to 1, and enter step c2; 步骤c2、根据周期数据变量,计算第一周期所有数据之和,以及第二周期所有数据之和,进入步骤c3;Step c2, according to the period data variable, calculate the sum of all data in the first period and the sum of all data in the second period, and enter step c3; 步骤c3、判断第一周期所有数据之和是否等于第二周期所有数据之和,如果:Step c3, determine whether the sum of all data in the first cycle is equal to the sum of all data in the second cycle, if: 是,进入步骤c4;Yes, go to step c4; 否,周期数据变量+1,返回步骤c2;No, cycle data variable +1, return to step c2; 步骤c4、周期数据变量即为多周期呼吸信号的周期。In step c4, the period data variable is the period of the multi-period breathing signal. 2.根据权利要求1所述的结合双周期判断的胸腹表面区域呼吸信号周期预测方法,其特征在于,步骤a所述的构建理想单周期呼吸信号,包括以下步骤:2. the thoracic and abdominal surface area breathing signal cycle prediction method according to claim 1 is characterized in that, the described construction of step a ideal single-cycle breathing signal, comprises the following steps: 步骤a1、根据人类正常呼吸的过程,将呼吸运动分为吸气过程、呼气过程和类暂停过程三个阶段;Step a1, according to the process of normal human breathing, the breathing movement is divided into three stages: an inhalation process, an exhalation process and a pause-like process; 步骤a2、分别确定吸气过程的持续时间T1、呼气过程的持续时间T2和类暂停过程的持续时间T3Step a2, respectively determining the duration T 1 of the inhalation process, the duration T 2 of the exhalation process and the duration T 3 of the pause-like process; 步骤a3、用周期为4T1相位为[0,π/2]的正弦函数模拟吸气过程,得到模拟吸气信号;用周期为4T2相位为[π/2,π]的正弦函数模拟呼气过程,得到模拟呼气信号;用周期为4T3相位为[π,3π/2]的正弦函数模拟类暂停过程,得到模拟类暂停信号;Step a3, simulate the inhalation process with a sine function with a period of 4T 1 and a phase of [0, π/2] to obtain a simulated inhalation signal; use a sine function with a period of 4T and a phase of [π/ 2 , π] to simulate exhalation. Breathing process, get the simulated exhalation signal; use the sine function whose period is 4T and 3 phase is [π, 3π/2] to simulate the pause process, and get the simulation pause signal; 步骤a4、将模拟吸气信号的相对幅值调整到[a,b]之间,将模拟呼气信号的相对幅值调整到[c,b]之间,将模拟类暂停信号的相对幅值调整到[a,c]之间;Step a4: Adjust the relative amplitude of the analog inspiratory signal to be between [a, b], adjust the relative amplitude of the analog expiratory signal to be between [c, b], and adjust the relative amplitude of the analog pause signal Adjust to between [a, c]; 步骤a5、将得到的模拟吸气信号、模拟呼气信号和模拟类暂停信号首尾相接,得到理想单周期呼吸信号。Step a5: Connect the obtained simulated inhalation signal, simulated expiratory signal and simulated pause signal end to end to obtain an ideal single-cycle breathing signal. 3.结合双周期判断的胸腹表面区域呼吸信号周期预测方法,其特征在于,由以下步骤组成:3. The method for predicting the breathing signal cycle of the thoracic-abdominal surface area in conjunction with the double-cycle judgment is characterized in that, it is made up of the following steps: 步骤a、构建理想单周期呼吸信号;Step a, constructing an ideal single-cycle breathing signal; 步骤a1、根据人类正常呼吸的过程,将呼吸运动分为吸气过程、呼气过程和类暂停过程三个阶段;Step a1, according to the process of normal human breathing, the breathing movement is divided into three stages: an inhalation process, an exhalation process and a pause-like process; 步骤a2、分别确定吸气过程的持续时间T1、呼气过程的持续时间T2和类暂停过程的持续时间T3Step a2, respectively determining the duration T 1 of the inhalation process, the duration T 2 of the exhalation process and the duration T 3 of the pause-like process; 步骤a3、用周期为4T1相位为[0,π/2]的正弦函数模拟吸气过程,得到模拟吸气信号;用周期为4T2相位为[π/2,π]的正弦函数模拟呼气过程,得到模拟呼气信号;用周期为4T3相位为[π,3π/2]的正弦函数模拟类暂停过程,得到模拟类暂停信号;Step a3, simulate the inhalation process with a sine function with a period of 4T 1 and a phase of [0, π/2] to obtain a simulated inhalation signal; use a sine function with a period of 4T and a phase of [π/ 2 , π] to simulate exhalation. Breathing process, get the simulated exhalation signal; use the sine function whose period is 4T and 3 phase is [π, 3π/2] to simulate the pause process, and get the simulation pause signal; 步骤a4、将模拟吸气信号的相对幅值调整到[a,b]之间,将模拟呼气信号的相对幅值调整到[c,b]之间,将模拟类暂停信号的相对幅值调整到[a,c]之间;Step a4: Adjust the relative amplitude of the analog inspiratory signal to be between [a, b], adjust the relative amplitude of the analog expiratory signal to be between [c, b], and adjust the relative amplitude of the analog pause signal Adjust to between [a, c]; 步骤a5、将得到的模拟吸气信号、模拟呼气信号和模拟类暂停信号首尾相接,得到理想单周期呼吸信号;Step a5, connecting the obtained simulated inhalation signal, simulated expiratory signal and simulated pause signal end to end to obtain an ideal single-cycle breathing signal; 步骤b、对步骤a得到的理想单周期呼吸信号进行周期延拓,得到理想多周期呼吸信号;Step b, performing periodic extension on the ideal single-cycle breathing signal obtained in step a, to obtain an ideal multi-cycle breathing signal; 步骤c、对步骤b得到的理想多周期呼吸信号进行周期提取,包括:Step c. Periodically extract the ideal multi-period breathing signal obtained in step b, including: 步骤c1、设定周期数据变量为1,进入步骤c2;Step c1, set the period data variable to 1, and enter step c2; 步骤c2、根据周期数据变量,计算第一周期所有数据之和,以及第二周期所有数据之和,进入步骤c3;Step c2, according to the period data variable, calculate the sum of all data in the first period and the sum of all data in the second period, and enter step c3; 步骤c3、判断第一周期所有数据之和是否等于第二周期所有数据之和,如果:Step c3, determine whether the sum of all data in the first cycle is equal to the sum of all data in the second cycle, if: 是,进入步骤c4;Yes, go to step c4; 否,周期数据变量+1,返回步骤c2;No, cycle data variable +1, return to step c2; 步骤c4、周期数据变量即为多周期呼吸信号的周期。In step c4, the period data variable is the period of the multi-period breathing signal. 4.一种理想单周期呼吸信号构建方法,其特征在于,包括以下步骤:4. an ideal single-cycle breathing signal construction method, is characterized in that, comprises the following steps: 步骤a1、根据人类正常呼吸的过程,将呼吸运动分为吸气过程、呼气过程和类暂停过程三个阶段;Step a1, according to the process of normal human breathing, the breathing movement is divided into three stages: an inhalation process, an exhalation process and a pause-like process; 步骤a2、分别确定吸气过程的持续时间T1、呼气过程的持续时间T2和类暂停过程的持续时间T3Step a2, respectively determining the duration T 1 of the inhalation process, the duration T 2 of the exhalation process and the duration T 3 of the pause-like process; 步骤a3、用周期为4T1相位为[0,π/2]的正弦函数模拟吸气过程,得到模拟吸气信号;用周期为4T2相位为[π/2,π]的正弦函数模拟呼气过程,得到模拟呼气信号;用周期为4T3相位为[π,3π/2]的正弦函数模拟类暂停过程,得到模拟类暂停信号;Step a3, simulate the inhalation process with a sine function with a period of 4T 1 and a phase of [0, π/2] to obtain a simulated inhalation signal; use a sine function with a period of 4T and a phase of [π/ 2 , π] to simulate exhalation. Breathing process, get the simulated exhalation signal; use the sine function whose period is 4T and 3 phase is [π, 3π/2] to simulate the pause process, and get the simulation pause signal; 步骤a4、将模拟吸气信号的相对幅值调整到[a,b]之间,将模拟呼气信号的相对幅值调整到[c,b]之间,将模拟类暂停信号的相对幅值调整到[a,c]之间;Step a4: Adjust the relative amplitude of the analog inspiratory signal to be between [a, b], adjust the relative amplitude of the analog expiratory signal to be between [c, b], and adjust the relative amplitude of the analog pause signal Adjust to between [a, c]; 步骤a5、将得到的模拟吸气信号、模拟呼气信号和模拟类暂停信号首尾相接,得到理想单周期呼吸信号。Step a5: Connect the obtained simulated inhalation signal, simulated expiratory signal and simulated pause signal end to end to obtain an ideal single-cycle breathing signal. 5.一种理想多周期呼吸信号周期提取方法,其特征在于,包括以下步骤:5. an ideal multi-period breathing signal period extraction method, is characterized in that, comprises the following steps: 步骤c1、设定周期数据变量为1,进入步骤c2;Step c1, set the period data variable to 1, and enter step c2; 步骤c2、根据周期数据变量,计算第一周期所有数据之和,以及第二周期所有数据之和,进入步骤c3;Step c2, according to the period data variable, calculate the sum of all data in the first period and the sum of all data in the second period, and enter step c3; 步骤c3、判断第一周期所有数据之和是否等于第二周期所有数据之和,如果:Step c3, determine whether the sum of all data in the first cycle is equal to the sum of all data in the second cycle, if: 是,进入步骤c4;Yes, go to step c4; 否,周期数据变量+1,返回步骤c2;No, cycle data variable +1, return to step c2; 步骤c4、周期数据变量即为多周期呼吸信号的周期。In step c4, the period data variable is the period of the multi-period breathing signal.
CN201910018163.5A 2019-01-09 2019-01-09 Chest and abdomen surface area respiratory signal period prediction method combining double period judgment Expired - Fee Related CN109741827B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910018163.5A CN109741827B (en) 2019-01-09 2019-01-09 Chest and abdomen surface area respiratory signal period prediction method combining double period judgment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910018163.5A CN109741827B (en) 2019-01-09 2019-01-09 Chest and abdomen surface area respiratory signal period prediction method combining double period judgment

Publications (2)

Publication Number Publication Date
CN109741827A true CN109741827A (en) 2019-05-10
CN109741827B CN109741827B (en) 2022-11-01

Family

ID=66364039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910018163.5A Expired - Fee Related CN109741827B (en) 2019-01-09 2019-01-09 Chest and abdomen surface area respiratory signal period prediction method combining double period judgment

Country Status (1)

Country Link
CN (1) CN109741827B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109741829A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101628154A (en) * 2008-07-16 2010-01-20 深圳市海博科技有限公司 Image guiding and tracking method based on prediction
CN101637388A (en) * 2008-07-29 2010-02-03 深圳市海博科技有限公司 Respiratory movement predicting method
CN102481127A (en) * 2009-08-13 2012-05-30 帝人制药株式会社 Arithmetic device for respiratory waveform information and medical equipment using respiratory waveform information
CN102613964A (en) * 2012-03-12 2012-08-01 深圳市视聆科技开发有限公司 Method and system for acquiring physiological signal cycle
TW201544074A (en) * 2014-05-22 2015-12-01 Apex Medical Corp Breathing waveform recognition method and system thereof
CN105286998A (en) * 2015-11-03 2016-02-03 苏州大学 Human body pleuroperitoneal cavity motion simulation device caused by breathing
CN109741828A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 A method for predicting the breathing signal cycle of the human thoracic and abdominal surface area during radiotherapy
CN109741830A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 A method for predicting the cycle of breathing signal in the thoracic and abdominal surface area based on the mixed judgment of single and double cycles
CN109741831A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of breathing signal cycle in human thoracic and abdominal surface area based on variance calculation
CN109741829A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101628154A (en) * 2008-07-16 2010-01-20 深圳市海博科技有限公司 Image guiding and tracking method based on prediction
CN101637388A (en) * 2008-07-29 2010-02-03 深圳市海博科技有限公司 Respiratory movement predicting method
CN102481127A (en) * 2009-08-13 2012-05-30 帝人制药株式会社 Arithmetic device for respiratory waveform information and medical equipment using respiratory waveform information
CN102613964A (en) * 2012-03-12 2012-08-01 深圳市视聆科技开发有限公司 Method and system for acquiring physiological signal cycle
TW201544074A (en) * 2014-05-22 2015-12-01 Apex Medical Corp Breathing waveform recognition method and system thereof
CN105286998A (en) * 2015-11-03 2016-02-03 苏州大学 Human body pleuroperitoneal cavity motion simulation device caused by breathing
CN109741828A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 A method for predicting the breathing signal cycle of the human thoracic and abdominal surface area during radiotherapy
CN109741830A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 A method for predicting the cycle of breathing signal in the thoracic and abdominal surface area based on the mixed judgment of single and double cycles
CN109741831A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of breathing signal cycle in human thoracic and abdominal surface area based on variance calculation
CN109741829A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周寿军: "基于后验概率的呼吸信号预测", 《中国医学生物工程学报》 *
孟晓亮: "基于三维傅里叶变换的胸腹表面测量", 《光学精密工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109741829A (en) * 2019-01-09 2019-05-10 哈尔滨理工大学 Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment
CN109741829B (en) * 2019-01-09 2022-10-28 哈尔滨理工大学 Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment

Also Published As

Publication number Publication date
CN109741827B (en) 2022-11-01

Similar Documents

Publication Publication Date Title
Lee et al. Local control after stereotactic body radiation therapy for stage I non-small cell lung cancer
EP3695882A1 (en) Computing radiotherapy dose distribution
Pollock et al. Breathing guidance in radiation oncology and radiology: a systematic review of patient and healthy volunteer studies
CN109741830A (en) A method for predicting the cycle of breathing signal in the thoracic and abdominal surface area based on the mixed judgment of single and double cycles
US20240000381A1 (en) System, device and method for automated auscultation
CN109741828A (en) A method for predicting the breathing signal cycle of the human thoracic and abdominal surface area during radiotherapy
Chang et al. Real-time respiratory tumor motion prediction based on a temporal convolutional neural network: Prediction model development study
CN109741831A (en) Prediction method of breathing signal cycle in human thoracic and abdominal surface area based on variance calculation
CN109741827A (en) Prediction method of respiratory signal cycle in thoracic-abdominal surface area combined with double-cycle judgment
Geneser et al. Quantifying variability in radiation dose due to respiratory-induced tumor motion
CN109741829A (en) Prediction method of respiratory signal cycle in thoracic and abdominal surface area combined with three-cycle judgment
Kim et al. Development of deformable moving lung phantom to simulate respiratory motion in radiotherapy
Dutta et al. Comparison of flow and gas washout characteristics between pressure control and high-frequency percussive ventilation using a test lung
CN109887573A (en) A method for predicting the motion of human thoracic and abdominal surfaces for radiotherapy
Yoganathan et al. Predicting respiratory motion using a novel patient specific dual deep recurrent neural networks
Cui et al. Commissioning and quality assurance for a respiratory training system based on audiovisual biofeedback
WO2022170607A1 (en) Positioning image conversion system
CN106570335A (en) Colorless transformation based on correlation model between tumor and marker point in stereoscopic radiotherapy
Waghorn et al. A comparison of the dosimetric effects of intrafraction motion on step‐and‐shoot, compensator, and helical tomotherapy‐based IMRT
WO2022165567A1 (en) Method for estimating the muscle pressure of a patient being ventilated by a pulmonary ventilator and auxiliary device for a pulmonary ventilator
Ichiji et al. Respiratory motion prediction for tumor following radiotherapy by using time-variant seasonal autoregressive techniques
Guidi et al. Real-time lung tumour motion modeling for adaptive radiation therapy using lego mindstorms
Schranz et al. A time-continuous model of respiratory mechanics of ARDS patients
Djoumessi et al. Proposed cardio-pulmonary model to investigate the effects of COVID-19 on the cardiovascular system
Shan et al. A comparative study on the generalization ability of back propagation neural network and support vector machine for tracking tumor motion in radiotherapy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20221101