CN104873221A - Ultrasonic-based liver fat quantification method and system - Google Patents

Ultrasonic-based liver fat quantification method and system Download PDF

Info

Publication number
CN104873221A
CN104873221A CN201510306622.1A CN201510306622A CN104873221A CN 104873221 A CN104873221 A CN 104873221A CN 201510306622 A CN201510306622 A CN 201510306622A CN 104873221 A CN104873221 A CN 104873221A
Authority
CN
China
Prior art keywords
liver
field region
ultrasonic
far
signal
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
Application number
CN201510306622.1A
Other languages
Chinese (zh)
Other versions
CN104873221B (en
Inventor
邵金华
孙锦
段后利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuxi Hisky Medical Technologies Co Ltd
Original Assignee
Wuxi Hisky Medical Technologies Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuxi Hisky Medical Technologies Co Ltd filed Critical Wuxi Hisky Medical Technologies Co Ltd
Priority to CN201510306622.1A priority Critical patent/CN104873221B/en
Publication of CN104873221A publication Critical patent/CN104873221A/en
Priority to HK16101627.1A priority patent/HK1213460A1/en
Application granted granted Critical
Publication of CN104873221B publication Critical patent/CN104873221B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

The invention provides an ultrasonic-based liver fat quantification method and system. According to recommendations, the system comprises an ultrasonic transmitter, an ultrasonic receiver, an ultrasonic signal processing module, a liver near-field positioning module, a liver far-field positioning module and an ultrasonic attenuation coefficient calculating module; the ultrasonic transmitter is used for transmitting pulse ultrasonic to liver tissues; the ultrasonic receiver is used for receiving ultrasonic echo signals of the liver tissues; the ultrasonic signal processing module is used for compensating, sampling, demodulating, envelope detecting or filtering of the ultrasonic echo signals of the liver tissues, so as to obtain processed ultrasonic echo corrected signals; the liver near-field positioning module is used for positioning a near-field region RN of a liver; the liver far-field positioning module is used for positioning a far-region of the liver; the ultrasonic attenuation coefficient calculating module is used for calculating an ultrasonic attenuation coefficient Epsilon of the liver.

Description

Based on hyperacoustic liver fat quantitative approach and system
Technical field
The present invention relates to a kind of based on hyperacoustic liver fat quantitative approach and system.
Background technology
At present, fatty liver is a kind of common DHD.Diagnosis morning of fatty liver, early treatment can effectively stop the progress of chronic hepatopathy and improve its prognosis.
Clinically, doctor can use ultrasonoscopy to carry out etiologic diagnosis to fatty liver, the method has that sensitivity is high, simple operation, low cost advantage, but its specificity is poor, comparatively large by doctor's subjective impact, therefore can only as a kind of etiologic diagnosis method.
CT (Computed Tomography, computed tomography) measures liver density and accurately can detect Liver fatty deposition, if be aided with correcting CT testing tube, then can accurate quantification liver fat, but there is the risk of radiation in it.
Nuclear magnetic resonance (Magnetic Resonance Imaging, and Magnetic Resonance Spectrum (MagneticResonance Spectrum MRI), MRS) a small amount of lipidosis in hepatocyte can be found, but it is the same with CT, it is expensive and the review time is long, not easily extensively promotes the use of.
Summary of the invention
In order to solve at least one problem mentioned in background technology, the invention provides a kind of based on hyperacoustic liver fat quantitative approach and system.
The present invention for achieving the above object, adopts following technical scheme: a kind of based on hyperacoustic liver fat quantitative approach, it to liver organization emission pulse ultrasonic, utilizes ultrasonic receiver to receive liver near-field region R respectively by ultrasonic transmitter nultrasound echo signal S 1with liver far-field region R fultrasound echo signal S 2, to ultrasound echo signal S 1and S 2carry out signal processing respectively and obtain ultrasonic echo corrected signal S p1and S p2, according to ultrasonic echo corrected signal S p1and S p2calculate the ultrasonic wave attenuation coefficient ε of liver, its computing formula is:
ϵ = ln ( I R N ) - ln ( I R F ) k * Δd * f ;
Wherein, for near-field region R nultrasonic echo corrected signal S p1the value of intensity or meansigma methods, for far-field region R fultrasonic echo corrected signal S p2the value of intensity or meansigma methods, Δ d is near-field region R nwith far-field region R fbetween distance, f is the mid frequency of ultrasonic echo corrected signal, and k is distance coefficient.
Further, the near-field region R of described liver n, and far-field region R fautomatically identify location, its localization method is as follows:
1), the ultrasonic echo corrected signal after process is divided into multiple detection subregion;
2), the eigenvalue of the ultrasonic echo corrected signal of liver in each detection subregion is calculated, and according to the near-field region R of this eigenvalue determination liver n, and far-field region R f.
Further, to the near-field region R of liver n, and far-field region R fwhen carrying out automatically identifying location, define a ultrasound echo signal and comprise n sampled point, the corresponding hepatic scan degree of depth is d, then every 1mm deep packet is containing n/d point, take z as spacing, n sampled point is divided into { d/z} section, { d/z}, for rounding up, is expressed as S i, i is the integer from 1 to t, t={d/z}-1, and the 1st section comprises { zn/d} point respectively to t segment signal;
Calculate each segment signal S respectively inakagami to distribute m value m i, the probability density function of ultrasound echo signal R can be expressed as:
f ( r ) = 2 m m r 2 m - 1 Γ ( m ) Ω m exp ( - m Ω r 2 ) U ( r ) ,
Wherein, Γ (m) and U (r) represents gamma function and unit-step function respectively; M is Nakagami distribution m value, can be calculated by following formula:
m = [ E ( R 2 ) ] 2 E [ R 2 - E ( R 2 ) ] 2 ,
Wherein, E (R 2), E [R 2-E (R 2)] 2for mean value function,
When m value is in 0 to 1 scope, ultrasound echo signal obeys pre-Rayleigh distribution; When m value equals 1, ultrasound echo signal obeys Rayleigh distribution; When m value is greater than 1, ultrasound echo signal obeys post-Rayleigh distribution,
Calculate each segment signal S respectively iaverage M iwith standard deviation SD i,
Ergodic signals S i, i=[1, t], if a certain segment signal S jmeet M j∈ [M n1, M n2], SD j∈ [SD n1, SD n2], d j∈ [d nP1, d nP2] and m j∈ [m n1, m n2], then by signal S jcorresponding region is set to liver near-field region R n, wherein, M n1and M n2be respectively the upper and lower threshold value of liver near-field region average, SD n1and SD n2be respectively the upper and lower threshold value of liver near-field region standard deviation, d nP1and d nP2be respectively the upper and lower threshold value of the liver near-field region degree of depth, m n1and m n2be respectively the upper and lower threshold value of liver near-field region Nakagami distribution m value;
Ergodic signals S i, i=[1, t], if a certain segment signal S kmeet M k∈ [M f1, M f2], SD k∈ [SD f1, SD f2], d k∈ [d fP1, d fP2] and m k∈ [m f1, m f2], then by S kcorresponding region is set to liver far-field region R f, wherein, M f1and M f2be respectively the upper and lower threshold value of liver far-field region average, SD f1and SD f2be respectively the upper and lower threshold value of liver far-field region standard deviation, d fP1and d fP2be respectively the upper and lower threshold value of the liver far-field region degree of depth, m f1and m f2be respectively the upper and lower threshold value of liver far-field region Nakagami distribution m value.
Further, the near-field region R of described liver n, and far-field region R falso constant depth can be adopted to locate, and its localization method is as follows:
The degree of depth is set at [d n1, d n2] region in scope is the near-field region R of liver n, the degree of depth is set at [d f1, d f2] region in scope is the far-field region R of liver f; Wherein, d n1and d n2be respectively the upper and lower depth threshold of near-field region; d f1and d f2be respectively the upper and lower depth threshold of far-field region.
Further, when constant depth is located, if those who are investigated need the near-field region R more accurately locating liver nwith far-field region R ftime, can according to the described near-field region R of body-mass index BMI correspondence setting of those who are investigated nupper and lower depth threshold d n1, d n2, and far-field region R fupper and lower depth threshold d f1, d f2.
Except above two kinds of locate modes, the near-field region R of described liver n, and far-field region R falso can carry out manual depth adjustment location, its localization method is as follows: user arranges liver near-field region R voluntarily n, and far-field region R f.
Further, the mode of described signal processing comprises following at least one or multiple: compensation, sampling, demodulation, envelope detected or Filtering Processing.
Further, the fat content C of liver is calculated by the ultrasonic wave attenuation coefficient ε of liver fAT, its computing formula is: C fAT=μ * ε+η, wherein, μ is the coefficient of ε, and η is constant.
Further, utilize above-mentioned ultrasound source, receive the ultrasound echo signal of kidney simultaneously, calculated the fat content C of liver by the ratio of the ultrasonic echo intensity of liver and kidney and the ultrasonic wave attenuation coefficient ε of liver fAT, its computing formula is: C fAT=α * r hR+ β * ε+γ, wherein, r hRfor the ratio of the ultrasonic echo intensity of liver and kidney, α is r hRcoefficient, β is the coefficient of ε, and γ is constant.
To achieve these goals, the present invention also provides a kind of liver fat quantitative system based on said method, comprising:
Ultrasonic transmitter, it is for liver organization emission pulse ultrasonic;
Ultrasonic receiver, it is for receiving the ultrasound echo signal of liver organization;
Ultrasonic signal processing module, it is for carrying out signal processing to the ultrasound echo signal of liver organization, to obtain the ultrasonic echo corrected signal after processing;
Liver near field locating module, it is for locating the near-field region R of liver n;
Liver far field locating module, it is for locating the far-field region R of liver f;
Ultrasonic attenuation coefficient computing module, it is according to described near-field region R nultrasonic echo corrected signal S p1with far-field region R fultrasonic echo corrected signal S p2calculate ultrasonic attenuation coefficient ε;
Display module, it is for showing ultrasonic attenuation coefficient ε;
Input module, it is for receiving the input information of user;
Main control module, it is for controlling and coordinating the work between above-mentioned modules.
Beneficial effect of the present invention: the present invention is based on hyperacoustic liver fat quantitative approach and system adds liver near field locating module and liver far field locating module, the ultrasonic attenuation coefficient of liver region can be obtained more accurately, thus improve the accuracy of liver fat detection; The present invention allows operator automatically and/or manually to set liver near-field region and liver far-field region, and setting and manual setting can operate alternately automatically, improve the convenience of operation under the prerequisite ensureing accuracy; The present invention compensates ultrasonic signal, sample, demodulation, at least one in envelope detected and Filtering Processing, to improve the accuracy and reliability that liver fat detects.
Accompanying drawing explanation
Fig. 1 is the system architecture diagram that the present invention is based on hyperacoustic liver fat quantitative approach and system.
Fig. 2 is the method flow diagram that the present invention is based on hyperacoustic liver fat quantitative approach and system.
Detailed description of the invention
Shown in Fig. 1, relate to a kind of based on hyperacoustic liver fat quantitative system,
Comprise ultrasonic transmitter, it is for liver organization emission pulse ultrasonic;
Ultrasonic receiver, it is for receiving the ultrasound echo signal of liver organization;
Ultrasonic signal processing module, it is for carrying out signal processing to the ultrasound echo signal of liver organization, to obtain the ultrasonic echo corrected signal after processing;
Liver near field locating module, it is for locating the near-field region R of liver n;
Liver far field locating module, it is for locating the far-field region R of liver f;
Ultrasonic attenuation coefficient computing module, it is according to described near-field region R nultrasonic echo corrected signal S p1with far-field region R fultrasonic echo corrected signal S p2calculate ultrasonic attenuation coefficient ε.
Wherein, described ultrasonic signal processing module to the signal processing mode that the ultrasound echo signal of liver organization carries out at least comprise following one or more: compensation, sampling, demodulation, envelope detected or Filtering Processing.Wherein, the computing formula of described ultrasonic attenuation coefficient ε is: wherein, for near-field region R nultrasonic echo corrected signal S p1the value of intensity, for far-field region R fultrasonic echo corrected signal S p2the value of intensity, Δ d is near-field region R nwith far-field region R fbetween distance, f is the mid frequency of ultrasonic echo corrected signal, and k is distance coefficient.
In addition, can also be near-field region R nultrasonic echo corrected signal S p1the meansigma methods of intensity. can also be far-field region R fultrasonic echo corrected signal S p2the meansigma methods of intensity.
After the ultrasonic wave attenuation coefficient obtaining liver, calculate hepatic fat content.Calculating hepatic fat content by ultrasonic wave attenuation coefficient is common practise, in this no longer superfluous words.
Ultrasonic attenuation coefficient ε in the present invention is also by ultrasonic echo corrected signal S pcarry out Autoregressive Spectrum Analysis (Girault J, Ossant F, Ouahabi A, et al.Time-varying autoregressive spectralestimation for ultrasound attenuation in tissue characterization.IEEE Transactionson Ultrasonics, Ferroelectrics and Frequency Control, 1998,45 (3): 650-659) obtain.
In native system, as k=1, two groups are adopted at the independently ultrasonic transducer organizing heteropleural to be measured for ultrasonic transmitter and ultrasonic receiver; As k=2, same group of ultrasonic transducer is adopted for ultrasonic transmitter and ultrasonic receiver.Visible, the supersonic generator that the present invention comprises and receptor can adopt same group of ultrasonic transducer, also two groups of independently ultrasonic transducers can be adopted, the advantage of same group of ultrasonic transducer is adopted to be cost-saving, reduce operation complexity, adopt the advantage of two groups of independently ultrasonic transducers to be to provide mode of operation more flexibly.
Also comprise fat content computing module in native system, it is for calculating the fat content C of liver fAT;
Display module, it is for showing ultrasonic attenuation coefficient ε;
Input module, it is for receiving the input information of user;
Main control module, it is for controlling and coordinating the work between above-mentioned modules.
In addition, in the present system, described fat content C fATcomputing formula be: C fAT=μ * ε+η, wherein, μ is the coefficient of ε, and η is constant.
In addition, in native system, the fat content C of liver fATalso by Liver and kidney echo when ultrasonic attenuation coefficient ε to the fat content C of liver fATcalculate, its computing formula is as follows:
Described fat content C fATcomputing formula be: C fAT=α * r hR+ β * ε+γ, wherein, r hRfor the ratio of the ultrasonic echo intensity of liver and kidney, α is r hRcoefficient, β is the coefficient of ε, and γ is constant.
Wherein, in above-mentioned parameter, μ, η, α, beta, gamma is demarcated by fatty liver calibration body mould.Can also be obtained by a certain amount of clinical sample statistics.Such as: α, beta, gamma can be set as α=62.6, β=168.1, γ=-27.9.
At setting μ, η, α, during beta, gamma, can set adaptively according to the BMI of those who are investigated.
In native system, liver near field locating module and liver far field locating module are by cooperatively interacting with to liver near-field region R n, and far-field region R fautomatically identify location, its localization method is as follows:
1), the ultrasonic echo corrected signal after process is divided into multiple detection subregion;
2), the eigenvalue of the ultrasonic echo corrected signal of liver in each detection subregion is calculated, and according to the near-field region R of this eigenvalue determination liver n, and far-field region R f.
In native system, liver near field locating module, liver far field locating module are also by cooperatively interacting with to liver near-field region R n, and far-field region R fbe fixed depth localization, its localization method is as follows:
The degree of depth is set at [d n1, d n2] region in scope is liver near-field region R n, the degree of depth is set at [d f1, d f2] region in scope is liver far-field region R f; Wherein, d n1and d n2be respectively the upper and lower depth threshold of near-field region; d f1and d f2be respectively the upper and lower depth threshold of far-field region.During use, can according to the upper and lower depth threshold d of the described near-field region of body-mass index BMI correspondence setting of those who are investigated n1, d n2, and the upper and lower depth threshold d of far-field region f1, d f2.
In native system, liver near field locating module, liver far field locating module are also by matching with to liver near-field region R n, and far-field region R fcarry out manual depth adjustment location, its localization method is as follows: user arranges liver near-field region R by described input module n, and far-field region R fdirty.
In native system, automatically identify liver near-field region R by such as under type n, far-field region R f:
If one ultrasound echo signal comprises n sampled point, corresponding scan depths is d (unit: mm), then every 1mm deep packet is containing n/d point.With z (unit: mm) for spacing, n sampled point is divided into d/z} section (d/z} is for rounding up), be expressed as S i, i is the integer from 1 to t, t={d/z}-1.1st section comprises { zn/d} point respectively to t segment signal.
Calculate each segment signal S respectively inakagami to distribute m value m i.Nakagami statistical model is the one of Ultrasonic tissue characterization technology, based on Nakagami statistical model, the probability density function of ultrasound echo signal R can be expressed as (P.M.Shankar, " A general statistical model for ultrasonicbackscattering from tissues; " IEEE Trans.Ultrason.Ferroelec.Freq.Contr., 2000 (47): 727-736):
f ( r ) = 2 m m r 2 m - 1 Γ ( m ) Ω m exp ( - m Ω r 2 ) U ( r ) ,
Wherein, Γ (m) and U (r) represents gamma function and unit-step function respectively; M is Nakagami distribution m value, can be calculated by following formula:
m = [ E ( R 2 ) ] 2 E [ R 2 - E ( R 2 ) ] 2 ,
Wherein, E (R 2), E [R 2-E (R 2)] 2for mean value function.When m value is in (0,1) scope, ultrasound echo signal obeys pre-Rayleigh distribution; When m value equals 1, ultrasound echo signal obeys Rayleigh distribution; When m value is greater than 1, ultrasound echo signal obeys post-Rayleigh distribution.
Calculate each segment signal S respectively iaverage M iwith standard deviation SD i.
Ergodic signals S i, i=[1, t], if a certain segment signal S jmeet M j∈ [M n1, M n2], SD j∈ [SD n1, SD n2], d j∈ [d nP1, d nP2] and m j∈ [m n1, m n2], then by signal S jcorresponding region is set to liver near-field region R n, wherein, M n1and M n2be respectively the upper and lower threshold value of liver near-field region average, SD n1and SD n2be respectively the upper and lower threshold value of liver near-field region standard deviation, d nP1and d nP2be respectively the upper and lower threshold value of the liver near-field region degree of depth, m n1and m n2be respectively the upper and lower threshold value of liver near-field region Nakagami distribution m value.
Ergodic signals S i, i=[1, t], if a certain segment signal S kmeet M k∈ [M f1, M f2], SD k∈ [SD f1, SD f2], d k∈ [d fP1, d fP2] and m k∈ [m f1, m f2], then by S kcorresponding region is set to liver far-field region R f, wherein, M f1and M f2be respectively the upper and lower threshold value of liver far-field region average, SD f1and SD f2be respectively the upper and lower threshold value of liver far-field region standard deviation, d fP1and d fP2be respectively the upper and lower threshold value of the liver far-field region degree of depth, m f1and m f2be respectively the upper and lower threshold value of liver far-field region Nakagami distribution m value.
It should be noted that, " comprising ", " comprising " of recording in literary composition or its any other variant are intended to contain comprising of nonexcludability, thus make to comprise the process of a series of key element, method, article or equipment and not only comprise those key elements, but also comprise other key elements clearly do not listed, or also comprise by the intrinsic key element of this process, method, article or equipment.When not more restrictions, the key element limited by statement " comprising ... ", and be not precluded within process, method, article or the equipment comprising described key element and also there is other identical element.
To the above-mentioned explanation of the disclosed embodiments, professional and technical personnel in the field are realized or uses the present invention.To be apparent for those skilled in the art to the multiple amendment of these embodiments, General Principle as defined herein can without departing from the spirit or scope of the present invention, realize in other embodiments.Therefore, the present invention can not be restricted to these embodiments shown in this article, but will meet the widest scope consistent with principle disclosed herein and features of novelty.

Claims (10)

1., based on a hyperacoustic liver fat quantitative approach, it is characterized in that, comprise the following steps:
To liver organization emission pulse ultrasonic;
Receive liver near-field region R respectively nultrasound echo signal S 1with liver far-field region R fultrasound echo signal S 2;
To S 1and S 2carry out signal processing respectively and obtain ultrasonic echo corrected signal S p1and S p2;
Calculate the ultrasonic wave attenuation coefficient ε of liver, its computing formula is:
ϵ = ln ( I R N ) - ln ( I R F ) k * Δd * f ;
Wherein, for near-field region R nultrasonic echo corrected signal S p1the value of intensity or meansigma methods, for far-field region R fultrasonic echo corrected signal S p2the value of intensity or meansigma methods, Δ d is near-field region R nwith far-field region R fbetween distance, f is the mid frequency of ultrasonic echo corrected signal, and k is distance coefficient.
2. as claimed in claim 1 based on hyperacoustic liver fat quantitative approach, it is characterized in that, described liver near-field region R nwith described liver far-field region R fautomatically identify location, its localization method is as follows:
1), by the ultrasonic echo corrected signal S after process p1and S p2be divided into multiple detection subregion;
2), calculate the eigenvalue of the ultrasonic echo corrected signal of liver in each detection subregion, and determine described liver near-field region R according to this eigenvalue nwith described liver far-field region R f.
3. one as claimed in claim 2 is based on hyperacoustic liver fat quantitative approach, it is characterized in that,
Define a ultrasound echo signal and comprise n sampled point, the corresponding hepatic scan degree of depth is d, then every 1mm deep packet is containing n/d point, take z as spacing, n sampled point is divided into { d/z} section, { d/z}, for rounding up, is expressed as S i, i is the integer from 1 to t, t={d/z}-1, and the 1st section comprises { zn/d} point respectively to t segment signal;
Calculate each segment signal S respectively inakagami to distribute m value m i, the probability density function of ultrasound echo signal R can be expressed as:
f ( r ) = 2 m m r 2 m - 1 Γ ( m ) Ω m exp ( - m Ω r 2 ) U ( r ) ,
Wherein, Γ (m) and U (r) represents gamma function and unit-step function respectively; M is Nakagami distribution m value, can be calculated by following formula:
m = [ E ( R 2 ) ] 2 E [ R 2 - E ( R 2 ) ] 2 ,
Wherein, E (R 2), E [R 2-E (R 2)] 2for mean value function,
When m value is in 0 to 1 scope, ultrasound echo signal obeys pre-Rayleigh distribution; When m value equals 1, ultrasound echo signal obeys Rayleigh distribution; When m value is greater than 1, ultrasound echo signal obeys post-Rayleigh distribution,
Calculate each segment signal S respectively iaverage M iwith standard deviation SD i,
Ergodic signals S i, i=[1, t], if a certain segment signal S jmeet M j∈ [M n1, M n2], SD j∈ [SD n1, SD n2], d j∈ [d nP1, d nP2] and m j∈ [m n1, m n2], then by signal S jcorresponding region is set to liver near-field region R n, wherein, M n1and M n2be respectively the upper and lower threshold value of liver near-field region average, SD n1and SD n2be respectively the upper and lower threshold value of liver near-field region standard deviation, d nP1and d nP2be respectively the upper and lower threshold value of the liver near-field region degree of depth, m n1and m n2be respectively the upper and lower threshold value of liver near-field region Nakagami distribution m value;
Ergodic signals S i, i=[1, t], if a certain segment signal S kmeet M k∈ [M f1, M f2], SD k∈ [SD f1, SD f2], d k∈ [d fP1, d fP2] and m k∈ [m f1, m f2], then by S kcorresponding region is set to liver far-field region R f, wherein, M f1and M f2be respectively the upper and lower threshold value of liver far-field region average, SD f1and SD f2be respectively the upper and lower threshold value of liver far-field region standard deviation, d fP1and d fP2be respectively the upper and lower threshold value of the liver far-field region degree of depth, m f1and m f2be respectively the upper and lower threshold value of liver far-field region Nakagami distribution m value.
4. as claimed in claim 1 a kind of based on hyperacoustic liver fat quantitative approach, it is characterized in that, described liver near-field region R nwith described liver far-field region R fbe fixed depth localization, its localization method is as follows:
The degree of depth is set at [d n1, d n2] region in scope is described liver near-field region R n, the degree of depth is set at [d f1, d f2] region in scope is described liver far-field region R f; Wherein, d n1and d n2be respectively the upper and lower depth threshold of near-field region; d f1and d f2be respectively the upper and lower depth threshold of far-field region.
5. one as claimed in claim 4 is based on hyperacoustic liver fat quantitative approach, it is characterized in that, according to the described liver near-field region R of body-mass index BMI correspondence setting of those who are investigated nupper and lower depth threshold d n1, d n2with described liver far-field region R fupper and lower depth threshold d f1, d f2.
6. as claimed in claim 1 a kind of based on hyperacoustic liver fat quantitative approach, it is characterized in that, described liver near-field region R nwith described liver far-field region R fcarry out manual depth adjustment location, its localization method is as follows: user arranges described liver near-field region R voluntarily nwith described liver far-field region R f.
7. as claimed in claim 1 based on hyperacoustic liver fat quantitative approach, it is characterized in that, the mode of described signal processing comprises at least one in compensation, sampling, demodulation, envelope detected and Filtering Processing.
8. as claimed in claim 1 based on hyperacoustic liver fat quantitative approach, it is characterized in that, calculated the fat content C of liver by the ultrasonic wave attenuation coefficient ε of liver fAT, its computing formula is: C fAT=μ * ε+η, wherein, μ is the coefficient of ε, and η is constant.
9. as claimed in claim 1 based on hyperacoustic liver fat quantitative approach, it is characterized in that, utilize above-mentioned ultrasound source, receive the ultrasound echo signal of kidney simultaneously, calculated the fat content C of liver by the ratio of the ultrasonic echo intensity of liver and kidney and the ultrasonic wave attenuation coefficient ε of liver fAT, its computing formula is: C fAT=α * r hR+ β * ε+γ, wherein, r hRfor the ratio of the ultrasonic echo intensity of liver and kidney, α is r hRcoefficient, β is the coefficient of ε, and γ is constant.
10., based on a liver fat quantitative system for method described in the arbitrary claim of claim 1-9, it is characterized in that, comprising:
Ultrasonic transmitter, it is for liver organization emission pulse ultrasonic;
Ultrasonic receiver, it is for receiving the ultrasound echo signal of liver organization;
Ultrasonic signal processing module, it is for carrying out signal processing to the ultrasound echo signal of liver organization, to obtain the ultrasonic echo corrected signal after processing;
Liver near field locating module, it is for locating the near-field region R of liver n;
Liver far field locating module, it is for locating the far-field region R of liver f;
Ultrasonic attenuation coefficient computing module, it is according to described near-field region R nultrasonic echo corrected signal S p1with far-field region R fultrasonic echo corrected signal S p2calculate ultrasonic attenuation coefficient ε;
Display module, it is for showing described ultrasonic attenuation coefficient ε;
Input module, it is for receiving the input information of user;
Main control module, it is for controlling and coordinating the work between above-mentioned modules.
CN201510306622.1A 2015-06-05 2015-06-05 Liver fat quantitative approach and system based on ultrasonic wave Active CN104873221B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201510306622.1A CN104873221B (en) 2015-06-05 2015-06-05 Liver fat quantitative approach and system based on ultrasonic wave
HK16101627.1A HK1213460A1 (en) 2015-06-05 2016-02-15 Ultrasonic based method and system for quantifying liver fat

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510306622.1A CN104873221B (en) 2015-06-05 2015-06-05 Liver fat quantitative approach and system based on ultrasonic wave

Publications (2)

Publication Number Publication Date
CN104873221A true CN104873221A (en) 2015-09-02
CN104873221B CN104873221B (en) 2018-03-13

Family

ID=53941096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510306622.1A Active CN104873221B (en) 2015-06-05 2015-06-05 Liver fat quantitative approach and system based on ultrasonic wave

Country Status (2)

Country Link
CN (1) CN104873221B (en)
HK (1) HK1213460A1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105930665B (en) * 2016-04-26 2018-06-29 北京工业大学 A kind of liver fat method for quantitative measuring based on ultrasonic attenuation coefficient
CN108366784A (en) * 2015-11-30 2018-08-03 奥林巴斯株式会社 The working procedure of ultrasound observation apparatus, the working method of ultrasound observation apparatus and ultrasound observation apparatus
CN110604595A (en) * 2019-05-21 2019-12-24 深圳迈瑞生物医疗电子股份有限公司 Fatty liver quantitative analysis method and fatty liver quantitative analysis system
CN111513762A (en) * 2019-10-30 2020-08-11 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic data processing method, ultrasonic data processing system and computer storage medium
WO2021018103A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and apparatus, device, and storage medium
CN112469338A (en) * 2018-08-29 2021-03-09 深圳迈瑞生物医疗电子股份有限公司 Liver detection device based on ultrasonic, ultrasonic equipment and ultrasonic imaging method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040125995A1 (en) * 2002-09-26 2004-07-01 Jeong Ji Wook System and method of measuring fat content in target organ and recording medium of recording fat content measuring program
WO2012158877A2 (en) * 2011-05-17 2012-11-22 University Of Rochester Non-invasive assessment of liver fat by crawling wave dispersion
CN102895003A (en) * 2011-07-27 2013-01-30 高鑫 Ultrasonic image quantitative diagnosing system and signal processing method thereof
CN103750861A (en) * 2014-01-21 2014-04-30 深圳市一体医疗科技股份有限公司 Liver fat detection system based on ultrasound
CN103784163A (en) * 2014-01-21 2014-05-14 深圳市一体医疗科技股份有限公司 Liver fat quantitative system based on ultrasonic

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040125995A1 (en) * 2002-09-26 2004-07-01 Jeong Ji Wook System and method of measuring fat content in target organ and recording medium of recording fat content measuring program
WO2012158877A2 (en) * 2011-05-17 2012-11-22 University Of Rochester Non-invasive assessment of liver fat by crawling wave dispersion
CN102895003A (en) * 2011-07-27 2013-01-30 高鑫 Ultrasonic image quantitative diagnosing system and signal processing method thereof
CN103750861A (en) * 2014-01-21 2014-04-30 深圳市一体医疗科技股份有限公司 Liver fat detection system based on ultrasound
CN103784163A (en) * 2014-01-21 2014-05-14 深圳市一体医疗科技股份有限公司 Liver fat quantitative system based on ultrasonic

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108366784A (en) * 2015-11-30 2018-08-03 奥林巴斯株式会社 The working procedure of ultrasound observation apparatus, the working method of ultrasound observation apparatus and ultrasound observation apparatus
CN105930665B (en) * 2016-04-26 2018-06-29 北京工业大学 A kind of liver fat method for quantitative measuring based on ultrasonic attenuation coefficient
CN112469338A (en) * 2018-08-29 2021-03-09 深圳迈瑞生物医疗电子股份有限公司 Liver detection device based on ultrasonic, ultrasonic equipment and ultrasonic imaging method
CN112469338B (en) * 2018-08-29 2023-11-03 深圳迈瑞生物医疗电子股份有限公司 Device for detecting liver based on ultrasound, ultrasound equipment and ultrasound imaging method
US11890133B2 (en) 2018-08-29 2024-02-06 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Ultrasound-based liver examination device, ultrasound apparatus, and ultrasound imaging method
CN110604595A (en) * 2019-05-21 2019-12-24 深圳迈瑞生物医疗电子股份有限公司 Fatty liver quantitative analysis method and fatty liver quantitative analysis system
WO2021018103A1 (en) * 2019-08-01 2021-02-04 无锡海斯凯尔医学技术有限公司 Ultrasonic signal processing method and apparatus, device, and storage medium
CN111513762A (en) * 2019-10-30 2020-08-11 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic data processing method, ultrasonic data processing system and computer storage medium

Also Published As

Publication number Publication date
CN104873221B (en) 2018-03-13
HK1213460A1 (en) 2016-07-08

Similar Documents

Publication Publication Date Title
CN104873221A (en) Ultrasonic-based liver fat quantification method and system
EP0066343B2 (en) Method and apparatus for measuring ultrasonic attenuation characteristics
Rickey et al. A wall-less vessel phantom for Doppler ultrasound studies
Nam et al. Simultaneous backscatter and attenuation estimation using a least squares method with constraints
WO2017107660A1 (en) Method and apparatus for measuring viscoelastic parameter of viscoelastic medium
EP3240484B1 (en) Method for quantifying the elasticity of a material by ultrasounds
CN104622509A (en) Ultrasonic diagnostic apparatus and elastic evaluation method
WO1991011146A1 (en) Ultrasonic method and apparatus for determination of vessel location and size
EP3111206A1 (en) Methods and systems for measuring properties with ultrasound
EP0038331B1 (en) Signal processing equipment for ultrasonic visualisation
US20180143165A1 (en) Local Speed Of Sound Estimation Method For Medical Ultrasound
KR20210069644A (en) Method for determining the speed of sound in a medium, an ultrasound imaging system embodying the method
US11737725B2 (en) Intelligent guided wave elastography
CN104116524B (en) A kind of ultrasonic attenuation coefficient bucking-out system and liver fat detection system
Daugschies et al. The preliminary evaluation of a 1 MHz ultrasound probe for measuring the elastic anisotropy of human cortical bone
Kimme-Smith et al. Assurance of consistent peak-velocity measurements with a variety of duplex Doppler instruments.
Branca et al. Comparative evaluation of ultrasound scanner accuracy in distance measurement
US11559285B2 (en) Reflection ultrasound tomographic imaging using full-waveform inversion
CN204765705U (en) Quantitative system of liver fat based on ultrasonic wave
Chenevert et al. Ultrasonic computed tomography of the breast. Improvement of image quality by use of cross-correlation time-of-flight and phase-insensitive attenuation measurements.
EP3698721A2 (en) Human tissue quasi-elastic coefficient and elasticity measurement methods and device
Shi et al. Deep learning for TOF extraction in bone ultrasound tomography
Lagalla et al. Image quality control in breast ultrasound
Wen et al. Design of ultrasonic tomography system for biomedical imaging
US20230380810A1 (en) Ultrasonic imaging system and ultrasonic imaging method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 1213460

Country of ref document: HK

CB02 Change of applicant information
CB02 Change of applicant information

Address after: 214000, Jiangsu, Wuxi new area, Taihu international science and Technology Park University Science Park, building B401, No. 530, No.

Applicant after: Wuxi Heskell Medical Technology Co., Ltd.

Address before: 214000 room B401, science and Technology Park, Taihu international science and Technology Park, Wuxi national hi tech Industrial Development Zone, Wuxi, Jiangsu, China, 530

Applicant before: Wuxi Heskell Medical Technology Co., Ltd.

GR01 Patent grant
GR01 Patent grant
REG Reference to a national code

Ref country code: HK

Ref legal event code: GR

Ref document number: 1213460

Country of ref document: HK

PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Liver fat quantitative method and system based on ultrasound

Effective date of registration: 20210330

Granted publication date: 20180313

Pledgee: Bank of China Co.,Ltd. Wuxi High tech Industrial Development Zone Branch

Pledgor: Wuxi Hisky Medical Technologies Co.,Ltd.

Registration number: Y2021980002279

PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20220114

Granted publication date: 20180313

Pledgee: Bank of China Co.,Ltd. Wuxi High tech Industrial Development Zone Branch

Pledgor: WUXI HISKY MEDICAL TECHNOLOGIES Co.,Ltd.

Registration number: Y2021980002279

PC01 Cancellation of the registration of the contract for pledge of patent right