Phantom for ultrasound measurement, and ultrasound CT device
10959706 ยท 2021-03-30
Assignee
Inventors
- Atsuro SUZUKI (Tokyo, JP)
- Kenichi Kawabata (Tokyo, JP)
- Takahide Terada (Tokyo, JP)
- Yushi Tsubota (Tokyo, JP)
- Wenjing Wu (Tokyo, JP)
- Kazuhiro Yamanaka (Tokyo, JP)
Cpc classification
A61B8/546
HUMAN NECESSITIES
A61B8/085
HUMAN NECESSITIES
G09B23/286
PHYSICS
A61B8/4416
HUMAN NECESSITIES
International classification
A61B8/00
HUMAN NECESSITIES
Abstract
Provided herein is a phantom capable of mimicking both a dense breast and a fatty breast. A phantom for ultrasound measurement includes: a first member that mimics an object of interest for measurement; and a second member having provided therein the first member. The second member has the property to decrease its sound speed with a temperature increase brought by external temperature control. The sound speed of the second member at a predetermined temperature is equal to the sound speed of a third member surrounding the second member. The first member and the second member are immiscible with each other.
Claims
1. An ultrasound CT device, comprising: a phantom for ultrasound measurement, comprising: a first member that mimics an object of interest for measurement; and a second member having provided therein the first member, the second member configured to decrease its sound speed with a temperature increase brought by external temperature control, the sound speed of the second member at a predetermined temperature being equal to a sound speed of a third member surrounding the second member, the first member and the second member being immiscible with each other; an input-output section capable of setting a preset temperature, a relative refractive index of the second member for a wave against the third member, or a breast type; and a controller configured to control a temperature of the third member surrounding the second member, according to at least one of the preset temperature, the relative refractive index of the second member for a wave against the third member, and the breast type, and wherein the controller calculates correction data by creating an image at varying temperatures and at varying solution concentrations for more than one size of the first member, and the controller uses the correction data to correct a quantitative value of a tumor obtained by clinical imaging.
2. The ultrasound CT device of claim 1, wherein the second member is one of an oil gel, urethane, and polyvinyl chloride plastisol.
3. The ultrasound CT device of claim 1, wherein the first member is liquid.
4. The ultrasound CT device of claim 3, wherein the liquid is one of water and salty water.
5. The ultrasound CT device of claim 4, wherein the salty water has a concentration that is higher than 0% and not more than 8%.
6. The ultrasound CT device of claim 1, wherein the third member is water.
7. The ultrasound CT device of claim 1, wherein the first member has a columnar shape or a spherical shape.
8. The ultrasound CT device of claim 1, wherein the first member has a surface with irregularities.
9. The ultrasound CT device of claim 3, wherein the phantom further comprises a plug or a film configured to seal the liquid.
10. The ultrasound CT device of claim 3, wherein the first member is embedded in the second member without being exposed at a surface of the second member.
11. The ultrasound CT device of claim 1, wherein the second member has one of a columnar shape, a tapered cup shape, and a bowl shape.
12. The ultrasound CT device of claim 1, wherein the controller controls the temperature of the third member, and sets a value of 0.9 to 1.05 for the relative refractive index of the second member for a wave against the third member.
13. The ultrasound CT device of claim 1, wherein the controller calculates at least one statistical information selected from a mean value, a standard deviation, a maximum value, a minimum value, and a median value across a plurality of facilities with regard to sound speeds and attenuations of the first member and the second member, and wherein the input-output section displays said at least one statistical information on a screen of the ultrasound CT device.
Description
BRIEF DESCRIPTION OF DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
DESCRIPTION OF EMBODIMENTS
(20) Embodiments of the present invention are described below with reference to the accompanying drawings. The drawings represent specific embodiments based on the principle of the present invention. However, the embodiments are intended to help understand the present invention, and are not to be construed to limit the present invention.
(21) The following embodiments are related to a phantom used for performance evaluation or correction of an ultrasound CT device. A typical configuration of an ultrasound CT device is described first, followed by how the relative refractive index of an object against the background water affects collection in imaging by an ultrasound CT device. A phantom for ultrasound measurement capable of mimicking both a dense breast and a fatty breast will be described last.
(22)
(23) The controller 6, the signal processor 7, and the memory section 8 may be provided by a general-purpose computer. The computer includes a processor such as a CPU (Central Processing Unit), a memory, and an auxiliary memory device such as hard disc. The processes of the controller 6 and the signal processor 7 may be implemented as functions of a program executed on the computer. The input-output section 9 is configured from, for example, an input section such as a keyboard and a pointing device (e.g., a mouse), and a display section such as a display.
(24)
(25)
(26) The following describes how the relative refractive index of the object against the background water affects collection in imaging by the ultrasound CT device 100. Because the object is in water in the ultrasound CT device 100, the transmitted ultrasound is received by the ultrasound transmitter and receiver units by first passing through water, and then the object, and again through water. Here, the sound speed difference between the object and the background water affects the position at which the transmitted waves are received. This is because the refraction of waves is determined by the sound speed difference between different media.
(27)
(28) In
(29)
(30) The following describes a phantom for ultrasound measurement capable of mimicking both a dense breast and a fatty breast. The property of the phantom 1 according to the embodiment is such that the first member 11 that mimics an object of interest for measurement is immiscible with the second member 12, and vice versa. The second member 12 has the property to decrease its sound speed with a temperature increase brought by external temperature control, and the sound speed of the second member 12 at a predetermined temperature is equal to the sound speed of the third member 13 (for example, water in the water tank 4) surrounding the second member 12.
(31) In the present embodiment, an oil gel is used as the material of the second member 12 of the phantom 1. The oil gel is produced by mixing a paraffin oil with a polymer called SEBS (Styrene-Ethylene/Butylenes-Styrene). In the present embodiment, the oil gel was produced by dissolving SEBS in a paraffin oil being stirred on a 175 C. hot plate, and cooling the mixture at ordinary temperature after deaeration.
(32)
(33)
(34) A feature of the phantom is that the sound speed and the attenuation rate of the inclusion region can be varied by using different solutions. Unlike the polyacrylamide gel, the oil gel has the desirable characteristic that it is not miscible with the solution in the inclusion region. The oil gel region is called a background region, a region mimicking normal tissue such as the mammary gland and the fat of a breast.
(35)
(36) Water has a sound speed of 1,520 [m/s] when it has a set temperature of 35 degrees Celsius in clinical imaging of a breast. A dense breast, which has more mammary gland and less fat, has a sound speed of 1,560 [m/s], and a fatty breast, which has more fat and less mammary gland, has a sound speed of 1,380 [m/s] (NPL 2). In imaging of a dense breast using the ultrasound CT device, the breast has a higher sound speed than the background water. In the case of a fatty breast, the breast has a lower sound speed than the background water. That is, by imaging the oil gel phantom 1 at a temperature below T.sub.0, the relative refractive index of the background oil gel region against the background water can be set to the value of when imaging a dense breast. Conversely, with the phantom 1 imaged at a temperature higher than T.sub.0, the relative refractive index of the background oil gel region against the background water can be set to the value of when imaging a fatty breast.
(37) The phantom 1 is disposed in the water tank 4, and the water surrounds the phantom 1. At the predetermined temperature T.sub.0, the oil gel (second member 12) has the same sound speed as the water (third member 13) surrounding the oil gel. As described above, the phantom 1 can mimic both a dense breast and a fatty breast by increasing and decreasing the temperature above and below T.sub.0. The third member 13 surrounding the oil gel is water in the present embodiment. However, the invention is not limited to this. A material that takes the same sound speed as the oil gel at the predetermined temperature also may be used instead of water.
(38)
(39)
(40) The following descriptions are given from the standpoint of refraction of waves, specifically, how the phantom imaged at different temperatures mimics the refraction of waves for breasts having different sound speeds. Here, the sound speed of water, and the sound speeds of the background region and the inclusion region of the phantom 1 being imaged are represented by S.sub.w(t), S.sub.BG(t), and S.sub.inc(t), respectively, where t is the temperature. By assuming that the temperature of water for clinical imaging is 35 C., the sound speed of water is S.sub.w(35)=1,520 [m/s], and the sound speeds of a breast and a tumor are S.sub.br and S.sub.c, respectively. In the phantom being imaged, the relative refractive index of the background region for ultrasound against water is n.sub.1 (=S.sub.BG(t)/S.sub.w(t)), and the relative refractive index of the inclusion region for ultrasound against the background region is n.sub.z (=S.sub.inc(t)/S.sub.BG(t)). The relative refractive index of a breast against water in clinical imaging is n.sub.3 (=S.sub.br/S.sub.w(35)), and the relative refractive index of a tumor against the breast is n.sub.4 (=S.sub.c/S.sub.br). The following relation is established by assuming that n.sub.1 and n.sub.3 are the same.
(41)
(42) It follows from this that S.sub.br can be determined as follows.
(43)
(44) Similarly, the following relation is established by assuming that n.sub.2 and n.sub.4 are the same.
(45)
(46) It follows from this that S.sub.c can be determined as follows.
(47)
(48)
(49)
(50) Method of Production of Background Region and Inclusion Region
(51) A method for producing the background region and the inclusion region is described below.
(52) The phantom 1 may have a different shape. A background region (oil gel) of a tapered shape with a top surface larger in diameter than the bottom surface can be produced by pouring a solution of SEBS in paraffin oil into a tapered cup-like container as shown in
(53) Suited as the container material is polypropylene, which can be cut. The inclusion region can be produced by inserting a round-rod mold into the paraffin oil solution before it gels, and cooling the solution. A cylindrical inclusion region is produced upon removing the round rod from the gel. To remove the mold, a narrow spatula is inserted at the area of contact between the mold and the gel with the oil gel being immersed in 30 to 40 degree Celsius water, and the round rod is pulled out while detaching the gel from the mold with the spatula. Examples of materials suited as the round rod include polypropylene, and fluororesin. The diameter of the columnar inclusion region is preferably about 1 to 50 mm, a diameter for mimicking breast cancer.
(54)
(55) The spherical space 11a and the cylindrical hole 11b are produced using a spherical mold having a thin round rod attached to the sphere. The diameter of the spherical space 11a is preferably about 1 to 50 mm. The round rod used to make the hole 11b has a diameter that is preferably no larger than a half of the sphere diameter. As with the case of the cylindrical inclusion region, the spherical space 11a and the hole 11b are produced by inserting the mold in a solution that has not gelled, and removing the mold after the solution has gelled. For removal, the spherical mold is pulled out by stretching the hole 11b with a spatula. Examples of materials suited for the mold include polypropylene, and fluororesin. Because the removal of the mold becomes more difficult as the hole becomes narrower, the mold material may use a water-soluble support material used for 3D printers, or a support material that can be dissolved with limonene. The spherical space and the hole can be produced by inserting the mold in the solution, and dissolving the mold after the solution has gelled. The inclusion region is sealed by disposing a plug 14 in the area (first surface 12a) connecting the hole 11b to outside.
(56)
(57) To exchange the solution, the inclusion region requires a hole for injecting the solution and removing the solution to outside. However, a phantom 1 may be produced that does not involve an exchange of solution.
(58)
(59)
(60) Phantom Material
(61) In the present embodiment, an oil gel is used for the material of the second member 12 (background region) of the phantom 1. Examples of materials other than oil gel include urethane, polyvinyl chloride plastisol (NPL 4, NPL 5).
(62) Performance Evaluation of Device Using Phantom
(63)
(64) First, a user sets a temperature for the measurement of the phantom 1. In this embodiment, the temperature is set in three patterns.
(65) (1) A temperature is directly entered through an interface on a display (S1801)
(66) (2) The relative refractive index or the sound speed difference of the phantom 1 against water is entered through an interface on a display (S1802-1)
(67) (3) A type of breast to be mimicked is selected (S1803-1)
(68) The memory section 8 contains a first table storing the relationship between temperature and the relative refractive index or the sound speed difference of the phantom 1 against water. In the case of (2), the controller 6 outputs a temperature that corresponds to the entered relative refractive index or sound speed difference, using the first table (S1802-2).
(69) The memory section 8 contains a second table storing the relationship between the breast type (for example, four types of extremely dense, heterogeneously dense, scattered fibroglandular, and fatty) and the corresponding relative refractive index or sound speed difference of the phantom 1 corresponding to each type against water. In the case of (3), the controller 6 obtains a relative refractive index or a sound speed difference corresponding to the selected breast type, using the second table, and outputs a temperature that corresponds to the acquired relative refractive index or sound speed difference, using the first table (S1803-2). As another example, the memory section 8 may contain a table storing the relationship between breast types and temperatures corresponding to each type.
(70) After the input made according to any of (1), (2), and (3), the controller 6 executes the temperature control process (S1804). Specifically, the controller 6 controls a temperature adjuster (not illustrated), and heats or cools the water inside the auxiliary tank 5 (S1805). The controller 6 sends the heated or cooled water to the water tank 4.
(71) The controller 6 puts itself in standby for a time period T.sub.1 (min) after sending water to the water tank 4 (S1806). After an elapsed time of T.sub.1 (min), the controller 6 determines whether the temperature of the water in the water tank 4 is at the preset temperature (S1807). The controller 6 repeats the steps S1806 and S1807 until the temperature reaches the preset temperature.
(72) Upon the temperature of the water in the water tank 4 reaching the preset temperature, the controller 6 controls the signal processor 7 to create an image of data for water to be used as correction data (S1808). The signal processor 7 outputs the image of water from data collected by the ultrasound transmitter and receiver unit array module 3.
(73) After the imaging of water is finished, the phantom 1 is installed in the water tank 4 (S1809). After an elapsed time T.sub.2 (min) (S1810), the controller 6 controls the signal processor 7 to create an image of the phantom 1 (S1811). The signal processor 7 outputs the sound speed value in the background region of the phantom 1 (S1812).
(74) The signal processor 7 determines whether the sound speed value in the background region is in the normal range (S1813). Here, the signal processor 7 determines that the sound speed value is in the normal range when the measured sound speed value is in a range of from VdV to V+dV, where V is the sound speed value in the background region of the phantom 1 at the preset temperature, and dV is the acceptable range.
(75) If the sound speed value is in the normal range in S1813, the signal processor 7 outputs the image of the phantom 1 (S1814). The signal processor 7 then outputs a quantitative value for the region of interest (inclusion region) (S1815). The controller 6 determines whether the quantitative value is in the normal range (S1816). When the quantitative value is in the normal range, the device is determined as being capable of producing an image as normal, and the process is finished. The device is inspected when the quantitative value is outside of the normal range in S1816 (S1819).
(76) When the sound speed value is outside of the normal range in S1813, the controller 6 determines whether a time period T.sub.3 (min) has elapsed since the start of phantom imaging (S1817). The steps from S1810 to S1813 are repeated when time T.sub.3 (min) has not elapsed. In the event where the sound speed value does not fall in the normal range even after a lapse of T.sub.3 (min) since imaging of the phantom started, it means that the phantom 1 has not reached the preset temperature even after a certain length of time. In this case, the device is inspected (S1818).
(77)
(78) The type of phantom entered through the first input section 1901a may be based on, for example, (1) the size and shape of the inclusion region, (2) the type of solution, and (3) the size, shape, and material of the background region. The first input section 1901a is configured so that these settings can be selected.
(79) Pushing the Start button on the GUI starts the procedures of
(80) Standardization Using Phantom
(81) Different facilities use different protocols for clinical imaging of breasts using the ultrasound CT device 100. In order to ensure the accuracy of diagnosis by clinical imaging, standardization is needed for parameters such as imaging conditions, and image processing methods. Ideally, a certain level of diagnosis accuracy should be ensured regardless of the facility conducting imaging.
(82) For standardization, the performance of a device needs to be quantitatively evaluated for different facilities, using the same phantom. The phantom 1 of the embodiment undergoes only small acoustic property changes over a long time period, and can be used for performance evaluations performed for standardization.
(83) As an example, the following conditions are used as the standardization protocol. The inclusion region in the oil gel/phantom is filled with salty water (3.5%). Water temperature is set to 25 C. The mode of ultrasound transmission is a fan beam. Ultrasound is projected in 256 directions. FBP (filtered back projection) is used for image reconstruction. The pixel size of image is 1 mm or less. Image filter is not used.
(84) The phantom 1 is imaged at different facilities using the ultrasound CT device 100. The ultrasound CT devices 100 at different facilities are connected to one another via a network. For example, image data may be collected by the ultrasound CT devices 100 at different facilities, and at least one statistical information selected from a mean value, a standard deviation, a maximum value, a minimum value, and a median value across the facilities may be calculated with regard to the sound speeds and the attenuations of the background region and the inclusion region. At each facility, the ultrasound CT device 100 displays the content of the standardization protocol, and the statistical information across the facilities. This makes it possible to compare the performance of the device at a given facility with the performance at other facilities. A user is thus able to confirm whether the performance of the ultrasound CT device 100 at his or her facility is at the same level of performance as compared to other facilities. A user may compare the numerical values (the sound speeds and the attenuations of the background region and the inclusion region) of the ultrasound CT device 100 at his or her facility with the statistical information, and inspect the device when the numerical values at his or her facility greatly differ from the statistical information. This contributes to the maintenance of the device, in addition to ensuring diagnosis accuracy.
(85) The standardization protocol is subject to change as the performance of the ultrasound CT device 100 improves, and the statistical information of quantitative values between facilities is updated every time a change is made to the standardization protocol.
(86) Correction of Quantitative Value for Object of Interest for Measurement
(87) The accuracy of the quantitative values of the sound speed and the attenuation imaged by the ultrasound CT device depends on the device's space resolution, which is determined by factors such as the frequency of the transmitted ultrasound, the size of the transmitter and receiver units, and the image reconstruction method. The accuracy of quantitative values decreases as the object of interest for measurement, specifically the tumor region or the inclusion region becomes smaller. The phantom 1 of the present embodiment has the inclusion region filled with a solution, and the sound speed of the inclusion region can be found with accuracy using a premeasured sound speed of the solution. The quantitative values of a tumor obtained by clinical imaging can be corrected in the manner described below.
(88) First, the ultrasound CT device 100 images the phantom 1 having an inclusion region of a spherical shape. The sound speed values of water, the background region, and the inclusion region in the imaged phantom are S.sub.w(t.sub.p), S.sub.BG(t.sub.p), and S.sub.inc(t.sub.p, d.sub.inc, c), respectively. Here, t.sub.p is the temperature of the phantom 1 being measured, d.sub.inc is the diameter of the inclusion region, and c is the salty water concentration. The quantitative value has a correction coefficient f(t.sub.p, d.sub.inc, c)=S.sub.inc(t.sub.p, d.sub.inc, c)/S.sub.BG(t.sub.p). A ROI (region of interest) is set for the reconstruction image of the phantom, with water, the background region, and the inclusion region having ROI values of ROI.sub.w(t.sub.p), ROI.sub.BG(t.sub.p), and ROI.sub.inc(t.sub.p, d.sub.inc, c), respectively. Because water and the background region can take a wide range of ROI, it can be assumed that ROI.sub.w(t.sub.p)=S.sub.w(t.sub.p), and ROI.sub.BG(t.sub.p)=S.sub.BG(t.sub.p). The proportion of the background region relative to water is r.sub.1.sup.p(t.sub.p)=ROI.sub.BG(t.sub.p)/ROI.sub.w(t.sub.p), and the proportion of the inclusion region relative to the background region is r.sub.2.sup.p(t.sub.p, d.sub.inc, c)=ROI.sub.inc(t.sub.p, d.sub.inc, c)/ROI.sub.BG(t.sub.p).
(89) For imaging, for example, a phantom having inclusion regions measuring 3, 5, 7, 10, 15, and 20 mm in diameter is measured at varying temperatures of 15, 17.5, 20, 22.5, 25, and 27.5 C., and at varying salty water concentrations of 0, 1, 2, 3, 4, 5, 6, 7, 8% in the imaged phantom. The ultrasound CT device 100 determines f(t.sub.p, d.sub.inc, c), r.sub.1.sup.p(t.sub.p), and r.sub.2.sup.p(t.sub.p, d.sub.inc, c) for all temperatures, diameters, and salty water concentrations. In this manner, the ultrasound CT device 100 calculates correction data as it creates an image at different temperatures and different solution concentrations for the inclusion regions of different sizes. The correction data so determined are stored in the memory section 8.
(90) For clinical imaging, sound speed values S.sub.w(t.sub.h), S.sub.bre, and S.sub.can(d.sub.can) are set for water, a breast, and a tumor region. Here, t.sub.h is the temperature at the time of clinical imaging, and d.sub.can is the diameter of when the tumor is approximated to a sphere. A ROI is set for the reconstruction image produced by clinical imaging, with water, the breast, and the tumor having ROI values of ROI.sub.w(t.sub.h), ROI.sub.bre, and ROI.sub.can(d.sub.can), respectively. Because water and breast can take a wide range of ROI, it can be assumed that ROI.sub.w(t.sub.h)=S.sub.w(t.sub.h), and ROI.sub.bre=S.sub.bre. The proportion of the breast region relative to water is r.sub.1.sup.h(t.sub.h)=ROI.sub.bre/ROI.sub.w(t.sub.h), and the proportion of the tumor region relative to the breast region is r.sub.2.sup.h(d.sub.can)=ROI.sub.can(d.sub.can)/ROI.sub.bre. The temperature at which r.sub.1.sup.p(t.sub.p) is closest to r.sub.1.sup.h(t.sup.h) is t.sub.p. The d.sub.inc value that is closest to the tumor diameter d.sub.can is d.sub.inc. The salty water concentration at which r.sub.2.sup.p(t.sub.p, d.sub.inc, c) is closest to r.sub.2.sup.h(d.sub.can) is c. The ultrasound CT device 100 determines the corrected quantitative value for the tumor region by ROI.sub.bref(t.sub.p, d.sub.inc, c). With this configuration, the ultrasound CT device 100 can correct the quantitative value of the clinically imaged tumor using the previously determined correction data.
(91) The phantom 1 for ultrasound measurement of the present embodiment includes the first member (inclusion region) 11 that mimics an object of interest for measurement, and the second member (background region) 12 having provided therein the first member 11. The properties of the first member 11 and the second member 12 are such that these are not miscible with each other. The second member 12 has the property to decrease its sound speed with a temperature increase brought by external temperature control, and the sound speed of the second member 12 at the predetermined temperature T.sub.0 is the same as the sound speed of the third member 13 (the solution in the water tank 4) surrounding the second member 12. With such a configuration, the second member 12 of the phantom 1 has a higher sound speed than the third member 13 at a temperature below T.sub.0, enabling the phantom 1 to mimic a dense breast. At a temperature higher than T.sub.0, the second member 12 has a lower sound speed than the third member 13, enabling the phantom 1 to mimic a fatty breast.
(92) The present invention is not limited to the embodiments described above, and includes many variations. The foregoing embodiments were described to help illustrate the present invention, and the invention is not necessarily required to include all of the configurations described above. A part of the configuration of a certain embodiment may be replaced with the configuration of some other embodiment. It is also possible to add the configuration of a certain embodiment to the configuration of some other embodiment. It is also possible to add other configuration to a part of the configurations of the embodiments, or delete and/or replace a part of the configurations of the embodiments.
(93) The functions of the controller 6 and the signal processor 7 may be implemented by software, specifically a processor interpreting and executing programs provided to implement the functions of these members. Information including programs and files for implementing functions may reside in memory, storage devices such as hard disc and SSD (Solid State Drive), or storage media such as IC cards, SD cards, and DVD. The foregoing functions may be implemented, either in part or as a whole, by hardware, for example, hardware designed with integrated circuits.
(94) The control lines and information lines used in the embodiments above merely represent lines that are considered to be necessary for the purpose of explanation, and do not necessarily represent the all control lines and information lines of a product. All configurations may be interconnected to one another.
REFERENCE SIGNS LIST
(95) 1: Phantom for ultrasound measurement 2: Bed 3: Ultrasound transmitter and receiver unit array module 4: Water tank 5: Auxiliary tank 6: Controller 7: Signal processor 8: Memory section 9: Input-output section 11: First member of phantom 12: Second member of phantom 13: Third member surrounding phantom 14: Plug 15: Film 100: Ultrasound CT device