A kind of method and device obtaining the respiratory cycle
Technical field
The invention belongs to field of medical technology, particularly relate to a kind of method and the device that utilize the PET device acquisition respiratory cycle.
Background technology
To when being carried out PET-CT scanning by sweep object, owing to certain artifact can be produced on scan image by the respiratory movement of sweep object. In order to the artifact cut down on scan image, it usually needs adopt respiratory gating technology to gather breath signal to obtain the respiratory cycle.
When conventional respiratory gating technology gathers breath signal, it is necessary to utilize respiration gate control equipment. But, and there is following shortcoming in respiration gate control equipment: 1) expensive, increase cost; 2) complicated operation, and scanning process can be impacted by the bad operation operating doctor; 3) need place with it patient or bind certain device, patient is impacted.
Therefore, it is necessary to provide a kind of respiration gate control equipment that substitutes to obtain method and the device of breath signal.
Summary of the invention
In view of this, the present invention provides a kind of method and the device that obtain the respiratory cycle based on PET device, avoids adopting respiration gate control equipment to obtain breath signal.
In order to realize foregoing invention object, present invention employs following technical scheme:
Obtain the method for respiratory cycle, comprising:
To carrying out when being carried out PET scan by sweep object of respiratory movement, gathering particular data, described particular data is the corresponding data that distribute in by sweep object tissue with radiopharmaceuticals obtained in PET scan;
From scanning to scanning terminate every preset duration statistics once described particular data quantity and; Wherein, a PET scan process can obtain multiple quantity and data; Described quantity has the cycle identical with respiratory movement with data; Described preset duration is less than the described lower value by the respiratory cycle scope of sweep object;
Analyze the quantity of described particular data and the cycle of data, to obtain by the respiratory cycle of sweep object.
Can selection of land, the quantity of the described particular data of described analysis and the cycle of data, to obtain by the respiratory cycle of sweep object, specifically comprise:
Choosing quantity and data from the quantity of described particular data and data, the quantity chosen and data are as benchmark quantity and data;
The quantity with described benchmark quantity and data match and data are searched from the quantity of described particular data and data;
The time point corresponding according to described benchmark quantity and data and the quantity matched and time point corresponding to data, obtain described benchmark quantity and cycle that data occur, and the cycle of described benchmark quantity and data was by the respiratory cycle of sweep object.
Can selection of land, described benchmark quantity and data are N number of continuous print quantity and data, described N be greater than 1 integer, using the region at described benchmark quantity and data place as reference data district, the quantity and the data that are positioned at described reference data district respectively are the 1st data, the 2nd data according to the sequencing of described quantity and data gathering ..., the i-th data ... N data; The described quantity from described particular data and data are searched the quantity with described benchmark quantity and data match and data, specifically comprise:
Number using N as the quantity and data that are evaluated data field, according to quantity and data gathering by first to after order, the quantity of the 1st data being later than in described reference data district and data are belonged to successively and are evaluated data field accordingly; Wherein in adjacent two data district, the 1st data that the 2nd data to N data in a front data field are respectively in a rear data field are to N-1 data; Or, the number using N as the quantity and data that are evaluated data field, according to the order from back to front of quantity and data gathering, is evaluated data field accordingly by belonging to successively early than the quantity of the N data in described reference data district and data; Wherein in adjacent two data district, the 2nd data that the 1st data to N-1 data in a rear data field are respectively in a front data field are to N data;
According to apart from described reference data district by closely judging the quantity in described reference data district and data to order far away successively and be evaluated the quantity in data field and whether the relation data meets pre-conditioned, until judged result is yes, it is determined that meet the pre-conditioned quantity being evaluated in data field and data are the quantity with described benchmark quantity and data match and data.
Can selection of land, the described quantity from described particular data and search the quantity with described benchmark quantity and data match and data, be specially:
According to searching the quantity with described benchmark quantity and data match and data successively at a distance of described benchmark quantity and data order from the close-by examples to those far off in predetermined interval; Described predetermined interval is grown up in the described higher limit by the respiratory cycle scope of sweep object when crossing over;
Time be positioned at described benchmark quantity and data when described predetermined interval before, the quantity that the statistics of time point the latest in described predetermined interval obtains and data are front data of described benchmark quantity and data;
Time be positioned at described benchmark quantity and data when described predetermined interval after, the quantity that the statistics of time point the earliest in described predetermined interval obtains and data are rear data of described benchmark quantity and data.
Can selection of land, described particular data is the counting rate inciding on the crystal of PET device.
Can selection of land, the crystal of PET device is divided into some groups, and described particular data is the counting rate inciding the crystal belonging to preset group.
Obtain the device of respiratory cycle, comprising:
Data acquisition unit, for carrying out when being carried out PET scan by sweep object of respiratory movement, gathering particular data, described particular data is the corresponding data that distribute in by sweep object tissue with radiopharmaceuticals obtained in PET scan;
Data statistics unit, for from scanning to scanning terminate every preset duration statistics once described particular data quantity and; Wherein, a PET scan process can obtain multiple quantity and data; Described quantity has the cycle identical with respiratory movement with data; Described preset duration is less than the described lower value by the respiratory cycle scope of sweep object;
Cycle analysis unit, for the cycle of the quantity and data of analyzing described particular data, to obtain by the respiratory cycle of sweep object.
Can selection of land, described cycle analysis unit, specifically comprises:
Choosing subelement, for choosing quantity and data from the quantity of described particular data and data, the quantity chosen and data are as benchmark quantity and data;
Search subelement, for searching the quantity with described benchmark quantity and data match and data from the quantity of described particular data and data;
Obtain subelement, for the time point that the time point corresponding according to described benchmark quantity and data and the quantity matched and data are corresponding, obtaining described benchmark quantity and the cycle of data appearance, the cycle of described benchmark quantity and data was by the respiratory cycle of sweep object.
Can selection of land, described benchmark quantity and data are N number of continuous print quantity and data, described N be greater than 1 integer, using the region at described benchmark quantity and data place as reference data district, the quantity and the data that are positioned at described reference data district respectively are the 1st data, the 2nd data according to the sequencing of described quantity and data gathering ..., the i-th data ... N data; Described subelement of searching specifically comprises:
It is evaluated data field ownership subelement, for the number using N as the quantity and data that are evaluated data field, according to quantity and data gathering by first to rear order, the quantity of the 1st data being later than in described reference data district and data being belonged to successively and is evaluated data field accordingly; Wherein in adjacent two data district, the 1st data that the 2nd data to N data in a front data field are respectively in a rear data field are to N-1 data; Or, the number using N as the quantity and data that are evaluated data field, according to the order from back to front of quantity and data gathering, is evaluated data field accordingly by belonging to successively early than the quantity of the N data in described reference data district and data; Wherein in adjacent two data district, the 2nd data that the 1st data to N-1 data in a rear data field are respectively in a front data field are to N data;
Judgment sub-unit, for according to apart from described reference data district by closely judging the quantity in described reference data district and data to order far away successively and be evaluated the quantity in data field and whether the relation data meets pre-conditioned, until judged result is yes, it is determined that meet the pre-conditioned quantity being evaluated in data field and data are the quantity with described benchmark quantity and data match and data.
Can selection of land, described particular data is the counting rate inciding on the crystal of PET device.
Compared to prior art, the present invention has following useful effect:
In respiratory, different sizes can be presented under different breathing states by the lobe of the lung of sweep object, meanwhile, cyclical movement can be done along with respiratory by tissues such as the livers of sweep object. Therefore, the distribution of radiopharmaceuticals in being organized by sweep object can change along with the stage of the respiratory cycle residing for sweep object, the number of its particular data sent also can change, so, to carrying out when being carried out PET scan by sweep object of respiratory movement, the distribution of quantity of the particular data gathered also has certain periodicity, and the cycle of this distributed number with by the respiratory cycle of sweep object identical. Based on this, provided by the invention acquisition the respiratory cycle method by every preset duration add up the particular data in a PET scan process quantity and, then the quantity of this number of series sum that analytic statistics obtains and the cycle of data appearance, these quantity obtained and the cycle of data were by the respiratory cycle of sweep object.
Therefore, the method for acquisition respiratory cycle provided by the invention can get the respiratory cycle by means of PET device and PET scan, it is not necessary to uses special respiration gate control equipment, therefore overcomes the various defects using respiration gate control equipment. Thus, compared to prior art, the method for acquisition respiratory cycle provided by the invention saves equipment cost, and simplify the operation of operator, in addition, due to without the need to the health of sweep object being placed or binds certain device, so improve by the comfort level of sweep object.
Accompanying drawing explanation
In order to be expressly understood the technical scheme of the present invention, the accompanying drawing used when describing the specific embodiment of the present invention below does a brief description. It is appreciated that these accompanying drawings are only the part embodiments of the present invention, those of ordinary skill in the art are under the prerequisite not paying creative work, it is also possible to obtain other accompanying drawing.
Fig. 1 is the method flow schematic diagram of the acquisition respiratory cycle that the embodiment of the present invention one provides;
Fig. 2 is the quantity of the described particular data of analysis and the method flow schematic diagram in the cycle of data of the embodiment of the present invention one offer;
Fig. 3 is the method flow schematic diagram utilizing mutuality counting method to search the quantity mated mutually and data;
Fig. 4 is the method flow schematic diagram of the acquisition respiratory cycle that the embodiment of the present invention two provides;
Fig. 5 is the quantity of the described particular data of analysis and the method flow schematic diagram in the cycle of data of the embodiment of the present invention two offer;
Fig. 6 is the structural representation of the device of the acquisition respiratory cycle that the embodiment of the present invention three provides;
Fig. 7 is the structural representation of the cycle analysis unit that the embodiment of the present invention three provides;
Fig. 8 is the structural representation searching subelement that the embodiment of the present invention three provides.
Embodiment
For make the object of the present invention, technical scheme, useful effect clearly, complete, below in conjunction with accompanying drawing, the specific embodiment of the present invention is described.
Embodiment one
Fig. 1 is the method flow schematic diagram of the acquisition respiratory cycle that the embodiment of the present invention one provides. As shown in Figure 1, the method comprises the following steps:
S101, to carrying out when being carried out PET scan by sweep object of respiratory movement, gather particular data:
It should be noted that, described particular data can be any data obtained in PET scan process. Specifically, described particular data be obtain in PET scan with the radiopharmaceuticals data that distribution is corresponding in by sweep object tissue. Such as, it is possible to be the counting rate of the �� photon on the crystal inciding PET device, it is possible to thinking the counting rate of the electric impulse signal that the photomultiplier of PET device exports, in addition, described particular data can also for meeting event or event.
Further, if particular data is the counting rate of the �� photon inciding on crystal, in order to reduce data noise, it is preferable that gather the counting rate of the �� photon inciding on multiple crystal. The counting rate of the �� photon that these collect is relatively big, is conducive to reducing data noise. Owing to a ring comprising multiple crystal, so, in order to convenient, described particular data can the counting rate of �� photon for inciding on a certain ring crystal.
In addition, during owing to being carried out respiratory movement by sweep object, the comparatively significant position of its motion is lung. So, in embodiments of the present invention, it is preferable that carry out PET scan to by the lung areas of sweep object and base of lung region. Further, it is possible to determined by the scanning area scope of sweep object by the mode of the flat sheet of scanning computed tomography.
S102, from scanning to scanning terminate every preset duration statistics once described particular data quantity and:
It should be noted that, under normal circumstances, those skilled in the art are according to the common practise of biology or medical field or search relevant data and can know such as, to by the respiratory cycle scope of sweep object, the respiratory cycle of people is usually about the 3-6 second. That is, can be known in advance by the higher limit of the respiratory cycle of sweep object and lower value. In order to guarantee at least can to add up within a respiratory cycle quantity and, described preset duration is less than the described lower value by the respiratory cycle scope of sweep object. Further, described preset duration is more little, and the data obtained within a respiratory cycle are more many, and the thus obtained respiratory cycle is more accurate. By this step, a PET scan process can obtain multiple quantity and the data of described particular data.
Due to by sweep object in respiratory, the different steps of the lobe of the lung in the respiratory cycle presents different sizes, cyclical movement can be done along with respiratory by tissues such as the livers of sweep object simultaneously, so being positioned over the radiopharmaceuticals by breathing in object tissue, by the different steps of the respiratory cycle of sweep object, be there is different distributions, therefore, the quantity of the particular data produced in the different steps of respiratory cycle is also different, so, to when being carried out PET scan by sweep object, also there is certain periodicity in the quantity of the particular data collected and data, and this quantity has the cycle identical with respiratory movement with data.
S103, analyze the quantity of described particular data and the cycle of data, to obtain by the respiratory cycle of sweep object:
As mentioned above, it is necessary, quantity has the cycle identical with the respiratory movement by sweep object with data, therefore, the cycle analyzing quantity and the data drawn was by the respiratory cycle of sweep object.
As an example of the present invention, the quantity of the described particular data of described analysis and the cycle of data, to obtain by the respiratory cycle of sweep object, it can realize by analyzing the cycle of wherein one or more quantity and data.
The quantity of the described particular data of described analysis and the cycle of data, to obtain by the process of the respiratory cycle of sweep object as shown in Figure 2, it specifically can comprise the following steps:
S201, choosing quantity and data from the quantity and data of described particular data, the quantity chosen and data are as benchmark quantity and data:
It should be noted that, the quantity chosen from the quantity of described particular data and data and data can be one, it is possible to think multiple. When for multiple quantity and data, it is preferable to multiple continuous print quantity and data.
When benchmark quantity and data are multiple continuous print quantity and during data, be conducive to improving the accuracy obtaining the respiratory cycle. This is because: within the same cycle, likely there will be two or more identical quantity and data, so, when the benchmark quantity chosen and data are one, that likely occur finding with benchmark quantity and data match quantity and data and described benchmark quantity and data be in fact positioned at the same cycle, and mate to think it is two different cycles by mistake mutually due to both, the cycle thus calculated is shorter than real actual cycle, thus causes the respiratory cycle obtained inaccurate. But, when benchmark quantity and data are multiple continuous print quantity and during data, it is necessary to when finding quantity and the data of multiple continuous print quantity and data match with this, just think that these quantity and data are the data in the another one cycle. When only having part data match, the data in the another one cycle can not be thought. Therefore, mate compared to data, select multiple continuous print quantity and data to reduce, as benchmark quantity and data, the probability that the acquisition cycle makes mistakes. And, the benchmark quantity chosen and the number of data are more many, and the probability of erroneous judgement is more low.
S202, from the quantity and data of described particular data, search the quantity with described benchmark quantity and data match and data:
Described is the quantity identical or close with data with described benchmark quantity and data with the quantity of data match with data with benchmark quantity. Described identical quantity is namely completely identical with data with described benchmark quantity with data, and the difference that described close quantity and data differ with described benchmark quantity and data is within the specific limits. The scope of this difference can be the numerical value that user determines according to working experience. In addition, it is possible to represent both differences by both variances.
In addition, due to the described particular data that counts in a PET scan process quantity and the time length crossed over of data be likely multiple different cycles. If all and benchmark quantity and data match quantity and data are all found out by the quantity of described particular data and data, can take a long time, and in order to find the cycle of quantity and data, there is no need the data finding all phases to mate, as long as finding the data mated mutually belonged in two adjacent periods to extrapolate the cycle. So, in order to improve search efficiency, it is possible to only search the quantity with described benchmark quantity and data match and data from the partial amt and data of described particular data.
Further, as long as finding the quantity with described benchmark quantity and data match and data within the cycle adjacent with the cycle of described benchmark quantity and data place, the cycle of data then can be got according to the difference of time point corresponding to these two data. So, as a preferred embodiment of the present invention, can by searching the quantity with described benchmark quantity and data match and data successively at a distance of described benchmark quantity and data order from the close-by examples to those far off in predetermined interval, until finding the quantity with described benchmark quantity and data match and data. Once after finding the quantity and data mated mutually, stop searching process, then calculating these two time differences between quantity and data mated mutually, this time difference is the cycle of quantity and data, also namely by the respiratory cycle of sweep object.
It should be noted that, when benchmark quantity and data are add up the quantity and data that obtain the earliest, then after described predetermined interval is positioned at described benchmark quantity and data. When benchmark quantity and data are add up the quantity and data that obtain the latest, then before described predetermined interval is positioned at described benchmark quantity and data. When all there are a lot of quantity and data in the front and back of benchmark quantity and data, and when the quantity of front and back and data are all many, then before predetermined interval can be positioned at described benchmark quantity and data, it is also possible to after being positioned at described benchmark quantity and data.
In embodiments of the present invention, in order to guarantee to find the quantity mated mutually and data, described predetermined interval cross over time grow up in the described higher limit by the respiratory cycle scope of sweep object. Further, in order to reduce searching work amount, this pre-set interval is preferably less than the higher limit of the suction periodic regime of 2 times. Such as, the higher limit of the respiratory cycle of people is 6 seconds, then predetermined interval cross over time grown up in 6 seconds, can be such as 7 seconds.
And, in order to as much as possible with benchmark quantity and the adjacent periods in data place cycle in find the data with described benchmark quantity and data match, time be positioned at described benchmark quantity and data when described predetermined interval before, the quantity that the statistics of time point the latest in described predetermined interval obtains and data are front data of described benchmark quantity and data. Time be positioned at described benchmark quantity and data when described predetermined interval after, the quantity that the statistics of time point the earliest in described predetermined interval obtains and data are rear data of described benchmark quantity and data.
As an example, when described benchmark quantity and data are data, described from the quantity of described particular data with data are searched the quantity with described benchmark quantity and data match and data specifically can by variance method from except quantity and the data of the described particular data of described benchmark quantity except data being searched and described benchmark quantity and data match and identical or close quantity and data. When searching the quantity mated mutually and data by variance method, benchmark quantity and data and the variance between the quantity mated with its phase and data are in preset threshold range.
As the another one example of the present invention, when described benchmark quantity and data are multiple data, described benchmark quantity data can be regarded as a reference data district, then with all data in reference data district integrally, the data field that is evaluated identical with data amount check respectively compares, with the data field found with described reference data district mates mutually. It should be noted that, the time length that described reference data district is crossed over is not limited by the embodiment of the present invention. So, the time length that the data in described reference data district are crossed over can be less than, be equal to or greater than by the length of the respiratory cycle of sweep object.
Setting benchmark quantity and data are N number of continuous print quantity and data, wherein, N be greater than 1 integer, then, quantity and data amount check in reference data district are N, definition is positioned at the quantity in described reference data district and data respectively are the 1st data, the 2nd data according to the sequencing of described quantity and data gathering ..., the i-th data ... N data. In this case, search and divide two kinds of situations with described benchmark quantity and the quantity of data match and the method for data, the first: search from front to back according to the sequencing of data, the 2nd kind of situation: search from back to front according to the sequencing of data. Lookup method corresponding to these two kinds of situations comprises the following steps:
First, by quantity and attribution data to being evaluated data field accordingly:
For two kinds of different situations, by also different, specific as follows to the method being evaluated data field accordingly to quantity and attribution data:
The first situation: the method searched from front to back according to the sequencing of data, is specially quantity and attribution data to being evaluated data field accordingly:
Number using N as the quantity and data that are evaluated data field, according to quantity and data gathering by first to after order, the quantity of the 1st data being later than in described reference data district and data are belonged to successively and are evaluated data field accordingly; Wherein in adjacent two data district, the 1st data that the 2nd data to N data in a front data field are respectively in a rear data field are to N-1 data:
Illustrate: assume that the quantity of the particular data gathered and data have 100, respectively be according to the sequencing gathered: the 1st data, the 2nd data ..., the i-th data ... 100th data.
Comprising 10 continuous print quantity and data in reference data district, and 10 data in this reference data district are 10 data gathered at first, it respectively is the 1st data, the 2nd data according to the sequencing gathered ..., the i-th data ... 10th data.
Then above-mentioned steps can example be: using 10 as the number of the quantity and data that are evaluated data field, from the 2nd data according to data gathering by first to after order, 2nd data are evaluated data field to the 11st attribution data to first, 3rd data are evaluated data field to the 12nd attribution data to the 2nd, analogize successively, until all quantity and data all belong to is evaluated data field accordingly. This kind of affiliation method makes the data in adjacent two data district have intersection overlapping, and namely the 2nd data to the 10th data in a front data field are that the 1st data in a rear data field are to the 9th data.
2nd kind of situation: the method searched from back to front according to the sequencing of data, is specially quantity and attribution data to being evaluated data field accordingly:
Number using N as the quantity and data that are evaluated data field, according to the order from back to front of quantity and data gathering, is evaluated data field accordingly by belonging to successively early than the quantity of the N data in described reference data district and data; Wherein in adjacent two data district, the 2nd data that the 1st data to N-1 data in a rear data field are respectively in a front data field are to N data:
Illustrate: assume that the quantity of the particular data gathered and data have 100, respectively be according to the sequencing gathered: the 1st data, the 2nd data ..., the i-th data ... 100th data.
10 continuous print quantity and data are comprised in reference data district, and 10 data in this reference data district are last 10 data gathered, it respectively is the 91st data, the 92nd data according to the sequencing gathered ..., the i-th data ... 100th data.
Then above-mentioned steps can example be: using 10 as the number of the quantity and data that are evaluated data field, according to data gathering by backward front order from the 99th data, 99th data are evaluated data field to the 90th attribution data to first, 98th data are evaluated data field to the 89th attribution data to the 2nd, analogize successively, until all quantity and data all belong to is evaluated data field accordingly. Equally, this kind of affiliation method makes the data in adjacent two data district have intersection, and the 2nd data that namely the 1st data to the 9th data in a rear data field are respectively in a front data field are to the 10th data.
It should be noted that, the cycle that quantity and attribution data can accurately be found data to occur to the method being evaluated data field accordingly by this kind.
Then, by apart from described reference data district by closely judging the quantity in described reference data district and data to order far away successively and be evaluated the quantity in data field and whether the relation data meets pre-conditioned, until judged result is yes, determine that meeting the pre-conditioned data field that is evaluated is the data field mated mutually with reference data district, meet the pre-conditioned quantity being evaluated in data field and data are the quantity with described benchmark quantity and data match and data:
According to for example above-mentioned, first judge whether the quantity in the first quantity being evaluated in data field and data and reference data district and the relation between data meet pre-conditioned, if, then determining that first to be evaluated data field be the data field mated mutually with reference data district for this, these first the 1st data being evaluated in data field are the quantity and data mated mutually with the 1st data in reference data district; If not, judge whether the quantity in the 2nd quantity being evaluated in data field and data and reference data district and the relation between data meet pre-conditioned, if, then determining the 2nd, to be evaluated data field be the data field mated mutually with reference data district, and the 2nd the 1st data being evaluated in data field are the quantity and data mated mutually with the 1st data in reference data district. Analogize successively, until above-mentioned judged result is yes.
It should be noted that, when adopting different methods carry out the quantity and the data that judge in described reference data district and be evaluated the quantity in data field and the relation between data, described pre-conditioned also thereupon different.
Specifically, when the difference of the average by quantity and data represents the quantity in described reference data district and data and is evaluated the quantity in data field and the relation between data, described pre-conditioned for whether the difference of the quantity in described reference data district and data and the average of the described quantity that is evaluated in data field and data is less than the first threshold value.
When representing the quantity in described reference data district and data by the variance of quantity and data and be evaluated the quantity in data field and the relation between data, described pre-conditioned for whether the variance of the quantity in described reference data district and data and the described quantity that is evaluated in data field and data is less than Second Threshold.
When representing the quantity in described reference data district and data with the cross-correlation coefficient of quantity and data and be evaluated the quantity in data field and the relation between data, whether the described pre-conditioned cross-correlation coefficient for the quantity in described reference data district and data and the described quantity that is evaluated in data field and data is in predetermined threshold value.
Wherein, the calculation formula of described cross-correlation coefficient is as follows:
Wherein: xiFor i-th data in reference data district, the span of i is all data amount checks in reference data district; X is the expectation of all data values in reference data district; yiFor being evaluated the i-th data of data field, Y is the expectation being evaluated all data values in data field. Wherein, cc value and 1 is more close, illustrates that being evaluated data field more mates with reference data district.
For the method for cross-correlation coefficient, the above-mentioned process searching quantity and the data mated mutually is described below.
The quantity that setting statistics obtains and data have 100, benchmark quantity and data are 5, so, 5 data are had in reference data district, and the data in this reference data district are add up 5 data obtained the earliest, now, as shown in Figure 3, it comprises the following steps the detailed process adopting above-mentioned lookup method:
A, using 5 as the number of the quantity being evaluated in data field and data, successively the 2nd to the 100th data are belonged to and different it is evaluated in data field, wherein, the sequence of data is the sequencing of timing statistics, 1st data are the data added up the earliest, and the 100th data are the data of last statistics. Like this, the 2-6 data belongs to first and is evaluated data field, and the 3-7 data belongs to the 2nd and be evaluated data field, analogizes successively.
First cross-correlation coefficient of b, the quantity calculated in the first quantity being evaluated in data field and data and reference data district and data.
C, judging described first cross-correlation coefficient whether in predetermined threshold value, if it does, then determine first, to be evaluated data field with reference data district be the data field mated mutually, and quantity in its two region and data are the data of coupling mutually; If not, steps d is performed.
D, calculate the 2nd and be evaluated in data field quantity and data and reference data district in quantity and the 2nd cross-correlation coefficient of data.
E, judging described 2nd cross-correlation coefficient whether in predetermined threshold value, if it does, then determine the 2nd, to be evaluated data field with reference data district be the data field mated mutually.
Fig. 3 illustrate only the multilevel iudge process that two are evaluated data field and reference data district, in fact, when being not less than the 3rd threshold value when the cross-correlation coefficient obtained, need to be evaluated data field and the cross-correlation coefficient in reference data district by the order from the close-by examples to those far off apart from described reference data district successively multilevel iudge, until described in be evaluated in data field quantity and data and reference data district in quantity and the cross-correlation coefficient of data be less than the 3rd threshold value.
It should be noted that, above-mentioned lookup method is described for cross-correlation coefficient, and in fact, above-mentioned cross-correlation coefficient can replace the difference into variance or average.
S203, the time point corresponding according to time point corresponding to described benchmark quantity and data and the quantity matched and data, calculate described benchmark quantity and the shortest time between data and the quantity matched and data poor, it is by the respiratory cycle of sweep object with the cycle that the cycle of the described benchmark quantity and data that obtain described particular data, described benchmark quantity and data occur:
It should be noted that, two time differences between quantity and data mated mutually that described benchmark quantity and the difference of the shortest time between data and the quantity matched and data are in adjacent two cycles, this time difference is the cycle of benchmark quantity and data, also being exactly a series of multiple quantity of described particular data and the cycle of data, this cycle was by the respiratory cycle of sweep object.
Specifically, when benchmark quantity and data only one time, if according to when searching the quantity and data that obtain with described benchmark quantity and data match successively apart from described benchmark quantity and data order from the close-by examples to those far off, first found is adjacent with the cycle of described benchmark quantity and data place with the cycle at data place with described benchmark quantity and the quantity of data match, then now, calculate this quantity mated mutually found and the difference of time point corresponding to the data time point corresponding with benchmark quantity and data, the difference so obtained be described benchmark quantity and the shortest time between data and the quantity matched and data poor, it is the cycle of quantity and data, the cycle of this quantity and data was by the respiratory cycle of sweep object.
When finding the quantity with the quantity in described reference data district and data match and data by entirety of reference data district, poor with the time point the earliest being evaluated place, data field with the time point the earliest at place, reference data district, the time length obtained is the cycle in reference data district, and the cycle in this reference data district was by the respiratory cycle of sweep object.
The embodiment of the method for the acquisition respiratory cycle provided for the embodiment of the present invention one above. In the method, based on the distribution rule of the quantity sum of the particular data in PET device and PET scan process, the cycle of quantity and data is drawn. Due to radiopharmaceuticals in the different steps of respiratory cycle different by the distribution in the tissue of sweep object, the size of the quantity sum of the particular data thus obtained based on the distribution of radiopharmaceuticals is also different along with the difference in stage respiratory cycle, and it has the period profile identical with the respiratory cycle. Therefore, cycle based on quantity and data can obtain by the respiratory cycle of sweep object. And the cycle of quantity and data was by the respiratory cycle of sweep object.
The method of acquisition respiratory cycle provided by the invention can get the respiratory cycle by means of PET device and PET scan, it is not necessary to uses special respiration gate control equipment, therefore overcomes the various defects using respiration gate control equipment. Thus, compared to prior art, the method for acquisition respiratory cycle provided by the invention saves equipment cost, and simplify the operation of operator, in addition, due to without the need to the health of sweep object being placed or binds certain device, so improve by the comfort level of sweep object.
In order to clearly understand the embodiment that the present invention obtains the method for respiratory cycle, below to be described the embodiment of the method for acquisition respiratory cycle described above for human body by sweep object. Referring specifically to embodiment two.
Embodiment two
As shown in Figure 4, this embodiment comprises the following steps:
S401, human body is carried out PET scan, gathers the �� photon inciding on crystal:
Radiopharmaceuticals occurs to produce �� photon, meet necessarily required �� photon and can incide on the crystal on PET device in the process of decay in human body. So, when human body is carried out PET scan, it is possible to gather the counting rate of the �� photon inciding on crystal.
It should be noted that, in embodiments of the present invention, in advance the crystal of PET device is divided into groups, such as using ring as unit, crystal is divided into groups by the ring at the crystal place according to PET device, and such as, PET device has 33 ring crystal, when human body is carried out PET scan, it is possible to gather the counting rate of the �� photon inciding on every ring crystal respectively. The counting rate often organizing the �� photon on crystal collected is the particular data described in the embodiment of the present invention. In addition, described particular data can also be the counting rate of the crystal that belongs to multiple groups. In general, the counting rate of the crystal that described particular data can belong to preset group for inciding. Described preset group can be one or more groups. Further, the position at described preset group place can be determined according to practical situation by operator.
S402, the counting rate terminating to add up every 100ms the �� photon on crystal from scanning to scanning, obtain the counting rate sum of multiple every 100ms of �� photon:
Due under normal circumstances, the time of a PET scan may the several seconds, so, in a PET scan, it is possible to obtain tens or more counting rate sum. Distribution and the respiratory cycle of this series of counting rate sum obtained have the identical cycle.
It should be noted that, the counting rate of the �� photon on described statistics crystal can be the counting rate on a crystal, it is possible to think the counting rate on multiple crystal. Owing to the numerical value of a counting rate sum is more big, data noise is more little, and the cycle obtained is more accurate, so, it is preferable that add up the counting rate on multiple crystal.
Owing to every ring crystal comprising multiple single crystal, so in order to convenient statistics, it is possible to taking the ring at crystal place as benchmark, add up the counting rate sum of the �� photon on a certain ring crystal.
Further, it is also possible to add up the counting rate sum of the �� photon on every ring crystal, with current PET device, have 33 ring crystal, then in each 100ms, just obtain 33 counting rate sums.
S403, time point according to a series of counting rate sum data of obtaining of statistics and each counting rate sum data, analyze the cycle of counting rate sum data, to obtain the respiratory cycle of human body:
Owing to the distribution of radiopharmaceuticals in tissue is subject to the impact of periodic breathing motion, so, the counting rate sum data obtained in a PET scan process also have certain periodicity. According to the time point that two same count rate sum data occur, analyze the cycle obtaining counting rate sum data, and then obtain by the respiratory cycle of sweep object.
In the present embodiment, analyzing the cycle obtaining counting rate sum data can adopt the arbitrary analytical procedure described in above-described embodiment one. The process how analyzing the cycle obtaining counting rate sum data it is described as benchmark quantity and data instance below using counting rate sum data. As shown in Figure 5, this process comprises the following steps:
S501, from a series of counting rate sum data of obtaining of statistics, choose counting rate sum data as benchmark quantity and data:
It should be noted that, when when statistical counting rate, within the same time period, the counting rate of the �� photon on the crystal on multiple different rings is all added up, like this, section just can obtain multiple counting rate sum at one time, and these counting rate sums that section obtains at one time are counting rate sum group. Now, benchmark quantity and data can be a counting rate sum group, and it is included in multiple counting rate sums that statistics of same time period obtains.
S502, from described counting rate sum data, search the counting rate sum data with described benchmark quantity and data match:
Owing to a PET scan process may comprise multiple respiratory cycle, each respiratory cycle includes the counting rate sum data with described benchmark quantity and data match, in order to obtain the cycle of counting rate sum data, by the difference of the corresponding time point of the counting rate sum data mated mutually that calculate in adjacent two cycles, so, search without the need in whole counting rate sum data that the PET scan process complete from obtains, only need to search from the data of the counting rate sum in a certain predetermined interval. As mentioned above, this predetermined interval is grown up in the described higher limit by the respiratory cycle scope of sweep object when crossing over, and timing statistics point in described predetermined interval quantity the earliest and data are upper data of described benchmark quantity and data, or, the quantity the latest of the timing statistics point in described predetermined interval and data are the next data of described benchmark quantity and data. The scope 3-6 second according to the human body respiration cycle, the time length that this predetermined interval is crossed over can be 7 seconds.
S503, the time point corresponding according to time point corresponding to described benchmark quantity and data and the counting rate sum data that match, calculate the shortest time between the two poor, to obtain the cycle of counting rate sum data, the cycle of described counting rate sum data is the respiratory cycle of human body.
The cycle of counting rate sum data can be analyzed by step S501 to S503, and then the respiratory cycle of human body can be got.
It should be noted that, above-described embodiment searches the counting rate sum data with a counting rate sum data match usually to obtain the cycle of counting rate sum data. Owing to may there are two or more identical counting rate sum data in one-period, so, the cycle adopting the method for counting rate sum data to obtain may be less than actual cycle, thus, and the possibility of result mistake drawn. It is data in the another one cycle to avoid thinking the data mated mutually within the same cycle by mistake, the embodiment of the present invention is preferably using multiple continuous print counting rate sum as benchmark quantity and data, such as, it is possible to using each counting rate sum 20 data in length when every 2 seconds as benchmark quantity and data. Like this, only finding the data with these 20 counting rate sum data match and just think in different cycles, just corresponding according to these data time point obtains the cycle. Thus this kind of method improves the accuracy rate in acquisition cycle compared to former approach.
The embodiment of the method in the acquisition human body respiration cycle provided for the embodiment of the present invention two above. The method, without the need to adopting respiration gate control equipment, thus overcomes the various defects using respiration gate control equipment. Thus, compared to prior art, the method for acquisition respiratory cycle provided by the invention saves equipment cost, and simplifies the operation of operator, in addition, due to without the need to placing on human body or bind certain device, so improve the comfort level of people.
Based on the method for the acquisition respiratory cycle described in above-described embodiment one and embodiment two, present invention also offers the device embodiment obtaining the respiratory cycle. See embodiment three.
Embodiment three
As shown in Figure 6, the device of this acquisition respiratory cycle comprises with lower unit:
Data acquisition unit 61, for carrying out when being carried out PET scan by sweep object of respiratory movement, gathering particular data, described particular data is the corresponding data that distribute in by sweep object tissue with radiopharmaceuticals obtained in PET scan;
Data statistics unit 62, for from scanning to scanning terminate every preset duration statistics once described particular data quantity and; Wherein, a PET scan process can obtain multiple quantity and data; Described quantity has the cycle identical with respiratory movement with data; Described preset duration is less than the described lower value by the respiratory cycle scope of sweep object;
Cycle analysis unit 63, for the cycle of the quantity and data of analyzing described particular data, to obtain by the respiratory cycle of sweep object.
By the device of acquisition respiratory cycle provided by the invention, based on the distribution rule of the quantity sum of the particular data in PET device and PET scan process, draw the cycle of quantity and data. Due to radiopharmaceuticals in the different steps of respiratory cycle different by the distribution in the tissue of sweep object, the size of the quantity sum of the particular data thus obtained based on the distribution of radiopharmaceuticals is also different along with the difference in stage respiratory cycle, and it has the period profile identical with the respiratory cycle. Therefore, cycle based on quantity and data can obtain by the respiratory cycle of sweep object. And the cycle of quantity and data was by the respiratory cycle of sweep object.
The device of acquisition respiratory cycle provided by the invention can get the respiratory cycle by means of PET device and PET scan, it is not necessary to uses special respiration gate control equipment, therefore overcomes the various defects using respiration gate control equipment. Thus, compared to prior art, the method for acquisition respiratory cycle provided by the invention saves equipment cost, and simplify the operation of operator, in addition, due to without the need to the health of sweep object being placed or binds certain device, so improve by the comfort level of sweep object.
Further, as shown in Figure 7, described cycle analysis unit 63, specifically comprises:
Choosing subelement 631, for choosing quantity and data from the quantity of described particular data and data, the quantity chosen and data are as benchmark quantity and data;
Search subelement 632, for searching the quantity with described benchmark quantity and data match and data from the quantity of described particular data and data;
Obtain subelement 633, for the time point that the time point corresponding according to described benchmark quantity and data and the quantity matched and data are corresponding, obtaining described benchmark quantity and the cycle of data appearance, the cycle of described benchmark quantity and data was by the respiratory cycle of sweep object.
Further, described benchmark quantity and data are N number of continuous print quantity and data, described N be greater than 1 integer, using the region at described benchmark quantity and data place as reference data district, the quantity and the data that are positioned at described reference data district respectively are the 1st data, the 2nd data according to the sequencing of described quantity and data gathering ..., the i-th data ... N data; As shown in Figure 8, search subelement 632 described in specifically to comprise:
It is evaluated data field ownership subelement 6321, for the number using N as the quantity and data that are evaluated data field, according to quantity and data gathering by first to rear order, the quantity of the 1st data being later than in described reference data district and data being belonged to successively and is evaluated data field accordingly; Wherein in adjacent two data district, the 1st data that the 2nd data to N data in a front data field are respectively in a rear data field are to N-1 data; Or, the number using N as the quantity and data that are evaluated data field, according to the order from back to front of quantity and data gathering, is evaluated data field accordingly by belonging to successively early than the quantity of the N data in described reference data district and data; Wherein in adjacent two data district, the 2nd data that the 1st data to N-1 data in a rear data field are respectively in a front data field are to N data:
Judgment sub-unit 6322, for according to apart from described reference data district by closely judging the quantity in described reference data district and data to order far away successively and be evaluated the quantity in data field and whether the relation data meets pre-conditioned, until judged result is yes, it is determined that meet the pre-conditioned quantity being evaluated in data field and data are the quantity with described benchmark quantity and data match and data.
It should be noted that, particular data described above is the counting rate inciding on the crystal of PET device. Further, the crystal of PET device is divided into some groups, and described particular data is the counting rate inciding the crystal belonging to preset group.
The above is only the preferred embodiment of the present invention; it is noted that for those skilled in the art, under the premise without departing from the principles of the invention; can also making some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.