CN105919624A - Encoding and imaging methods and devices of high-pulse repeated scanning frequency signals - Google Patents
Encoding and imaging methods and devices of high-pulse repeated scanning frequency signals Download PDFInfo
- Publication number
- CN105919624A CN105919624A CN201610351777.1A CN201610351777A CN105919624A CN 105919624 A CN105919624 A CN 105919624A CN 201610351777 A CN201610351777 A CN 201610351777A CN 105919624 A CN105919624 A CN 105919624A
- Authority
- CN
- China
- Prior art keywords
- signal
- scanning
- frequency
- scanning signal
- anglec
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention discloses encoding and imaging methods and devices of high-pulse repeated scanning frequency signals. The encoding method comprises steps as follows: a scanning signal set is determined according to scanning depth and comprises at least two signals; a rotation angle factor of fractional order Fourier transform in the best signal convergence area in the scanning signal set is determined; the starting frequency of each scanning signal is calculated according to the rotation angle factor; the scanning signals corresponding to the starting frequency are encoded according to the starting frequency of the scanning signals. The corresponding relation between a received signal and a transmitted signal can be determined, limitation of the scanning distance to the pulse repeated scanning frequency can be eliminated, and a high-speed target is scanned.
Description
Technical field
The present invention relates to medical ultrasonic technical field, particularly relate to a kind of high impulse and repeat scanning frequency signal
Coding, formation method and device.
Background technology
Ultrasound wave, when human body is propagated, runs into moving histoorgan or blood flow cell, in reflection and scattering
Time can cause the frequency shift (FS) of echo, i.e. Doppler frequency shift.By extracting and analyzing Doppler shift information,
Can obtain the movable information of histoorgan, i.e. D type (Doppler) is ultrasonic.Ultrasonic Doppler imaging is used for
Obtain the information in terms of the function information of people's intracorporeal organ, anatomic information and hematodinamics, as blood flow to,
Speed and cardiac motion state etc..
The detectable velocity interval of doppler imaging, depends on the pulse of target in sampling gate is repeated scanning frequency
Rate (Pulse-Recurrence-Frequency, PRF), for high-speed target, needs higher pulse weight
Multiple scanning frequency.But owing to ultrasound wave spread speed in tissue is certain, so in scanning distance
During deeper destination organization, due to propagation distance farther out, so the time of sound wave round trip is longer, this is just
Constrain pulse and repeat scanning frequency, so there is physics lance in principle between the sampling gate degree of depth and high PRF
Shield.
At present, when the PRF needed is beyond the restriction of the current scanning degree of depth, the used method of current techniques is
Reducing the PRF integral multiple needed, this has resulted in several sampling gate, and last Doppler results will be this
The mixed information of the target travel in several sampling gates, this cannot truly reflect the speed of required sampled targets
Degree information.
Summary of the invention
In view of this, the embodiment of the present invention provides a kind of high impulse to repeat the coding of scanning frequency signal, imaging
Method and device, to solve the technical problem that target pulse repetition scanning frequency is limited by propagation distance.
First aspect, embodiments provides a kind of high impulse and repeats the coded method of scanning frequency signal
Including:
Determine that scanning signal set, described scanning signal set include at least two signal according to the scanning degree of depth；
Determine the anglec of rotation of the Fourier Transform of Fractional Order of the optimal zone of convergency of the signal in scanning signal set
The factor；
The initial frequency of described each scanning signal is calculated according to the described anglec of rotation factor；
Initial frequency pair according to the described scanning signal scanning signal corresponding with described initial frequency is compiled
Code.
Second aspect, the embodiment of the present invention additionally provides method for ultrasonic imaging in a kind of Doppler Imaging systems,
Including:
Receive ultrasound wave coded echo signal, and described ultrasound wave coded echo signal is converted to digital signal；
Described digital signal is carried out Beam synthesis, generates radiofrequency signal；
To described radiofrequency signal anti-aliasing and decoding process, and by inphase-quadrature modem, generate inphase quadrature
Signal；
Carry out imaging processing according to the signal after inphase-quadrature modem, generate ultrasonography.
The third aspect, the embodiment of the present invention additionally provides a kind of high impulse and repeats the coding dress of scanning frequency signal
Put, including:
Set determines module, for determining scanning signal set, described scanning signal set according to the scanning degree of depth
Including at least two signal；
Angular factors determines module, for determining the mark of the optimal zone of convergency of the signal in scanning signal set
The anglec of rotation factor of rank Fourier transformation；
Initial frequency computing module, for calculating described each scanning signal according to the described anglec of rotation factor
Initial frequency；
Coding module, for sweep corresponding with described initial frequency of initial frequency pair according to described scanning signal
Look into signal to encode.
Fourth aspect, the embodiment of the present invention additionally provides a kind of ultrasonic imaging apparatus, including:
Digital signal conversion module, is used for receiving ultrasound wave coded echo signal, and is encoded by described ultrasound wave
Echo-signal is converted to digital signal；
Radiofrequency signal generation module, for described digital signal carries out Beam synthesis, generates radiofrequency signal；
Inphase quadrature signal generation module, is used for described radiofrequency signal anti-aliasing and decoding process, and passes through
Inphase-quadrature modem, generates inphase quadrature signal；
Image generation module, for carrying out imaging processing according to the signal after inphase-quadrature modem, generates ultrasonic
Ripple image.
The high impulse that the embodiment of the present invention provides repeats the coding of scanning frequency signal, formation method and device,
By different transmitting signals is encoded, and the signal received is decoded accordingly, can be true
Surely the signal received and the corresponding relation launched between signal.Can solve the problem that scanning distance repeats for pulse
The restriction of scanning frequency, it is achieved that the scanning to high-speed target.
Accompanying drawing explanation
The detailed description that non-limiting example is made made with reference to the following drawings by reading, the present invention
Other features, objects and advantages will become more apparent upon:
Fig. 1 is the flow process that the high impulse that the embodiment of the present invention one provides repeats the coded method of scanning frequency signal
Schematic diagram；
Fig. 2 is the flow process that the high impulse that the embodiment of the present invention two provides repeats the coded method of scanning frequency signal
Schematic diagram；
Fig. 3 is the schematic flow sheet of the method for ultrasonic imaging that the embodiment of the present invention three provides；
Fig. 4 is the schematic flow sheet of the method for ultrasonic imaging that the embodiment of the present invention four provides；
Fig. 5 is the structure that the high impulse that the embodiment of the present invention five provides repeats the code device of scanning frequency signal
Schematic diagram；
Fig. 6 is the structural representation of the ultrasonic imaging apparatus that the embodiment of the present invention six provides.
Detailed description of the invention
The present invention is described in further detail with embodiment below in conjunction with the accompanying drawings.It is understood that this
Specific embodiment described by place is used only for explaining the present invention, rather than limitation of the invention.The most also need
Be noted that for the ease of describing, accompanying drawing illustrate only part related to the present invention and not all in
Hold.
Embodiment one
Fig. 1 repeats the flow process of the coded method of scanning frequency signal for the high impulse that the embodiment of the present invention one provides
Schematic diagram, the situation that the scanning signal that the present embodiment is applicable to launch ultrasonic system encodes, should
Method can be performed by the code device of pulse repetition scanning frequency signal, and this device can be by software/hardware side
Formula realizes, and can be integrated in corresponding ultrasonic wave detecting system.
Seeing Fig. 1, described high impulse repeats the coded method of scanning frequency signal, including:
According to the scanning degree of depth, S110, determines that scanning signal set, described scanning signal set include at least two
Signal.
Modern ultrasound detection technique uses pulse reflection method, i.e. to one the short pulse punching of tested object emission
Ultrasound wave, then receive it.And so forth.In unit interval, the number of times driving probe to launch ultrasound wave is exactly
Pulse recurrence frequency.
For high-speed moving object, such as when histoorgan or blood flow cell carry out scanning, due to the speed of blood flow
Degree relative motion is very fast, needs to improve the scanning frequency of response, i.e. launches multiple scanning arteries and veins within the unit interval
Rush signal.Exemplary, these scanning pulse signals composition scanning signal set, each scanning signal set
Include several scanning signals.Scanning signal set includes at least two scanning signal.
Scanning signal set can be determined by the scanning degree of depth, is D to the current scanning degree of depth, then the pulse of its maximum
Repeating scanning frequency is Prfmax=C/ (2*D), and wherein C is SVEL.In the present embodiment, permissible
The all scanning pulse signal seats sent in maximum pulse is repeated the cycle that scanning frequency Prfmax is corresponding
Scanning pulse collection.
S120, determines the rotation of the Fourier Transform of Fractional Order of the optimal zone of convergency of the signal in scanning signal set
The gyration factor.
Fourier Transform of Fractional Order (Fractional Fourier Transform, FrFT) is Fourier transformation
A kind of generalized form, signal contains the signal letter at time-domain and frequency-domain on fractional number order Fourier simultaneously
Breath, makes FrFT have the character that many conventional Fourier transform do not possess, can be used to process non-stationary signal.
The basic meaning of Fourier space is, arbitrary periodic function (signal) can be decomposed into infinite multiple frequency and be
The sine wave of its fundamental frequency integer (including zero) times and cosine wave sum.Fourier transformation is then by its week
Phase expands to infinite formation.Can be understood as signal coordinate axes in time-frequency plane to rotate counterclockwise around initial point
Method for expressing after arbitrarily angled.
The multiple pulse signals launched for ultrasonic system, can regard as periodic signal, for this week
Phase signal, has the property that signal energy is gradually assembled along with the fractional order territory of signal is gradually increased,
It is impulse function to frequency domain representation.Based on this characteristic, can determine unique according to Fourier Transform of Fractional Order
The factor of influence anglec of rotation factor.
Concrete, the anglec of rotation factor can be calculated by equation below:
Aifa=1-2/pi*atan (Tx_fs/BandWidthOfFM)；Wherein, aifa be the anglec of rotation because of
Son；Tx_fs represents the sample frequency that system is launched；The modulation band of BandWidthOfFM presentation code signal
Wide.
S130, calculates the initial frequency of described each scanning signal according to the described anglec of rotation factor.
Concrete, may include steps of:
The relative bandwidth of each scanning signal in described scanning signal set is determined according to the described anglec of rotation factor
Frequency；
Mid frequency according to described relative bandwidth frequency and the pulse of ultrasonic system transmitting calculates described each
The initial frequency of scanning signal.
Several line chart pictures crossing at 0 can be shown as in time-domain and frequency-domain plane for scanning signal set,
The line segment of the most each mistake 0 represents the relative frequency of each scanning signal in scanning signal set.
Concrete, can be calculated by equation below:
Bandwidth_Sig=4*nZ*abs (sin (aifa*pi/2))/TimeOfFM+D*f_Atten；
Wherein, the modulation duration of TimeOfFM presentation code signal, D represents the maximum scanning degree of depth, f_Atten
Represent the frequency attenuation quotient with change in depth.NZ is constant coefficient, typically can take 1, for every n
Scanning signal, can be calculated the relative frequency of this scanning signal according to above-mentioned formula.
For the probe of any one ultrasonic system, its driving frequency can set according to different diagnostic purposes
It is set to different frequencies.And driving frequency is when the mid frequency set, and work efficiency is the highest.Scanning is believed
Number set in each scanning signal frequency should probe mid frequency certain limit in.And can basis
In scanning signal set, each relative frequency of scanning signal and the mid frequency of probe are calculated scanning signal
The frequency of each scanning signal in set.Concrete, can be calculated by equation below:
Fstart=f0 BandWidthOfFM/2+n*bandwidth_Sig.
Wherein, f0 represents the mid frequency of impulse ejection.
S140, enters according to the scanning signal that the initial frequency pair of described scanning signal is corresponding with described initial frequency
Row coding.
Owing to scanning signal does not has the change of phase place, correspondence scanning can be generated according to the frequency of scanning signal
The waveshape signal of signal, encodes scanning signal according to described waveshape signal.Concrete, can pass through
Equation below obtains the coding signal that scanning signal is corresponding:
Base_excitation (i)=cos (2*pi* (fstart*t+BandWidthOfFM/2/TimeOfFM
*t.^2))；Wherein t=(0,1/Tx_fs, 2/Tx_fs ... TimeOfFM).
The high impulse that the present embodiment provides repeats the coded method of scanning frequency signal, utilizes ultrasonic system institute
The scanning signal set launched sets the time-domain and frequency-domain plane of the anglec of rotation factor in Fourier Transform of Fractional Order
The feature of interior gathering, by calculating the frequency of each scanning signal, determines the corresponding of each scanning signal then
Waveform, and according to described waveform, scanning signal is encoded, so that ultrasonic system is receiving scanning
Reception parameter can be determined according to the waveform after scanning Signal coding during signal echo, and then realize obtaining at a high speed
Motion sample target velocity information accurately.
Embodiment two
Fig. 2 is the flow process that the high impulse that the embodiment of the present invention two provides repeats the coded method of scanning frequency signal
Schematic diagram.The present embodiment, based on above-described embodiment, increases following steps: determine scanning signal set number,
And determine each scanning signal scanning sequence number in corresponding scanning set；Accordingly, by described according to institute
State the anglec of rotation factor and calculate the initial frequency of described each scanning signal, be specifically optimized for: according to described rotation
The gyration factor and scanning signal scanning sequence number in corresponding scanning set calculate described each scanning signal
Initial frequency.
Accordingly, the high impulse that the present embodiment is provided repeats the coded method of scanning frequency signal, specifically wraps
Include:
S210, determines scanning signal set number, and determines that each scanning signal is in corresponding scanning set
Scanning sequence number.
Scanning signal set can be determined by the scanning degree of depth, is D to the current scanning degree of depth, then the pulse of its maximum
Repeating scanning frequency is Prfmax=C/ (2*D), and wherein C is SVEL.I.e. in ultrasonic system,
Within a transmitting cycle, at most can launch Prfmax ultrasonic signal, and high-speed motion is swept
Looking into target, it is Prf_desire that its required pulse repeats scanning frequency, if the pulse required repeats to sweep
Looking into frequency is that Prf_desire is more than Prfmax, then in 1 transmitting cycle, it is impossible to complete to send requirement
Pulse repeat scanning number of times, need within multiple cycles, complete required pulse and repeat scanning signal
Sending, each cycle can send a scanning signal set.Concrete, can be calculated by equation below
The scanning signal set number N of transmitting needed for obtaining:
CEIL (Prf_desire/Prfmax), CEIL (.) expression rounds up.
The signal launched for i & lt, it is the transmitting of the n-th sequence in some scanning signal set
Signal, concrete, scanning signal scanning sequence in corresponding scanning set can be determined by equation below
Number:
N=mod (i, N).
S220, determines the rotation of the Fourier Transform of Fractional Order of the optimal zone of convergency of the signal in scanning signal set
The gyration factor.
S230, the scanning sequence number in corresponding scanning set according to the described anglec of rotation factor and scanning signal
Calculate the initial frequency of described each scanning signal.
Concrete, may include steps of:
According to the described anglec of rotation factor and scanning signal, the scanning sequence number in corresponding scanning set determines institute
State the relative bandwidth frequency of each scanning signal in scanning signal set；Concrete, equation below can be passed through
It is calculated:
Bandwidth_Sig=4*nZ*abs (sin (aifa*pi/2))/TimeOfFM+D*f_Atten；
NZ typically takes 1, the modulation duration of TimeOfFM presentation code signal, and D represents the maximum scanning degree of depth,
F_Atten represents the frequency attenuation quotient with change in depth；
Exist according to described relative bandwidth frequency, the mid frequency of the pulse of ultrasonic system transmitting and scanning signal
The corresponding scanning sequence number in scanning set calculates the initial frequency of described each scanning signal.Concrete, can
To be calculated in the following way:
Fstart=f0 BandWidthOfFM/2+ (n (N-1)/2) * bandwidth_Sig；
F0 represents the mid frequency of impulse ejection, and N is the scanning signal set quantity needing to launch.
S240, enters according to the scanning signal that the initial frequency pair of described scanning signal is corresponding with described initial frequency
Row coding.
The present embodiment is by increasing following steps: determines scanning signal set number, and determines each scanning signal
Scanning sequence number in corresponding scanning set；Accordingly, calculate described according to the described anglec of rotation factor
The initial frequency of described each scanning signal, is specifically optimized for: believe according to the described anglec of rotation factor and scanning
Scanning sequence number number in corresponding scanning set calculates the initial frequency of described each scanning signal.Can root
Repeat scanning frequency according to maximum pulse and the pulse repetition scanning frequency required by sampled targets is carried out different
Set packet, and the different scanning signals in set are encoded.
Embodiment three
The schematic flow sheet of the method for ultrasonic imaging that Fig. 3 provides for the embodiment of the present invention three, the present embodiment can
Utilizing the situation of ultrasonic signal imaging be applicable to Doppler Imaging systems, the method can be become by ultrasound wave
As device performs, this device can be realized by software/hardware mode, and can be integrated in corresponding ultrasound examination
In system.
See Fig. 3, described method for ultrasonic imaging, including:
S310, receives ultrasound wave coded echo signal, and described ultrasound wave coded echo signal is converted to number
Word signal.
Ultrasound probe in ultrasonic system is used for launching ultrasound wave and receiving ultrasonic echo.And pass through ripple
The amplification of echo, filtering and Analog-digital Converter, obtain digitalized ultrasonic ripple echo；Described reception device can
To include the low noise amplifier circuit, Anti-aliasing Filter Circuits and the D/A conversion circuit that are linked in sequence.Institute
State low noise amplifier circuit to be connected with described ultrasound probe.
S320, carries out Beam synthesis by described digital signal, generates radiofrequency signal.
Beam synthesis refer to the multispot array arrange certain geometries respectively shakes the output of unit through postponing,
After weighted sum summation etc. out, output is made to have space directivity.Exemplary, delay stack can be used
Method, it is thus achieved that ultrasound RF data.
S330, to described radiofrequency signal anti-aliasing and decoding process, and by quadrature demodulation, is just generating homophase
Hand over signal.
Exemplary, anti-aliasing filtering can be used to remove device, on the premise of sample frequency is constant, by setting
The low pass filter of good frequency filters the frequency content being not required to occur.Owing to the signal of input includes I/Q two
Point, need by quadrature demodulation, to generate the orthogonal signalling of homophase.
S340, carries out imaging processing according to the signal after inphase-quadrature modem, generates ultrasonography.
The present embodiment is by carrying out anti-aliasing and decoding process to the ultrasonic signal received, by receive
Coding signal processing, generates inphase quadrature signal, it is possible to achieve obtain high-speed motion sampled targets the fastest
Degree information.
Embodiment four
Fig. 4 is the schematic flow sheet of the method for ultrasonic imaging that the embodiment of the present invention four provides.The present embodiment with
Based on above-described embodiment, by radiofrequency signal anti-aliasing, specifically it is optimized for: according to default length samples institute
Stating radiofrequency signal, described default length is according to the width of sampling gate and receives ultrasound wave coded echo signal
Sample frequency determines；Determine the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal；According to
The described radiofrequency signal that the sampling of described anglec of rotation factor pair obtains carries out Fourier Transform of Fractional Order；Set band
The cut-off frequency of bandpass filter, and be filtered processing to the radiofrequency signal after conversion；To filtered radio frequency
Signal carries out inverse Fourier Transform of Fractional Order.
See Fig. 4, described method for ultrasonic imaging, including:
S410, receives ultrasound wave coded echo signal, and described ultrasound wave coded echo signal is converted to number
Word signal.
S420, carries out Beam synthesis by described digital signal, generates radiofrequency signal.
S430, according to radiofrequency signal described in default length samples, described default length is according to sampling gate
The sample frequency of width and reception ultrasound wave coded echo signal determines.
Exemplary, can determine sampling length in the following way, sampling length L=w_sig*2/C*Rx_fs,
Wherein, w_sig is the width of sampling gate, and C is acoustic wave propagation velocity, and Rx_fs represents the sampling receiving signal
Frequency.
S440, determines the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal.
Exemplary, the anglec of rotation factor of Fourier Transform of Fractional Order can be determined in the following way, rotation
Gyration factor aifa_rx=1-2/pi*atan (Rx_fs/BandWidthOfFM), wherein
The modulation bandwidth of BandWidthOfFM presentation code signal；Rx_fs represents the sample frequency receiving signal.
S450, is carried out in fractional order Fu according to the described radiofrequency signal that the sampling of described anglec of rotation factor pair obtains
Leaf transformation.
Radiofrequency signal is the compound of multiple signal, there is serious cross term interference.Fourier Transform of Fractional Order
Video distribution cross term can be suppressed.Exemplary, in the following way radiofrequency signal can be carried out mark
Rank Fourier transformation, the signal Sig_frft=frft (sig, aifa_rx) after Fourier Transform of Fractional Order,
Frft (.) represents Fourier Transform of Fractional Order function.
S460, sets the cut-off frequency of band filter, and is filtered processing to the radiofrequency signal after conversion.
For i & lt scanning signal, after Fourier Transform of Fractional Order, n-th Z excessively 0 band is a width of:
Bandwidth_Sig_rx=4*nZ*abs (sin (aifa_rx*pi/2))/TimeOfFM+
D*f_Atten；Wherein, the modulation duration of TimeOfFM presentation code signal, D represents the maximum scanning degree of depth,
F_Atten represents the frequency attenuation quotient with change in depth.Z is constant coefficient, typically can take 1, for
Each scanning signal, can be calculated the relative frequency of this scanning signal according to above-mentioned formula.So then filtering
The mid frequency of ripple device is:
F0_bandfilter_i=f0+ (n (N-1)/2) * bandwidth_Sig_rx
Then the upper lower limiting frequency of this band filter is respectively as follows:
Bandfilter_lowFreq=f0_bandfilter_i-bandwidth_Sig_rx/2；
Bandfilter_highFreq=f0_bandfilter_i+bandwidth_Sig_rx/2.
S470, carries out inverse Fourier Transform of Fractional Order to filtered radiofrequency signal.
Filtered radiofrequency signal is carried out inverse Fourier Transform of Fractional Order, makes described filtered radiofrequency signal
Recover original frequency.
S480, to the signal decoding process after anti-aliasing, and by quadrature demodulation, generates inphase quadrature signal.
S490, carries out imaging processing according to the signal after inphase-quadrature modem, generates ultrasonography.
The present embodiment, by by radiofrequency signal anti-aliasing, is specifically optimized for: described in default length samples
Radiofrequency signal, described default length is according to the width of sampling gate and adopting of reception ultrasound wave coded echo signal
Sample frequency determines；Determine the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal；According to institute
State the described radiofrequency signal that obtains of anglec of rotation factor pair sampling and carry out Fourier Transform of Fractional Order；Set band logical
The cut-off frequency of wave filter, and be filtered processing to the radiofrequency signal after conversion；Filtered radio frequency is believed
Number carry out inverse Fourier Transform of Fractional Order.In less frequency range, radiofrequency signal can be filtered,
Remove noise wave removing, improve the precision of ultrasonic scanning.
Embodiment five
Fig. 5 is the structure that the high impulse that the embodiment of the present invention five provides repeats the code device of scanning frequency signal
Schematic diagram, as it is shown in figure 5, described device includes:
Set determines module 510, for determining scanning signal set, described scanning signal collection according to the scanning degree of depth
Conjunction includes at least two signal；
Angular factors determines module 520, for determining dividing of the optimal zone of convergency of the signal in scanning signal set
The anglec of rotation factor of number rank Fourier transformation；
Initial frequency computing module 530, for calculating described each scanning signal according to the described anglec of rotation factor
Initial frequency；
Coding module 540 is corresponding with described initial frequency for the initial frequency pair according to described scanning signal
Scanning signal encodes.
The high impulse that the present embodiment provides repeats the code device of scanning frequency signal, by different transmittings
Signal encodes, and decodes the signal received accordingly, it may be determined that the signal received with
Launch the corresponding relation between signal.Scanning distance can be released and pulse is repeated to the restriction of scanning frequency,
Achieve the scanning to high-speed target.
On the basis of the various embodiments described above, described initial frequency computing module, including:
Relative bandwidth determines unit, for determining in described scanning signal set according to the described anglec of rotation factor
The relative bandwidth of each scanning signal；
Frequency computing unit, in the pulse launched according to described relative bandwidth frequency and ultrasonic system
Frequency of heart calculates the initial frequency of described each scanning signal.
On the basis of the various embodiments described above, described device also includes:
Scanning sequence number determines module, is used for determining scanning signal set number, and determines that each scanning signal is right
Scanning sequence number in the scanning set answered；
Accordingly, described relative bandwidth determines unit, is used for:
According to the described anglec of rotation factor and scanning signal, the scanning sequence number in corresponding scanning set determines institute
State the relative bandwidth of each scanning signal in scanning signal set.
The high impulse that the embodiment of the present invention is provided repeats the code device of scanning frequency signal and can be used for performing basis
The high impulse that invention any embodiment provides repeats the coded method of scanning frequency signal, possesses corresponding function
Module, it is achieved identical beneficial effect.
Embodiment six
Fig. 6 is the structural representation of the ultrasonic imaging apparatus that the embodiment of the present invention six provides, as shown in Figure 6,
Described device includes:
Digital signal conversion module 610, is used for receiving ultrasound wave coded echo signal, and is compiled by described ultrasound wave
Code echo-signal is converted to digital signal；
Radiofrequency signal generation module 620, for described digital signal carries out Beam synthesis, generates radiofrequency signal；
Inphase quadrature signal generation module 630, is used for described radiofrequency signal anti-aliasing and decoding process, and leads to
Cross inphase-quadrature modem, generate inphase quadrature signal；
Image generation module 640, for carrying out imaging processing according to the signal after inphase-quadrature modem, generates super
Sonogram picture.
On the basis of the various embodiments described above, described inphase quadrature signal generation module, it is used for:
According to radiofrequency signal described in default length samples, described default length according to the width of sampling gate and
The sample frequency receiving ultrasound wave coded echo signal determines；
Determine the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal；
Fractional order Fourier change is carried out according to the described radiofrequency signal that the sampling of described anglec of rotation factor pair obtains
Change；
Set the cut-off frequency of band filter, and be filtered processing to the radiofrequency signal after conversion；
Filtered radiofrequency signal is carried out inverse Fourier Transform of Fractional Order.
The ultrasonic imaging apparatus that the present embodiment provides, by carrying out anti-aliasing to the ultrasonic signal received
And decoding process, the coding signal processing that will receive, generate inphase quadrature signal, it is possible to achieve obtain height
Speed motion sample target velocity information accurately.
The ultrasonic imaging apparatus that the embodiment of the present invention is provided can be used for performing any embodiment of the present invention to be provided
Method for ultrasonic imaging, possess corresponding functional module, it is achieved identical beneficial effect.
Obviously, it will be understood by those skilled in the art that each module or each operation of the above-mentioned present invention can be led to
Cross terminal unit as above to implement.Alternatively, the embodiment of the present invention can be able to perform with computer installation
Program realize, perform such that it is able to be stored in storing in device by processor, described journey
Sequence can be stored in a kind of computer-readable recording medium, and storage medium mentioned above can be read-only depositing
Reservoir, disk or CD etc.；Or they to be fabricated to respectively each integrated circuit modules, or by them
In multiple modules or operation be fabricated to single integrated circuit module and realize.So, the present invention is not restricted to
The combination of any specific hardware and software.
Note, above are only presently preferred embodiments of the present invention and institute's application technology principle.Those skilled in the art
It will be appreciated that the invention is not restricted to specific embodiment described here, can enter for a person skilled in the art
Row various obvious changes, readjust and substitute without departing from protection scope of the present invention.Therefore, though
So by above example, the present invention is described in further detail, but the present invention be not limited only to
Upper embodiment, without departing from the inventive concept, it is also possible to include other Equivalent embodiments more,
And the scope of the present invention is determined by scope of the appended claims.
Claims (10)
1. the coded method of a high impulse repetition scanning frequency signal, it is characterised in that including:
Determine that scanning signal set, described scanning signal set include at least two signal according to the scanning degree of depth；
Determine the anglec of rotation of the Fourier Transform of Fractional Order of the optimal zone of convergency of the signal in scanning signal set
The factor；
The initial frequency of described each scanning signal is calculated according to the described anglec of rotation factor；
Initial frequency pair according to the described scanning signal scanning signal corresponding with described initial frequency is compiled
Code.
Method the most according to claim 1, it is characterised in that described according to the described anglec of rotation factor
Calculate the initial frequency of described each scanning signal, including:
The opposite band of each scanning signal in described scanning signal set is determined according to the described anglec of rotation factor
Wide；
Mid frequency according to described relative bandwidth frequency and the pulse of ultrasonic system transmitting calculates described each
The initial frequency of scanning signal.
Method the most according to claim 2, it is characterised in that described method also includes:
Determine scanning signal set number, and determine each scanning signal scanning sequence in corresponding scanning set
Number；
Accordingly, described each scanning letter is determined in described scanning signal set according to the described anglec of rotation factor
Number relative bandwidth, including:
According to the described anglec of rotation factor and scanning signal, the scanning sequence number in corresponding scanning set determines institute
State the relative bandwidth of each scanning signal in scanning signal set.
4. a method for ultrasonic imaging, it is characterised in that including:
Receive ultrasound wave coded echo signal, and described ultrasound wave coded echo signal is converted to digital signal；
Described digital signal is carried out Beam synthesis, generates radiofrequency signal；
To described radiofrequency signal anti-aliasing and decoding process, and by inphase-quadrature modem, generate inphase quadrature
Signal；
Carry out imaging processing according to the signal after inphase-quadrature modem, generate ultrasonography.
Method the most according to claim 4, it is characterised in that described to described radiofrequency signal anti-aliasing,
Including:
According to radiofrequency signal described in default length samples, described default length according to the width of sampling gate and
The sample frequency receiving ultrasound wave coded echo signal determines；
Determine the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal；
Fractional order Fourier change is carried out according to the described radiofrequency signal that the sampling of described anglec of rotation factor pair obtains
Change；
Set the cut-off frequency of band filter, and be filtered processing to the radiofrequency signal after conversion；
Filtered radiofrequency signal is carried out inverse Fourier Transform of Fractional Order.
6. the code device of a high impulse repetition scanning frequency signal, it is characterised in that including:
Set determines module, for determining scanning signal set, described scanning signal set according to the scanning degree of depth
Including at least two signal；
Angular factors determines module, for determining the mark of the optimal zone of convergency of the signal in scanning signal set
The anglec of rotation factor of rank Fourier transformation；
Initial frequency computing module, for calculating described each scanning signal according to the described anglec of rotation factor
Initial frequency；
Coding module, for sweep corresponding with described initial frequency of initial frequency pair according to described scanning signal
Look into signal to encode.
Device the most according to claim 1, it is characterised in that described initial frequency computing module, bag
Include:
Relative bandwidth determines unit, for determining in described scanning signal set according to the described anglec of rotation factor
The relative bandwidth of each scanning signal；
Frequency computing unit, in the pulse launched according to described relative bandwidth frequency and ultrasonic system
Frequency of heart calculates the initial frequency of described each scanning signal.
Device the most according to claim 7, it is characterised in that described device also includes:
Scanning sequence number determines module, is used for determining scanning signal set number, and determines that each scanning signal is right
Scanning sequence number in the scanning set answered；
Accordingly, described relative bandwidth determines unit, is used for:
According to the described anglec of rotation factor and scanning signal, the scanning sequence number in corresponding scanning set determines institute
State the relative bandwidth of each scanning signal in scanning signal set.
9. a ultrasonic imaging apparatus, it is characterised in that including:
Digital signal conversion module, is used for receiving ultrasound wave coded echo signal, and is encoded by described ultrasound wave
Echo-signal is converted to digital signal；
Radiofrequency signal generation module, for described digital signal carries out Beam synthesis, generates radiofrequency signal；
Inphase quadrature signal generation module, is used for described radiofrequency signal anti-aliasing and decoding process, and passes through
Inphase-quadrature modem, generates inphase quadrature signal；
Image generation module, for carrying out imaging processing according to the signal after inphase-quadrature modem, generates ultrasonic
Ripple image.
Device the most according to claim 9, it is characterised in that described inphase quadrature signal generates mould
Block, is used for:
According to radiofrequency signal described in default length samples, described default length according to the width of sampling gate and
The sample frequency receiving ultrasound wave coded echo signal determines；
Determine the anglec of rotation factor of the Fourier Transform of Fractional Order of described radiofrequency signal；
Fractional order Fourier change is carried out according to the described radiofrequency signal that the sampling of described anglec of rotation factor pair obtains
Change；
Set the cut-off frequency of band filter, and be filtered processing to the radiofrequency signal after conversion；
Filtered radiofrequency signal is carried out inverse Fourier Transform of Fractional Order.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610351777.1A CN105919624B (en) | 2016-05-25 | 2016-05-25 | A kind of high impulse repeats the coding and device of scanning frequency signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610351777.1A CN105919624B (en) | 2016-05-25 | 2016-05-25 | A kind of high impulse repeats the coding and device of scanning frequency signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105919624A true CN105919624A (en) | 2016-09-07 |
| CN105919624B CN105919624B (en) | 2019-09-10 |
Family
ID=56842269
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610351777.1A Active CN105919624B (en) | 2016-05-25 | 2016-05-25 | A kind of high impulse repeats the coding and device of scanning frequency signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105919624B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109199448A (en) * | 2018-09-04 | 2019-01-15 | 深圳市理邦精密仪器股份有限公司 | Pulsed-Wave Doppler imaging method and device based on HPRF |
| CN110559014A (en) * | 2019-08-28 | 2019-12-13 | 华南理工大学 | fractional order Fourier transform echo imaging method and system based on probe compensation |
| CN111093515A (en) * | 2018-03-30 | 2020-05-01 | 深圳迈瑞生物医疗电子股份有限公司 | Blood flow imaging processing method and ultrasonic imaging equipment |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6350240B1 (en) * | 2000-02-01 | 2002-02-26 | Medison Co., Ltd | Ultrasound imaging method and apparatus based on pulse compression technique using modified golay codes |
| US20090209858A1 (en) * | 2008-02-18 | 2009-08-20 | The Board Of Trustees Of The University Of Illinois | System and method for ultrasonic image processing |
| CN101785684A (en) * | 2010-01-27 | 2010-07-28 | 天津大学 | Ultrasonic endoscopic synthetic aperture imaging system and synthetic method of larger aperture |
| CN102133089A (en) * | 2011-03-28 | 2011-07-27 | 李鹏 | Blood flow speed estimation method and method for generating color blood flow chart with same |
| CN102724155A (en) * | 2012-05-17 | 2012-10-10 | 哈尔滨工程大学 | High frequency domain energy concentration synchronization method based on fractional Fourier transform |
| CN102999473A (en) * | 2012-10-18 | 2013-03-27 | 中国人民解放军电子工程学院 | Detection and parameter estimation method for linear frequency modulation signals |
| CN103070701A (en) * | 2012-12-28 | 2013-05-01 | 深圳先进技术研究院 | Ultrasonic modulation imaging system and method |
| CN104980390A (en) * | 2015-05-31 | 2015-10-14 | 厦门大学 | Demodulation method of double-rate combined chirp signal |
-
2016
- 2016-05-25 CN CN201610351777.1A patent/CN105919624B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6350240B1 (en) * | 2000-02-01 | 2002-02-26 | Medison Co., Ltd | Ultrasound imaging method and apparatus based on pulse compression technique using modified golay codes |
| US20090209858A1 (en) * | 2008-02-18 | 2009-08-20 | The Board Of Trustees Of The University Of Illinois | System and method for ultrasonic image processing |
| CN101785684A (en) * | 2010-01-27 | 2010-07-28 | 天津大学 | Ultrasonic endoscopic synthetic aperture imaging system and synthetic method of larger aperture |
| CN102133089A (en) * | 2011-03-28 | 2011-07-27 | 李鹏 | Blood flow speed estimation method and method for generating color blood flow chart with same |
| CN102724155A (en) * | 2012-05-17 | 2012-10-10 | 哈尔滨工程大学 | High frequency domain energy concentration synchronization method based on fractional Fourier transform |
| CN102999473A (en) * | 2012-10-18 | 2013-03-27 | 中国人民解放军电子工程学院 | Detection and parameter estimation method for linear frequency modulation signals |
| CN103070701A (en) * | 2012-12-28 | 2013-05-01 | 深圳先进技术研究院 | Ultrasonic modulation imaging system and method |
| CN104980390A (en) * | 2015-05-31 | 2015-10-14 | 厦门大学 | Demodulation method of double-rate combined chirp signal |
Non-Patent Citations (10)
| Title |
|---|
| 刘宝华: "多分量LFM信号在分数阶Fourier域的参数估计与分离", 《中国科技论文》 * |
| 刘贵栋，沈毅: "基于分数阶余弦变换的医学超声编码激励技术", 《中国医学物理学杂志》 * |
| 周奇: "高频超声编码激励系统研究与设计", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
| 周浩等: "医学超声编码激励技术研究进展", 《浙江大学学报（工学版）》 * |
| 张月云: "基于时频分析的Chirp 信号的参数估计方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 彭旗宇，高上凯: "医学超声成像中伪Chirp编码激励技术的仿真研究", 《中国超声诊断杂志》 * |
| 李家强等: "基于高斯短时分数阶傅里叶变换的多分量LFM信号检测与参数估计", 《电子与信息学报》 * |
| 李进: "分数阶 Fourier 变换在超声信号处理中的应用", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 李靖等: "基于分数阶傅里叶变换的 chirp 信号时频分析", 《系统工程与电子技术》 * |
| 王靓等: "一种基于chirp编码的高帧速率超声成像方法", 《现代生物医学进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111093515A (en) * | 2018-03-30 | 2020-05-01 | 深圳迈瑞生物医疗电子股份有限公司 | Blood flow imaging processing method and ultrasonic imaging equipment |
| US11602324B2 (en) | 2018-03-30 | 2023-03-14 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Flow imaging processing method and ultrasound imaging device |
| CN111093515B (en) * | 2018-03-30 | 2023-10-20 | 深圳迈瑞生物医疗电子股份有限公司 | Blood flow imaging processing method and ultrasonic imaging equipment |
| CN109199448A (en) * | 2018-09-04 | 2019-01-15 | 深圳市理邦精密仪器股份有限公司 | Pulsed-Wave Doppler imaging method and device based on HPRF |
| CN109199448B (en) * | 2018-09-04 | 2021-08-27 | 深圳市理邦精密仪器股份有限公司 | Pulse wave Doppler imaging method and device based on HPRF |
| CN110559014A (en) * | 2019-08-28 | 2019-12-13 | 华南理工大学 | fractional order Fourier transform echo imaging method and system based on probe compensation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105919624B (en) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN204495996U (en) | broadband sonar receiver | |
| O'Donnell | Coded excitation system for improving the penetration of real-time phased-array imaging systems | |
| CN101595396A (en) | The apparatus and method of high-resolution ultrasound displacement measurement | |
| JPH11309147A (en) | System and method for imaging flow of ultrasonic scatterer | |
| CN105997142A (en) | Ultrasonic system transmitting signal compositing and imaging method and device | |
| CN105919624A (en) | Encoding and imaging methods and devices of high-pulse repeated scanning frequency signals | |
| CN105997147B (en) | A kind of ultrasonic pulse Doppler imaging method and device | |
| CN107817297A (en) | A kind of ultrasonic imaging processing method and processing system based on ultrasonic echo RF data | |
| CN105167802A (en) | Doppler imaging method and device | |
| US10107901B2 (en) | Acoustic signal phase rotation processor | |
| US6758815B2 (en) | Apparatus and method for indicating mechanical stiffness properties of body tissue | |
| CN108186050B (en) | Doppler blood flow velocity imaging method and system based on ultrasonic channel data | |
| CN103654861A (en) | Ultrasonic umbilical cord blood flow measuring system and ultrasonic umbilical cord blood flow measuring method both based on pseudorandom codes | |
| CN1886094B (en) | Doppler velocity detection device and ultrasonographic device using the same | |
| Win et al. | Identification and removal of reverberation in ultrasound imaging | |
| JP2010063829A (en) | Ultrasonic imaging apparatus | |
| JP5801956B2 (en) | Propagation velocity measuring device, propagation velocity measuring program, and propagation velocity measuring method | |
| CN101149433A (en) | Supersonic movement detection method coding using reversing Golay dode | |
| CN107870034B (en) | Underwater acoustic velocity measurement method based on phase difference | |
| JPH11309146A (en) | System and method for imaging flow of ultrasonic scatterer | |
| CN108652666A (en) | A kind of generation method and device of doppler blood flow imaging | |
| WO2018172882A1 (en) | Doppler ultrasound using sparse emission pattern | |
| Jin et al. | A new scheme of coded ultrasound using Golay codes | |
| Wagner et al. | Compressed beamforming with applications to ultrasound imaging | |
| CN104297756B (en) | Underwater rotating target scale sonar imaging method and system |
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 | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology Co.,Ltd. Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd. |
|
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd. Address before: 215000 5th floor, building a, 4th floor, building C, 27 Xinfa Road, Suzhou Industrial Park, Suzhou City, Jiangsu Province Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd. |