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.