CN106097384A - A kind of method and apparatus determining gate-control signal - Google Patents

Abstract

The invention discloses a kind of method and apparatus determining gate-control signal, belong to medical domain.Described method includes: detect in the sub-period of the acquisition positron emission tomography PET device each preset duration within the detection period respectively buries in oblivion dot position information；According to described each sub-period corresponding respectively bury in oblivion dot position information, determine the direction of motion of lesion tissue, determine in lesion tissue motion amplitude information in the described direction of motion described in each sub-period；According in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, determine gate-control signal, carry out PET image reconstruction based on described gate-control signal.Use the present invention, diagnostic system cost can be reduced.

Description

A kind of method and apparatus determining gate-control signal

Technical field

The present invention relates to medical domain, particularly to a kind of method and apparatus determining gate-control signal.

Background technology

In medical domain, PET/CT (positron emission tomography/computed can be passed through Tomography, positron emission tomography/computed tomography) equipment is to human body or animal body tumor, cardiac system Disease and nervous system disease carry out early diagnosis, the image that i.e. can be obtained scanned object scan by PET/CT equipment Determine the position that focus occurs.But, when PET imaging, owing to imaging process has certain duration, exhaling of scanned object Inhale motion or heartbeat makes the tissue of focus (such as lung) spatially periodic shift, cause the PET of gained to scheme As fuzzy, making it difficult to maybe can not identify the little pathological changes in PET image.

The PET image caused for eliminating respiratory movement obscures, and often uses gating technology in PET, i.e. by specific Equipment gathers gate-control signal corresponding to respiratory movement, and (wherein, gate-control signal could be for reflecting that scanned object is when scanning Respiratory movement process in long (such as 12 minutes) or the signal of heartbeat process), at present, conventional collection gate-control signal Method generally requires dependence external equipment, such as, is scanned object body surface by infrared collection device or camera acquisition Motor message, wherein, motor message can be that scanned object is under lying status, and the height of body surface a part is over time The height value being changed, this motor message collected is gate-control signal.It is then possible to according to the gate-control signal obtained, Carry out PET image reconstruction, the breathing cycle that i.e. can reflect, the same motion that will repeat many times according to gate-control signal In PET data be grouped together, the PET data being grouped together is redeveloped into the image of a motion stage.

During realizing the present invention, inventor finds that prior art at least there is problems in that

When gathering gate-control signal based on above-mentioned processing mode, need to increase extra external equipment, and the valency of external equipment Lattice the most costly, thus, cause diagnostic system cost to increase.

Summary of the invention

In order to solve problem of the prior art, embodiments provide a kind of method determining gate-control signal and dress Put.Described technical scheme is as follows:

First aspect, it is provided that a kind of method determining gate-control signal, described method includes:

Obtain in the sub-period of the positron emission tomography PET device each preset duration within the detection period and detect To respectively bury in oblivion dot position information；

According to described each sub-period corresponding respectively bury in oblivion dot position information, determine the direction of motion of lesion tissue, determine In lesion tissue motion amplitude information in the described direction of motion described in each sub-period；

According in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, determine gate letter Number, carry out PET image reconstruction based on described gate-control signal.

Optionally, described according to described each sub-period corresponding respectively bury in oblivion dot position information, determine the fortune of lesion tissue Dynamic direction, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

For each sub-period, according to described sub-period corresponding respectively bury in oblivion dot position information, determine in described sub-period Difference is buried in oblivion corresponding the burying in oblivion of dot position information and is counted out；

Bury in oblivion corresponding the burying in oblivion of dot position information according to difference in described each sub-period to count out, determine lesion tissue The direction of motion, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period.

Optionally, described for each sub-period, according to described sub-period corresponding respectively bury in oblivion dot position information, determine institute State in sub-period and different to bury in oblivion corresponding the burying in oblivion of dot position information and count out, including:

For each sub-period, according to the radial position respectively buried in oblivion in dot position information corresponding to described sub-period, angle Position and axial location mark, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted Mesh, wherein, described axial location mark is intended to indicate that the information burying in oblivion axial location a little；

Described bury in oblivion corresponding the burying in oblivion of dot position information count out according to different in described each sub-period, determine focus group The direction of motion knitted, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

For each sub-period, by radial position-angle position-axial location mark group each described in described sub-period A new radial position-angle is merged in the multiple radial position-angle positions meeting preset proximity condition in conjunction Position grouping, and corresponding for multiple radial positions-angle position combination is buried in oblivion sum of counting out, as described new radial direction position Put-angle position corresponding the burying in oblivion of combination count out, and obtains position, radial position-angle position new in described sub-period-axially Put corresponding the burying in oblivion of identifier combination to count out；

Corresponding according to radial position-angle position-axial location identifier combination new described in described each sub-period Bury in oblivion and count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period in the described direction of motion On motion amplitude information.

Optionally, described according to radial position-angle position new described in described each sub-period-axial location mark Corresponding the burying in oblivion of combination is counted out, and determines the direction of motion of lesion tissue, determines in lesion tissue described in each sub-period in institute State the motion amplitude information in the direction of motion, including:

Radial position-angle position-axial location identifier combination new described in described each sub-period is carried out data Restructuring processes, and obtains the radial position-angle position-axial location identifier combination after data recombination, according to described each period of the day from 11 p.m. to 1 a.m Radial position-angle position new described in Duan-corresponding burying in oblivion of axial location identifier combination is counted out, and determines described every height Radial position-angle position-corresponding burying in oblivion of axial location identifier combination after data recombination described in period is counted out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding burying in oblivion to count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period described Motion amplitude information in the direction of motion.

Optionally, described according to the radial position-angle position after data recombination described in described each sub-period-axially Station location marker corresponding the burying in oblivion of combination is counted out, and determines the direction of motion of lesion tissue, determines in focus described in each sub-period It is organized in the motion amplitude information in the described direction of motion, including:

According to radial position-angle position after data recombination described in nucleic decay formula and described each sub-period- Corresponding the burying in oblivion of axial location identifier combination is counted out, and determines the radial direction position after data recombination described in described each sub-period Put-revised the burying in oblivion of decay corresponding to angle position-axial location identifier combination count out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described according to the radial position-angle position after data recombination described in described each sub-period-axially Revised the burying in oblivion of decay of station location marker combination correspondence is counted out, and determines the direction of motion of lesion tissue, determines at every height The motion amplitude information in the described direction of motion of the lesion tissue described in period, including:

For the radial position-angle position-axial location mark after described data recombination each in described each sub-period Know the radial position in combination-angle position combination, after the data recombination belonging to described radial position-angle position combination Revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination count out, determine described radial direction Revised the burying in oblivion of decay of position-angle position combination correspondence is counted out；

By the revised radial position-angle buried in oblivion and count out not less than preset number threshold value that decays in each sub-period Position grouping, the decay of the radial position-angle position-axial location identifier combination after affiliated data recombination is revised falls into oblivion Vanishing point number is set to zero, obtains the radial position-angle position-axial location mark after eliminating noise in described each sub-period Revised the burying in oblivion of decay knowing combination corresponding is counted out；

According to the radial position-angle position-axial location mark group after elimination noise described in described each sub-period Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described according to eliminating the radial position-angle position after noise-axially described in described each sub-period Revised the burying in oblivion of decay of station location marker combination correspondence is counted out, and determines the direction of motion of lesion tissue, determines at every height The motion amplitude information in the described direction of motion of the lesion tissue described in period, including:

According to default principal component analysis PCA algorithm, to by the radial direction after elimination noise described in described each sub-period Position-angle position-axial location identifier combination corresponding decay revised bury in oblivion count out constitute matrix carry out main one-tenth Analyze, obtain at least one eigenvalue, and at least one eigenvalue characteristic of correspondence vector described, wherein, described feature It is the biggest that intensity of variation on the direction that the characteristic vector of described eigenvalue is corresponding of counting out is buried in oblivion in the biggest expression of value；

Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by described First eigenvector or direction corresponding to described second feature vector are as the direction of motion of lesion tissue；

For each sub-period, according to eliminating the radial position-angle position after noise-axially described in described sub-period Revised the burying in oblivion of decay corresponding to station location marker combination is counted out, and described first eigenvector or described second feature to Amount, determines in lesion tissue motion amplitude information in the described direction of motion described in described sub-period.

Optionally, described basis lesion tissue described in each sub-period in the described direction of motion motion amplitude letter Breath, determines gate-control signal, including:

To described each sub-period with at lesion tissue motion amplitude in the described direction of motion described in each sub-period Motion amplitude information-the time function of information structure carries out interpolation processing, obtains the motion amplitude information-time after interpolation processing Function；

Using the motion amplitude information-time function after described interpolation processing as gate-control signal.

Optionally, described basis lesion tissue described in each sub-period in the described direction of motion motion amplitude letter Breath, determines gate-control signal, including:

To described each sub-period with at lesion tissue motion amplitude in the described direction of motion described in each sub-period Motion amplitude information-the time function of information structure carries out frequency domain transform process, obtains described motion amplitude information-time function Corresponding frequency spectrum；

Based on having the band filter of default center frequency and pre-set bandwidths, described frequency spectrum is filtered, is filtered Frequency spectrum after ripple；

Described filtered frequency spectrum is carried out the inverse transformation of described frequency domain transform, obtains gate-control signal.

Second aspect, it is provided that a kind of device determining gate-control signal, described device includes:

Acquisition module, for obtaining the positron emission tomography PET device each preset duration within the detection period Sub-period in detect respectively bury in oblivion dot position information；

First determines module, for according to described each sub-period corresponding respectively bury in oblivion dot position information, determine focus group The direction of motion knitted, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period；

Second determines module, for according at the motion width in the described direction of motion of the lesion tissue described in each sub-period Degree information, determines gate-control signal, carries out PET image reconstruction based on described gate-control signal.

Optionally, described first determines module, including:

First determines submodule, for for each sub-period, believing according to the position of respectively burying in oblivion that described sub-period is corresponding Breath, determines in described sub-period and different to bury in oblivion corresponding the burying in oblivion of dot position information and count out；

Second determines submodule, for burying in oblivion corresponding the burying in oblivion a little of dot position information according to different in described each sub-period Number, determines the direction of motion of lesion tissue, determines in lesion tissue fortune in the described direction of motion described in each sub-period Dynamic amplitude information.

Optionally, described first determines submodule, is used for:

For each sub-period, according to the radial position respectively buried in oblivion in dot position information corresponding to described sub-period, angle Position and axial location mark, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted Mesh, wherein, described axial location mark is intended to indicate that the information burying in oblivion axial location a little；

Described second determines submodule, is used for:

For each sub-period, by radial position-angle position-axial location mark group each described in described sub-period A new radial position-angle is merged in the multiple radial position-angle positions meeting preset proximity condition in conjunction Position grouping, and corresponding for multiple radial positions-angle position combination is buried in oblivion sum of counting out, as described new radial direction position Put-angle position corresponding the burying in oblivion of combination count out, and obtains position, radial position-angle position new in described sub-period-axially Put corresponding the burying in oblivion of identifier combination to count out；

Corresponding according to radial position-angle position-axial location identifier combination new described in described each sub-period Bury in oblivion and count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period in the described direction of motion On motion amplitude information.

Optionally, described second determines submodule, is used for:

Radial position-angle position-axial location identifier combination new described in described each sub-period is carried out data Restructuring processes, and obtains the radial position-angle position-axial location identifier combination after data recombination, according to described each period of the day from 11 p.m. to 1 a.m Radial position-angle position new described in Duan-corresponding burying in oblivion of axial location identifier combination is counted out, and determines described every height Radial position-angle position-corresponding burying in oblivion of axial location identifier combination after data recombination described in period is counted out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding burying in oblivion to count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period described Motion amplitude information in the direction of motion.

Optionally, described second determines submodule, is used for:

According to radial position-angle position after data recombination described in nucleic decay formula and described each sub-period- Corresponding the burying in oblivion of axial location identifier combination is counted out, and determines the radial direction position after data recombination described in described each sub-period Put-revised the burying in oblivion of decay corresponding to angle position-axial location identifier combination count out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described second determines submodule, is used for:

For the radial position-angle position-axial location mark after described data recombination each in described each sub-period Know the radial position in combination-angle position combination, after the data recombination belonging to described radial position-angle position combination Revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination count out, determine described radial direction Revised the burying in oblivion of decay of position-angle position combination correspondence is counted out；

By the revised radial position-angle buried in oblivion and count out not less than preset number threshold value that decays in each sub-period Position grouping, the decay of the radial position-angle position-axial location identifier combination after affiliated data recombination is revised falls into oblivion Vanishing point number is set to zero, obtains the radial position-angle position-axial location mark after eliminating noise in described each sub-period Revised the burying in oblivion of decay knowing combination corresponding is counted out；

According to the radial position-angle position-axial location mark group after elimination noise described in described each sub-period Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described second determines submodule, is used for:

According to default principal component analysis PCA algorithm, to by the radial direction after elimination noise described in described each sub-period Position-angle position-axial location identifier combination corresponding decay revised bury in oblivion count out constitute matrix carry out main one-tenth Analyze, obtain at least one eigenvalue, and at least one eigenvalue characteristic of correspondence vector described, wherein, described feature It is the biggest that intensity of variation on the direction that the characteristic vector of described eigenvalue is corresponding of counting out is buried in oblivion in the biggest expression of value；

Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by described First eigenvector or direction corresponding to described second feature vector are as the direction of motion of lesion tissue；

For each sub-period, according to eliminating the radial position-angle position after noise-axially described in described sub-period Revised the burying in oblivion of decay corresponding to station location marker combination is counted out, and described first eigenvector or described second feature to Amount, determines in lesion tissue motion amplitude information in the described direction of motion described in described sub-period.

Optionally, described second determines module, is used for:

To described each sub-period with at lesion tissue motion amplitude in the described direction of motion described in each sub-period Motion amplitude information-the time function of information structure carries out interpolation processing, obtains the motion amplitude information-time after interpolation processing Function；

Using the motion amplitude information-time function after described interpolation processing as gate-control signal.

Optionally, described second determines module, including:

Frequency domain transform submodule, for described each sub-period and in lesion tissue described in each sub-period in described fortune Motion amplitude information-time function that motion amplitude information on dynamic direction is constituted carries out frequency domain transform process, obtains described fortune The frequency spectrum that dynamic amplitude information-time function is corresponding；

Filtering submodule, is used for based on the band filter with default center frequency and pre-set bandwidths, to described frequency spectrum It is filtered, obtains filtered frequency spectrum；

Frequency domain inverse transformation submodule, for described filtered frequency spectrum carries out the inverse transformation of described frequency domain transform, obtains Gate-control signal.

The technical scheme that the embodiment of the present invention provides has the benefit that

In the embodiment of the present invention, gate-control signal, nothing can be determined according to the dot position information of respectively burying in oblivion that PET device detects Extra external equipment need to be increased and gathers gate-control signal PET device detection respectively buries in oblivion dot position information when, thus, can To reduce diagnostic system cost.

Accompanying drawing explanation

For the technical scheme being illustrated more clearly that in the embodiment of the present invention, in embodiment being described below required for make Accompanying drawing be briefly described, it should be apparent that, below describe in accompanying drawing be only some embodiments of the present invention, for From the point of view of those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to obtain other according to these accompanying drawings Accompanying drawing.

Fig. 1 is a kind of method flow diagram determining gate-control signal that the embodiment of the present invention provides；

Fig. 2 is the structural representation of a kind of terminal that the embodiment of the present invention provides；

Fig. 3 is the schematic diagram of a kind of crystal detection pair that the embodiment of the present invention provides；

Fig. 4 is a kind of apparatus structure schematic diagram determining gate-control signal that the embodiment of the present invention provides；

Fig. 5 is a kind of apparatus structure schematic diagram determining gate-control signal that the embodiment of the present invention provides；

Fig. 6 is a kind of apparatus structure schematic diagram determining gate-control signal that the embodiment of the present invention provides.

Detailed description of the invention

For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing to embodiment party of the present invention Formula is described in further detail.

Embodiments provide a kind of method determining gate-control signal, as it is shown in figure 1, the executive agent of the method is Terminal.Wherein, this terminal can be the terminal having and determining gate-control signal function, such as, can be and PET (positron Emission tomography, positron emission tomography) equipment is connected, can determine according to the data of PET device collection The computer equipment of gate-control signal.

Terminal can include processor 210, transceiver 220, memorizer 230, input block 240, display unit 250, sound The parts such as frequency circuit 260 and power supply 270, as in figure 2 it is shown, it will be understood by those skilled in the art that the terminal knot shown in Fig. 2 Structure is not intended that the restriction to terminal, can include that ratio illustrates more or less of parts or combines some parts, or not Same parts are arranged.Wherein:

Processor 210 can be the control centre of terminal, utilizes each of various interface and the whole terminal unit of connection Individual part, such as transceiver 220 and memorizer 230 etc., by run or perform the software program that is stored in memorizer 230 and/ Or module, and call the data being stored in memorizer 230, perform the various functions of terminal and process data, thus to end End carries out integral monitoring.Optionally, processor 210 can include one or more process core.In the present invention, processor 210 It is determined for the relevant treatment of gate-control signal.Transceiver 220 may be used for receiving and sending data, and terminal can be by receiving Sending out device 220 and receive and send data, terminal can pass through the Internet transceiving data, and transceiver can be network interface card.

Memorizer 230 can be used for storing software program and module, and processor 210 is stored in memorizer 230 by operation Software program and module, thus perform various function application and data process.Memorizer 230 can mainly include storing journey Sequence district and storage data field, wherein, storage program area can store the application program (ratio needed for operating system, at least one function As determined gate-control signal function etc.) etc.；Storage data field can store the data that the use according to terminal created and (such as bury in oblivion a little Positional information etc.) etc..Additionally, memorizer 230 can include high-speed random access memory, it is also possible to include non-volatile memories Device, for example, at least one disk memory, flush memory device or other volatile solid-state parts.Input block 240 is permissible For receiving numeral or the character information of input, and produce the keyboard relevant with user setup and function control, mouse, behaviour Make bar, optics or the input of trace ball signal.Display unit 250 is displayed for the information that inputted 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, video Constitute with its combination in any.Display unit 250 can include display floater 251, optionally, can use LCD (Liquid Crystal Display, liquid crystal display), OLED (Organic Light-Emitting Diode, Organic Light Emitting Diode) Display floater 251 is configured etc. form.Voicefrequency circuit 260, speaker 261, microphone 262 can provide between user and terminal Audio interface, the voice data received can be converted to the signal of telecommunication by voicefrequency circuit 260.Power supply 270 can pass through power management System is logically contiguous with processor 210, thus realizes 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 alternating current power supply, recharging system, power failure monitor Circuit, power supply changeover device or the random component such as inverter, power supply status indicator.Terminal also includes memorizer, Yi Jiyi Individual or more than one program, one of them or more than one program is stored in memorizer, and is configured to by one Or more than one processor performs to state one or more than one program to perform the determination door described in each embodiment following The method of control signal.

Below in conjunction with detailed description of the invention, being described in detail the handling process shown in Fig. 1, content can be as Under:

Step 101, obtains the period of the day from 11 p.m. to 1 a.m of the positron emission tomography PET device each preset duration within the detection period Detect in Duan respectively buries in oblivion dot position information.

Wherein, burying in oblivion can be a little that positron buries in oblivion the location point producing a pair photon with internal negatron, falls into oblivion Vanishing point positional information could be for representing the information burying in oblivion a position, such as, can be almost to be simultaneously received a pair photon The numbering of two crystal detections and axial location mark, wherein, axial location mark could be for representing the axle buried in oblivion a little To the information of position, axial location mark can represent by the numbering of detector rings each in PET device, such as, and each detection Device ring includes that 60 crystal detections, numbering are followed successively by 1,2 ... 60, PET device includes that 15 detector rings, numbering are successively altogether 1,2 ... 15, certain buries in oblivion dot position information can be (40,50,2), i.e. represents and buries in oblivion a pair photon produced, quilt at this No. 40 of 2nd detector rings and No. 50 crystal detections are almost detected simultaneously by.

Bury in oblivion dot position information and can also is that radial position, angle position and axial location mark, i.e. by radial position, angle Degree position and axial location three parameter group cooperations of mark are one and bury in oblivion dot position information, wherein, by detector rings numbering and What axial location mark formed buries in oblivion dot position information and by burying in oblivion that radial position, angle position and axial location mark form Dot position information has one-to-one relationship, can turn mutually between detector rings numbering and radial position-angle position combination Changing, angle position can be the angle between line and vertical curve between the crystal detection of two detector rings numbering correspondences, Value is between 0-180, and radial position can be the distance of the line between the center of circle to two crystal detection of detector rings.

It should be noted is that, the axial location mark that what PET device detected bury in oblivion in dot position information can be used The numbering of detector rings represents.Specifically, it is positioned at same detection when the crystal detection being almost detected simultaneously by a pair photon During device ring, this buries in oblivion the numbering that a corresponding axial location mark is i.e. this detector rings.But, almost it is detected simultaneously by a pair The crystal detection of photon is also possible to be positioned at different detector rings, as it is shown on figure 3, now, this buries in oblivion a corresponding axial location Mark can combine by the numbering of two detector rings that crystal detection is positioned at and represent, it is also possible to visit by the two pre-set The mark surveying the combination of device ring corresponding represents, wherein, can be previously provided with the axle that the combination of each two detector rings is corresponding in this locality To station location marker, may be greater than the numeral of the detector rings number that PET device is comprised, such as, PET device has 15 spies Surveying device ring, numbering is 1-15 successively, can pre-set two crystal detections of 16 expression and be positioned at the 1st, 2 detector rings, 17 represents Two crystal detections are positioned at the 1st, 3 detector rings, 18 represent two crystal detections be positioned at the 1st, 4 detector rings, by that analogy, often Two detector rings combine the most in advance to there being a numbering.So, PET device detects when burying in oblivion, can be according to detection Bury in oblivion the detector rings at crystal detection place a little to this, bury in oblivion for this and distribute corresponding axial location mark.

In force, at medical domain, can by PET device to human body or animal body tumor, cardiac system disease and Nervous system disease carries out early diagnosis, i.e. can determine generation by the image that scanned object scan is obtained by PET device The position of focus.Concrete, first can to scanned target injection radionuclide (as¹⁸F、¹¹C), wherein, radioactive nucleus Element is positron radionuclide, and radionuclide (wherein, can occur the tissue of focus can absorb greatly along with metabolism enters internal tissue The radionuclide of amount), then, positron will be buried in oblivion with internal negatron, produces a pair heterodromous photon, This pair heterodromous photon will be detected by two crystal detections being positioned at symmetric position, wherein, and each spy of PET device Survey the crystal detection on device ring with same number, and each crystal detection in each detector rings have unique numbering, PET device often detects a pair smooth period of the day from 11 p.m. to 1 a.m, can be to representing that this information (i.e. burying in oblivion dot position information) burying in oblivion a position is carried out Storage, i.e. can deposit the numbering of two crystal detections and axial location mark that are almost detected simultaneously by a pair photon Storage, finally, terminal can bury in oblivion dot position information according to what PET device detected, forms PET image.And then, doctor can root Determining whether this sweep object has the tissue of generation focus according to PET image, wherein, PET device respectively buries in oblivion a position letter in detection During breath, need certain detection duration (general detection duration is about at about 12 minutes), at PET device, scanned object is swept Retouching period, scanned object can carry out repeatedly respiratory movement or heartbeat, in such cases, if scanned object has generation disease During the tissue of stove, due to the motion of lesion tissue, it will cause the PET image obtained to occur fuzzy.

At present, the problem that gating technology solves to cause the PET image obtained to obscure due to lesion tissue motion can be passed through, To specifically introduce terminal and determine the detailed process of gate-control signal below.

The whole detection Time segments division of PET device can be multiple sub-periods with identical preset duration by terminal, its In, preset duration can be 0.5 second, it is also possible to is 1 second.Terminal can obtain PET device and detect many within the detection period Individual burying in oblivion dot position information, wherein, what PET device detected respectively buries in oblivion dot position information can be according to respectively burying in oblivion a position letter The list type data that the detection time that breath is corresponding is arranged in order, and then, terminal can be corresponding based on respectively burying in oblivion dot position information The detection time, determine that detect in the sub-period of each preset duration within the detection period respectively buries in oblivion dot position information.

Step 102, according to each sub-period corresponding respectively bury in oblivion dot position information, determine the direction of motion of lesion tissue, Determine in each sub-period lesion tissue motion amplitude information in the movement direction.

Wherein, motion amplitude information can be the information of reflection lesion tissue motion amplitude in the movement direction.

In force, for each sub-period, terminal can according to each sub-period corresponding respectively bury in oblivion dot position information, Determine the direction of motion of lesion tissue, and according to each sub-period corresponding respectively bury in oblivion dot position information, determine lesion tissue The trend over time of motion amplitude information in the movement direction, i.e. determines the lesion tissue motion width at each sub-period Degree information.It is to say, when burying in oblivion dot position information and comprising crystal detection numbering and axial location mark, can be according to respectively falling into oblivion Crystal detection numbering and axial location that vanishing point is corresponding identify, and determine the direction of motion of lesion tissue, and determine at every height Period lesion tissue motion amplitude information in the movement direction.Radial position, angle position is comprised when burying in oblivion dot position information When identifying with axial location, disease can be determined according to respectively burying in oblivion a corresponding radial position, angle position and axial location mark The direction of motion of stove tissue, and determine in each sub-period lesion tissue motion amplitude information in the movement direction.

Optionally, (can bury in oblivion according to the different occurrence number burying in oblivion dot position information in each sub-period and count Mesh), determine the direction of motion of lesion tissue, and determine in lesion tissue described in each sub-period in the described direction of motion Motion amplitude information, accordingly, the processing procedure of step 102 can be such that for each sub-period, corresponding according to sub-period Respectively bury in oblivion dot position information, determine in sub-period and different to bury in oblivion corresponding the burying in oblivion of dot position information and count out；According to each period of the day from 11 p.m. to 1 a.m Different in Duan bury in oblivion corresponding the burying in oblivion of dot position information and count out, determine the direction of motion of lesion tissue, determine at each period of the day from 11 p.m. to 1 a.m Duan Suoshu lesion tissue motion amplitude information in the movement direction.

In force, terminal get each sub-period corresponding respectively bury in oblivion dot position information (wherein, each sub-period pair What respectively burying in oblivion of answering comprised repetition in dot position information buries in oblivion dot position information) after, different can bury in oblivion a position to each The occurrence number of information is added up, and obtains each different occurrence number burying in oblivion dot position information (in the most each sub-period Difference is buried in oblivion corresponding the burying in oblivion of dot position information and is counted out), i.e. can get the sinogram Sinogram number that each sub-period is corresponding According to, wherein, Sinogram data can include that each different occurrence number burying in oblivion dot position information and correspondence thereof (is buried in oblivion a little Number), Sinogram data can be used to reflect the distribution situation buried in oblivion a little.Such as, a position of burying in oblivion of the first sub-period is believed 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 occurrence number that this sub-period is buried in oblivion dot position information is added up, the information that the Sinogram packet obtained contains For: 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 to record in table form each to bury in oblivion dot position information and burying in oblivion of correspondence is counted out, and needs explanation Being that, in practical operation, each sub-period can include substantial amounts of dot position information of burying in oblivion, above-mentioned example is simply to illustrate that by arranging List-directed data is to the processing procedure of Sinogram data.After terminal obtains the Sinogram data of each sub-period, can basis The Sinogram data of each sub-period, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period Motion amplitude information in the movement direction.

Optionally, terminal can be according to the radial position respectively buried in oblivion in dot position information corresponding to each sub-period, angle Position and axial location mark, determine that difference is buried in oblivion corresponding the burying in oblivion of dot position information and counted out, and accordingly, processing procedure is permissible As follows: for each sub-period, according to the radial position respectively buried in oblivion in dot position information, angle position and axle that sub-period is corresponding To station location marker, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, wherein, axle The information burying in oblivion axial location a little it is intended to indicate that to station location marker.Accordingly, processing procedure can be such that for every height Period, many by radial position-angle position-axial location identifier combination each in sub-period meets preset proximity condition A new radial position-angle position combination is merged in individual radial position-angle position, and by multiple radial positions- What angle position combination was corresponding buries in oblivion sum of counting out, and counts as new radial position-angle position corresponding burying in oblivion of combination Mesh, obtains radial position-angle position-corresponding burying in oblivion of axial location identifier combination new in sub-period and counts out；According to often Radial position-angle position new in individual sub-period-corresponding burying in oblivion of axial location identifier combination is counted out, and determines lesion tissue The direction of motion, determine in each sub-period lesion tissue motion amplitude information in the described direction of motion.

In force, terminal obtain each sub-period corresponding bury in oblivion dot position information after, wherein, bury in oblivion dot position information Can be radial position, angle position and axial location mark, each radial position-angle position-axial location can be added up The occurrence number of identifier combination, i.e. can add up each radial position-angle position-corresponding burying in oblivion of axial location identifier combination Count out, namely can be according to the corresponding radial position respectively buried in oblivion in dot position information, angle position and the axle of each sub-period To station location marker, determining the Sinogram data that each sub-period is corresponding, wherein, the form of Sinogram data can be in advance Storage, after often determining that a radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, can be by This buries in oblivion the relevant position added in Sinogram data of counting out, and comprises in the Sinogram data of the most each sub-period Each radial position-angle position-axial location identifier combination is fixing, each radial position-angle position-axial location It is change that corresponding the burying in oblivion of identifier combination is counted out, and is that the dot position information of respectively burying in oblivion according to each sub-period is determined, Sinogram data can be the form of three-dimensional table.Such as, for the sub-Sinogram data that certain axial location mark is corresponding, As shown in table 1, wherein, according to from left to right, the order that numerical value increases successively from top to bottom arranges, each radial position and angle The number that position grouping is corresponding is this and buries in oblivion corresponding the burying in oblivion of dot position information and count out.Such as, axial location is identified 18, position (0,0) corresponding 5 i.e. represent radial position be 0, angle position be 0, axial location be designated 18 bury in oblivion a position It is 5 that corresponding the burying in oblivion of information is counted out.

Table 1

After terminal obtains the Sinogram data that each sub-period is corresponding, for the son that each axial location mark is corresponding Sinogram data, can make sub-Sinogram data according to default proximity condition, along radial position direction and angle Locality is compressed.Concrete, terminal can by the preset number of the most adjacent (i.e. meet and preset close to condition) radially A new radial position-angle position combination, wherein, new radial position-angle position are merged in position-angle position Put the average that the radial position in combination can be the radial position in multiple radial position-angle position combination, new radial direction Angle position in position-angle position combination can be the equal of the angle position in multiple radial position-angle position combination Value, and can burying in oblivion preset number radial position-angle position combination correspondence to count out and be added, obtain falls into oblivion Vanishing point number sum is counted out, for each axial location mark as new radial position-angle position corresponding burying in oblivion of combination Knowing, terminal all can carry out data compression according to above-mentioned data compression method along radial position direction and direction, angle position, Radial position-angle position new in sub-period-corresponding burying in oblivion of axial location identifier combination is counted out.Such as, such as table 1 institute 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 Individual new radial position-angle position combination-(0.5,2.5), (2,0), (3,0), (2,5), (3,5) are merged into one new Radial position-combination-(2.5,2.5), angle position, merge into a new footpath by (0,10), (1,10), (0,12), (1,12) To position-combination-(0.5,11), angle position, a new radial direction is merged in (2,10), (3,10), (2,12), (3,12) Position-combination-(2.5,11), angle position.As such, it is possible to make each radial direction position in the Sinogram data after data compression Put-angle position-corresponding burying in oblivion of axial location identifier combination count out increase, i.e. can improve each radial position-angle The data statistics amount of position-axial location identifier combination, and then, the interference of noise can be reduced.

After terminal carries out data compression to the Sinogram data of each sub-period, can be according to new in each sub-period Radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, and determines the direction of motion of lesion tissue, really It is scheduled on lesion tissue motion amplitude information in the described direction of motion described in each sub-period.Concrete, can be according to presetting Principal component analysis PCA algorithm, to by radial position-angle position-axial location identifier combination pair new in each sub-period The matrix constituted of counting out of burying in oblivion answered carries out principal component analysis, obtains at least one eigenvalue, and at least one eigenvalue Characteristic of correspondence vector, wherein, the biggest expression of eigenvalue is buried in oblivion to count out and is become on the direction that the characteristic vector of eigenvalue is corresponding Change degree is the biggest；Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, will First eigenvector or direction corresponding to second feature vector are as the direction of motion of lesion tissue；For each sub-period, root Count out and fisrt feature according to radial position-angle position new in sub-period-corresponding burying in oblivion of axial location identifier combination Vector or second feature vector, determine in sub-period lesion tissue motion amplitude information in the movement direction.

Optionally, terminal can also be to radial position-angle position-axial location identifier combination new in each sub-period Carry out data recombination process, accordingly, processing procedure can be such that to radial position-angle position new in each sub-period- Axial location identifier combination carries out data recombination process, obtains the radial position-angle position-axial location mark after data recombination Know combination, count according to radial position-angle position new in each sub-period-corresponding burying in oblivion of axial location identifier combination Mesh, determines radial position-angle position-corresponding the burying in oblivion a little of axial location identifier combination after data recombination in each sub-period Number；According to the radial position-angle position after data recombination in each sub-period-corresponding burying in oblivion of axial location identifier combination Count out, determine the direction of motion of lesion tissue, determine in the motion in the described direction of motion of each sub-period lesion tissue Amplitude information.

In force, for each sub-period, after terminal carries out data compression to the Sinogram data of sub-period, first Can be according to default data recombination algorithm (such as SSRB (Single Slice Rebin, monolayer is recombinated) algorithm, FORE (Fourier Rebin, Fourier recombinate) algorithm), represent almost be simultaneously received the spy of photon pair for burying in oblivion dot position information Survey device crystal to be positioned at the point of burying in oblivion of different detector rings (i.e. axial location mark buries in oblivion a position letter more than detector rings number Cease corresponding burying in oblivion a little, or axial location mark buries in oblivion dot position information correspondence with what two detector rings numbering combinations represented Bury in oblivion a little), determine the radial position-angle position-axial location identifier combination after data recombination, say, that Ke Yigen According to default data recombination algorithm, determine these bury in oblivion a corresponding data recombination after bury in oblivion dot position information, such as, PET Equipment includes 15 detector rings altogether, and certain buries in oblivion a corresponding dot position information of burying in oblivion is (1,10,17), wherein, pre-sets Axial location mark 17 expression two different detector rings be 1,3, can be determined this according to default data recombination algorithm Bury in oblivion a corresponding new axial location and be designated 2, i.e. through data recombination process after, this bury in oblivion a little bury in oblivion dot position information For (10,20,2), the axial location mark that i.e. will be greater than detector rings number is converted to the axial location mark less than detector number of rings Know.Then, terminal can be to each radial position-angle position-axial location identifier combination of data recombination with without data weight Radial position-angle position-corresponding the burying in oblivion of axial location identifier combination of group is counted out and is reintegrated, and determines every height In period, radial position-angle position-corresponding the burying in oblivion of axial location identifier combination after data recombination is counted out, and is i.e. counted According to the Sinogram data after restructuring.

After terminal obtains the Sinogram data after data recombination, can be according to the footpath after data recombination in each sub-period To position-angle position-corresponding the burying in oblivion of axial location identifier combination count out, determine the direction of motion of lesion tissue, determine In each sub-period lesion tissue motion amplitude information in the movement direction.Concrete, can divide according to default main constituent Analysis PCA algorithm, to corresponding by the radial position-angle position after data recombination in each sub-period-axial location identifier combination Bury in oblivion count out constitute matrix carry out principal component analysis, obtain at least one eigenvalue, and at least one eigenvalue pair The characteristic vector answered, wherein, the biggest expression of eigenvalue is buried in oblivion to count out and is changed on the direction that the characteristic vector of eigenvalue is corresponding Degree is the biggest；Choose first eigenvector corresponding to eigenvalue of maximum or second feature vector corresponding to secondary big eigenvalue, by the One characteristic vector or direction corresponding to second feature vector are as the direction of motion of lesion tissue；For each sub-period, according to In sub-period, radial position-angle position-corresponding the burying in oblivion of axial location identifier combination after data recombination is counted out and One characteristic vector or second feature vector, determine in sub-period lesion tissue motion amplitude information in the movement direction.

Optionally, it is contemplated that the nucleic decay impact on statistical data, the footpath after data recombination in each sub-period is obtained After position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, can be to corresponding respectively the falling into oblivion of each sub-period Vanishing point number carries out decay and revises, and accordingly, processing procedure can be such that according to number in nucleic decay formula and each sub-period Count out according to radial position-angle position-corresponding the burying in oblivion of axial location identifier combination after restructuring, determine in each sub-period Revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination after data recombination is counted out；Root According to after the decay correction that the radial position-angle position-axial location identifier combination after data recombination in each sub-period is corresponding Bury in oblivion and count out, determine the direction of motion of lesion tissue, determine in each sub-period lesion tissue in the described direction of motion Motion amplitude information.

In force, for each sub-period, obtain under this sub-period the radial position-angle position after data recombination- After corresponding the burying in oblivion of axial location identifier combination is counted out, can be according to nucleic decay formula, to data recombination under this sub-period After radial position-angle position-corresponding burying in oblivion of axial location identifier combination count out carry out decay revise, obtain this period of the day from 11 p.m. to 1 a.m Revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination after the lower data recombination of section is counted Mesh, concrete, position, the radial position-angle position after data recombination-axially can be obtained in each sub-period according to formula (1) Put identifier combination revised the burying in oblivion of corresponding decay to count out,

$S^{\′}\left(t\right)=S\left(t\right)\·e^{t\·ln\left(2\right)/T_{1/2}}---\left(1\right)$

Wherein, radial position-angle position-axial location identifier combination pair after data recombination under S ' (t) is sub-period t Revised the burying in oblivion of decay answered is counted out, and S (t) is radial position-angle position after data recombination under sub-period t-axially Burying in oblivion before the decay correction that station location marker combination is corresponding is counted out, T_1/2For the half-life of nucleic,For nucleic Decay formula.

Terminal obtains the radial position-angle position-axial location identifier combination pair in each sub-period after data recombination The decay answered is revised bury in oblivion and count out after, can be according to the radial position after data recombination in each sub-period-angle position Put-revised the burying in oblivion of decay corresponding to axial location identifier combination count out, determine the direction of motion of lesion tissue, determine Each sub-period lesion tissue motion amplitude information in the described direction of motion.Concrete, can be according to default main constituent Analyze PCA algorithm, to by the radial position-angle position-axial location identifier combination pair after data recombination in each sub-period The revised matrix constituted of counting out of burying in oblivion of decay answered carries out principal component analysis, obtains at least one eigenvalue, Yi Jizhi Few eigenvalue characteristic of correspondence vector, wherein, the characteristic vector pair counted out at eigenvalue is buried in oblivion in the biggest expression of eigenvalue On the direction answered, intensity of variation is the biggest；Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second Characteristic vector, using direction corresponding to first eigenvector or second feature vector as the direction of motion of lesion tissue；For often Individual sub-period, according to the decay that the radial position-angle position-axial location identifier combination after data recombination in sub-period is corresponding Revised bury in oblivion count out and first eigenvector or second feature vector, determine in sub-period lesion tissue in motion side Motion amplitude information upwards.

Optionally, it is contemplated that the noise impact on determining the direction of motion, the footpath after data recombination in each sub-period is obtained After position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, can be to corresponding respectively the falling into oblivion of each sub-period Vanishing point number carries out denoising Processing, accordingly, after processing procedure can be such that for data recombination each in each sub-period Radial position in radial position-angle position-axial location identifier combination-angle position combination, according to radial position-angle After the decay correction that radial position-angle position-axial location identifier combination after data recombination belonging to position grouping is corresponding Bury in oblivion and count out, determine that revised the burying in oblivion of decay corresponding to radial position-angle position combination is counted out；By each period of the day from 11 p.m. to 1 a.m Decay in Duan the revised radial position-angle position combination buried in oblivion and count out not less than preset number threshold value, affiliated number According to the decay of radial position-angle position-axial location identifier combination after restructuring revised bury in oblivion to count out be set to zero, Obtain the decay correction that the radial position-angle position-axial location identifier combination after eliminating noise in each sub-period is corresponding After bury in oblivion and count out；According to the radial position-angle position-axial location identifier combination after elimination noise in each sub-period Corresponding revised the burying in oblivion of decay is counted out, and determines the direction of motion of lesion tissue, determines in each sub-period lesion tissue Motion amplitude information in the movement direction.

In force, can be previously stored with quantity threshold in terminal, wherein, each sub-time period can be corresponding to having Quantity threshold, concrete, terminal can be previously stored with percentage ratio, the quantity threshold that each sub-period is corresponding is i.e. this son Maximum in period buries in oblivion the product counted out with preset percentage.The burying in oblivion dot position information and may comprise of PET device collection Noise information, it is contemplated that the noise information impact on determining the direction of motion, can be to data recombination each in each sub-period after Attenuation correction corresponding to radial position-angle position-axial location identifier combination after bury in oblivion to count out and carry out denoising Processing, Concrete, firstly for the radial position-angle position-axial location identifier combination after data recombination each in each sub-period In radial position-angle position combination, can be to the radial direction position after the data recombination belonging to radial position-angle position combination Put-decay corresponding to angle position-axial location identifier combination revised bury in oblivion to count out be added, as radially Revised the burying in oblivion of decay of position-angle position combination correspondence is counted out, then, for each sub-period, it is judged that radially position Put-decay corresponding to angle position combination is revised buries in oblivion the magnitude relationship counted out with preset number threshold value, determines this son Decayed in period the revised radial position-angle position combination buried in oblivion and count out not less than preset number threshold value, and then, can To determine the revised number buried in oblivion and count out not less than belonging to radial position-angle position combination of preset number threshold value that decays According to the radial position-angle position-axial location identifier combination after restructuring, can be by the radial direction position after its affiliated data recombination Put-decay corresponding to angle position-axial location identifier combination revised bury in oblivion to count out be set to zero, decay revised Burying in oblivion counts out exceedes the radial position of preset number threshold value-angle position combination, and the affiliated radial position after data recombination- The decay of angle position-axial location identifier combination is revised, and to bury in oblivion holding of counting out constant, obtains disappearing in each sub-period Decay revised bury in oblivion corresponding except the radial position-angle position-axial location identifier combination after noise is counted out.? After, terminal can be corresponding according to the radial position-angle position-axial location identifier combination after elimination noise in each sub-period Revised the burying in oblivion of decay count out, determine the direction of motion of lesion tissue, determine each sub-period lesion tissue fortune Motion amplitude information on dynamic direction.

Optionally, it is possible to use Principal Component Analysis Algorithm, the direction of motion of lesion tissue is determined and in this direction of motion On motion amplitude information, accordingly, processing procedure can be such that according to preset principal component analysis PCA algorithm, to by each Sub-period eliminates decay corresponding to the radial position-angle position-axial location identifier combination after noise is revised buries in oblivion The matrix constituted of counting out carries out principal component analysis, obtains at least one eigenvalue, and the spy that at least one eigenvalue is corresponding Levying vector, wherein, the biggest expression of eigenvalue is buried in oblivion and is counted out on the direction that the characteristic vector of eigenvalue is corresponding intensity of variation more Greatly；Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by fisrt feature Vector or direction corresponding to second feature vector are as the direction of motion of lesion tissue；For each sub-period, according to sub-period Described in eliminate that decay corresponding to the radial position-angle position-axial location identifier combination after noise is revised buries in oblivion a little Number, and first eigenvector or second feature vector, determine at sub-period lesion tissue motion width in the movement direction Degree information.

In force, corresponding each sub-period, terminal obtains the radial position-angle position after eliminating noise in this sub-period Put-decay corresponding to axial location identifier combination revised bury in oblivion count out after, it can be put in order according to certain, The vector of composition n × 1 dimension, can use X (t_i) represent, wherein, n is positive integer, and n can be the footpath in Sinogram data form To the product of the quantity of the quantity of position, the quantity of angle position and axial location mark, t_iExpression i-th sub-period, i=1, 2 ... M, wherein, M is the number of sub-period, and the most one-dimensional numerical value is i.e. the radial position-angle after the elimination noise that this dimension is corresponding Revised the burying in oblivion of decay spending position-axial location identifier combination corresponding is counted out, and then, can be by the X of each sub-period (t_i) composition matrix, can represent with matrix S, i.e. S=[X (t₁), X (t₂), X (t₃)…X(t_M)].And then, terminal can utilize PCA (Principal Component Analysis, principal component analysis) algorithm, carries out principal component analysis to matrix S.Specifically , after terminal obtains matrix S, the average of each line number value in matrix S can be calculated, each data of matrix S often row are deducted The average that this row is corresponding, obtains matrix S ', and then, the covariance matrix Σ of matrix S ' can be calculated according to formula (2),

$\mathrm{\Σ}=\frac{1}{M}{S}^{\′}\·S^{\′T}---\left(2\right)$

Wherein, S '^TThe transposed matrix of representing matrix S ', and then, eigenvalue and the feature of covariance matrix Σ can be calculated Vector, wherein, it is the biggest that intensity of variation on the direction that the characteristic vector of eigenvalue is corresponding of counting out is buried in oblivion in the biggest expression of eigenvalue, In this programme, the direction that eigenvalue of maximum characteristic of correspondence vector is corresponding can be the corresponding fortune of moving that exercise intensity is maximum Dynamic direction is (due to the intensity of variation of the size response data of eigenvalue, owing to the motion such as respiratory movement or heartbeat causes disease Stove histokinesis, and then, cause the distribution buried in oblivion a little to change, i.e. the data of each position in matrix S change, because of This, can be using the direction of eigenvalue of maximum characteristic of correspondence vector that obtains according to matrix S as the maximum motion of exercise intensity The corresponding direction of motion), by that analogy, certain direction of motion moved can be determined, such as, typically according to the characteristic vector obtained In the case of, respirometric exercise intensity more than the exercise intensity of heartbeat, can by eigenvalue of maximum characteristic of correspondence to The direction of amount correspondence is as the respirometric direction of motion, using direction corresponding for secondary big eigenvalue characteristic of correspondence vector as the heart The dirty direction of motion beaten, in particular cases, if the exercise intensity of heartbeat is more than respirometric exercise intensity, then may be used Using by direction corresponding for eigenvalue of maximum characteristic of correspondence vector as the direction of motion of heartbeat, by corresponding for time big eigenvalue Direction corresponding to characteristic vector as the respirometric direction of motion.And then, it is corresponding that terminal can obtain eigenvalue of maximum First eigenvector or secondary big eigenvalue corresponding second feature vector, by corresponding to first eigenvector or second feature vector Direction, as the direction of motion of lesion tissue, i.e. can be chosen direction corresponding to one of them characteristic vector and make as required The direction of motion for lesion tissue, it is preferred that can be using direction corresponding for first eigenvector as the motion 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, the X after i.e. can calculating first eigenvector or second feature vector and deducting average (t_i) dot product, wherein it is possible to by projection corresponding for each sub-period of obtaining, as in sub-period lesion tissue in motion side Motion amplitude information upwards, can be designated as W (t), i.e. W (t) and can be used to indicate that the lesion tissue motion amplitude at t sub-period Information.

Step 103, according in each sub-period lesion tissue motion amplitude information in the movement direction, determines gate letter Number, carry out PET image reconstruction based on gate-control signal.

Wherein, gate-control signal can be the motor message determined by the motion amplitude information that each sub-period is corresponding.

In force, terminal determines the lesion tissue motion amplitude information in the movement direction at each sub-period After, can be by motion amplitude information corresponding for each sub-period, as gate-control signal.Can extract and swept from gate-control signal Retouch each moment that when object is exhaled maximum, (or during air-breathing maximum) is corresponding, by the duration decile between the most adjacent two moment For the sub-duration of preset number, the dot position information of burying in oblivion detected in period of the day from 11 p.m. to 1 a.m length corresponding in each duration is divided in one Rise, wherein it is possible to think that the dot position information of burying in oblivion being grouped together is to be scanned detection when object is in same motion Arrive, and then, 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.Such as, the moment that scanned object is corresponding when exhaling maximum is t1, t2, t3, t4, t5 respectively, respectively by t2- The duration of t1, t3-t2, t4-t3, t5-t4 is divided into 3 deciles, and numbering is respectively 1-3, can be by all period of the day from 11 p.m. to 1 a.m of numbered 1 The dot position information of respectively burying in oblivion detected in length is grouped together, and buries in oblivion a position according to these and carries out PET image reconstruction, obtains To should the PET image of motion stage, the dot position information of respectively burying in oblivion detected in the period of the day from 11 p.m. to 1 a.m length of all numbered 2 is divided Together, bury in oblivion a position according to these and carry out PET image reconstruction, obtain to should the PET image of motion stage, in all volumes Number it is that the dot position information of respectively burying in oblivion detected in the period of the day from 11 p.m. to 1 a.m length of 3 is grouped together, buries in oblivion a position according to these and carry out PET figure As rebuilding, obtain to should the PET image of motion stage, and then, the PET image of multiple motion stage can be obtained.

Optionally, terminal can be by the motion amplitude information through interpolation processing, as gate-control signal, accordingly, and step The processing procedure of 103 can be such that to each sub-period with at each sub-period lesion tissue motion width in the movement direction Degree information structure motion amplitude information-time function carry out interpolation processing, obtain after interpolation processing motion amplitude information-time Between function；Using the motion amplitude information-time function after interpolation processing as gate-control signal.

In force, terminal obtains terminal and determines the lesion tissue motion width in the movement direction at each sub-period After degree information W (t), based on default interpolation algorithm, W (t) can be carried out interpolation processing, obtain the motion width after interpolation processing Degree information-time function W ' (t), concrete, W (t) can be carried out linear interpolation, by the motion amplitude that adjacent moment is corresponding Information, determines the motion amplitude information that the intermediate time of this adjacent moment is corresponding, can also use non-linear interpolation in this programme Algorithm, is not limited thereof.Obtain motion amplitude information-time function W ' (t) after interpolation processing, W ' (t) can be made For gate-control signal.As such, it is possible to obtain the moment that maximum range of motion information is corresponding more accurately, in order to carry out PET image weight Build.

Optionally, it is also possible to using motion amplitude information corresponding for each sub-period after after filtering as gate-control signal, Accordingly, the processing procedure of step 103 can be such that to each sub-period and in each sub-period lesion tissue in the direction of motion On motion amplitude information constitute motion amplitude information-time function carry out frequency domain transform process, obtain motion amplitude information- The frequency spectrum that time function is corresponding；Based on having the band filter of default center frequency and pre-set bandwidths, frequency spectrum is filtered, Obtain filtered frequency spectrum；Filtered frequency spectrum is carried out the inverse transformation of frequency domain transform, obtains gate-control signal.

In force, after terminal obtains each motion amplitude information W (t) corresponding for sub-period t, W (t) can be carried out soon Speed Fourier transformation (FFT), and then, obtain the frequency spectrum of W (t), wherein, the frequency spectrum of W (t) can be calculated according to formula (3),

H (f)=| FFT (W (t)) | (3)

Wherein, H (f) is i.e. the frequency spectrum of W (t), and f represents frequency, and then, terminal can be to the frequency of motion amplitude information W (t) Spectrum is filtered, wherein it is possible to by band filter (the such as gaussian filtering with default center frequency and pre-set bandwidths Device) it is filtered, obtain filtered frequency spectrum, afterwards, filtered frequency spectrum is carried out inversefouriertransform process, obtains Gate-control signal.

It addition, when motion amplitude information W (t) is motor message corresponding to respiratory movement, default center frequency can be When presetting respiratory frequency value (such as, presetting respiratory frequency value can be 0.5Hz), the gate-control signal obtained is that respiratory movement is corresponding Gate-control signal, when motion amplitude information W (t) is motor message corresponding to heartbeat, default center frequency is for presetting the heart When jumping frequency values (such as, presetting palmic rate value can be 1Hz), the gate-control signal obtained is the gate letter that heartbeat is corresponding Number.

In the embodiment of the present invention, gate-control signal, nothing can be determined according to the dot position information of respectively burying in oblivion that PET device detects Extra external equipment need to be increased and gathers gate-control signal PET device detection respectively buries in oblivion dot position information when, thus, can To reduce diagnostic system cost.

Based on identical technology design, the embodiment of the present invention additionally provides a kind of device determining gate-control signal, such as Fig. 4 institute Showing, this device includes:

Acquisition module 410, for obtain positron emission tomography PET device detect in the period each default time Detect in long sub-period respectively buries in oblivion dot position information；

First determines module 420, for according to described each sub-period corresponding respectively bury in oblivion dot position information, determine focus The direction of motion of tissue, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period；

Second determines module 430, for according in lesion tissue fortune in the described direction of motion described in each sub-period Dynamic amplitude information, determines gate-control signal, carries out PET image reconstruction based on described gate-control signal.

Optionally, as it is shown in figure 5, described first determines module 420, including:

First determines submodule 4201, for for each sub-period, according to described sub-period corresponding respectively bury in oblivion a position Confidence ceases, and determines in described sub-period and different to bury in oblivion corresponding the burying in oblivion of dot position information and count out；

Second determines submodule 4202, for burying in oblivion corresponding the falling into oblivion of dot position information according to different in described each sub-period Vanishing point number, determines the direction of motion of lesion tissue, determines in lesion tissue described in each sub-period in the described direction of motion Motion amplitude information.

Optionally, described first determines submodule 4201, is used for:

For each sub-period, according to the radial position respectively buried in oblivion in dot position information corresponding to described sub-period, angle Position and axial location mark, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted Mesh, wherein, described axial location mark is intended to indicate that the information burying in oblivion axial location a little；

Described second determines submodule 4202, is used for:

For each sub-period, by radial position-angle position-axial location mark group each described in described sub-period A new radial position-angle is merged in the multiple radial position-angle positions meeting preset proximity condition in conjunction Position grouping, and corresponding for multiple radial positions-angle position combination is buried in oblivion sum of counting out, as described new radial direction position Put-angle position corresponding the burying in oblivion of combination count out, and obtains position, radial position-angle position new in described sub-period-axially Put corresponding the burying in oblivion of identifier combination to count out；

Corresponding according to radial position-angle position-axial location identifier combination new described in described each sub-period Bury in oblivion and count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period in the described direction of motion On motion amplitude information.

Optionally, described second determines submodule 4202, is used for:

Radial position-angle position-axial location identifier combination new described in described each sub-period is carried out data Restructuring processes, and obtains the radial position-angle position-axial location identifier combination after data recombination, according to described each period of the day from 11 p.m. to 1 a.m Radial position-angle position new described in Duan-corresponding burying in oblivion of axial location identifier combination is counted out, and determines described every height Radial position-angle position-corresponding burying in oblivion of axial location identifier combination after data recombination described in period is counted out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding burying in oblivion to count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period described Motion amplitude information in the direction of motion.

Optionally, described second determines submodule 4202, is used for:

According to radial position-angle position after data recombination described in nucleic decay formula and described each sub-period- Corresponding the burying in oblivion of axial location identifier combination is counted out, and determines the radial direction position after data recombination described in described each sub-period Put-revised the burying in oblivion of decay corresponding to angle position-axial location identifier combination count out；

According to the radial position-angle position after data recombination described in described each sub-period-axial location mark group Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described second determines submodule 4202, is used for:

For the radial position-angle position-axial location mark after described data recombination each in described each sub-period Know the radial position in combination-angle position combination, after the data recombination belonging to described radial position-angle position combination Revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination count out, determine described radial direction Revised the burying in oblivion of decay of position-angle position combination correspondence is counted out；

By the revised radial position-angle buried in oblivion and count out not less than preset number threshold value that decays in each sub-period Position grouping, the decay of the radial position-angle position-axial location identifier combination after affiliated data recombination is revised falls into oblivion Vanishing point number is set to zero, obtains the radial position-angle position-axial location mark after eliminating noise in described each sub-period Revised the burying in oblivion of decay knowing combination corresponding is counted out；

According to the radial position-angle position-axial location mark group after elimination noise described in described each sub-period Close corresponding revised the burying in oblivion of decay to count out, determine the direction of motion of lesion tissue, determine sick described in each sub-period Stove is organized in the motion amplitude information in the described direction of motion.

Optionally, described second determines submodule 4202, is used for:

According to default principal component analysis PCA algorithm, to by the radial direction after elimination noise described in described each sub-period Position-angle position-axial location identifier combination corresponding decay revised bury in oblivion count out constitute matrix carry out main one-tenth Analyze, obtain at least one eigenvalue, and at least one eigenvalue characteristic of correspondence vector described, wherein, described feature It is the biggest that intensity of variation on the direction that the characteristic vector of described eigenvalue is corresponding of counting out is buried in oblivion in the biggest expression of value；

Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by described First eigenvector or direction corresponding to described second feature vector are as the direction of motion of lesion tissue；

For each sub-period, according to eliminating the radial position-angle position after noise-axially described in described sub-period Revised the burying in oblivion of decay corresponding to station location marker combination is counted out, and described first eigenvector or described second feature to Amount, determines in lesion tissue motion amplitude information in the described direction of motion described in described sub-period.

Optionally, described second determines module 430, is used for:

To described each sub-period with at lesion tissue motion amplitude in the described direction of motion described in each sub-period Motion amplitude information-the time function of information structure carries out interpolation processing, obtains the motion amplitude information-time after interpolation processing Function；

Using the motion amplitude information-time function after described interpolation processing as gate-control signal.

Optionally, as shown in Figure 6, described second determines module 430, including:

Frequency domain transform submodule 4301, for described each sub-period and in lesion tissue described in each sub-period in institute Motion amplitude information-the time function stating the composition of the motion amplitude information in the direction of motion carries out frequency domain transform process, obtains institute State the frequency spectrum that motion amplitude information-time function is corresponding；

Filtering submodule 4302, is used for based on the band 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 described filtered frequency spectrum being carried out the inverse transformation of described frequency domain transform, Obtain gate-control signal.

In the embodiment of the present invention, gate-control signal, nothing can be determined according to the dot position information of respectively burying in oblivion that PET device detects Extra external equipment need to be increased and gathers gate-control signal PET device detection respectively buries in oblivion dot position information when, thus, can To reduce diagnostic system cost.

It should be understood that the device of determination gate-control signal that above-described embodiment provides is when determining gate-control signal, only with The division of above-mentioned each functional module is illustrated, and in actual application, can distribute above-mentioned functions by not as desired With functional module complete, the internal structure of terminal will be divided into different functional modules, with complete described above all Or partial function.It addition, above-described embodiment provide determination gate-control signal device and determine gate-control signal method implement Example belongs to same design, and it implements process and refers to embodiment of the method, repeats no more here.

One of ordinary skill in the art will appreciate that all or part of step realizing above-described embodiment can pass through hardware Completing, it is also possible to instruct relevant hardware by program and complete, described program can be stored 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 only presently preferred embodiments of the present invention, not in order to limit the present invention, all spirit in the present invention and Within principle, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.

Claims (18)

1. the method determining gate-control signal, it is characterised in that described method includes:

The sub-period of the acquisition positron emission tomography PET device each preset duration within the detection period detects Respectively bury in oblivion dot position information；

According to described each sub-period corresponding respectively bury in oblivion dot position information, determine the direction of motion of lesion tissue, determine often Lesion tissue motion amplitude information in the described direction of motion described in individual sub-period；

According in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, determine gate-control signal, PET image reconstruction is carried out based on described gate-control signal.

Method the most according to claim 1, it is characterised in that described according to described each sub-period corresponding respectively bury in oblivion a little Positional information, determines the direction of motion of lesion tissue, determines in lesion tissue described in each sub-period in the described direction of motion Motion amplitude information, including:

For each sub-period, according to described sub-period corresponding respectively bury in oblivion dot position information, determine in described sub-period different Bury in oblivion corresponding the burying in oblivion of dot position information to count out；

Bury in oblivion corresponding the burying in oblivion of dot position information according to difference in described each sub-period to count out, determine the motion of lesion tissue Direction, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period.

Method the most according to claim 2, it is characterised in that described for each sub-period, according to described sub-period pair That answers respectively buries in oblivion dot position information, determines in described sub-period and different to bury in oblivion corresponding the burying in oblivion of dot position information and count out, including:

For each sub-period, according to the radial position respectively buried in oblivion in dot position information corresponding to described sub-period, angle position Identify with axial location, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, its In, described axial location mark is intended to indicate that the information burying in oblivion axial location a little；

Described bury in oblivion corresponding the burying in oblivion of dot position information count out according to different in described each sub-period, determine lesion tissue The direction of motion, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

For each sub-period, by radial position-angle position-axial location identifier combination each described in described sub-period A new radial position-angle position is merged in the multiple radial position-angle positions meeting preset proximity condition Combination, and corresponding for multiple radial positions-angle position combination is buried in oblivion sum of counting out, as described new radial position- Angle position corresponding the burying in oblivion of combination is counted out, and obtains radial position-angle position-axial location mark new in described sub-period Know corresponding the burying in oblivion of combination to count out；

According to radial position-angle position new described in described each sub-period-corresponding burying in oblivion of axial location identifier combination Count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period in the described direction of motion Motion amplitude information.

Method the most according to claim 3, it is characterised in that described according to radial direction new described in described each sub-period Position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, and determines the direction of motion of lesion tissue, determines Lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

Radial position-angle position-axial location identifier combination new described in described each sub-period is carried out data recombination Process, obtain the radial position-angle position-axial location identifier combination after data recombination, according in described each sub-period Described new radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, and determines described each sub-period Described in radial position-angle position-corresponding burying in oblivion of axial location identifier combination after data recombination count out；

According to the radial position-angle position-axial location identifier combination pair after data recombination described in described each sub-period Burying in oblivion of answering is counted out, and determines the direction of motion of lesion tissue, determines in lesion tissue described in each sub-period in described motion Motion amplitude information on direction.

Method the most according to claim 4, it is characterised in that described according to data recombination described in described each sub-period After radial position-angle position-corresponding burying in oblivion of axial location identifier combination count out, determine the motion side of lesion tissue To, determine in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

According to radial position-angle position after data recombination described in nucleic decay formula and described each sub-period-axially Station location marker corresponding the burying in oblivion of combination is counted out, and determines the radial position-angle after data recombination described in described each sub-period Revised the burying in oblivion of decay spending position-axial location identifier combination corresponding is counted out；

According to the radial position-angle position-axial location identifier combination pair after data recombination described in described each sub-period Revised the burying in oblivion of decay answered is counted out, and determines the direction of motion of lesion tissue, determines in focus group described in each sub-period It is woven in the motion amplitude information in the described direction of motion.

Method the most according to claim 5, it is characterised in that described according to data recombination described in described each sub-period After revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination count out, determine focus group The direction of motion knitted, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

For the radial position-angle position after described data recombination each in described each sub-period-axial location mark group Radial position in conjunction-angle position combination, according to the footpath after the data recombination belonging to described radial position-angle position combination Decay revised bury in oblivion corresponding to position-angle position-axial location identifier combination is counted out, and determines described radial direction position Put-revised the burying in oblivion of decay corresponding to angle position combination count out；

By the revised radial position-angle position buried in oblivion and count out not less than preset number threshold value of decaying in each sub-period Combination, the decay of the radial position-angle position-axial location identifier combination after affiliated data recombination is revised buries in oblivion a little Number is set to zero, obtains the radial position-angle position-axial location mark group after eliminating noise in described each sub-period Close corresponding revised the burying in oblivion of decay to count out；

According to the radial position-angle position-axial location identifier combination pair after elimination noise described in described each sub-period Revised the burying in oblivion of decay answered is counted out, and determines the direction of motion of lesion tissue, determines in focus group described in each sub-period It is woven in the motion amplitude information in the described direction of motion.

Method the most according to claim 6, it is characterised in that described according to eliminating noise described in described each sub-period After revised the burying in oblivion of decay corresponding to radial position-angle position-axial location identifier combination count out, determine focus group The direction of motion knitted, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period, including:

According to default principal component analysis PCA algorithm, to by eliminate described in described each sub-period the radial position after noise- The revised matrix constituted of counting out of burying in oblivion of decay corresponding to angle position-axial location identifier combination carries out main constituent and divides Analysis, obtains at least one eigenvalue, and at least one eigenvalue characteristic of correspondence vector described, and wherein, described eigenvalue is more It is the biggest that intensity of variation on the direction that the characteristic vector of described eigenvalue is corresponding of counting out is buried in oblivion in big expression；

Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by described first Characteristic vector or direction corresponding to described second feature vector are as the direction of motion of lesion tissue；

For each sub-period, according to the radial position-angle position-axial location after elimination noise described in described sub-period Revised the burying in oblivion of decay corresponding to identifier combination is counted out, and described first eigenvector or described second feature vector, Determine in lesion tissue motion amplitude information in the described direction of motion described in described sub-period.

Method the most according to claim 1, it is characterised in that described basis in lesion tissue described in each sub-period in institute State the motion amplitude information in the direction of motion, determine gate-control signal, including:

To described each sub-period with in lesion tissue motion amplitude information in the described direction of motion described in each sub-period Motion amplitude information-the time function constituted carries out interpolation processing, obtains the motion amplitude information after interpolation processing-time letter Number；

Using the motion amplitude information-time function after described interpolation processing as gate-control signal.

Method the most according to claim 1, it is characterised in that described basis in lesion tissue described in each sub-period in institute State the motion amplitude information in the direction of motion, determine gate-control signal, including:

To described each sub-period with in lesion tissue motion amplitude information in the described direction of motion described in each sub-period Motion amplitude information-the time function constituted carries out frequency domain transform process, obtains described motion amplitude information-time function corresponding Frequency spectrum；

Based on having the band filter of default center frequency and pre-set bandwidths, described frequency spectrum is filtered, after being filtered Frequency spectrum；

Described filtered frequency spectrum is carried out the inverse transformation of described frequency domain transform, obtains gate-control signal.

10. the device determining gate-control signal, it is characterised in that described device includes:

Acquisition module, for obtaining the son of the positron emission tomography PET device each preset duration within the detection period Detect in period respectively buries in oblivion dot position information；

First determines module, for according to described each sub-period corresponding respectively bury in oblivion dot position information, determine lesion tissue The direction of motion, determines in lesion tissue motion amplitude information in the described direction of motion described in each sub-period；

Second determines module, for according to lesion tissue described in each sub-period in the described direction of motion motion amplitude letter Breath, determines gate-control signal, carries out PET image reconstruction based on described gate-control signal.

11. devices according to claim 10, it is characterised in that described first determines module, including:

First determines submodule, for for each sub-period, according to described sub-period corresponding respectively bury in oblivion dot position information, really Different in fixed described sub-period bury in oblivion corresponding the burying in oblivion of dot position information and count out；

Second determines submodule, counts for burying in oblivion corresponding the burying in oblivion of dot position information according to difference in described each sub-period Mesh, determines the direction of motion of lesion tissue, determines in lesion tissue motion in the described direction of motion described in each sub-period Amplitude information.

12. devices according to claim 11, it is characterised in that described first determines submodule, is used for:

For each sub-period, according to the radial position respectively buried in oblivion in dot position information corresponding to described sub-period, angle position Identify with axial location, determine that each radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, its In, described axial location mark is intended to indicate that the information burying in oblivion axial location a little；

Described second determines submodule, is used for:

For each sub-period, by radial position-angle position-axial location identifier combination each described in described sub-period A new radial position-angle position is merged in the multiple radial position-angle positions meeting preset proximity condition Combination, and corresponding for multiple radial positions-angle position combination is buried in oblivion sum of counting out, as described new radial position- Angle position corresponding the burying in oblivion of combination is counted out, and obtains radial position-angle position-axial location mark new in described sub-period Know corresponding the burying in oblivion of combination to count out；

According to radial position-angle position new described in described each sub-period-corresponding burying in oblivion of axial location identifier combination Count out, determine the direction of motion of lesion tissue, determine in lesion tissue described in each sub-period in the described direction of motion Motion amplitude information.

13. devices according to claim 12, it is characterised in that described second determines submodule, is used for:

Radial position-angle position-axial location identifier combination new described in described each sub-period is carried out data recombination Process, obtain the radial position-angle position-axial location identifier combination after data recombination, according in described each sub-period Described new radial position-angle position-corresponding burying in oblivion of axial location identifier combination is counted out, and determines described each sub-period Described in radial position-angle position-corresponding burying in oblivion of axial location identifier combination after data recombination count out；

According to the radial position-angle position-axial location identifier combination pair after data recombination described in described each sub-period Burying in oblivion of answering is counted out, and determines the direction of motion of lesion tissue, determines in lesion tissue described in each sub-period in described motion Motion amplitude information on direction.

14. devices according to claim 13, it is characterised in that described second determines submodule, is used for:

According to radial position-angle position after data recombination described in nucleic decay formula and described each sub-period-axially Station location marker corresponding the burying in oblivion of combination is counted out, and determines the radial position-angle after data recombination described in described each sub-period Revised the burying in oblivion of decay spending position-axial location identifier combination corresponding is counted out；

According to the radial position-angle position-axial location identifier combination pair after data recombination described in described each sub-period Revised the burying in oblivion of decay answered is counted out, and determines the direction of motion of lesion tissue, determines in focus group described in each sub-period It is woven in the motion amplitude information in the described direction of motion.

15. devices according to claim 14, it is characterised in that described second determines submodule, is used for:

For the radial position-angle position after described data recombination each in described each sub-period-axial location mark group Radial position in conjunction-angle position combination, according to the footpath after the data recombination belonging to described radial position-angle position combination Decay revised bury in oblivion corresponding to position-angle position-axial location identifier combination is counted out, and determines described radial direction position Put-revised the burying in oblivion of decay corresponding to angle position combination count out；

By the revised radial position-angle position buried in oblivion and count out not less than preset number threshold value of decaying in each sub-period Combination, the decay of the radial position-angle position-axial location identifier combination after affiliated data recombination is revised buries in oblivion a little Number is set to zero, obtains the radial position-angle position-axial location mark group after eliminating noise in described each sub-period Close corresponding revised the burying in oblivion of decay to count out；

According to the radial position-angle position-axial location identifier combination pair after elimination noise described in described each sub-period Revised the burying in oblivion of decay answered is counted out, and determines the direction of motion of lesion tissue, determines in focus group described in each sub-period It is woven in the motion amplitude information in the described direction of motion.

16. devices according to claim 15, it is characterised in that described second determines submodule, is used for:

According to default principal component analysis PCA algorithm, to by eliminate described in described each sub-period the radial position after noise- The revised matrix constituted of counting out of burying in oblivion of decay corresponding to angle position-axial location identifier combination carries out main constituent and divides Analysis, obtains at least one eigenvalue, and at least one eigenvalue characteristic of correspondence vector described, and wherein, described eigenvalue is more It is the biggest that intensity of variation on the direction that the characteristic vector of described eigenvalue is corresponding of counting out is buried in oblivion in big expression；

Choose first eigenvector corresponding to eigenvalue of maximum or secondary big eigenvalue corresponding second feature vector, by described first Characteristic vector or direction corresponding to described second feature vector are as the direction of motion of lesion tissue；

For each sub-period, according to the radial position-angle position-axial location after elimination noise described in described sub-period Revised the burying in oblivion of decay corresponding to identifier combination is counted out, and described first eigenvector or described second feature vector, Determine in lesion tissue motion amplitude information in the described direction of motion described in described sub-period.

17. devices according to claim 10, it is characterised in that described second determines module, is used for:

To described each sub-period with in lesion tissue motion amplitude information in the described direction of motion described in each sub-period Motion amplitude information-the time function constituted carries out interpolation processing, obtains the motion amplitude information after interpolation processing-time letter Number；

Using the motion amplitude information-time function after described interpolation processing as gate-control signal.

18. devices according to claim 10, it is characterised in that described second determines module, including:

Frequency domain transform submodule, for described each sub-period and in lesion tissue described in each sub-period in described motion side Motion amplitude information-time function that motion amplitude information upwards is constituted carries out frequency domain transform process, obtains described motion width The frequency spectrum that degree information-time function is corresponding；

Filtering submodule, is used for, based on the band filter with default center frequency and pre-set bandwidths, carrying out described frequency spectrum Filtering, obtains filtered frequency spectrum；

Frequency domain inverse transformation submodule, for described filtered frequency spectrum carries out the inverse transformation of described frequency domain transform, is gated Signal.