CN105919624B - A kind of high impulse repeats the coding and device of scanning frequency signal - Google Patents
A kind of high impulse repeats the coding and device of scanning frequency signal Download PDFInfo
- Publication number
- CN105919624B CN105919624B CN201610351777.1A CN201610351777A CN105919624B CN 105919624 B CN105919624 B CN 105919624B CN 201610351777 A CN201610351777 A CN 201610351777A CN 105919624 B CN105919624 B CN 105919624B
- Authority
- CN
- China
- Prior art keywords
- scanning
- signal
- frequency
- scanning signal
- rotation angular
- 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
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 coding, imaging method and devices that a kind of high impulse repeats scanning frequency signal, wherein the coding method includes: to determine that scanning signal set, the scanning signal set include at least two signals according to scanning depth；Determine the rotation angular factors of the Fourier Transform of Fractional Order of the best zone of convergency of signal in scanning signal set；The initial frequency of each scanning signal is calculated according to the rotation angular factors；It is encoded according to the initial frequency pair of scanning signal scanning signal corresponding with the initial frequency.It can determine the signal received and emit the corresponding relationship between signal.It is able to solve the limitation that scanning distance repeats pulse scanning frequency, realizes the scanning to high-speed target.
Description
Technical field
The present invention relates to medical ultrasonic technical field more particularly to a kind of high impulse repeat scanning frequency signal coding,
Imaging method and device.
Background technique
Ultrasonic wave encounters moving histoorgan or blood flow cell, can make in reflection and scattering when human body is propagated
At the frequency shift (FS) of echo, i.e. Doppler frequency shift.By extracting and analyzing Doppler shift information, available histoorgan
Motion information, i.e. D type (Doppler) ultrasound.Ultrasonic Doppler imaging is used to obtain the functional information of people's intracorporeal organ, dissection letter
Information in terms of breath and Hemodynamics, such as blood flow direction, speed and cardiac motion state.
The detectable velocity interval of doppler imaging repeats scanning frequency depending on the pulse to target in sampling gate
(Pulse-Recurrence-Frequency, PRF) needs higher pulse to repeat scanning frequency high-speed target.But by
In ultrasonic wave, spread speed is certain in tissue, so in the scanning deeper destination organization of distance, due to propagating
Farther out, so the time of sound wave round trip is longer, this just constrains pulse and repeats scanning frequency distance, so sampling gate depth
Contradiction between high PRF there are physics in principle.
Currently, the used method of current techniques is to needs when the PRF of needs exceeds the limitation of current scanning depth
PRF integral multiple reduces, this has resulted in several sampling gates, and last Doppler results will be that target in this several sampling gate is transported
Dynamic mixed information, this cannot really reflect the velocity information of required sampled targets.
Summary of the invention
In view of this, the embodiment of the present invention provide the coding that a kind of high impulse repeats scanning frequency signal, imaging method and
Device repeats the technical issues of scanning frequency is limited by propagation distance to solve target pulse.
In a first aspect, the embodiment of the invention provides a kind of high impulse repeat scanning frequency signal coding method include:
Determine that scanning signal set, the scanning signal set include at least two signals according to scanning depth；
Determine the rotation angular factors of the Fourier Transform of Fractional Order of the best zone of convergency of signal in scanning signal set；
The initial frequency of each scanning signal is calculated according to the rotation angular factors；
It is encoded according to the initial frequency pair of scanning signal scanning signal corresponding with the initial frequency.
Second aspect, the embodiment of the invention also provides method for ultrasonic imaging in a kind of Doppler Imaging systems, comprising:
Ultrasonic wave coded echo signal is received, and the ultrasonic wave coded echo signal is converted into digital signal；
The digital signal is subjected to Beam synthesis, generates radiofrequency signal；
To the radiofrequency signal anti-aliasing and decoding process, and by inphase-quadrature modem, inphase quadrature signal is generated；
Imaging is carried out according to the signal after inphase-quadrature modem, generates ultrasonography.
The third aspect, the embodiment of the invention also provides the code device that a kind of high impulse repeats scanning frequency signal, packets
It includes:
Gather determining module, for determining that scanning signal set, the scanning signal set include extremely according to scanning depth
Few two signals；
Angular factors determining module, in fractional order Fu for determining the best zone of convergency of signal in scanning signal set
The rotation angular factors of leaf transformation；
Initial frequency computing module, for calculating the starting frequency of each scanning signal according to the rotation angular factors
Rate；
Coding module, for the initial frequency pair scanning signal corresponding with the initial frequency according to the scanning signal
It is encoded.
Fourth aspect, the embodiment of the invention also provides a kind of ultrasonic imaging apparatus, comprising:
Digital signal conversion module is believed for receiving ultrasonic wave coded echo signal, and by the ultrasonic wave coded echo
Number be converted to digital signal；
Radiofrequency signal generation module generates radiofrequency signal for the digital signal to be carried out Beam synthesis；
Inphase quadrature signal generation module is used for the radiofrequency signal anti-aliasing and decoding process, and passes through positive intermodulation
System demodulation, generates inphase quadrature signal；
Image generation module generates ultrasonography for carrying out imaging according to the signal after inphase-quadrature modem.
High impulse provided in an embodiment of the present invention repeats coding, imaging method and the device of scanning frequency signal, by right
Different transmitting signals are encoded, and are decoded accordingly to the signal received, can determine the signal that receives with
Emit the corresponding relationship between signal.It is able to solve the limitation that scanning distance repeats pulse scanning frequency, is realized to height
The scanning of fast target.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, of the invention other
Feature, objects and advantages will become more apparent upon:
Fig. 1 is the process signal for the coding method that the high impulse that the embodiment of the present invention one provides repeats scanning frequency signal
Figure；
Fig. 2 is the process signal for the coding method that high impulse provided by Embodiment 2 of the present invention repeats scanning frequency signal
Figure；
Fig. 3 is the flow diagram for the method for ultrasonic imaging that the embodiment of the present invention three provides；
Fig. 4 is the flow diagram for the method for ultrasonic imaging that the embodiment of the present invention four provides；
Fig. 5 is the structural representation for the code device that the high impulse that the embodiment of the present invention five provides repeats scanning frequency signal
Figure；
Fig. 6 is the structural schematic diagram for the ultrasonic imaging apparatus that the embodiment of the present invention six provides.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It also should be noted that in order to just
In description, only some but not all contents related to the present invention are shown in the drawings.
Embodiment one
Fig. 1 is the process signal for the coding method that the high impulse that the embodiment of the present invention one provides repeats scanning frequency signal
Figure, the present embodiment are applicable to the case where encoding to the scanning signal that ultrasonic system emits, and this method can be by pulse
The code device of scanning frequency signal is repeated to execute, which can be realized by software/hardware mode, and can be integrated in corresponding
In ultrasonic wave detecting system.
Referring to Fig. 1, the high impulse repeats the coding method of scanning frequency signal, comprising:
S110 determines that scanning signal set, the scanning signal set include at least two signals according to scanning depth.
Modern ultrasound detection technique uses pulse reflection method, i.e., to one short pulse ultrasound of tested object emission
Wave, then receive it.And so on.In unit time, the number of driving probe transmitting ultrasonic wave is exactly pulse recurrence frequency.
For high-speed moving object, such as when histoorgan or the progress scanning of blood flow cell, since the speed of blood flow is opposite
Movement is very fast, needs to improve the scanning frequency of response, i.e., emits multiple scanning pulse signals within the unit time.Illustratively,
These scanning pulse signals form scanning signal set, include several scanning signals in each scanning signal set.Scanning letter
Number set include at least two scanning signals.
Scanning signal set can be determined by scanning depth, be D to current scanning depth, then its maximum pulse repeats to sweep
Looking into frequency is Prfmax=C/ (2*D), and wherein C is acoustic speed.In the present embodiment, maximum pulse can be repeated into scanning
All scanning pulse signals seat scanning pulse collection sent in the frequency Prfmax corresponding period.
S120 determines the rotation angle of the Fourier Transform of Fractional Order of the best zone of convergency of signal in scanning signal set
The factor.
Fourier Transform of Fractional Order (Fractional Fourier Transform, FrFT) is one kind of Fourier transformation
Generalized form, signal contain signal in the information of time-domain and frequency-domain simultaneously on fractional number order Fourier, have FrFT and permitted
The property that more conventional Fourier transforms do not have can be used to handle non-stationary signal.The basic meaning of Fourier space is, arbitrarily
Periodic function (signal) sine wave that infinite multiple frequencies are its basic frequency integer (including zero) times and remaining can be decomposed into
The sum of string wave.And Fourier transformation is expanded its period to infinite formation.Signal be can be understood as in time-frequency plane
Representation method of the reference axis after origin any rotation counterclockwise.
For multiple pulse signals that ultrasonic system is emitted, can regard as periodic signal, for the periodic signal,
It has the property that signal energy is gradually assembled as the fractional order domain of signal is gradually increased, is impact letter to frequency domain representation
Number.Based on this characteristic, unique impact factor rotation angular factors can be determined according to Fourier Transform of Fractional Order.
Specifically, rotation angular factors can be calculated by following formula:
Aifa=1-2/pi*atan (Tx_fs/BandWidthOfFM)；Wherein, aifa is rotation angular factors；Tx_fs
The sample frequency of expression system transmitting；The modulation bandwidth of BandWidthOfFM presentation code signal.
S130 calculates the initial frequency of each scanning signal according to the rotation angular factors.
Specifically, may include steps of:
The relative bandwidth frequency of each scanning signal in the scanning signal set is determined according to the rotation angular factors；
The centre frequency of the pulse emitted according to the relative bandwidth frequency and ultrasonic system calculates each scanning
The initial frequency of signal.
Several line segment images for crossing at 0 point can be shown as in time-domain and frequency-domain plane for scanning signal set, wherein often
The line segment of 0 point of a mistake all represents the relative frequency of each scanning signal in scanning signal set.
Specifically, can be calculated by following formula:
Bandwidth_Sig=4*nZ*abs (sin (aifa*pi/2))/TimeOfFM+D*f_Atten；
Wherein, the modulation duration of TimeOfFM presentation code signal, D indicate maximum scanning depth, and f_Atten indicates frequency
With the attenuation coefficient of change in depth.NZ is constant coefficient, can generally take 1, can be according to above-mentioned for every n scanning signal
The relative frequency of the scanning signal is calculated in formula.
For the probe of any one ultrasonic system, driving frequency can be set as not according to different diagnostic purposes
Same frequency.And driving frequency is in the centre frequency of setting, working efficiency highest.Each scanning letter in scanning signal set
Number frequency should be in a certain range of the centre frequency of probe.And it can be according to scanning signal each in scanning signal set
The frequency of each scanning signal in scanning signal set is calculated in relative frequency and the centre frequency of probe.Specifically, can be with
It is calculated by following formula:
Fstart=f0-BandWidthOfFM/2+n*bandwidth_Sig.
Wherein, f0 indicates the centre frequency of impulse ejection.
S140 is compiled according to the initial frequency pair of scanning signal scanning signal corresponding with the initial frequency
Code.
Since there is no the variations of phase for scanning signal, it can be corresponded to according to the frequency of scanning signal and generate scanning signal
Waveform signal encodes scanning signal according to the waveform signal.Specifically, scanning letter can be obtained by following formula
Number corresponding encoded signal:
Base_excitation (i)=cos (2*pi* (fstart*t+BandWidthOfFM/2/TimeOfFM*t.^
2))；Wherein t=(0,1/Tx_fs, 2/Tx_fs ... ... TimeOfFM).
High impulse provided in this embodiment repeats the coding method of scanning frequency signal, is emitted using ultrasonic system
Scanning signal set sets the characteristics of assembling in the time-domain and frequency-domain plane for rotating angular factors in Fourier Transform of Fractional Order, leads to
The frequency for calculating each scanning signal is crossed, then determines the respective waveforms of each scanning signal, and according to the waveform to scanning
Signal is encoded, so that ultrasonic system can be according to the wave after scanning Signal coding when receiving scanning signal echo
Shape, which determines, receives parameter, and then realizes and obtain the accurate velocity information of high-speed motion sampled targets.
Embodiment two
Fig. 2 is the process signal for the coding method that high impulse provided by Embodiment 2 of the present invention repeats scanning frequency signal
Figure.The present embodiment increases following steps based on above-described embodiment: determining scanning signal set number, and determines each scanning
Scanning serial number of the signal in corresponding scanning set；Correspondingly, by described described every according to rotation angular factors calculating
The initial frequency of a scanning signal is specific to optimize are as follows: according to the rotation angular factors and scanning signal in corresponding scanning collection
Scanning serial number in conjunction calculates the initial frequency of each scanning signal.
Correspondingly, high impulse provided by the present embodiment repeats the coding method of scanning frequency signal, specifically include:
S210 determines scanning signal set number, and determines scanning sequence of each scanning signal in corresponding scanning set
Number.
Scanning signal set can be determined by scanning depth, be D to current scanning depth, then its maximum pulse repeats to sweep
Looking into frequency is Prfmax=C/ (2*D), and wherein C is acoustic speed.I.e. in ultrasonic system, in a transmit cycle,
Prfmax ultrasonic signal can at most be emitted, and for the scanning target of high-speed motion, required pulse repeats to sweep
Looking into frequency is Prf_desire, if it is that Prf_desire is greater than Prfmax that desired pulse, which repeats scanning frequency, at 1
In transmit cycle, it cannot complete to send desired pulse repetition scanning number, need to complete to be required within multiple periods
Pulse repeats the transmission of scanning signal, and each period can send a scanning signal set.Specifically, can be by following public
The scanning signal set number N that formula emits needed for being calculated:
CEIL (Prf_desire/Prfmax), CEIL () expression round up.
It is the transmitting signal of n-th of sequence in some scanning signal set, tool for the signal of i-th transmitting
Body, scanning serial number of the scanning signal in corresponding scanning set can be determined by following formula:
N=mod (i, N).
S220 determines the rotation angle of the Fourier Transform of Fractional Order of the best zone of convergency of signal in scanning signal set
The factor.
S230 calculates institute according to the scanning serial number of the rotation angular factors and scanning signal in corresponding scanning set
State the initial frequency of each scanning signal.
Specifically, may include steps of:
According to it is described rotation angular factors and scanning signal in corresponding scanning set scanning serial number determine described in sweep
Look into the relative bandwidth frequency of each scanning signal in signal set；Specifically, can be calculated by following formula:
Bandwidth_Sig=4*nZ*abs (sin (aifa*pi/2))/TimeOfFM+D*f_Atten；
NZ generally takes the modulation duration of 1, TimeOfFM presentation code signal, and D indicates maximum scanning depth, f_Atten table
Show frequency with the attenuation coefficient of change in depth；
The centre frequency and scanning signal of the pulse emitted according to the relative bandwidth frequency, ultrasonic system are corresponding
Scanning serial number in scanning set calculates the initial frequency of each scanning signal.Specifically, can count in the following way
It obtains:
Fstart=f0-BandWidthOfFM/2+ (n-(N-1)/2) * bandwidth_Sig；
F0 indicates the centre frequency of impulse ejection, and N is the scanning signal set quantity for needing to emit.
S240 is compiled according to the initial frequency pair of scanning signal scanning signal corresponding with the initial frequency
Code.
The present embodiment is by increasing following steps: determining scanning signal set number, and determines each scanning signal in correspondence
Scanning set in scanning serial number；Correspondingly, calculating each scanning signal according to the rotation angular factors for described
Initial frequency, it is specific to optimize are as follows: according to the scanning of the rotation angular factors and scanning signal in corresponding scanning set
Serial number calculates the initial frequency of each scanning signal.Scanning frequency can be repeated according to maximum pulse to sampled targets institute
It is required that pulse repeat scanning frequency and carry out different set groupings, and the different scanning signals in set are encoded.
Embodiment three
Fig. 3 is the flow diagram for the method for ultrasonic imaging that the embodiment of the present invention three provides, and the present embodiment is applicable to
The case where being imaged in Doppler Imaging systems using ultrasonic signal, this method can be executed by ultrasonic imaging apparatus, should
Device can be realized by software/hardware mode, and can be integrated in corresponding ultrasonic wave detecting system.
Referring to Fig. 3, the method for ultrasonic imaging, comprising:
S310 receives ultrasonic wave coded echo signal, and the ultrasonic wave coded echo signal is converted to digital signal.
Ultrasonic probe in ultrasonic system is for emitting ultrasonic wave and receiving ultrasonic echo.And pass through wave echo
Amplification, filtering and Analog-digital Converter, obtain digitalized ultrasonic wave echo；The reception device may include the low of sequential connection
Noise amplifier circuit, Anti-aliasing Filter Circuits and D/A conversion circuit.The low noise amplifier circuit and the ultrasonic wave
Probe connection.
The digital signal is carried out Beam synthesis, generates radiofrequency signal by S320.
Beam synthesis refers to the output by member of respectively shaking in multispot array that certain geometries arrange by delay, weighted sum
After summation etc. comes out, make output that there is space directivity.Illustratively, the method that can use delay stack, obtains ultrasonic wave
Rf data.
S330 generates inphase quadrature signal to the radiofrequency signal anti-aliasing and decoding process, and by quadrature demodulation.
Illustratively, device can be gone using anti-aliasing filtering, under the premise of sample frequency is constant, by setting frequency
Low-pass filter filters the frequency content for being not required to occur.Since the signal of input includes I/Q two parts, need through orthogonal solution
It adjusts, to generate the orthogonal signalling with phase.
S340 carries out imaging according to the signal after inphase-quadrature modem, generates ultrasonography.
The present embodiment is believed the coding received by carrying out anti-aliasing and decoding process to the ultrasonic signal received
Number processing, generate inphase quadrature signal, may be implemented obtain the accurate velocity information of high-speed motion sampled targets.
Example IV
Fig. 4 is the flow diagram for the method for ultrasonic imaging that the embodiment of the present invention four provides.The present embodiment is with above-mentioned reality
It is specific to optimize are as follows: described pre- according to radiofrequency signal described in preset length samples by radiofrequency signal anti-aliasing based on applying example
If length according to the width of sampling gate and receive the sample frequency of ultrasonic wave coded echo signal and determine；Determine the radio frequency letter
Number Fourier Transform of Fractional Order rotation angular factors；The radio frequency letter sampled according to the rotation angle factor pair
Number carry out Fourier Transform of Fractional Order；The cutoff frequency of bandpass filter is set, and transformed radiofrequency signal is filtered
Processing；Inverse Fourier Transform of Fractional Order is carried out to filtered radiofrequency signal.
Referring to fig. 4, the method for ultrasonic imaging, comprising:
S410 receives ultrasonic wave coded echo signal, and the ultrasonic wave coded echo signal is converted to digital signal.
The digital signal is carried out Beam synthesis, generates radiofrequency signal by S420.
S430, according to radiofrequency signal described in preset length samples, the preset length according to the width of sampling gate and
The sample frequency for receiving ultrasonic wave coded echo signal determines.
Illustratively, sampling length can be determined in the following way, sampling length L=w_sig*2/C*Rx_fs,
In, w_sig is the width of sampling gate, and C is acoustic wave propagation velocity, and Rx_fs indicates to receive the sample frequency of signal.
S440 determines the rotation angular factors of the Fourier Transform of Fractional Order of the radiofrequency signal.
Illustratively, the rotation angular factors of Fourier Transform of Fractional Order can be determined in the following way, rotate angle
Factor aifa_rx=1-2/pi*atan (Rx_fs/BandWidthOfFM), wherein BandWidthOfFM presentation code signal
Modulation bandwidth；Rx_fs indicates to receive the sample frequency of signal.
S450 carries out fractional order Fourier change according to the radiofrequency signal that the rotation angle factor pair samples
It changes.
Radiofrequency signal is the compound of multiple signals, and there are serious cross term interferences.Fourier Transform of Fractional Order can press down
Video distribution cross term processed.Illustratively, Fourier Transform of Fractional Order, score can be carried out to radiofrequency signal in the following way
Signal Sig_frft=frft (sig, aifa_rx) after rank Fourier transformation, frft () indicate Fourier Transform of Fractional Order letter
Number.
S460 sets the cutoff frequency of bandpass filter, and is filtered to transformed radiofrequency signal.
For i-th scanning signal, 0 bandwidth is crossed for n-th Z after Fourier Transform of Fractional Order are as follows:
Bandwidth_Sig_rx=4*nZ*abs (sin (aifa_rx*pi/2))/TimeOfFM+D*f_Atten；Its
In, the modulation duration of TimeOfFM presentation code signal, D indicates maximum scanning depth, and f_Atten indicates frequency with change in depth
Attenuation coefficient.Z is constant coefficient, can generally take 1, for each scanning signal, can be calculated according to above-mentioned formula
The relative frequency of the scanning signal.So the then centre frequency of filter are as follows:
F0_bandfilter_i=f0+ (n-(N-1)/2) * bandwidth_Sig_rx
Then the upper lower limiting frequency of the bandpass 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.
Inverse Fourier Transform of Fractional Order is carried out to filtered radiofrequency signal, the filtered radiofrequency signal is made to restore former
Some frequencies.
S480 generates inphase quadrature signal to the signal decoding process after anti-aliasing, and by quadrature demodulation.
S490 carries out imaging according to the signal after inphase-quadrature modem, generates ultrasonography.
The present embodiment is by specifically optimizing radiofrequency signal anti-aliasing are as follows: believes according to radio frequency described in preset length samples
Number, the preset length is determined according to the width of sampling gate and the sample frequency of reception ultrasonic wave coded echo signal；It determines
The rotation angular factors of the Fourier Transform of Fractional Order of the radiofrequency signal；It is sampled according to the rotation angle factor pair
The radiofrequency signal carries out Fourier Transform of Fractional Order；The cutoff frequency of bandpass filter is set, and transformed radio frequency is believed
It number is filtered；Inverse Fourier Transform of Fractional Order is carried out to filtered radiofrequency signal.It can be in lesser frequency range
It is interior that radiofrequency signal is filtered, noise wave removing is removed, the precision of ultrasonic scanning is improved.
Embodiment five
Fig. 5 is the structural representation for the code device that the high impulse that the embodiment of the present invention five provides repeats scanning frequency signal
Figure, as shown in figure 5, described device includes:
Gather determining module 510, for determining that scanning signal set, the scanning signal set include according to scanning depth
At least two signals；
Angular factors determining module 520, for determining the fractional order of the best zone of convergency of signal in scanning signal set
The rotation angular factors of Fourier transformation；
Initial frequency computing module 530, for calculating rising for each scanning signal according to the rotation angular factors
Beginning frequency；
Coding module 540, for the initial frequency pair scanning corresponding with the initial frequency according to the scanning signal
Signal is encoded.
High impulse provided in this embodiment repeats the code device of scanning frequency signal, by different transmitting signals into
Row coding, and the signal received is decoded accordingly, can determine the signal received and emits pair between signal
It should be related to.The limitation that scanning distance repeats pulse scanning frequency can be released, the scanning to high-speed target is realized.
On the basis of the various embodiments described above, the initial frequency computing module, comprising:
Relative bandwidth determination unit is each swept for being determined in the scanning signal set according to the rotation angular factors
Look into the relative bandwidth of signal；
Frequency computing unit, for the centre frequency according to the relative bandwidth frequency and the pulse of ultrasonic system transmitting
Calculate the initial frequency of each scanning signal.
On the basis of the various embodiments described above, described device further include:
Scanning serial number determining module for determining scanning signal set number, and determines that each scanning signal is swept corresponding
Look into the scanning serial number in set；
Correspondingly, the relative bandwidth determination unit, is used for:
According to it is described rotation angular factors and scanning signal in corresponding scanning set scanning serial number determine described in sweep
Look into the relative bandwidth of each scanning signal in signal set.
The code device that high impulse provided by the embodiment of the present invention repeats scanning frequency signal can be used for executing the present invention
The high impulse that any embodiment provides repeats the coding method of scanning frequency signal, has corresponding functional module, realizes identical
Beneficial effect.
Embodiment six
Fig. 6 is the structural schematic diagram for the ultrasonic imaging apparatus that the embodiment of the present invention six provides, as shown in fig. 6, the dress
It sets and includes:
Digital signal conversion module 610, for receiving ultrasonic wave coded echo signal, and by the ultrasonic wave coded echo
Signal is converted to digital signal；
Radiofrequency signal generation module 620 generates radiofrequency signal for the digital signal to be carried out Beam synthesis；
Inphase quadrature signal generation module 630 is used for the radiofrequency signal anti-aliasing and decoding process, and by orthogonal
Modulation /demodulation generates inphase quadrature signal；
Image generation module 640 generates ultrasonogram for carrying out imaging according to the signal after inphase-quadrature modem
Picture.
On the basis of the various embodiments described above, the inphase quadrature signal generation module is used for:
According to radiofrequency signal described in preset length samples, the preset length is super according to the width of sampling gate and reception
The sample frequency of acoustic coded echo-signal determines；
Determine the rotation angular factors of the Fourier Transform of Fractional Order of the radiofrequency signal；
Fourier Transform of Fractional Order is carried out according to the radiofrequency signal that the rotation angle factor pair samples；
The cutoff frequency of bandpass filter is set, and transformed radiofrequency signal is filtered；
Inverse Fourier Transform of Fractional Order is carried out to filtered radiofrequency signal.
Ultrasonic imaging apparatus provided in this embodiment, by carrying out anti-aliasing and decoding to the ultrasonic signal received
The encoded signal received processing is generated inphase quadrature signal by processing, and it is accurate that acquisition high-speed motion sampled targets may be implemented
Velocity information.
Ultrasonic imaging apparatus provided by the embodiment of the present invention can be used for executing the super of any embodiment of that present invention offer
Sound wave imaging method has corresponding functional module, realizes identical beneficial effect.
Obviously, it will be understood by those skilled in the art that each module of the above invention or each operation can be by as above
The terminal device is implemented.Optionally, the embodiment of the present invention can be realized with the program that computer installation can be performed, thus
It can be stored in storage device and be executed by a processor, the program can store computer-readable deposits in a kind of
In storage media, storage medium mentioned above can be read-only memory, disk or CD etc.；Or they are fabricated to
Each integrated circuit modules, or by them multiple modules or operation be fabricated to single integrated circuit module and realize.This
Sample, the present invention is not limited to the combinations of any specific hardware and software.
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that
The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation,
It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention
It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also
It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.
Claims (2)
1. the coding method that a kind of high impulse repeats scanning frequency signal characterized by comprising
Determine that scanning signal set, the scanning signal set include at least two signals according to scanning depth；
Determine the rotation angular factors of the Fourier Transform of Fractional Order of the best zone of convergency of signal in scanning signal set；
The initial frequency of each scanning signal is calculated according to the rotation angular factors；
It is encoded according to the initial frequency pair of scanning signal scanning signal corresponding with the initial frequency；
Wherein, the initial frequency that each scanning signal is calculated according to the rotation angular factors, comprising:
It determines scanning signal set number, and determines scanning serial number of each scanning signal in corresponding scanning set；
The scanning letter is determined according to the scanning serial number of the rotation angular factors and scanning signal in corresponding scanning set
Number set in each scanning signal relative bandwidth；
The centre frequency of the pulse emitted according to the relative bandwidth frequency and ultrasonic system calculates rising for each scanning signal
Beginning frequency.
2. the code device that a kind of high impulse repeats scanning frequency signal characterized by comprising
Gather determining module, for determining that scanning signal set, the scanning signal set include at least two according to scanning depth
A signal；
Angular factors determining module, for determining that the fractional order Fourier of the best zone of convergency of signal in scanning signal set becomes
The rotation angular factors changed；
Initial frequency computing module, for calculating the initial frequency of each scanning signal according to the rotation angular factors；
Coding module is carried out for the initial frequency pair scanning signal corresponding with the initial frequency according to the scanning signal
Coding；
Wherein, the angular factors determining module includes: scanning serial number determining module, for determining scanning signal set number, and
Determine scanning serial number of each scanning signal in corresponding scanning set；
Relative bandwidth determination unit, is used for: according to the rotation angular factors and scanning signal in corresponding scanning set
Scanning serial number determines the relative bandwidth of each scanning signal in the scanning signal set；
Frequency computing unit, for being calculated according to the centre frequency of the relative bandwidth frequency and the pulse of ultrasonic system transmitting
The initial frequency of each scanning signal.
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 CN105919624A (en) | 2016-09-07 |
| CN105919624B true 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) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019183990A1 (en) * | 2018-03-30 | 2019-10-03 | 深圳迈瑞生物医疗电子股份有限公司 | Blood flow mapping processing method and ultrasonic imaging device |
| CN109199448B (en) * | 2018-09-04 | 2021-08-27 | 深圳市理邦精密仪器股份有限公司 | Pulse wave Doppler imaging method and device based on HPRF |
| CN110559014B (en) * | 2019-08-28 | 2021-10-26 | 华南理工大学 | Fractional order Fourier transform echo imaging method and system based on probe compensation |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100362000B1 (en) * | 2000-02-01 | 2002-11-22 | 주식회사 메디슨 | Ultrasound imaging method and apparatus based on pulse compression technique using modified golay code |
| US20090209858A1 (en) * | 2008-02-18 | 2009-08-20 | The Board Of Trustees Of The University Of Illinois | System and method for ultrasonic image processing |
| CN101785684B (en) * | 2010-01-27 | 2012-10-10 | 天津大学 | 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 |
| CN102724155B (en) * | 2012-05-17 | 2014-09-17 | 哈尔滨工程大学 | 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 |
| CN104980390B (en) * | 2015-05-31 | 2018-01-09 | 厦门大学 | A kind of demodulation method of diclinic rate combination chirp signals |
-
2016
- 2016-05-25 CN CN201610351777.1A patent/CN105919624B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105919624A (en) | 2016-09-07 |
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 | |
| US8864671B2 (en) | Methods and systems for color flow imaging | |
| JP2014006066A (en) | Interference detection device, interference remover, rader device, interference detection method and interference detection program | |
| CN105919624B (en) | A kind of high impulse repeats the coding and device of scanning frequency signal | |
| CN106529561A (en) | Sparkle artifact detection in Ultrasound color flow | |
| CN105997142A (en) | Ultrasonic system transmitting signal compositing and imaging method and device | |
| CN105997147B (en) | A kind of ultrasonic pulse Doppler imaging method and device | |
| JPH11309145A (en) | System and method for estimating and displaying doppler frequency shift | |
| CN107817297A (en) | A kind of ultrasonic imaging processing method and processing system based on ultrasonic echo RF data | |
| US10107901B2 (en) | Acoustic signal phase rotation processor | |
| US6758815B2 (en) | Apparatus and method for indicating mechanical stiffness properties of body tissue | |
| CN103654861A (en) | Ultrasonic umbilical cord blood flow measuring system and ultrasonic umbilical cord blood flow measuring method both based on pseudorandom codes | |
| CN105745556B (en) | For the method and apparatus by sonic transducer measurement relative velocity | |
| CN1886094B (en) | Doppler velocity detection device and ultrasonographic device using the same | |
| Tarbox et al. | Low-power ultrasound imaging systems using time delay spectrometry | |
| CN108652666B (en) | Doppler blood flow imaging generation method and device | |
| Fu et al. | Barker coded excitation using LFM carrier for improving axial resolution in ultrasound imaging | |
| JPH11309146A (en) | System and method for imaging flow of ultrasonic scatterer | |
| Jin et al. | A new scheme of coded ultrasound using Golay codes | |
| Wagner et al. | Compressed beamforming with applications to ultrasound imaging | |
| Andria et al. | Digital Methods for very accurate Ultrasonic sensor measurements | |
| JP2003250802A (en) | Ultrasonic diagnostic apparatus | |
| Alam et al. | Reduction of computational complexity in the butterfly search technique [ultrasonic blood flow imaging] | |
| CN1360877A (en) | Coding exciting method and equipment for liquid flow image |
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. |