CN105266813B - A kind of motor message extracting method and device from gated three-dimensional cardiac imaging - Google Patents
A kind of motor message extracting method and device from gated three-dimensional cardiac imaging Download PDFInfo
- Publication number
- CN105266813B CN105266813B CN201510855793.XA CN201510855793A CN105266813B CN 105266813 B CN105266813 B CN 105266813B CN 201510855793 A CN201510855793 A CN 201510855793A CN 105266813 B CN105266813 B CN 105266813B
- Authority
- CN
- China
- Prior art keywords
- signal
- directions
- heartbeat
- dimensional
- cardiac imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The present invention provides a kind of motor message extracting methods and device from gated three-dimensional cardiac imaging, are related to three-dimensional cardiac imaging technical field.This method includes：According to TureFISP sequences, cine cardiac imaging data are obtained using the radial mixing acquisition mode of the three-dimensional of more echoes；Each acquisition layer in cine cardiac imaging data radially samples the signal value that multiple points are obtained on line, and determines the average value of the signal value of multiple points, as navigation signal；One-dimensional discrete Fourier transformation is carried out to image area according to navigation signal, and is weighted and acquired heartbeat and respiratory movement mixed signal；Heartbeat and respiratory movement mixed signal are filtered by bandpass filter, get heartbeat signal and breath signal；According to heartbeat signal and breath signal, cine cardiac imaging data permutation is synchronized, generates three-dimensional cardiac film image.The present invention solves the inaccuracy of the motor message acquired in the prior art, the problem of influencing to rebuild the quality of heart film image.
Description
Technical field
The present invention relates to three-dimensional cardiac imaging technical field more particularly to a kind of movement letters from gated three-dimensional cardiac imaging
Number extracting method and device.
Background technology
Currently, having benefited from splendid resolution ratio of the magnetic resonance imaging with soft tissue, any angle imaging and non-destructive testing etc.
Numerous advantages, cardiac magnetic resonance imaging (cardiac magnetic resonance imaging, CMRI) are evaluation heart
Structure and function provides powerful measure, has been applied to the science of cardiac shape, function, perfusion, myocardial ischemia and activity etc.
Research.However traditional magnetic resonance imaging speed is slower, and it is long the time required to imaging, it is easy by examinee's heartbeat and exhales
The influence for inhaling movement causes imaging process to introduce motion artifacts, limits CMRI technologies and examined in Disease Clinicals such as heart, angiocarpy
Effect in disconnected.Therefore, how to inhibit motion artifacts, ensure that picture quality has become the pass that CRMI further develops and applies
Key.
And it is current, magnetic resonance inhibits motion artifacts imaging method to be broadly divided into two classes.The first kind is to be based on quick sampling
Mr imaging technique, including：Fast imaging sequences (asymmetric blipped echoplanar single pulse technique EPI, steady state free precession based on physics
Sequence SSFP, single excitation spin-echo sequence HASTE etc.)；Multi-channel parallel imaging method based on system hardware；Based on information
The sparse imaging method of theory.The appearance of these fast imaging techniques promotes the development of cine cardiac imaging, but this kind of
Requirement of the three-dimensional cardiac film imaging for temporal resolution is still not achieved in method, it is still necessary to overcome technical difficult point and bottle
Neck.Second class is to apply electrocardio and respiration gate control in imaging process based on gate airmanship, eliminates movement for CMRI's
It influences.Currently, the clinical CMRI technologies used need examinee repeatedly to hold the breath or to electrocardiogram
(Electrocardiogram, ECG) and breath signal are carried out with CMRI data collection synchronous, this just extends sweep time, and
And increase the complexity of pulse train.In addition, ECG signal is generally obtained by heart leads, gradient being switched fast and penetrating
Frequency pulse meeting severe jamming ECG signal, field strength is higher, and interference is stronger, causes to be difficult to obtain stable ECG signal.
How obtaining reliable ECG signal, to further apply clinical research to CMRI most important.In recent years, CMRI is led
There is a kind of new gate film imaging technique certainly in domain --- and from gated imaging technique, this method utilizes acquisition for heart and breathing
CMRI data acquisitions heart and respiratory movement signal, then according to the heartbeat breath signal of extraction by CMRI data permutations
Extremely in different movement phases, to obtain heart film image.This imaging technique from gate had not both needed electrocardio triggering,
It holds the breath without examinee, greatly improves scan efficiency.In addition, someone has further developed the technology, it is proposed that
Heart based on three-dimensional radial mixing acquisition is sampled in layer choosing direction using Descartes, phase and frequency from imaging method is gated
Direction is encoded to sample using radial, using layer choosing direction navigation signal weighted value as after motor message filtering with obtain heartbeat and
Breath signal.This method is successfully realized three-dimensional whole-heartedly dynamic film and is imaged, however due to being influenced by eddy current effect, each layer
Navigation signal position shifts, and the motor message obtained is caused to be inaccurate reliably, the final heart film figure for influencing to rebuild
Picture.
Invention content
The embodiment of the present invention provides a kind of motor message extracting method and device from gated three-dimensional cardiac imaging, with solution
Certainly current three-dimensional whole-heartedly dynamic film imaging is influenced by eddy current effect, and each layer of navigation signal position shifts, and causes
The motor message of acquisition is inaccurate, the problem of influencing the reconstruction of heart film image.
In order to achieve the above objectives, the present invention adopts the following technical scheme that：
A kind of motor message extracting method from gated three-dimensional cardiac imaging, including：
According to TureFISP sequences, cine cardiac imaging number is obtained using the radial mixing acquisition mode of the three-dimensional of more echoes
According to；The cine cardiac imaging data include each acquisition layer radially sampling line；
The signal value of multiple points is obtained on each acquisition layer radially sampling line, and the signal value of determining the multiple point is flat
Mean value, using the average value of the signal value of the multiple point as navigation signal；
One-dimensional discrete Fourier transformation is carried out to image area according to the navigation signal, and is weighted and acquired heartbeat and exhale
Inhale movement mixed signal；
The heartbeat and respiratory movement mixed signal are filtered by bandpass filter, get heartbeat signal and exhaled
Inhale signal；
According to the heartbeat signal and breath signal, the cine cardiac imaging data permutation is synchronized, generates three
Tie up heart film image.
Herein, described according to TureFISP sequences, cardiac electric is obtained using the radial mixing acquisition mode of the three-dimensional of more echoes
Shadow imaging data, including：
It obtains and relevant first coordinate system of cine cardiac imaging data；First coordinate system is mutually perpendicular to including three
Acquisition direction, respectively directions kz, the directions kx, the directions ky；
Descartes's sampling is carried out along the directions kz, acquires multiple acquisition layers radially sampling line；
The multiple acquisition layer is radially sampled to one predetermined angle of plane internal rotation that line forms in the directions kx and the directions ky,
And continue to carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data.
Specifically, carrying out one-dimensional discrete Fourier transformation to image area according to the navigation signal, and weights and acquire
Heartbeat and respiratory movement mixed signal, including：
According to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；
M (z) is amplitude of the navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
Specifically, the heartbeat and respiratory movement mixed signal are filtered by bandpass filter, heartbeat is got
Signal and breath signal, including：
According to human normal palmic rate range 0.1HZ to 0.5HZ, human normal respiratory rate 0.6HZ to 3HZ is to described
Heartbeat and respiratory movement mixed signal are filtered, and obtain the heartbeat signal and breath signal.
A kind of motor message extraction element from gated three-dimensional cardiac imaging, including：
Cine cardiac imaging data capture unit, for according to TureFISP sequences, being mixed using the three-dimensional of more echoes is radial
It closes acquisition mode and obtains cine cardiac imaging data；The cine cardiac imaging data include each acquisition layer radially sampling line；
Navigation signal determination unit for obtaining the signal value of multiple points on each acquisition layer radially sampling line, and determines
The average value of the signal value of the multiple point, using the average value of the signal value of the multiple point as navigation signal；
Mixed signal determination unit is used for according to navigation signal progress one-dimensional discrete Fourier transformation to image area,
And it weights and acquires heartbeat and respiratory movement mixed signal；
Filter processing unit, for the heartbeat and respiratory movement mixed signal to be filtered by bandpass filter,
Get heartbeat signal and breath signal；
Three-dimensional cardiac film image generation unit is used for according to the heartbeat signal and breath signal, to the cardiac electric
Shadow imaging data rearranges synchronization, generates three-dimensional cardiac film image.
Specifically, the cine cardiac imaging data capture unit, including：
Coordinate system acquisition module, for obtaining and relevant first coordinate system of cine cardiac imaging data；Described first sits
Mark system includes three orthogonal acquisition directions, the respectively directions kz, the directions kx, the directions ky；
Three-dimensional radial mixing sampling module acquires multiple acquisitions for carrying out Descartes's sampling along the directions kz
The radial sampling line of layer；The multiple acquisition layer plane internal rotation one that radially sampling line forms in the directions kx with the directions ky is preset
Angle, and continue to carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data.
In addition, the mixed signal determination unit, is specifically used for：
According to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；
M (z) is amplitude of the navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
In addition, the filter processing unit, is specifically used for：
According to human normal palmic rate range 0.1HZ to 0.5HZ, human normal respiratory rate 0.6HZ to 3HZ is to described
Heartbeat and respiratory movement mixed signal are filtered, and obtain the heartbeat signal and breath signal.
A kind of motor message extracting method and device from gated three-dimensional cardiac imaging provided in an embodiment of the present invention, according to
TureFISP sequences obtain cine cardiac imaging data using the radial mixing acquisition mode of the three-dimensional of more echoes；In each acquisition layer
The signal value of multiple points is obtained on radial sampling line, and determines the average value of the signal value of multiple points, by the signal value of multiple points
Average value as navigation signal；One-dimensional discrete Fourier transformation is carried out to image area according to navigation signal, and is weighted and obtained
To heartbeat and respiratory movement mixed signal；Heartbeat and respiratory movement mixed signal are filtered by bandpass filter, obtained
To heartbeat signal and breath signal；According to heartbeat signal and breath signal, cine cardiac imaging data permutation is synchronized, it is raw
At three-dimensional cardiac film image.Each acquisition layer is not only radially sampled the median of line by the present invention when obtaining navigation signal
As navigation signal, but the signal value of multiple points is obtained on each acquisition layer radially sampling line, and determines the signal of multiple points
The average value of value is as navigation signal, this is because being influenced by eddy current effect, each radially samples line signal maximum
Position may shift, and the median that each acquisition layer radially samples line may be differed with practical navigation signal, cause to filter
The motor message obtained is inaccurate, seriously affects and rebuilds heart film picture quality.As it can be seen that present invention institute compared with the existing technology
Acquisition motor message is more accurate, and the heart film picture quality rebuild is improved.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention without having to pay creative labor, may be used also for those of ordinary skill in the art
With obtain other attached drawings according to these attached drawings.
Fig. 1 is a kind of flow of motor message extracting method from gated three-dimensional cardiac imaging provided in an embodiment of the present invention
Figure one；
Fig. 2 is a kind of flow of motor message extracting method from gated three-dimensional cardiac imaging provided in an embodiment of the present invention
Figure two；
Fig. 3 is the schematic diagram of relevant first coordinate system of cine cardiac imaging data provided in an embodiment of the present invention；
Fig. 4 is the schematic diagram of the relevant coordinate system of cine cardiac imaging data in the prior art；
Fig. 5 is schematic diagram of the result compared with the result that the prior art generates caused by the embodiment of the present invention；
Fig. 6 provides a kind of structure of the motor message extraction element from gated three-dimensional cardiac imaging for the embodiment of the present invention and shows
It is intended to one；
Fig. 7 provides a kind of structure of the motor message extraction element from gated three-dimensional cardiac imaging for the embodiment of the present invention and shows
It is intended to two.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
As shown in Figure 1, the embodiment of the present invention provides a kind of motor message extracting method from gated three-dimensional cardiac imaging, packet
It includes：
Step 101, according to TureFISP sequences, heart film is obtained using the radial mixing acquisition mode of the three-dimensional of more echoes
Imaging data.
Wherein, the cine cardiac imaging data include each acquisition layer radially sampling line.The TureFISP sequences are true
Stable state precession fast imaging sequences.
Step 102 obtains the signal value of multiple points on each acquisition layer radially sampling line, and determines the letter of the multiple point
Number value average value, using the average value of the signal value of the multiple point as navigation signal.
Step 103 carries out one-dimensional discrete Fourier transformation to image area according to the navigation signal, and weights and acquire
Heartbeat and respiratory movement mixed signal.
The heartbeat and respiratory movement mixed signal are filtered by step 104 by bandpass filter, get heartbeat
Signal and breath signal.
Step 105, according to the heartbeat signal and breath signal, it is same to the cine cardiac imaging data permutation
Step generates three-dimensional cardiac film image.
A kind of motor message extracting method from gated three-dimensional cardiac imaging provided in an embodiment of the present invention, according to
TureFISP sequences obtain cine cardiac imaging data using the radial mixing acquisition mode of the three-dimensional of more echoes；In each acquisition layer
The signal value of multiple points is obtained on radial sampling line, and determines the average value of the signal value of multiple points, by the signal value of multiple points
Average value as navigation signal；One-dimensional discrete Fourier transformation is carried out to image area according to navigation signal, and is weighted and obtained
To heartbeat and respiratory movement mixed signal；Heartbeat and respiratory movement mixed signal are filtered by bandpass filter, obtained
To heartbeat signal and breath signal；According to heartbeat signal and breath signal, cine cardiac imaging data permutation is synchronized, it is raw
At three-dimensional cardiac film image.Each acquisition layer is not only radially sampled the median of line by the present invention when obtaining navigation signal
As navigation signal, but the signal value of multiple points is obtained on each acquisition layer radially sampling line, and determines the signal of multiple points
The average value of value is as navigation signal, this is because being influenced by eddy current effect, each radially samples line signal maximum
Position may shift, and the median that each acquisition layer radially samples line may be differed with practical navigation signal, cause to filter
The motor message obtained is inaccurate, seriously affects and rebuilds heart film picture quality.As it can be seen that present invention institute compared with the existing technology
Acquisition motor message is more accurate, and the heart film picture quality rebuild is improved.
In order to make those skilled in the art be better understood by the present invention, a more specifically embodiment is set forth below,
As shown in Fig. 2, the embodiment of the present invention provides a kind of motor message extracting method from gated three-dimensional cardiac imaging, including：
Step 201 obtains and relevant first coordinate system of cine cardiac imaging data；First coordinate system includes three
Orthogonal acquisition direction, the respectively directions kz, the directions kx, the directions ky.
Step 202 carries out Descartes's sampling along the directions kz, acquires multiple acquisition layers radially sampling line.
The multiple acquisition layer is radially sampled the plane internal rotation one that line forms in the directions kx and the directions ky by step 203
Predetermined angle, and continue to carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data.
Wherein, the cine cardiac imaging data include each acquisition layer radially sampling line.For example, it is desired in the directions heart kz
Upper 14 layers of acquisition, then radial direction, which samples line, just 14.But it is not only limited to this, in actual example, generally radially samples line
Item number it is numerous, such as 4800.
Step 204 obtains the signal value of multiple points on each acquisition layer radially sampling line, and determines the letter of the multiple point
Number value average value, using the average value of the signal value of the multiple point as navigation signal.
For example, the first coordinate system as shown in Figure 3, including the directions kz, the directions kx, the directions ky (constitute a virtual circle
Column), with time change, radially plane internal rotation one of the sampling line in the directions kx and the directions ky composition is default for multiple acquisition layers
Angle (such as 112.2 °), wherein 1DFT indicate one-dimensional discrete Fourier transformation.In the prior art, as shown in figure 4, generally will
Acquisition layer radially samples the midrange of line as navigation signal.And the present invention is unlike the prior art, as shown in figure 3, it is each
Acquisition layer radially samples the signal value (such as 5 in Fig. 3) that multiple points are obtained on line, in this way, in the letter for getting multiple points
After number value S1, S2, S3, S4, S5, navigation signal can be determined as：It obtains in this way
Navigation signal closer to navigation signal actual value.
Step 205 carries out one-dimensional discrete Fourier transformation to image area according to the navigation signal, and weights and acquire
Heartbeat and respiratory movement mixed signal.
It in one embodiment, herein can be according to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；
M (z) is amplitude of the navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
Step 206, according to human normal palmic rate range 0.1HZ to 0.5HZ, human normal respiratory rate 0.6HZ is extremely
3HZ is filtered the heartbeat and respiratory movement mixed signal, obtains the heartbeat signal and breath signal.
Step 207, according to the heartbeat signal and breath signal, it is same to the cine cardiac imaging data permutation
Step generates three-dimensional cardiac film image.
Due to being not depend on the triggering of any physiological signal from acquisition sequential in the scanning of gated cardiac magnetic resonant imaging examination,
Only the data of image and gate-control signal are marked plus the corresponding time.Therefore, the heartbeat that can be obtained herein according to filtering is believed
Number and the gomma of breath signal rearrange organization chart picture sampling line, and then reconstruct the three-dimensional cardiac film of different phases
Image.For example heart heart beat cycle can be divided into 22 sections by we, that is, be defined as having 22 different hearts in a heart beat cycle
Motion state.There are 22 periods between two wave crests in label heartbeat signal period in this way, then adopting according to sampling line
All sampling lines are rearranged the three-dimensional heart for being stored in the corresponding period, and then reconstructing 22 different phases by the collection time
Dirty film image.
Through the above steps 201 to the available result of step 207 compared with prior art, as shown in Figure 5.Wherein,
The parts a in Fig. 5 are heartbeat signal waveform in the prior art；The parts b are the getable heartbeat signal of the embodiment of the present invention
Waveform；The parts c are ECG electrocardiosignals；D partial left sides are the three-dimensional cardiac film image of the prior art, it can be seen that in arrow
There are apparent artifacts for directed section, and d partial right sides are three-dimensional cardiac film image caused by the embodiment of the present invention, are had
Effect inhibits the artifact that eddy current effect generates, and improves the accuracy of motor message, the three-dimensional cardiac film image matter rebuild
Amount is preferable.
A kind of motor message extracting method from gated three-dimensional cardiac imaging provided in an embodiment of the present invention, according to
TureFISP sequences obtain cine cardiac imaging data using the radial mixing acquisition mode of the three-dimensional of more echoes；In each acquisition layer
The signal value of multiple points is obtained on radial sampling line, and determines the average value of the signal value of multiple points, by the signal value of multiple points
Average value as navigation signal；One-dimensional discrete Fourier transformation is carried out to image area according to navigation signal, and is weighted and obtained
To heartbeat and respiratory movement mixed signal；Heartbeat and respiratory movement mixed signal are filtered by bandpass filter, obtained
To heartbeat signal and breath signal；According to heartbeat signal and breath signal, cine cardiac imaging data permutation is synchronized, it is raw
At three-dimensional cardiac film image.Each acquisition layer is not only radially sampled the median of line by the present invention when obtaining navigation signal
As navigation signal, but the signal value of multiple points is obtained on each acquisition layer radially sampling line, and determines the signal of multiple points
The average value of value is as navigation signal, this is because being influenced by eddy current effect, each radially samples line signal maximum
Position may shift, and the median that each acquisition layer radially samples line may be differed with practical navigation signal, cause to filter
The motor message obtained is inaccurate, seriously affects and rebuilds heart film picture quality.As it can be seen that present invention institute compared with the existing technology
Acquisition motor message is more accurate, and the heart film picture quality rebuild is improved.
Corresponding to above-mentioned Fig. 1 and embodiment of the method shown in Fig. 2, as shown in fig. 6, the embodiment of the present invention provides one kind from door
The motor message extraction element of three-dimensional cardiac imaging is controlled, including：
Cine cardiac imaging data capture unit 31, can be radial using the three-dimensional of more echoes according to TureFISP sequences
Mixing acquisition mode obtains cine cardiac imaging data；Cine cardiac imaging data include each acquisition layer radially sampling line.
Navigation signal determination unit 32 can obtain the signal value of multiple points on each acquisition layer radially sampling line, and really
The average value of the signal value of fixed multiple points, using the average value of the signal value of multiple points as navigation signal.
Mixed signal determination unit 33 can carry out one-dimensional discrete Fourier transformation to image area according to navigation signal, and
Weighting acquires heartbeat and respiratory movement mixed signal.
Heartbeat and respiratory movement mixed signal can be filtered by bandpass filter, be obtained by filter processing unit 34
Get heartbeat signal and breath signal.
Three-dimensional cardiac film image generation unit 35, can be according to heartbeat signal and breath signal, to cine cardiac imaging
Data permutation synchronizes, and generates three-dimensional cardiac film image.
In addition, as shown in fig. 7, the cine cardiac imaging data capture unit 31, may include：
Coordinate system acquisition module 311 can obtain and relevant first coordinate system of cine cardiac imaging data；First coordinate
System includes three orthogonal acquisition directions, the respectively directions kz, the directions kx, the directions ky.
Three-dimensional radial mixing sampling module 312, can carry out Descartes's sampling along the directions kz, acquire multiple acquisition layers
Radial sampling line；Multiple acquisition layers are radially sampled to one predetermined angle of plane internal rotation that line forms in the directions kx and the directions ky,
And continue to carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data.
In addition, mixed signal determination unit 33, it specifically can be according to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；m(z)
The amplitude for being navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
In addition, filter processing unit 34, it specifically can be according to human normal palmic rate range 0.1HZ to 0.5HZ, human body
Normal respiratory rate 0.6HZ to 3HZ is filtered heartbeat and respiratory movement mixed signal, obtains heartbeat signal and breathing is believed
Number.
It is worth noting that a kind of motor message from gated three-dimensional cardiac imaging provided in an embodiment of the present invention extracts dress
The specific implementation set may refer to the embodiment of the method corresponding to above-mentioned Fig. 1 and Fig. 2, and details are not described herein again.
A kind of motor message extraction element from gated three-dimensional cardiac imaging provided in an embodiment of the present invention, according to
TureFISP sequences obtain cine cardiac imaging data using the radial mixing acquisition mode of the three-dimensional of more echoes；In each acquisition layer
The signal value of multiple points is obtained on radial sampling line, and determines the average value of the signal value of multiple points, by the signal value of multiple points
Average value as navigation signal；One-dimensional discrete Fourier transformation is carried out to image area according to navigation signal, and is weighted and obtained
To heartbeat and respiratory movement mixed signal；Heartbeat and respiratory movement mixed signal are filtered by bandpass filter, obtained
To heartbeat signal and breath signal；According to heartbeat signal and breath signal, cine cardiac imaging data permutation is synchronized, it is raw
At three-dimensional cardiac film image.Each acquisition layer is not only radially sampled the median of line by the present invention when obtaining navigation signal
As navigation signal, but the signal value of multiple points is obtained on each acquisition layer radially sampling line, and determines the signal of multiple points
The average value of value is as navigation signal, this is because being influenced by eddy current effect, each radially samples line signal maximum
Position may shift, and the median that each acquisition layer radially samples line may be differed with practical navigation signal, cause to filter
The motor message obtained is inaccurate, seriously affects and rebuilds heart film picture quality.As it can be seen that present invention institute compared with the existing technology
Acquisition motor message is more accurate, and the heart film picture quality rebuild is improved.
It should be understood by those skilled in the art that, the embodiment of the present invention can be provided as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the present invention
Apply the form of example.Moreover, the present invention can be used in one or more wherein include computer usable program code computer
The computer program production implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The present invention be with reference to according to the method for the embodiment of the present invention, the flow of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram
The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided
Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine so that the instruction generation executed by computer or the processor of other programmable data processing devices can be real
The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to
Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or
The instruction offer executed on other programmable devices may be implemented in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Specific embodiment is applied in the present invention, and principle and implementation of the present invention are described, above example
Explanation be merely used to help understand the present invention method and its core concept；Meanwhile for those of ordinary skill in the art,
According to the thought of the present invention, there will be changes in the specific implementation manner and application range, in conclusion in this specification
Appearance should not be construed as limiting the invention.
Claims (6)
1. a kind of motor message extracting method from gated three-dimensional cardiac imaging, which is characterized in that including：
According to TureFISP sequences, cine cardiac imaging data are obtained using the radial mixing acquisition mode of the three-dimensional of more echoes；Institute
It includes each acquisition layer radially sampling line to state cine cardiac imaging data；
On each acquisition layer radially sampling line the signal values of multiple points is obtained, and determines that the signal value of the multiple point is averaged
Value, using the average value of the signal value of the multiple point as navigation signal；
One-dimensional discrete Fourier transformation is carried out to image area according to the navigation signal, and is weighted and acquired heartbeat and breathing fortune
Dynamic mixed signal；
The heartbeat and respiratory movement mixed signal are filtered by bandpass filter, heartbeat signal is got and breathing is believed
Number；
According to the heartbeat signal and breath signal, the cine cardiac imaging data permutation is synchronized, generates the three-dimensional heart
Dirty film image；
It is described according to TureFISP sequences, cine cardiac imaging number is obtained using the radial mixing acquisition mode of the three-dimensional of more echoes
According to, including：
It obtains and relevant first coordinate system of cine cardiac imaging data；First coordinate system includes three and orthogonal adopts
Collect direction, the respectively directions kz, the directions kx, the directions ky；
Descartes's sampling is carried out along the directions kz, acquires multiple acquisition layers radially sampling line；
By the multiple acquisition layer one predetermined angle of plane internal rotation that radially sampling line forms in the directions kx and the directions ky, and after
Continue and carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data；
On each acquisition layer radially sampling line the signal values of multiple points is obtained, and determines that the signal value of the multiple point is averaged
Value, using the average value of the signal value of the multiple point as navigation signal, including：
In the first coordinate system, including the directions kz, the directions kx, the directions ky, with time change, multiple acquisition layers radially sample
Line obtains 5 points in one predetermined angle of plane internal rotation in the directions kx and the directions ky composition on each acquisition layer radially sampling line
Signal value navigation signal is determined as after signal value S1, S2, S3, S4, the S5 for getting 5 points：Actual value of the obtained navigation signal closer to navigation signal.
2. the motor message extracting method according to claim 1 from gated three-dimensional cardiac imaging, which is characterized in that according to
The navigation signal carries out one-dimensional discrete Fourier transformation to image area, and weights and acquire heartbeat and respiratory movement mixing letter
Number, including：
According to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；m(z)
For amplitude of the navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
3. the motor message extracting method according to claim 2 from gated three-dimensional cardiac imaging, which is characterized in that by institute
It states heartbeat and respiratory movement mixed signal is filtered by bandpass filter, get heartbeat signal and breath signal, including：
According to human normal palmic rate range 0.1HZ to 0.5HZ, human normal respiratory rate 0.6HZ to 3HZ is to the heartbeat
And respiratory movement mixed signal is filtered, and obtains the heartbeat signal and breath signal.
4. a kind of motor message extraction element from gated three-dimensional cardiac imaging, which is characterized in that including：
Cine cardiac imaging data capture unit, for according to TureFISP sequences, being adopted using the radial mixing of the three-dimensional of more echoes
Mode set obtains cine cardiac imaging data；The cine cardiac imaging data include each acquisition layer radially sampling line；
Navigation signal determination unit, the signal value for obtaining multiple points on the radial sampling line of each acquisition layer, and described in determination
The average value of the signal value of multiple points, using the average value of the signal value of the multiple point as navigation signal；
Mixed signal determination unit for carrying out one-dimensional discrete Fourier transformation to image area according to the navigation signal, and adds
Power acquires heartbeat and respiratory movement mixed signal；
Filter processing unit is obtained for the heartbeat and respiratory movement mixed signal to be filtered by bandpass filter
To heartbeat signal and breath signal；
Three-dimensional cardiac film image generation unit, for according to the heartbeat signal and breath signal, to the heart film at
As data permutation synchronization, three-dimensional cardiac film image is generated；
The cine cardiac imaging data capture unit, including：
Coordinate system acquisition module, for obtaining and relevant first coordinate system of cine cardiac imaging data；First coordinate system
Including three orthogonal acquisition directions, the respectively directions kz, the directions kx, the directions ky；
Three-dimensional radial mixing sampling module acquires multiple acquisition layer diameters for carrying out Descartes's sampling along the directions kz
To sampling line；By the multiple acquisition layer one preset angle of plane internal rotation that radially sampling line forms in the directions kx and the directions ky
Degree, and continue to carry out Descartes's sampling in the directions kz, until obtaining whole cine cardiac imaging data；
In navigation signal determination unit, the signal value of multiple points is obtained on each acquisition layer radially sampling line, and determination is described more
The average value of the signal value of a point, using the average value of the signal value of the multiple point as navigation signal, including：
In the first coordinate system, including the directions kz, the directions kx, the directions ky, with time change, multiple acquisition layers radially sample
Line obtains 5 points in one predetermined angle of plane internal rotation in the directions kx and the directions ky composition on each acquisition layer radially sampling line
Signal value navigation signal is determined as after signal value S1, S2, S3, S4, the S5 for getting 5 points：Actual value of the obtained navigation signal closer to navigation signal.
5. the motor message extraction element according to claim 4 from gated three-dimensional cardiac imaging, which is characterized in that described
Mixed signal determination unit, is specifically used for：
According to formula：
Acquire heartbeat and respiratory movement mixed signal S；Wherein, z is coordinate of the navigation signal on the directions kz；m(z)
For amplitude of the navigation signal on image area；N is that Descartes acquires the number of plies on the directions kz.
6. the motor message extraction element according to claim 5 from gated three-dimensional cardiac imaging, which is characterized in that described
Filter processing unit is specifically used for：
According to human normal palmic rate range 0.1HZ to 0.5HZ, human normal respiratory rate 0.6HZ to 3HZ is to the heartbeat
And respiratory movement mixed signal is filtered, and obtains the heartbeat signal and breath signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510855793.XA CN105266813B (en) | 2015-11-30 | 2015-11-30 | A kind of motor message extracting method and device from gated three-dimensional cardiac imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510855793.XA CN105266813B (en) | 2015-11-30 | 2015-11-30 | A kind of motor message extracting method and device from gated three-dimensional cardiac imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105266813A CN105266813A (en) | 2016-01-27 |
| CN105266813B true CN105266813B (en) | 2018-10-12 |
Family
ID=55137257
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510855793.XA Active CN105266813B (en) | 2015-11-30 | 2015-11-30 | A kind of motor message extracting method and device from gated three-dimensional cardiac imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105266813B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107133549B (en) * | 2016-02-29 | 2020-11-24 | 上海联影医疗科技有限公司 | ECT motion gating signal acquisition method and ECT image reconstruction method |
| CA3096034A1 (en) * | 2018-04-05 | 2019-10-10 | Siemens Medical Solutions Usa, Inc. | Motion signal derived from imaging data |
| CN109949206B (en) * | 2019-01-30 | 2023-09-26 | 上海联影医疗科技股份有限公司 | Method, device, equipment and storage medium for generating motion artifact image |
| CN110327066B (en) * | 2019-03-26 | 2022-11-25 | 上海联影医疗科技股份有限公司 | Cardiac motion signal acquisition method, cardiac motion signal acquisition apparatus, computer device, and storage medium |
| WO2024037019A1 (en) * | 2023-04-21 | 2024-02-22 | 武汉迈瑞医疗技术研究院有限公司 | Ultrasonic imaging method and system for fetal heart |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102077108A (en) * | 2008-04-28 | 2011-05-25 | 康奈尔大学 | Tool for accurate quantification in molecular mri |
| CN102078196A (en) * | 2009-11-27 | 2011-06-01 | 西门子公司 | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10706592B2 (en) * | 2014-01-06 | 2020-07-07 | Cedars-Sinai Medical Center | Systems and methods for myocardial perfusion MRI without the need for ECG gating and additional systems and methods for improved cardiac imaging |
-
2015
- 2015-11-30 CN CN201510855793.XA patent/CN105266813B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102077108A (en) * | 2008-04-28 | 2011-05-25 | 康奈尔大学 | Tool for accurate quantification in molecular mri |
| CN102078196A (en) * | 2009-11-27 | 2011-06-01 | 西门子公司 | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance |
Non-Patent Citations (2)
| Title |
|---|
| Respiratory and Cardiac Self-Gated Free-Breathing Cardiac CINE imaging with multi-echo 3D Hybrid radial SSFP Acquisition;jingliu等;《magnetic resonacne medical》;20100531;第63卷(第5期);第1-3页，第10-19页，摘要 * |
| 基于TrueFisp序列的心脏自门控信号提取方法;王志康等;《中国医疗器械杂志》;20061231;第30卷(第1期);1-3页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105266813A (en) | 2016-01-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105266813B (en) | A kind of motor message extracting method and device from gated three-dimensional cardiac imaging | |
| US10058257B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US8417007B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US8332010B2 (en) | Method for non-contrast enhanced magnetic resonance angiography | |
| US9113810B2 (en) | System and method for ungated non-contrast enhanced magnetic resonance angiography | |
| JP6073570B2 (en) | Magnetic resonance imaging system and method | |
| US8086297B2 (en) | Dark blood delayed enhancement magnetic resonance viability imaging techniques for assessing subendocardial infarcts | |
| JP5854575B2 (en) | Magnetic resonance imaging system | |
| US20150157238A1 (en) | Method For Non-Contrast Enhanced Magnetic Resonance Angiography | |
| WO2014120953A1 (en) | Method for accurate and robust cardiac motion self-gating in magnetic resonance imaging | |
| US20100222666A1 (en) | Method and apparatus for breath-held mr data acquisition using interleaved acquisition | |
| KR101809213B1 (en) | Method and apparatus for multi-slice imaging of t2-relaxation time | |
| Barkauskas et al. | Quantification of left ventricular functional parameter values using 3D spiral bSSFP and through-time Non-Cartesian GRAPPA | |
| US8154287B2 (en) | System and method for ghost magnetic resonance imaging | |
| Syed et al. | Magnetic resonance imaging of congenital heart disease | |
| JP2004329614A (en) | Magnetic resonance imaging apparatus | |
| CN109800800A (en) | A kind of method and apparatus that magnetic resonance imaging obtains deep learning training set | |
| CN206209898U (en) | Three-dimensional cardiac image re-construction system | |
| US10162035B2 (en) | System and method for controlling motion effects in magnetic resonance imaging (MRI) images | |
| JP4443918B2 (en) | Magnetic resonance imaging apparatus and signal processing method of magnetic resonance imaging apparatus | |
| US20140303482A1 (en) | Magnetic resonance imaging method for imaging components with short transverse relaxation times (t2) in a human or an animal heart | |
| JP5738120B2 (en) | Non-contrast angiogram reconstruction method and magnetic resonance imaging apparatus | |
| JP4738540B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| Jin et al. | A study of in vitro and in vivo MR of free-breathing whole-heart 3D coronary angiography using parallel imaging | |
| Fogel | General Principles of Cardiac Magnetic Resonance Imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |