CN106097384B - A kind of method and apparatus of determining gate-control signal - Google Patents
A kind of method and apparatus of determining gate-control signal Download PDFInfo
- Publication number
- CN106097384B CN106097384B CN201610374708.2A CN201610374708A CN106097384B CN 106097384 B CN106097384 B CN 106097384B CN 201610374708 A CN201610374708 A CN 201610374708A CN 106097384 B CN106097384 B CN 106097384B
- Authority
- CN
- China
- Prior art keywords
- period
- sub
- oblivion
- angle position
- radial position
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000033001 locomotion Effects 0.000 claims abstract description 302
- 230000003902 lesion Effects 0.000 claims abstract description 151
- 238000002600 positron emission tomography Methods 0.000 claims abstract description 11
- 230000006798 recombination Effects 0.000 claims description 95
- 238000005215 recombination Methods 0.000 claims description 95
- 238000012545 processing Methods 0.000 claims description 37
- 238000001228 spectrum Methods 0.000 claims description 30
- 238000000513 principal component analysis Methods 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000008030 elimination Effects 0.000 claims description 15
- 238000003379 elimination reaction Methods 0.000 claims description 15
- 239000003550 marker Substances 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 29
- 239000013078 crystal Substances 0.000 description 20
- 230000008569 process Effects 0.000 description 16
- 230000015654 memory Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 12
- 201000010099 disease Diseases 0.000 description 11
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 11
- 230000000241 respiratory effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000013144 data compression Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 230000008521 reorganization Effects 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 208000012902 Nervous system disease Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000004903 cardiac system Anatomy 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000036387 respiratory rate Effects 0.000 description 2
- 230000009131 signaling function Effects 0.000 description 2
- 238000012879 PET imaging Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
The invention discloses a kind of method and apparatus of determining gate-control signal, belong to medical domain.The described method includes: obtain detected in the sub-period of each preset duration of the positron emission tomography PET device when detecting in section respectively bury in oblivion dot position information;Respectively bury in oblivion dot position information according to each sub-period is corresponding, determine the direction of motion of lesion tissue, determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement;According to the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement, gate-control signal is determined, PET image reconstruction is carried out based on the gate-control signal.Using the present invention, diagnostic system cost can be reduced.
Description
Technical field
The present invention relates to medical domain, in particular to a kind of method and apparatus of determining gate-control signal.
Background technique
In the field of medicine, PET/CT (positron emission tomography/computed can be passed through
Tomography, positron emission tomography/computed tomography) equipment is to human body or animal body tumour, cardiac system
Disease and the nervous system disease are early diagnosed, it can the image scanned by PET/CT equipment to scanned object
Determine the position that lesion occurs.However, since imaging process has certain duration, scanned object is exhaled in PET imaging
It inhales movement or heartbeat makes to be occurred the tissue (such as lung) of lesion spatially periodic shift, resulting PET is caused to scheme
As obscuring, the small lesion in PET image is made it difficult to or can not identified.
It is fuzzy to eliminate PET image caused by respiratory movement, gating technology is often used in PET, i.e., by specific
(wherein, gate-control signal can be for reflecting scanned object in scanning the corresponding gate-control signal of equipment acquisition respiratory movement
The signal of respiratory movement process or heartbeat process in long (such as 12 minutes)), currently, commonly acquiring gate-control signal
Method is generally required by external equipment, for example, acquiring scanned object body surface by infrared collection device or video camera
Motor message, wherein motor message can be scanned object and be under lying status, and the height of body surface a part is with the time
The height value being changed, the collected motor message is gate-control signal.It is then possible to according to obtained gate-control signal,
Carry out PET image reconstruction, it can the respiratory cycle reflected according to gate-control signal, will duplicate same motion many times
In PET data be grouped together, the PET data being grouped together is redeveloped into the image of a motion stage.
In the implementation of the present invention, the inventor finds that the existing technology has at least the following problems:
When acquiring gate-control signal based on above-mentioned processing mode, need to increase additional external equipment, and the valence of external equipment
Lattice are often more expensive, thus, lead to diagnostic system increased costs.
Summary of the invention
In order to solve problems in the prior art, the embodiment of the invention provides the methods and dress of a kind of determining gate-control signal
It sets.The technical solution is as follows:
In a first aspect, providing a kind of method of determining gate-control signal, which comprises
It obtains and is detected in the sub-period of each preset duration of the positron emission tomography PET device when detecting in section
To respectively bury in oblivion dot position information;
Respectively bury in oblivion dot position information according to each sub-period is corresponding, determine the direction of motion of lesion tissue, determines
The motion amplitude information of the lesion tissue described in each sub-period in the direction of movement;
According to the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement, gate letter is determined
Number, PET image reconstruction is carried out based on the gate-control signal.
Optionally, described respectively to bury in oblivion dot position information according to each sub-period is corresponding, determine the fortune of lesion tissue
Dynamic direction determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement, comprising:
For each sub-period, respectively buries in oblivion dot position information according to the sub-period is corresponding, determine in the sub-period
Difference buries in oblivion that dot position information is corresponding to bury in oblivion a number;
Bury in oblivion that dot position information is corresponding to bury in oblivion a number according to difference in each sub-period, determines lesion tissue
The direction of motion determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement.
Optionally, described for each sub-period, respectively bury in oblivion dot position information according to the sub-period is corresponding, determines institute
It states difference in sub-period and buries in oblivion that dot position information is corresponding to bury in oblivion a number, comprising:
For each sub-period, according to the corresponding radial position respectively buried in oblivion in dot position information of the sub-period, angle
Position and axial position mark, determine that each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion points
Mesh, wherein the axial position mark is intended to indicate that the information of the axial position for the vanishing point that falls into oblivion;
It is described to bury in oblivion that dot position information is corresponding to bury in oblivion a number according to difference in each sub-period, determine lesion group
The direction of motion knitted determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement, comprising:
For each sub-period, by each radial position-angle position-axial position mark group described in the sub-period
A new radial position-angle is merged into the multiple radial position-angle positions for meeting preset proximity condition in conjunction
Position grouping, and bury in oblivion the sum of number for multiple radial positions-angle position combination is corresponding, as the new radial position
Set-angle position combination it is corresponding bury in oblivion a number, obtain radial position-angle position-axial direction position new in the sub-period
Set that identifier combination is corresponding to bury in oblivion a number;
New radial position-angle position-axial position identifier combination according to each sub-period is corresponding
Bury in oblivion a number, determine the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the direction of motion
On motion amplitude information.
Optionally, the new radial position-angle position according to each sub-period-axial position mark
Combine it is corresponding bury in oblivion a number, determine the direction of motion of lesion tissue, determine the lesion tissue described in each sub-period in institute
State the motion amplitude information in the direction of motion, comprising:
Data are carried out to radial position-angle position new described in each sub-period-axial position identifier combination
Reorganization, radial position-angle position-axial position identifier combination after obtaining data recombination, according to each period of the day from 11 p.m. to 1 a.m
New radial position-angle position-axial position identifier combination is corresponding described in section buries in oblivion a number, determines every height
Radial position-angle position-axial position identifier combination after data recombination described in period is corresponding to bury in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
Close it is corresponding bury in oblivion a number, determine the direction of motion of lesion tissue, determine the lesion tissue described in each sub-period described
Motion amplitude information in the direction of motion.
Optionally, radial position-angle position-axial direction after the data recombination according to each sub-period
Station location marker combination is corresponding to bury in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion described in each sub-period
The motion amplitude information of tissue in the direction of movement, comprising:
Radial position-angle position-after the data recombination described in formula and each sub-period that decayed according to nucleic
Axial position identifier combination is corresponding to bury in oblivion a number, the radial position after determining data recombination described in each sub-period
It sets the corresponding decay of-angle position-axial position identifier combination and revised buries in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, radial position-angle position-axial direction after the data recombination according to each sub-period
The corresponding decay of station location marker combination is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines in every height
The motion amplitude information of lesion tissue described in period in the direction of movement, comprising:
For radial position-angle position-axial position mark after each data recombination in each sub-period
Radial position-angle position combination in combination is known, after data recombination belonging to the radial position-angle position combination
Radial position-angle position-axial position identifier combination it is corresponding decay it is revised bury in oblivion a number, determine the radial direction
Combination corresponding decay in position-angle position is revised to bury in oblivion a number;
To decay the revised radial position-angle buried in oblivion a number and be less than preset number threshold value in each sub-period
Position grouping, radial position-angle position-axial position identifier combination decay after affiliated data recombination is revised to fall into oblivion
Vanishing point number is set as zero, obtains eliminating radial position-angle position-axial position mark after noise in each sub-period
The corresponding decay of knowledge combination is revised to bury in oblivion a number;
Radial position-angle position-axial position mark group after the elimination noise according to each sub-period
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, radial position-angle position-axial direction after the elimination noise according to each sub-period
The corresponding decay of station location marker combination is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines in every height
The motion amplitude information of lesion tissue described in period in the direction of movement, comprising:
According to preset principal component analysis PCA algorithm, to the radial direction after the elimination noise as described in each sub-period
The corresponding decay of position-angle position-axial position identifier combination is revised bury in oblivion the matrix that number is constituted carry out it is main at
Analysis, obtains at least one characteristic value and the corresponding feature vector of at least one described characteristic value, wherein the feature
It is bigger that number variation degree on the corresponding direction of feature vector of the characteristic value is buried in oblivion in the bigger expression of value;
The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, it will be described
The direction of motion of first eigenvector or the corresponding direction of the second feature vector as lesion tissue;
Radial position-angle position-axial direction for each sub-period, after the elimination noise according to the sub-period
Station location marker combination it is corresponding decay it is revised bury in oblivion a number and the first eigenvector or the second feature to
Amount determines the motion amplitude information of the lesion tissue described in the sub-period in the direction of movement.
Optionally, the motion amplitude letter according to the lesion tissue described in each sub-period in the direction of movement
Breath, determines gate-control signal, comprising:
To the motion amplitude of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information that information is constituted-function of time carries out interpolation processing, motion amplitude information-time after obtaining interpolation processing
Function;
Using motion amplitude information-function of time after the interpolation processing as gate-control signal.
Optionally, the motion amplitude letter according to the lesion tissue described in each sub-period in the direction of movement
Breath, determines gate-control signal, comprising:
To the motion amplitude of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information that information is constituted-function of time carries out frequency-domain transform processing, obtains the motion amplitude information-function of time
Corresponding frequency spectrum;
Based on the bandpass filter with default center frequency and pre-set bandwidths, the frequency spectrum is filtered, is filtered
Frequency spectrum after wave;
The inverse transformation that the frequency-domain transform is carried out to the filtered frequency spectrum, obtains gate-control signal.
Second aspect, provides a kind of device of determining gate-control signal, and described device includes:
Module is obtained, for obtaining each preset duration of the positron emission tomography PET device when detecting in section
Sub-period in detect respectively bury in oblivion dot position information;
First determining module determines lesion group for respectively burying in oblivion dot position information according to each sub-period is corresponding
The direction of motion knitted determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement;
Second determining module, for the movement width according to the lesion tissue described in each sub-period in the direction of movement
Information is spent, determines gate-control signal, PET image reconstruction is carried out based on the gate-control signal.
Optionally, first determining module, comprising:
First determines submodule, for respectively burying in oblivion point confidence according to the sub-period is corresponding for each sub-period
Breath determines in the sub-period that difference buries in oblivion that dot position information is corresponding to bury in oblivion a number;
Second determines submodule, for burying in oblivion the corresponding vanishing point that falls into oblivion of dot position information according to difference in each sub-period
Number determines the direction of motion of lesion tissue, determines the fortune of the lesion tissue described in each sub-period in the direction of movement
Dynamic amplitude information.
Optionally, it described first determines submodule, is used for:
For each sub-period, according to the corresponding radial position respectively buried in oblivion in dot position information of the sub-period, angle
Position and axial position mark, determine that each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion points
Mesh, wherein the axial position mark is intended to indicate that the information of the axial position for the vanishing point that falls into oblivion;
Described second determines submodule, is used for:
For each sub-period, by each radial position-angle position-axial position mark group described in the sub-period
A new radial position-angle is merged into the multiple radial position-angle positions for meeting preset proximity condition in conjunction
Position grouping, and bury in oblivion the sum of number for multiple radial positions-angle position combination is corresponding, as the new radial position
Set-angle position combination it is corresponding bury in oblivion a number, obtain radial position-angle position-axial direction position new in the sub-period
Set that identifier combination is corresponding to bury in oblivion a number;
New radial position-angle position-axial position identifier combination according to each sub-period is corresponding
Bury in oblivion a number, determine the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the direction of motion
On motion amplitude information.
Optionally, it described second determines submodule, is used for:
Data are carried out to radial position-angle position new described in each sub-period-axial position identifier combination
Reorganization, radial position-angle position-axial position identifier combination after obtaining data recombination, according to each period of the day from 11 p.m. to 1 a.m
New radial position-angle position-axial position identifier combination is corresponding described in section buries in oblivion a number, determines every height
Radial position-angle position-axial position identifier combination after data recombination described in period is corresponding to bury in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
Close it is corresponding bury in oblivion a number, determine the direction of motion of lesion tissue, determine the lesion tissue described in each sub-period described
Motion amplitude information in the direction of motion.
Optionally, it described second determines submodule, is used for:
Radial position-angle position-after the data recombination described in formula and each sub-period that decayed according to nucleic
Axial position identifier combination is corresponding to bury in oblivion a number, the radial position after determining data recombination described in each sub-period
It sets the corresponding decay of-angle position-axial position identifier combination and revised buries in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, it described second determines submodule, is used for:
For radial position-angle position-axial position mark after each data recombination in each sub-period
Radial position-angle position combination in combination is known, after data recombination belonging to the radial position-angle position combination
Radial position-angle position-axial position identifier combination it is corresponding decay it is revised bury in oblivion a number, determine the radial direction
Combination corresponding decay in position-angle position is revised to bury in oblivion a number;
To decay the revised radial position-angle buried in oblivion a number and be less than preset number threshold value in each sub-period
Position grouping, radial position-angle position-axial position identifier combination decay after affiliated data recombination is revised to fall into oblivion
Vanishing point number is set as zero, obtains eliminating radial position-angle position-axial position mark after noise in each sub-period
The corresponding decay of knowledge combination is revised to bury in oblivion a number;
Radial position-angle position-axial position mark group after the elimination noise according to each sub-period
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, it described second determines submodule, is used for:
According to preset principal component analysis PCA algorithm, to the radial direction after the elimination noise as described in each sub-period
The corresponding decay of position-angle position-axial position identifier combination is revised bury in oblivion the matrix that number is constituted carry out it is main at
Analysis, obtains at least one characteristic value and the corresponding feature vector of at least one described characteristic value, wherein the feature
It is bigger that number variation degree on the corresponding direction of feature vector of the characteristic value is buried in oblivion in the bigger expression of value;
The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, it will be described
The direction of motion of first eigenvector or the corresponding direction of the second feature vector as lesion tissue;
Radial position-angle position-axial direction for each sub-period, after the elimination noise according to the sub-period
Station location marker combination it is corresponding decay it is revised bury in oblivion a number and the first eigenvector or the second feature to
Amount determines the motion amplitude information of the lesion tissue described in the sub-period in the direction of movement.
Optionally, second determining module, is used for:
To the motion amplitude of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information that information is constituted-function of time carries out interpolation processing, motion amplitude information-time after obtaining interpolation processing
Function;
Using motion amplitude information-function of time after the interpolation processing as gate-control signal.
Optionally, second determining module, comprising:
Frequency-domain transform submodule is used for each sub-period and the lesion tissue described in each sub-period in the fortune
Motion amplitude information-function of time that motion amplitude information on dynamic direction is constituted carries out frequency-domain transform processing, obtains the fortune
The dynamic corresponding frequency spectrum of amplitude information-function of time;
Submodule is filtered, for based on the bandpass filter with default center frequency and pre-set bandwidths, to the frequency spectrum
It is filtered, obtains filtered frequency spectrum;
Frequency domain inverse transformation submodule is obtained for carrying out the inverse transformation of the frequency-domain transform to the filtered frequency spectrum
Gate-control signal.
Technical solution provided in an embodiment of the present invention has the benefit that
In the embodiment of the present invention, the dot position information of respectively burying in oblivion that can be detected according to PET device determines gate-control signal, nothing
Additional external equipment acquisition gate-control signal need to be increased when PET device detects and respectively buries in oblivion dot position information, thus, it can
To reduce diagnostic system cost.
Detailed description of the invention
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is a kind of method flow diagram of determining gate-control signal provided in an embodiment of the present invention;
Fig. 2 is a kind of structural schematic diagram of terminal provided in an embodiment of the present invention;
Fig. 3 is a kind of schematic diagram of crystal detection pair provided in an embodiment of the present invention;
Fig. 4 is a kind of apparatus structure schematic diagram of determining gate-control signal provided in an embodiment of the present invention;
Fig. 5 is a kind of apparatus structure schematic diagram of determining gate-control signal provided in an embodiment of the present invention;
Fig. 6 is a kind of apparatus structure schematic diagram of determining gate-control signal provided in an embodiment of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to embodiment party of the present invention
Formula is described in further detail.
The embodiment of the invention provides a kind of methods of determining gate-control signal, as shown in Figure 1, the executing subject of this method is
Terminal.Wherein, which can be with the terminal of gate-control signal function is determined, for example, can be and PET (positron
Emission tomography, positron emission tomography) equipment be connected, can according to PET device acquire data determine
The computer equipment of gate-control signal.
Terminal may include processor 210, transceiver 220, memory 230, input unit 240, display unit 250, sound
The components such as frequency circuit 260 and power supply 270, as shown in Fig. 2, it will be understood by those skilled in the art that terminal knot shown in Figure 2
The restriction of structure not structure paired terminal may include perhaps combining certain components or not than illustrating more or fewer components
Same component layout.Wherein:
Processor 210 can be the control centre of terminal, utilize each of various interfaces and the entire terminal device of connection
A part, such as transceiver 220 and memory 230, by run or execute the software program being stored in memory 230 and/
Or module, and the data being stored in memory 230 are called, the various functions and processing data of terminal are executed, thus to end
End carries out integral monitoring.Optionally, processor 210 may include one or more processing cores.In the present invention, processor 210
It is determined for the relevant treatment of gate-control signal.Transceiver 220 can be used for sending and receiving data, and terminal can pass through receipts
Hair device 220 sends and receivees data, and terminal can be network interface card by internet sending and receiving data, transceiver.
Memory 230 can be used for storing software program and module, and processor 210 is stored in memory 230 by operation
Software program and module, thereby executing various function application and data processing.Memory 230 can mainly include storage journey
Sequence area and storage data area, wherein storing program area can the (ratio of application program needed for storage program area, at least one function
Such as determine gate-control signal function) etc.;Storage data area, which can be stored, uses created data (such as the vanishing point that falls into oblivion according to terminal
Location information etc.) etc..In addition, memory 230 may include high-speed random access memory, it can also include non-volatile memories
Device, for example, at least a disk memory, flush memory device or other volatile solid-state parts.Input unit 240 can be with
Number for receiving input or character information, and generate keyboard related with user setting and function control, mouse, behaviour
Make bar, optics or trackball signal input.Display unit 250 is displayed for information input by user or is supplied to use
The information at family and the various graphical user interface of terminal, these graphical user interface can be by figure, text, icon, videos
It is constituted with any combination thereof.Display unit 250 may include display panel 251, optionally, can use LCD (Liquid
Crystal Display, liquid crystal display), OLED (Organic Light-Emitting Diode, Organic Light Emitting Diode)
Etc. forms configure display panel 251.Voicefrequency circuit 260, loudspeaker 261, microphone 262 can provide between user and terminal
The audio data received can be converted to electric signal by audio interface, voicefrequency circuit 260.Power supply 270 can pass through power management
System and processor 210 are logically contiguous, to realize management charging, electric discharge and power managed etc. by power-supply management system
Function.Power supply 270 can also include one or more direct current or AC power source, recharging system, power failure monitor
The random components such as circuit, power adapter or inverter, power supply status indicator.Terminal further includes having memory, Yi Jiyi
A perhaps more than one program one of them or more than one program is stored in memory, and is configured to by one
Either the execution of more than one processor states one or more than one program to execute determination door described in following each embodiments
The method for controlling signal.
Below in conjunction with specific embodiment, process flow shown in FIG. 1 is described in detail, content can be as
Under:
Step 101, the period of the day from 11 p.m. to 1 a.m of each preset duration of the positron emission tomography PET device when detecting in section is obtained
What is detected in section respectively buries in oblivion dot position information.
Wherein, the vanishing point that falls into oblivion can be positive electron and intracorporal negative electron buries in oblivion the location point for generating a pair of of photon, fall into oblivion
Vanishing point location information can be intended to indicate that the information for burying in oblivion a position, for example, can be almost while receiving a pair of of photon
Two crystal detections number and axial position mark, wherein axial position, which identifies, can be intended to indicate that the axis of the vanishing point that falls into oblivion
To the information of position, axial position mark can be indicated with the number of detector rings each in PET device, for example, each detection
Device ring includes 60 crystal detections, numbers and be followed successively by 1,2 ... 60, and PET device includes 15 detector rings altogether, and number is successively
1,2 ... 15, some, which buries in oblivion dot position information, can be (40,50,2), that is, indicate a pair of of photon in the vanishing point generation of falling into oblivion, quilt
No. 40 of 2nd detector rings and No. 50 crystal detections are almost detected simultaneously by.
Burying in oblivion dot position information can also be radial position, angle position and axial position mark, i.e., by radial position, angle
Degree position and axial position identify three parameter combinations and as one bury in oblivion dot position information, wherein by detector rings number and
The dot position information of burying in oblivion of axial position mark composition is buried in oblivion with what is be made of radial position, angle position and axial position mark
Dot position information have one-to-one relationship, detector rings number combined with radial position-angle position between can turn mutually
It changes, angle position can be two detector rings and number angle between line and vertical line between corresponding crystal detection,
For value between 0-180, radial position can be the center of circle of detector rings to the distance of the line between two crystal detections.
It should be noted is that the axial position mark buried in oblivion in dot position information that PET device detects can be used
The number of detector rings indicates.Specifically, when the crystal detection for being almost detected simultaneously by a pair of of photon is located at the same detection
When device ring, the corresponding axial position mark of the vanishing point that falls into oblivion is the number of the detector rings.But almost it is detected simultaneously by a pair
The crystal detection of photon is also possible to be located at different detector rings, as shown in figure 3, the corresponding axial position of vanishing point at this point, this falls into oblivion
The number combination for identifying two detector rings that can be located at crystal detection indicates, can also use the two pre-set spies
It surveys device ring and combines corresponding mark expression, wherein every two detector rings can be previously provided in local and combine corresponding axis
To station location marker, the number for the detector rings number that PET device is included may be greater than, for example, PET device shares 15 spies
Device ring is surveyed, number is successively 1-15, and can preset 16 expression, two crystal detections indicates positioned at the 1st, 2 detector rings, 17
Two crystal detections are located at the 1st, 3 detector rings, 18 two crystal detections of expression are located at the 1st, 4 detector rings, and so on, often
Two detector rings combinations are all corresponding with a number in advance.In this way, PET device detects when falling into oblivion vanishing point, it can be according to detection
The detector rings to be fallen into oblivion where the crystal detection of vanishing point to this distribute corresponding axial position mark for the vanishing point that falls into oblivion.
In an implementation, in medical domain, can by PET device to human body or animal body tumour, cardiac system disease and
The nervous system disease is early diagnosed, it can is determined and is occurred to the image that scanned object scans by PET device
The position of lesion.Specifically, first radionuclide can be injected (such as to scanned object18F、11C), wherein radioactive nucleus
Element is positron radionuclide, and radionuclide can enter intracorporal tissue with metabolism, and (tissue that lesion wherein, occurs can absorb greatly
The radionuclide of amount), then, positive electron will be buried in oblivion with intracorporal negative electron, a pair of heterodromous photon is generated,
Two crystal detections that this pair of heterodromous photon will be located at symmetric position detect, wherein each spy of PET device
The crystal detection on device ring with same number is surveyed, and each crystal detection in each detector rings has unique number,
When PET device often detects a pair of of photon, it can bury in oblivion the information (burying in oblivion dot position information) of a position to this is indicated and carry out
Storage, it can number and axial position mark to two crystal detections for being almost detected simultaneously by a pair of of photon are deposited
Storage, finally, terminal can be detected according to PET device bury in oblivion dot position information, form PET image.In turn, doctor can root
Determine whether the sweep object has the tissue that lesion occurs according to PET image, wherein PET device respectively buries in oblivion point confidence in detection
When breath, certain detection duration (general detection duration is about at 12 minutes or so) are needed, scanned object is swept in PET device
During retouching, scanned object will do it multiple respiratory movement or heartbeat, in such cases, if scanned object has generation disease
When the tissue of stove, due to the movement of lesion tissue, it will the PET image caused occurs fuzzy.
Currently, can solve the problems, such as that the PET image caused is moved due to lesion tissue to be obscured by gating technology,
Terminal will be specifically introduced and determines the detailed process of gate-control signal below.
The entire detection Time segments division of PET device can be multiple sub-periods with identical preset duration by terminal,
In, preset duration can be 0.5 second, be also possible to 1 second.The available PET device of terminal detects in section more when detecting
It is a to bury in oblivion dot position information, wherein what PET device detected respectively bury in oblivion, and dot position information can be according to respectively burying in oblivion point confidence
The list type data that corresponding detection time is arranged successively are ceased, in turn, terminal can be corresponding based on dot position information is respectively buried in oblivion
Detection time determines that is detected in the sub-period of each preset duration in section when detecting respectively buries in oblivion dot position information.
Step 102, respectively bury in oblivion dot position information according to each sub-period is corresponding, determine the direction of motion of lesion tissue,
Determine the motion amplitude information in each sub-period lesion tissue in the movement direction.
Wherein, motion amplitude information can be the information of the motion amplitude of reflection lesion tissue in the movement direction.
In an implementation, for each sub-period, terminal can respectively bury in oblivion dot position information according to each sub-period is corresponding,
It determines the direction of motion of lesion tissue, and respectively buries in oblivion dot position information according to each sub-period is corresponding, determine lesion tissue
Motion amplitude information in the movement direction changes with time trend, that is, determines lesion tissue in the movement width of each sub-period
Spend information.That is, when burying in oblivion dot position information and including crystal detection number and axial position mark, it can be according to respectively falling into oblivion
The corresponding crystal detection number of vanishing point and axial position mark, determine the direction of motion of lesion tissue, and determine in every height
The motion amplitude information of period lesion tissue in the movement direction.It include radial position, angle position when burying in oblivion dot position information
When with axial position mark, it can be identified according to the corresponding radial position of the vanishing point that respectively falls into oblivion, angle position and axial position, determine disease
The direction of motion of stove tissue, and determine the motion amplitude information in each sub-period lesion tissue in the movement direction.
Optionally, (points can be buried in oblivion according to the different frequency of occurrence for burying in oblivion dot position information in each sub-period
Mesh), determine the direction of motion of lesion tissue, and determine the lesion tissue described in each sub-period in the direction of movement
Motion amplitude information, correspondingly, the treatment process of step 102 can be such that for each sub-period, it is corresponding according to sub-period
Respectively bury in oblivion dot position information, determines in sub-period that difference buries in oblivion that dot position information is corresponding to bury in oblivion a number;According to each period of the day from 11 p.m. to 1 a.m
Difference buries in oblivion that dot position information is corresponding to bury in oblivion a number in section, determines the direction of motion of lesion tissue, determines in each period of the day from 11 p.m. to 1 a.m
The motion amplitude information of Duan Suoshu lesion tissue in the movement direction.
In an implementation, terminal gets that each sub-period is corresponding respectively to bury in oblivion dot position information (wherein, each sub-period pair
Respectively burying in oblivion in dot position information for answering buries in oblivion dot position information comprising duplicate) after, different it can bury in oblivion a position to each
The frequency of occurrence of information is counted, and obtains each different frequency of occurrence for burying in oblivion dot position information (in i.e. each sub-period
Difference buries in oblivion that dot position information is corresponding to bury in oblivion a number), the corresponding sinogram Sinogram number of each sub-period can be obtained
It may include each different burying in oblivion dot position information and its corresponding frequency of occurrence (falling into oblivion vanishing point according to, wherein Sinogram data
Number), Sinogram data can be used to reflect the distribution situation for the vanishing point that falls into oblivion.For example, the first sub-period buries in oblivion point confidence
Breath include (1,10,1), (5,10,1), (1,10,3), (10,20,1), (1,10,3), (10,20,1), (5,10,1), (10,
20,1)、(1,10,3)、(1,10,3)、(1,10,2)、(1,10,1)、(1,10,2)、(5,10,1)、(1,10,1)、(5,10,
1), the frequency of occurrence for burying in oblivion dot position information of the sub-period is counted, the information that obtained Sinogram data include
Are as follows: 3 (1,10,1), 4 (5,10,1), 4 (1,10,3), 3 (10,20,1), 2 (1,10,2), wherein Sinogram
Data can be in table form record each bury in oblivion dot position information and its it is corresponding bury in oblivion a number, need to illustrate
It is that in actual operation, each sub-period, which will include, largely buries in oblivion dot position information, and above-mentioned example is simply to illustrate that by arranging
Treatment process of the list-directed data to Sinogram data.It, can basis after terminal obtains the Sinogram data of each sub-period
The Sinogram data of each sub-period determine the direction of motion of lesion tissue, determine the lesion tissue described in each sub-period
Motion amplitude information in the movement direction.
Optionally, terminal can be according to the corresponding radial position respectively buried in oblivion in dot position information of each sub-period, angle
Position and axial position mark, determine that difference buries in oblivion that dot position information is corresponding to bury in oblivion a number, correspondingly, treatment process can be with
It is as follows: for each sub-period, according to the corresponding radial position respectively buried in oblivion in dot position information of sub-period, angle position and axis
To station location marker, determines that each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion a number, wherein axis
The information of the axial position for the vanishing point that falls into oblivion is intended to indicate that station location marker.Correspondingly, treatment process can be such that for every height
Period will meet the more of preset proximity condition in radial position-angle position-axial position identifier combination each in sub-period
New radial position-angle position combination is merged into a radial position-angle position, and by multiple radial positions-
Angle position combination is corresponding to bury in oblivion the sum of number, buries in oblivion points as new radial position-angle position combination is corresponding
Mesh, obtains that radial position-angle position-axial position identifier combination new in sub-period is corresponding to bury in oblivion a number;According to every
New radial position-angle position-axial position identifier combination is corresponding in a sub-period buries in oblivion a number, determines lesion tissue
The direction of motion, determine motion amplitude information in each sub-period lesion tissue in the direction of movement.
In an implementation, after terminal obtains that each sub-period is corresponding and buries in oblivion dot position information, wherein bury in oblivion dot position information
It can be radial position, angle position and axial position mark, each radial position-angle position-axial position can be counted
The frequency of occurrence of identifier combination, it can count each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion
Point number, also can be according to the corresponding radial position respectively buried in oblivion in dot position information of each sub-period, angle position and axis
To station location marker, the corresponding Sinogram data of each sub-period are determined, wherein the format of Sinogram data can be in advance
Storage, one radial position-angle position-axial position identifier combination of every determination is corresponding bury in oblivion a number after, can will
This buries in oblivion the corresponding position that a number is added in Sinogram data, i.e., includes in the Sinogram data of each sub-period
Each radial position-angle position-axial position identifier combination is fixed, each radial position-angle position-axial position
The corresponding number of burying in oblivion of identifier combination is variation, is respectively to bury in oblivion what dot position information was determined according to each sub-period,
Sinogram data can be the format of three-dimensional table.For example, corresponding sub- Sinogram data are identified to Mr. Yu's axial position,
As shown in table 1, wherein arranged according to the sequence that from left to right, numerical value is sequentially increased from top to bottom, each radial position and angle
The corresponding number of position grouping is that this buries in oblivion that dot position information is corresponding to bury in oblivion a number.For example, being identified for axial position
18, position (0,0) corresponding 5 is to indicate that radial position is 0, angle position 0, axial position are identified as 18 and bury in oblivion a position
The corresponding number of burying in oblivion of information is 5.
Table 1
After terminal obtains the corresponding Sinogram data of each sub-period, corresponding son is identified for each axial position
Sinogram data can make sub- Sinogram data according to preset proximity condition, along radial position direction and angle
Locality is compressed.Specifically, terminal can be radial by the preset number per adjacent (meeting default close to condition)
A new radial position-angle position combination is merged into position-angle position, wherein new radial position-angle position
Set the mean value for the radial position that the radial position in combination can be in multiple radial positions-angle position combination, new radial direction
Angle position in position-angle position combination can be the equal of the angle position in multiple radial positions-angle position combination
Value, and can preset number radial position-angle position combine a corresponding number of burying in oblivion and be added, what is obtained falls into oblivion
The sum of vanishing point number buries in oblivion a number as new radial position-angle position combination is corresponding, for each axial position mark
Know, terminal all can carry out data compression along radial position direction and angle position direction according to above-mentioned data compression method, obtain
Into sub-period, new radial position-angle position-axial position identifier combination is corresponding buries in oblivion a number.For example, such as 1 institute of table
The sub- Sinogram data shown, preset number can be 4, and (0,0), (1,0), (0,5), (1,5) can be merged into one by terminal
(2,0), (3,0), (2,5), (3,5) are merged into one newly by a new radial position-angle position combination-(0.5,2.5)
Radial position-angle position combination-(2.5,2.5), merges into a new diameter for (0,10), (1,10), (0,12), (1,12)
To position-angle position combination-(0.5,11), a new radial direction is merged into (2,10), (3,10), (2,12), (3,12)
Position-angle position combination-(2.5,11).In this way, the radial position of each of Sinogram data after data compression can be made
It sets the corresponding number of burying in oblivion of-angle position-axial position identifier combination to increase, it can improve each radial position-angle
Position-axial position identifier combination data statistics amount can reduce the interference of noise in turn.
It, can be according to new in each sub-period after terminal carries out data compression to the Sinogram data of each sub-period
Radial position-angle position-axial position identifier combination is corresponding to bury in oblivion a number, determines the direction of motion of lesion tissue, really
It is scheduled on the motion amplitude information of lesion tissue described in each sub-period in the direction of movement.Specifically, can be according to default
Principal component analysis PCA algorithm, to by radial position-angle position-axial position identifier combination pair new in each sub-period
The matrix that a number is constituted of burying in oblivion answered carries out principal component analysis, obtains at least one characteristic value and at least one characteristic value
Corresponding feature vector, wherein the bigger expression of characteristic value is buried in oblivion a number and become on the corresponding direction of feature vector of characteristic value
Change degree is bigger;The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, it will
The direction of motion of first eigenvector or the corresponding direction of second feature vector as lesion tissue;For each sub-period, root
Bury in oblivion a number and fisrt feature according to radial position-angle position-axial position identifier combination new in sub-period is corresponding
Vector or second feature vector determine the motion amplitude information in sub-period lesion tissue in the movement direction.
Optionally, terminal can also be to radial position-angle position-axial position identifier combination new in each sub-period
Data recombination processing is carried out, correspondingly, treatment process can be such that radial position-angle position-new in each sub-period
Axial position identifier combination carries out data recombination processing, radial position-angle position-axial position mark after obtaining data recombination
Know combination, buries in oblivion points according to radial position-angle position-axial position identifier combination new in each sub-period is corresponding
Mesh determines the corresponding vanishing point that falls into oblivion of radial position-angle position-axial position identifier combination in each sub-period after data recombination
Number;Bury in oblivion according to axial position identifier combination is corresponding radial position-angle position-after data recombination in each sub-period
Point number, determines the direction of motion of lesion tissue, determines the movement in each sub-period lesion tissue in the direction of movement
Amplitude information.
In an implementation, for each sub-period, after terminal carries out data compression to the Sinogram data of sub-period, first
It can be according to preset data recombination algorithm (such as SSRB (Single Slice Rebin, single layer recombination) algorithm, FORE
(Fourier Rebin, Fourier's recombination) algorithm), the spy of photon pair is indicated almost while received to bury in oblivion dot position information
Device crystal is surveyed to be located at the vanishing points that falls into oblivion of different detector rings (i.e. axial position mark buries in oblivion point confidence greater than detector rings number
It ceases corresponding fall into oblivion vanishing point or axial position mark and buries in oblivion dot position information correspondence with what the number combination of two detector rings indicated
The vanishing point that falls into oblivion), radial position-angle position-axial position identifier combination after determining data recombination, that is to say, that Ke Yigen
According to preset data recombination algorithm, bury in oblivion dot position information after determining these corresponding data recombinations of vanishing point that falls into oblivion, for example, PET
Equipment includes 15 detector rings altogether, and the corresponding dot position information of burying in oblivion of certain vanishing point that falls into oblivion is (1,10,17), wherein is preset
Axial position mark 17 indicate that two different detector rings are 1,3, can be according to preset data recombination algorithm, determining should
The corresponding new axial position of the vanishing point that falls into oblivion is identified as 2, i.e., after data recombination is handled, which buries in oblivion dot position information
For (10,20,2), that is, the axial position that will be greater than detector rings number identifies the axial position mark be converted to less than detector number of rings
Know.Then, terminal can each radial position-angle position-axial position identifier combination to data recombination and without data weight
The corresponding number of burying in oblivion of radial position-angle position-axial position identifier combination of group is reintegrated, and determines every height
Radial position-angle position-axial position identifier combination in period after data recombination is corresponding to bury in oblivion a number to get to number
According to the Sinogram data after recombination.
It, can be according to the diameter after data recombination in each sub-period after terminal obtains the Sinogram data after data recombination
Bury in oblivion a number to position-angle position-axial position identifier combination is corresponding, determine the direction of motion of lesion tissue, determines
In the motion amplitude information of each sub-period lesion tissue in the movement direction.Specifically, can be according to preset principal component point
PCA algorithm is analysed, to corresponding by radial position-angle position after data recombination in each sub-period-axial position identifier combination
The matrix that number is constituted of burying in oblivion carry out principal component analysis, obtain at least one characteristic value and at least one characteristic value pair
The feature vector answered, wherein the bigger expression of characteristic value is buried in oblivion a number and changed on the corresponding direction of feature vector of characteristic value
Degree is bigger;The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, by the
The direction of motion of one feature vector or the corresponding direction of second feature vector as lesion tissue;For each sub-period, according to
Radial position-angle position-axial position identifier combination in sub-period after data recombination is corresponding to bury in oblivion a number and
One feature vector or second feature vector determine the motion amplitude information in sub-period lesion tissue in the movement direction.
Optionally, it is contemplated that the influence that nucleic decays to statistical data obtains the diameter in each sub-period after data recombination
To position-angle position-axial position identifier combination is corresponding bury in oblivion a number after, corresponding to each sub-period can respectively fall into oblivion
Vanishing point number carries out decay amendment, correspondingly, treatment process can be such that according to number in nucleic decay formula and each sub-period
Bury in oblivion a number according to radial position-angle position-axial position identifier combination after recombination is corresponding, determines in each sub-period
The corresponding decay of radial position-angle position-axial position identifier combination after data recombination is revised to bury in oblivion a number;Root
After the corresponding decay amendment of radial position-angle position-axial position identifier combination after data recombination in each sub-period
Bury in oblivion a number, determine the direction of motion of lesion tissue, determine each sub-period lesion tissue in the direction of movement
Motion amplitude information.
In an implementation, for each sub-period, radial position-angle position-under the sub-period after data recombination is obtained
Axial position identifier combination is corresponding bury in oblivion a number after, can according to nucleic decay formula, to data recombination under the sub-period
The corresponding number of burying in oblivion of radial position-angle position-axial position identifier combination afterwards carries out decay amendment, obtains the period of the day from 11 p.m. to 1 a.m
The corresponding decay of radial position-angle position-axial position identifier combination after the lower data recombination of section is revised to bury in oblivion points
Mesh, specifically, radial position-angle position-axial direction position in each sub-period after data recombination can be obtained according to formula (1)
Set identifier combination it is corresponding decay it is revised bury in oblivion a number,
Wherein, S ' (t) is radial position-angle position-axial position identifier combination pair under sub-period t after data recombination
The decay answered is revised to bury in oblivion a number, and S (t) is radial position-angle position-axial direction under sub-period t after data recombination
Station location marker, which combines, buries in oblivion a number, T before corresponding decay is corrected1/2For the half-life period of nucleic,For nucleic
Decay formula.
Terminal obtains radial position-angle position-axial position identifier combination pair in each sub-period after data recombination
The decay answered is revised bury in oblivion a number after, can be according to the radial position after data recombination in each sub-period-angle position
Set-the corresponding decay of axial position identifier combination is revised buries in oblivion a number, it determines the direction of motion of lesion tissue, determines and exist
The motion amplitude information of each sub-period lesion tissue in the direction of movement.Specifically, can be according to preset principal component
PCA algorithm is analyzed, to by radial position-angle position-axial position identifier combination pair after data recombination in each sub-period
The revised matrix progress principal component analysis burying in oblivion a number and constituting of the decay answered, obtains at least one characteristic value, Yi Jizhi
Few corresponding feature vector of a characteristic value, wherein the bigger expression of characteristic value buries in oblivion a number in the feature vector pair of characteristic value
Variation degree is bigger on the direction answered;Choose the corresponding first eigenvector of maximum eigenvalue or secondary big characteristic value corresponding second
Feature vector, using first eigenvector or the corresponding direction of second feature vector as the direction of motion of lesion tissue;For every
A sub-period, according to radial position-angle position after data recombination in sub-period-corresponding decay of axial position identifier combination
It is revised to bury in oblivion a number and first eigenvector or second feature vector, it determines in sub-period lesion tissue in movement side
Upward motion amplitude information.
Optionally, it is contemplated that influence of the noise to the direction of motion is determined obtains the diameter in each sub-period after data recombination
To position-angle position-axial position identifier combination is corresponding bury in oblivion a number after, corresponding to each sub-period can respectively fall into oblivion
Vanishing point number carries out denoising, correspondingly, treatment process can be such that for after data recombination each in each sub-period
Radial position-angle position combination in radial position-angle position-axial position identifier combination, according to radial position-angle
After the corresponding decay amendment of radial position-angle position-axial position identifier combination after data recombination belonging to position grouping
Bury in oblivion a number, determine that the corresponding decay of radial position-angle position combination is revised and bury in oblivion a number;By each period of the day from 11 p.m. to 1 a.m
Decay the revised radial position-angle position combination buried in oblivion a number and be less than preset number threshold value in section, affiliated number
It is set as zero according to the revised number of burying in oblivion of radial position-angle position-axial position identifier combination decay after recombination,
It obtains eliminating the corresponding decay amendment of radial position-angle position-axial position identifier combination after noise in each sub-period
Afterwards bury in oblivion a number;According to radial position-angle position-axial position identifier combination after elimination noise in each sub-period
Corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines in each sub-period lesion tissue
Motion amplitude information in the movement direction.
In an implementation, quantity threshold can be previously stored in terminal, wherein each sub- period can be corresponding with accordingly
Quantity threshold, specifically, can be previously stored with percentage in terminal, the corresponding quantity threshold of each sub-period is the son
Maximum in period buries in oblivion the product of a number and preset percentage.What PET device acquired, which buries in oblivion dot position information, to include
Noise information, it is contemplated that noise information to determine the direction of motion influence, can be to each data recombination in each sub-period after
Radial position-angle position-corresponding attenuation correction of axial position identifier combination after the number buried in oblivion carry out a denoising,
Specifically, firstly for radial position-angle position-axial position identifier combination after each data recombination in each sub-period
In radial position-angle position combination, can be to the radial position after data recombination belonging to radial position-angle position combination
It sets the corresponding revised number of burying in oblivion of decay of-angle position-axial position identifier combination to be added, as radial direction
Combination corresponding decay in position-angle position is revised to bury in oblivion a number, then, for each sub-period, judges radial position
Set-angle position combines the revised size relation for burying in oblivion a number Yu preset number threshold value of corresponding decay, determine the son
Decay the revised radial position-angle position combination buried in oblivion a number and be less than preset number threshold value in period, in turn, can
To determine that the revised number of burying in oblivion of decay is less than number belonging to radial position-angle position combination of preset number threshold value
It, can be by the radial position after the data recombination belonging to it according to radial position-angle position-axial position identifier combination after recombination
It sets the revised number of burying in oblivion of the corresponding decay of-angle position-axial position identifier combination and is set as zero, decay revised
Bury in oblivion radial position-angle position combination of the number more than preset number threshold value, the radial position-after affiliated data recombination
The revised number of burying in oblivion of angle position-axial position identifier combination decay remains unchanged, and obtains disappearing in each sub-period
Bury in oblivion a number except the corresponding decay of radial position-angle position-axial position identifier combination after noise is revised.Most
Afterwards, terminal can be corresponding according to radial position-angle position-axial position identifier combination after elimination noise in each sub-period
Decay it is revised bury in oblivion a number, determine the direction of motion of lesion tissue, determination is being transported in each sub-period lesion tissue
Motion amplitude information on dynamic direction.
Optionally, it can use Principal Component Analysis Algorithm, determine the direction of motion of lesion tissue and in the direction of motion
On motion amplitude information, correspondingly, treatment process can be such that according to preset principal component analysis PCA algorithm, to by each
Corresponding revised bury in oblivion that decay of radial position-angle position-axial position identifier combination after eliminating noise in sub-period
The matrix that point number is constituted carries out principal component analysis, obtains at least one characteristic value and the corresponding spy of at least one characteristic value
Levy vector, wherein a number is buried in oblivion in the bigger expression of characteristic value, and variation degree is more on the corresponding direction of feature vector of characteristic value
Greatly;The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, by fisrt feature
The direction of motion of vector or the corresponding direction of second feature vector as lesion tissue;For each sub-period, according to sub-period
Described in the revised vanishing point that falls into oblivion of the corresponding decay of radial position-angle position-axial position identifier combination after elimination noise
Number and first eigenvector or second feature vector determine the movement width in sub-period lesion tissue in the movement direction
Spend information.
In an implementation, corresponding each sub-period, terminal obtain eliminating radial position-angle position after noise in the sub-period
Set-axial position identifier combination it is corresponding decay it is revised bury in oblivion a number after, it can be put in order according to certain,
The vector that n × 1 is tieed up is formed, X (t can be usedi) indicate, wherein n is positive integer, and n can be the diameter in Sinogram data format
The product of the quantity identified to the quantity, the quantity of angle position and axial position of position, tiIndicate i-th of sub-period, i=1,
2 ... M, wherein M is the number of sub-period, is the corresponding radial position-angle eliminated after noise of the dimension per one-dimensional numerical value
The corresponding decay of degree position-axial position identifier combination is revised to bury in oblivion a number, in turn, can be by the X of each sub-period
(ti) composition matrix, it can be indicated with matrix S, i.e. S=[X (t1), X (t2), X (t3)…X(tM)].In turn, terminal can use
PCA (Principal Component Analysis, principal component analysis) algorithm carries out principal component analysis to matrix S.Specifically
, after terminal obtains matrix S, each data of the every row of matrix S can be subtracted with the mean value of line number value each in calculating matrix S
The corresponding mean value of the row, obtains matrix S ', in turn, can according to the covariance matrix Σ of formula (2) calculating matrix S ',
Wherein, S 'TThe transposed matrix of representing matrix S ' can calculate the characteristic value and feature of covariance matrix Σ in turn
Vector, wherein it is bigger that number variation degree on the corresponding direction of feature vector of characteristic value is buried in oblivion in the bigger expression of characteristic value,
In the present solution, the corresponding direction of the corresponding feature vector of maximum eigenvalue can be the corresponding fortune of the maximum movement of exercise intensity
Direction is moved (due to the variation degree of the size response data of characteristic value, since the movements such as respiratory movement or heartbeat lead to disease
Stove histokinesis, in turn, the distribution for the vanishing point that causes to fall into oblivion change, i.e. the data of each position in matrix S change, because
This, can be using the direction of the corresponding feature vector of the maximum eigenvalue obtained according to matrix S as the maximum movement of exercise intensity
The corresponding direction of motion), and so on, the direction of motion of certain movement can be determined according to the feature vector of acquisition, for example, generally
In the case of, respirometric exercise intensity be greater than heartbeat exercise intensity, can by the corresponding feature of maximum eigenvalue to
Corresponding direction is measured as the respirometric direction of motion, by time corresponding direction of the corresponding feature vector of big characteristic value as the heart
The direction of motion of dirty beating in special circumstances, can if the exercise intensity of heartbeat is greater than respirometric exercise intensity
Using the direction of motion by the corresponding direction of the corresponding feature vector of maximum eigenvalue as heartbeat, time big characteristic value is corresponded to
The corresponding direction of feature vector as the respirometric direction of motion.In turn, the available maximum eigenvalue of terminal is corresponding
First eigenvector or the corresponding second feature vector of secondary big characteristic value, first eigenvector or second feature vector is corresponding
The direction of motion of the direction as lesion tissue, it can as needed, choose the corresponding direction of one of feature vector and make
For the direction of motion of lesion tissue, it is preferred that can be using the corresponding direction of first eigenvector as the movement side of lesion tissue
To.
After terminal determines the direction of motion of lesion tissue, above-mentioned each of matrix S ' can be calculated and be listed in first eigenvector
Or the projection of second feature vector, it can calculate first eigenvector or second feature vector and subtract the X after mean value
(ti) dot product, wherein the corresponding projection of each sub-period that can will be obtained, as in sub-period lesion tissue in movement side
Upward motion amplitude information, can be denoted as W (t), i.e. W (t) can be used to indicate that lesion tissue in the motion amplitude of t sub-period
Information.
Step 103, according to the motion amplitude information in each sub-period lesion tissue in the movement direction, gate letter is determined
Number, PET image reconstruction is carried out based on gate-control signal.
Wherein, gate-control signal can be the motor message determined by the corresponding motion amplitude information of each sub-period.
In an implementation, terminal determines lesion tissue in the motion amplitude information in the movement direction of each sub-period
It afterwards, can be by the corresponding motion amplitude information of each sub-period, as gate-control signal.It can extract and be swept from gate-control signal
Corresponding each moment when (or when air-breathing maximum) retouch object expiration maximum, by the duration equal part between each adjacent two moment
For the sub- duration of preset number, the dot position information of burying in oblivion detected in the period of the day from 11 p.m. to 1 a.m corresponding in each duration length is divided in one
It rises, wherein it is considered that the dot position information of burying in oblivion being grouped together is detection when scanned object is in same motion
It arrives, in turn, PET image reconstruction can be carried out according to the dot position information of burying in oblivion being grouped together, obtain the corresponding sports stage
PET image.For example, scanned object is t1, t2, t3, t4, t5 respectively at the time of correspondence when exhaling maximum, respectively by t2-
The duration of t1, t3-t2, t4-t3, t5-t4 are divided into 3 equal parts, and number is respectively 1-3, can will number the period of the day from 11 p.m. to 1 a.m for being 1 all
The dot position information of respectively burying in oblivion detected in length is grouped together, and is buried in oblivion a position according to these and is carried out PET image reconstruction, obtains
The PET image of the corresponding motion stage divides the dot position information of respectively burying in oblivion detected in the period of the day from 11 p.m. to 1 a.m that all numbers are 2 is long
Together, bury in oblivion a position according to these and carry out PET image reconstruction, the PET image for corresponding to the motion stage is obtained, in all volumes
Number for 3 the period of the day from 11 p.m. to 1 a.m it is long in the dot position information of respectively burying in oblivion that detects be grouped together, bury in oblivion a position according to these and carry out a PET figure
As reconstruction, the PET image for corresponding to the motion stage is obtained, in turn, the PET image of available multiple motion stages.
Optionally, terminal can be by the motion amplitude information Jing Guo interpolation processing, as gate-control signal, correspondingly, step
103 treatment process can be such that the movement width to each sub-period and in each sub-period lesion tissue in the movement direction
Spend information constitute motion amplitude information-function of time carry out interpolation processing, after obtaining interpolation processing motion amplitude information-when
Between function;Using motion amplitude information-function of time after interpolation processing as gate-control signal.
In an implementation, terminal obtains terminal and determines lesion tissue in the movement width in the movement direction of each sub-period
After spending information W (t), interpolation processing can be carried out to W (t) based on preset interpolation algorithm, the movement width after obtaining interpolation processing
Information-function of time W ' (t) is spent, specifically, linear interpolation can be carried out to W (t), passes through the corresponding motion amplitude of adjacent moment
Information determines the corresponding motion amplitude information of the intermediate time of the adjacent moment, can also use non-linear interpolation in this programme
Algorithm is not limited thereof.Motion amplitude information-function of time W ' (t) after obtaining interpolation processing can make W ' (t)
For gate-control signal.In this way, at the time of can more accurately obtaining maximum range of motion information and correspond to, to carry out PET image weight
It builds.
Optionally, can also will by the corresponding motion amplitude information of filtered each sub-period as gate-control signal,
Correspondingly, the treatment process of step 103 can be such that each sub-period and in each sub-period lesion tissue in the direction of motion
On motion amplitude information constitute motion amplitude information-function of time carry out frequency-domain transform processing, obtain motion amplitude information-
The corresponding frequency spectrum of the function of time;Based on the bandpass filter with default center frequency and pre-set bandwidths, frequency spectrum is filtered,
Obtain filtered frequency spectrum;The inverse transformation that frequency-domain transform is carried out to filtered frequency spectrum, obtains gate-control signal.
In an implementation, after terminal obtains the corresponding motion amplitude information W (t) of each sub-period t, W (t) can be carried out fast
Fast Fourier transformation (FFT) obtains the frequency spectrum of W (t) in turn, wherein and the frequency spectrum of W (t) can be calculated according to formula (3),
H (f)=| FFT (W (t)) | (3)
Wherein, H (f) is the frequency spectrum of W (t), and f indicates frequency, and in turn, terminal can be to the frequency of motion amplitude information W (t)
Spectrum is filtered, wherein can pass through bandpass filter with default center frequency and pre-set bandwidths (such as gaussian filtering
Device) it is filtered, filtered frequency spectrum is obtained, later, inversefouriertransform processing is carried out to filtered frequency spectrum, is obtained
Gate-control signal.
In addition, default center frequency can be when motion amplitude information W (t) is respiratory movement corresponding motor message
When default respiratory rate value (for example, default respiratory rate value can be 0.5Hz), obtained gate-control signal is corresponding for respiratory movement
Gate-control signal, when motion amplitude information W (t) is the corresponding motor message of heartbeat, default center frequency is the default heart
When jumping frequency values (for example, default palmic rate value can be 1Hz), obtained gate-control signal is the corresponding gate letter of heartbeat
Number.
In the embodiment of the present invention, the dot position information of respectively burying in oblivion that can be detected according to PET device determines gate-control signal, nothing
Additional external equipment acquisition gate-control signal need to be increased when PET device detects and respectively buries in oblivion dot position information, thus, it can
To reduce diagnostic system cost.
Based on the same technical idea, the embodiment of the invention also provides a kind of devices of determining gate-control signal, such as Fig. 4 institute
Show, which includes:
Module 410 is obtained, when each presetting in section when detecting for obtaining positron emission tomography PET device
What is detected in long sub-period respectively buries in oblivion dot position information;
First determining module 420 determines lesion for respectively burying in oblivion dot position information according to each sub-period is corresponding
The direction of motion of tissue determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement;
Second determining module 430, for the fortune according to the lesion tissue described in each sub-period in the direction of movement
Dynamic amplitude information, determines gate-control signal, carries out PET image reconstruction based on the gate-control signal.
Optionally, as shown in figure 5, first determining module 420, comprising:
First determines submodule 4201, for respectively burying in oblivion point according to the sub-period is corresponding for each sub-period
Confidence breath determines in the sub-period that difference buries in oblivion that dot position information is corresponding to bury in oblivion a number;
Second determines submodule 4202, and for being buried in oblivion according to difference in each sub-period, dot position information is corresponding to fall into oblivion
Vanishing point number determines the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the direction of movement
Motion amplitude information.
Optionally, it described first determines submodule 4201, is used for:
For each sub-period, according to the corresponding radial position respectively buried in oblivion in dot position information of the sub-period, angle
Position and axial position mark, determine that each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion points
Mesh, wherein the axial position mark is intended to indicate that the information of the axial position for the vanishing point that falls into oblivion;
Described second determines submodule 4202, is used for:
For each sub-period, by each radial position-angle position-axial position mark group described in the sub-period
A new radial position-angle is merged into the multiple radial position-angle positions for meeting preset proximity condition in conjunction
Position grouping, and bury in oblivion the sum of number for multiple radial positions-angle position combination is corresponding, as the new radial position
Set-angle position combination it is corresponding bury in oblivion a number, obtain radial position-angle position-axial direction position new in the sub-period
Set that identifier combination is corresponding to bury in oblivion a number;
New radial position-angle position-axial position identifier combination according to each sub-period is corresponding
Bury in oblivion a number, determine the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the direction of motion
On motion amplitude information.
Optionally, it described second determines submodule 4202, is used for:
Data are carried out to radial position-angle position new described in each sub-period-axial position identifier combination
Reorganization, radial position-angle position-axial position identifier combination after obtaining data recombination, according to each period of the day from 11 p.m. to 1 a.m
New radial position-angle position-axial position identifier combination is corresponding described in section buries in oblivion a number, determines every height
Radial position-angle position-axial position identifier combination after data recombination described in period is corresponding to bury in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
Close it is corresponding bury in oblivion a number, determine the direction of motion of lesion tissue, determine the lesion tissue described in each sub-period described
Motion amplitude information in the direction of motion.
Optionally, it described second determines submodule 4202, is used for:
Radial position-angle position-after the data recombination described in formula and each sub-period that decayed according to nucleic
Axial position identifier combination is corresponding to bury in oblivion a number, the radial position after determining data recombination described in each sub-period
It sets the corresponding decay of-angle position-axial position identifier combination and revised buries in oblivion a number;
According to radial position-angle position after data recombination described in each sub-period-axial position mark group
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, it described second determines submodule 4202, is used for:
For radial position-angle position-axial position mark after each data recombination in each sub-period
Radial position-angle position combination in combination is known, after data recombination belonging to the radial position-angle position combination
Radial position-angle position-axial position identifier combination it is corresponding decay it is revised bury in oblivion a number, determine the radial direction
Combination corresponding decay in position-angle position is revised to bury in oblivion a number;
To decay the revised radial position-angle buried in oblivion a number and be less than preset number threshold value in each sub-period
Position grouping, radial position-angle position-axial position identifier combination decay after affiliated data recombination is revised to fall into oblivion
Vanishing point number is set as zero, obtains eliminating radial position-angle position-axial position mark after noise in each sub-period
The corresponding decay of knowledge combination is revised to bury in oblivion a number;
Radial position-angle position-axial position mark group after the elimination noise according to each sub-period
It closes that corresponding decay is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determine the disease described in each sub-period
The motion amplitude information of stove tissue in the direction of movement.
Optionally, it described second determines submodule 4202, is used for:
According to preset principal component analysis PCA algorithm, to the radial direction after the elimination noise as described in each sub-period
The corresponding decay of position-angle position-axial position identifier combination is revised bury in oblivion the matrix that number is constituted carry out it is main at
Analysis, obtains at least one characteristic value and the corresponding feature vector of at least one described characteristic value, wherein the feature
It is bigger that number variation degree on the corresponding direction of feature vector of the characteristic value is buried in oblivion in the bigger expression of value;
The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, it will be described
The direction of motion of first eigenvector or the corresponding direction of the second feature vector as lesion tissue;
Radial position-angle position-axial direction for each sub-period, after the elimination noise according to the sub-period
Station location marker combination it is corresponding decay it is revised bury in oblivion a number and the first eigenvector or the second feature to
Amount determines the motion amplitude information of the lesion tissue described in the sub-period in the direction of movement.
Optionally, second determining module 430, is used for:
To the motion amplitude of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information that information is constituted-function of time carries out interpolation processing, motion amplitude information-time after obtaining interpolation processing
Function;
Using motion amplitude information-function of time after the interpolation processing as gate-control signal.
Optionally, as shown in fig. 6, second determining module 430, comprising:
Frequency-domain transform submodule 4301 is used for each sub-period and the lesion tissue described in each sub-period in institute
It states motion amplitude information-function of time that the motion amplitude information in the direction of motion is constituted and carries out frequency-domain transform processing, obtain institute
State the corresponding frequency spectrum of motion amplitude information-function of time;
Submodule 4302 is filtered, for based on the bandpass filter with default center frequency and pre-set bandwidths, to described
Frequency spectrum is filtered, and obtains filtered frequency spectrum;
Frequency domain inverse transformation submodule 4303, for carrying out the inverse transformation of the frequency-domain transform to the filtered frequency spectrum,
Obtain gate-control signal.
In the embodiment of the present invention, the dot position information of respectively burying in oblivion that can be detected according to PET device determines gate-control signal, nothing
Additional external equipment acquisition gate-control signal need to be increased when PET device detects and respectively buries in oblivion dot position information, thus, it can
To reduce diagnostic system cost.
It should be understood that the device of determining gate-control signal provided by the above embodiment is when determining gate-control signal, only with
The division progress of above-mentioned each functional module can according to need and for example, in practical application by above-mentioned function distribution by not
Same functional module is completed, i.e., the internal structure of terminal is divided into different functional modules, to complete whole described above
Or partial function.In addition, the device of determining gate-control signal provided by the above embodiment and the method for determining gate-control signal are implemented
Example belongs to same design, and specific implementation process is detailed in embodiment of the method, and which is not described herein again.
Those of ordinary skill in the art will appreciate that realizing that all or part of the steps of above-described embodiment can pass through hardware
It completes, relevant hardware can also be instructed to complete by program, the program can store in a kind of computer-readable
In storage medium, storage medium mentioned above can be read-only memory, disk or CD etc..
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (12)
1. a kind of method of determining gate-control signal, which is characterized in that the described method includes:
Step 1, it obtains and is examined in the sub-period of each preset duration of the positron emission tomography PET device when detecting in section
What is measured respectively buries in oblivion dot position information;
Step 2, respectively bury in oblivion dot position information according to each sub-period is corresponding, determine the direction of motion of lesion tissue, really
It is scheduled on the motion amplitude information of lesion tissue described in each sub-period in the direction of movement;
Step 3, the motion amplitude information according to the lesion tissue described in each sub-period in the direction of movement determines gate
Signal carries out PET image reconstruction based on the gate-control signal;
Wherein, the step 2 specifically includes:
For each sub-period, according to the corresponding radial position respectively buried in oblivion in dot position information of the sub-period, angle position
It is identified with axial position, determines that each radial position-angle position-axial position identifier combination is corresponding and bury in oblivion a number,
In, the axial position mark is intended to indicate that the information of the axial position for the vanishing point that falls into oblivion;
For each sub-period, each axial position is identified, by each radial position-angle position described in the sub-period
Set-axial position identifier combination in meet multiple radial position-angle positions of preset proximity condition and merge into one
New radial position-angle position combination, and angle position combination is corresponding buries in oblivion the sum of number by multiple radial positions-,
Bury in oblivion a number as the new radial position-angle position combination is corresponding, obtains radial position new in the sub-period
Set that-angle position-axial position identifier combination is corresponding to bury in oblivion a number;
According to preset principal component analysis PCA algorithm, to by radial position-angle position-axial direction position new in each sub-period
The corresponding matrix progress principal component analysis burying in oblivion a number and constituting of identifier combination is set, at least one characteristic value, Yi Jisuo are obtained
State the corresponding feature vector of at least one characteristic value, wherein the bigger expression of characteristic value buries in oblivion a number in the characteristic value
The corresponding direction of feature vector on variation degree it is bigger;
The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, by described first
The direction of motion of feature vector or the corresponding direction of the second feature vector as lesion tissue;
It is corresponding according to radial position-angle position new in the sub-period-axial position identifier combination for each sub-period
Bury in oblivion a number and the first eigenvector or the second feature vector, determine the lesion described in the sub-period
The motion amplitude information of tissue in the direction of movement.
2. the method according to claim 1, wherein the step 2 further comprises:
Data recombination is carried out to radial position-angle position new described in each sub-period-axial position identifier combination
Processing, radial position-angle position-axial position identifier combination after obtaining data recombination, according in each sub-period
New radial position-angle position-axial position the identifier combination is corresponding to bury in oblivion a number, determines each sub-period
Described in radial position-angle position-axial position identifier combination after data recombination is corresponding buries in oblivion a number;
According to radial position-angle position-axial position identifier combination pair after data recombination described in each sub-period
That answers buries in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the movement
Motion amplitude information on direction.
3. according to the method described in claim 2, it is characterized in that, the data recombination according to each sub-period
Radial position-angle position-axial position identifier combination afterwards is corresponding to bury in oblivion a number, determines the movement side of lesion tissue
To the motion amplitude information of the determining lesion tissue described in each sub-period in the direction of movement, comprising:
Radial position-angle position-axial direction after the data recombination described in formula and each sub-period that decayed according to nucleic
Station location marker combination is corresponding to bury in oblivion a number, radial position-angle after determining data recombination described in each sub-period
The corresponding decay of degree position-axial position identifier combination is revised to bury in oblivion a number;
According to radial position-angle position-axial position identifier combination pair after data recombination described in each sub-period
The decay answered is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion group described in each sub-period
Knit motion amplitude information in the direction of movement.
4. according to the method described in claim 3, it is characterized in that, the data recombination according to each sub-period
The corresponding decay of radial position-angle position-axial position identifier combination afterwards is revised to bury in oblivion a number, determines lesion group
The direction of motion knitted determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement, comprising:
For radial position-angle position after each data recombination in each sub-period-axial position mark group
Radial position-angle position combination in conjunction, according to the diameter after data recombination belonging to the radial position-angle position combination
Bury in oblivion a number to the corresponding decay of position-angle position-axial position identifier combination is revised, determines the radial position
Set-combination corresponding decay in angle position revised buries in oblivion a number;
To decay the revised radial position-angle position buried in oblivion a number and be less than preset number threshold value in each sub-period
It combines, the revised vanishing point that falls into oblivion of radial position-angle position-axial position identifier combination decay after affiliated data recombination
Number is set as zero, obtains eliminating radial position-angle position-axial position mark group after noise in each sub-period
It closes corresponding decay and revised buries in oblivion a number;
Radial position-angle position-axial position identifier combination pair after the elimination noise according to each sub-period
The decay answered is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion group described in each sub-period
Knit motion amplitude information in the direction of movement.
5. the method according to claim 1, wherein the step 3 includes:
To the motion amplitude information of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information of composition-function of time carries out interpolation processing, motion amplitude information-time letter after obtaining interpolation processing
Number;
Using motion amplitude information-function of time after the interpolation processing as gate-control signal.
6. the method according to claim 1, wherein the step 3 includes:
To the motion amplitude information of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information of composition-function of time carries out frequency-domain transform processing, and it is corresponding to obtain the motion amplitude information-function of time
Frequency spectrum;
Based on the bandpass filter with default center frequency and pre-set bandwidths, the frequency spectrum is filtered, after obtaining filtering
Frequency spectrum;
The inverse transformation that the frequency-domain transform is carried out to the filtered frequency spectrum, obtains gate-control signal.
7. a kind of device of determining gate-control signal, which is characterized in that described device includes:
Module is obtained, for obtaining the son of each preset duration of the positron emission tomography PET device when detecting in section
What is detected in period respectively buries in oblivion dot position information;
First determining module determines lesion tissue for respectively burying in oblivion dot position information according to each sub-period is corresponding
The direction of motion determines the motion amplitude information of the lesion tissue described in each sub-period in the direction of movement;
Second determining module, for the motion amplitude letter according to the lesion tissue described in each sub-period in the direction of movement
Breath, determines gate-control signal, carries out PET image reconstruction based on the gate-control signal;
Wherein, first determining module further comprises:
First determines submodule, for respectively burying in oblivion in dot position information according to the sub-period is corresponding for each sub-period
Radial position, angle position and axial position mark, determine each radial position-angle position-axial position identifier combination
It is corresponding to bury in oblivion a number, wherein the axial position mark is intended to indicate that the information of the axial position for the vanishing point that falls into oblivion;
Second determines submodule, will be described in the sub-period for being identified for each axial position for each sub-period
Meet multiple radial position-angles of preset proximity condition in each radial position-angle position-axial position identifier combination
Position grouping merges into a new radial position-angle position combination, and multiple radial positions-angle position combination is corresponded to
Bury in oblivion the sum of number, bury in oblivion a number as the new radial position-angle position combination is corresponding, obtain the son
New radial position-angle position-axial position identifier combination is corresponding in period buries in oblivion a number;
According to preset principal component analysis PCA algorithm, to by radial position-angle position-axial direction position new in each sub-period
The corresponding matrix progress principal component analysis burying in oblivion a number and constituting of identifier combination is set, at least one characteristic value, Yi Jisuo are obtained
State the corresponding feature vector of at least one characteristic value, wherein the bigger expression of characteristic value buries in oblivion a number in the characteristic value
The corresponding direction of feature vector on variation degree it is bigger;
The corresponding first eigenvector of maximum eigenvalue or the corresponding second feature vector of secondary big characteristic value are chosen, by described first
The direction of motion of feature vector or the corresponding direction of the second feature vector as lesion tissue;
It is corresponding according to radial position-angle position new in the sub-period-axial position identifier combination for each sub-period
Bury in oblivion a number and the first eigenvector or the second feature vector, determine the lesion described in the sub-period
The motion amplitude information of tissue in the direction of movement.
8. device according to claim 7, which is characterized in that described second determines submodule, is also used to:
Data recombination is carried out to radial position-angle position new described in each sub-period-axial position identifier combination
Processing, radial position-angle position-axial position identifier combination after obtaining data recombination, according in each sub-period
New radial position-angle position-axial position the identifier combination is corresponding to bury in oblivion a number, determines each sub-period
Described in radial position-angle position-axial position identifier combination after data recombination is corresponding buries in oblivion a number;
According to radial position-angle position-axial position identifier combination pair after data recombination described in each sub-period
That answers buries in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion tissue described in each sub-period in the movement
Motion amplitude information on direction.
9. device according to claim 8, which is characterized in that described second determines submodule, is also used to:
Radial position-angle position-axial direction after the data recombination described in formula and each sub-period that decayed according to nucleic
Station location marker combination is corresponding to bury in oblivion a number, radial position-angle after determining data recombination described in each sub-period
The corresponding decay of degree position-axial position identifier combination is revised to bury in oblivion a number;
According to radial position-angle position-axial position identifier combination pair after data recombination described in each sub-period
The decay answered is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion group described in each sub-period
Knit motion amplitude information in the direction of movement.
10. device according to claim 9, which is characterized in that described second determines submodule, is also used to:
For radial position-angle position after each data recombination in each sub-period-axial position mark group
Radial position-angle position combination in conjunction, according to the diameter after data recombination belonging to the radial position-angle position combination
Bury in oblivion a number to the corresponding decay of position-angle position-axial position identifier combination is revised, determines the radial position
Set-combination corresponding decay in angle position revised buries in oblivion a number;
To decay the revised radial position-angle position buried in oblivion a number and be less than preset number threshold value in each sub-period
It combines, the revised vanishing point that falls into oblivion of radial position-angle position-axial position identifier combination decay after affiliated data recombination
Number is set as zero, obtains eliminating radial position-angle position-axial position mark group after noise in each sub-period
It closes corresponding decay and revised buries in oblivion a number;
Radial position-angle position-axial position identifier combination pair after the elimination noise according to each sub-period
The decay answered is revised to bury in oblivion a number, determines the direction of motion of lesion tissue, determines the lesion group described in each sub-period
Knit motion amplitude information in the direction of movement.
11. device according to claim 7, which is characterized in that second determining module is used for:
To the motion amplitude information of each sub-period and the lesion tissue described in each sub-period in the direction of movement
The motion amplitude information of composition-function of time carries out interpolation processing, motion amplitude information-time letter after obtaining interpolation processing
Number;
Using motion amplitude information-function of time after the interpolation processing as gate-control signal.
12. device according to claim 7, which is characterized in that second determining module, comprising:
Frequency-domain transform submodule is used for each sub-period and the lesion tissue described in each sub-period in the movement side
The motion amplitude information that upward motion amplitude information is constituted-function of time carries out frequency-domain transform processing, obtains the movement width
Spend the corresponding frequency spectrum of information-function of time;
Submodule is filtered, for being carried out to the frequency spectrum based on the bandpass filter with default center frequency and pre-set bandwidths
Filtering, obtains filtered frequency spectrum;
Frequency domain inverse transformation submodule is gated for carrying out the inverse transformation of the frequency-domain transform to the filtered frequency spectrum
Signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374708.2A CN106097384B (en) | 2016-05-30 | 2016-05-30 | A kind of method and apparatus of determining gate-control signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610374708.2A CN106097384B (en) | 2016-05-30 | 2016-05-30 | A kind of method and apparatus of determining gate-control signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106097384A CN106097384A (en) | 2016-11-09 |
CN106097384B true CN106097384B (en) | 2018-12-18 |
Family
ID=57229541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610374708.2A Active CN106097384B (en) | 2016-05-30 | 2016-05-30 | A kind of method and apparatus of determining gate-control signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106097384B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106073812B (en) * | 2016-05-30 | 2018-09-18 | 江苏赛诺格兰医疗科技有限公司 | A kind of method and apparatus of determining gate-control signal |
CN106236123B (en) * | 2016-09-05 | 2019-02-19 | 赛诺联合医疗科技(北京)有限公司 | A kind of method and apparatus of determining gate-control signal |
US10307116B2 (en) * | 2017-02-22 | 2019-06-04 | Uih America, Inc. | System and method for detecting organ motion |
CN108634974B (en) * | 2018-04-03 | 2022-03-04 | 东软医疗系统股份有限公司 | Gate control signal determining method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101454801A (en) * | 2006-02-28 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | Local motion compensation based on list mode data |
CN101702232A (en) * | 2009-10-16 | 2010-05-05 | 昆明理工大学 | Respiration correction technique in positron emission tomography |
CN104352245A (en) * | 2014-11-24 | 2015-02-18 | 沈阳东软医疗系统有限公司 | Method and device for acquiring breathing cycle |
CN105411617A (en) * | 2015-12-11 | 2016-03-23 | 沈阳东软医疗系统有限公司 | Respiratory gating equipment and method as well as MCU (Microcontroller Unit) |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8351671B2 (en) * | 2007-07-26 | 2013-01-08 | Koninklijke Philips Electronics N.V. | Motion correction in nuclear imaging |
-
2016
- 2016-05-30 CN CN201610374708.2A patent/CN106097384B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101454801A (en) * | 2006-02-28 | 2009-06-10 | 皇家飞利浦电子股份有限公司 | Local motion compensation based on list mode data |
CN101702232A (en) * | 2009-10-16 | 2010-05-05 | 昆明理工大学 | Respiration correction technique in positron emission tomography |
CN104352245A (en) * | 2014-11-24 | 2015-02-18 | 沈阳东软医疗系统有限公司 | Method and device for acquiring breathing cycle |
CN105411617A (en) * | 2015-12-11 | 2016-03-23 | 沈阳东软医疗系统有限公司 | Respiratory gating equipment and method as well as MCU (Microcontroller Unit) |
Also Published As
Publication number | Publication date |
---|---|
CN106097384A (en) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106097384B (en) | A kind of method and apparatus of determining gate-control signal | |
CN101765865B (en) | Motion correction in nuclear imaging | |
CN109961834B (en) | Image diagnosis report generation method and device | |
CN106073812B (en) | A kind of method and apparatus of determining gate-control signal | |
US9424644B2 (en) | Methods and systems for evaluating bone lesions | |
US10289953B2 (en) | Dose optimization based on outcome quality | |
JP2021013726A (en) | Medical apparatus | |
CN102067176B (en) | Radiological imaging incorporating local motion monitoring, correction, and assessment | |
CN106236123B (en) | A kind of method and apparatus of determining gate-control signal | |
CN106232009B (en) | Tomographic apparatus and the method that tomoscan image is reconstructed by tomographic apparatus | |
CN107133996A (en) | Produce the method and PET/CT systems for the decay pattern rebuild for PET data | |
DE102020118559A1 (en) | Deep learning for perfusion in medical imaging | |
CN103534730A (en) | List mode dynamic image reconstruction | |
CN106163405B (en) | Tomographic apparatus and method of displaying tomographic image by tomographic apparatus | |
CN107635469A (en) | The estimation of the decay pattern met based on the scattering in PET system | |
CN104114091B (en) | Free-air correction core image reconstruction | |
CN106999135A (en) | The outer activity estimations of the FOV of use plain film and previous patient data in Positron emission tomography | |
CN108292443A (en) | The PET image reconstruction acted on behalf of using lesion and processing | |
CN108271346A (en) | Dose of radiation reduces and mixes the improved consistency between the period in imaging research | |
CN110536640A (en) | From the extract real-time of the noise robustness of the respiratory movement signal in PET table data | |
CN106102581A (en) | The segmentation of the structure of movement in view data | |
EP3408771B1 (en) | Predictive model for optimizing clinical workflow | |
CN105103194B (en) | Reconstructed image data visualization | |
CN104541183B (en) | The virtual framework rebuild for the distributed column table schema flight time moved using continuous bed | |
US10929974B2 (en) | Determining calcium content from spectral CT data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |