FIELD OF THE INVENTION
[0001] The invention relates to a heat sink parameter determination apparatus, and computer
program for determining a parameter of a heat sink within an object. The invention
relates further to a temperature distribution determination apparatus, and computer
program for determining a temperature distribution within the object, which use the
determined heat sink parameter. The invention relates also to a heating system comprising
the temperature distribution determination apparatus.
BACKGROUND OF THE INVENTION
[0002] US 2007/106157 A1 discloses a system for selectively delivering a thermal therapy to a specific mass
of tissue. The system comprises a thermal source, an ultrasound imaging probe, a means
for identifying when the specific mass of tissue has reached a predetermined temperature,
and a controller being adapted to implement the functions of using the thermal source
to change a temperature of the specific mass of tissue and collecting a first set
of ultrasound data as the specific mass of tissue returns to a previous temperature.
The controller is further adapted to implement the functions of using the first set
of ultrasound data to determine a thermal diffusivity parameter corresponding to the
specific mass of tissue, using the thermal source to change a temperature of the specific
mass of tissue to the predetermined value, while collecting a second set of ultrasound
data, and determining a length of time required to change the temperature of the specific
mass of tissue to the predetermined value. The controller is further adapted to implement
the function of using the thermal diffusivity parameter, the second set of ultrasound
data, the length of time, and a biological tissue heat transfer model to determine
a thermal source magnitude parameter corresponding to the thermal source and the specific
mass of tissue, thereby calibrating the biological tissue heat transfer model to the
thermal source and the specific mass of tissue.
[0003] EP 2 387 963 A1 discloses a temperature distribution determining apparatus for determining a temperature
distribution within an object caused by applying energy to the object, wherein the
temperature distribution determining apparatus comprises a temperature distribution
measuring unit for measuring a spatially and temporally dependent first temperature
distribution in the object, while the energy is applied to the object such that the
object is heated to a temperature within a first temperature range, and a temperature
distribution estimating unit for estimating a spatially and temporally dependent second
temperature distribution in the object within a second temperature range, which is
different to the first temperature range, based on the spatial and temporal dependence
of the measured first temperature distribution.
[0004] US 2013/296743 A1 discloses an ultrasound scanner comprising transmitters operable to scan with beams
of ultrasound in a scan sequence and a processor configured to create the scan sequence
as a function of a therapy plan and to turn off the transmitters for at least a portion
of the scan sequence. The ultrasound scanner further comprises a trigger output configured
to output triggers to a therapy device, wherein the output triggers are responsive
to the scan sequence.
[0005] WO 2010/140125 A1 discloses a therapeutic system comprising an ultrasound therapy unit arranged to
insonify at least a portion of a body of a patient with high intensity ultrasound
and a magnetic resonance imaging unit arranged to acquire magnetic resonance signals
from the portion of the body and to reconstruct a thermographic magnetic resonance
image from the magnetic resonance signals. The therapeutic system further comprises
an ultrasound diagnostic unit arranged to acquire ultrasound signals from the portion
of the body and to derive at least one local temperature value from the ultrasound
signals.
[0006] US 2011/0060221 A1 discloses a temperature distribution determination apparatus for determining a temperature
distribution within a person, while a part of the person is ablated during, for instance,
a radio frequency (RF) ablation procedure. The temperature distribution determination
apparatus is adapted to measure temperature related ultrasound data representing the
temperature at different locations within the person and to use these ultrasound data
for modeling the temperature distribution within the person with a time-dependent
machine-trained model. In an embodiment anatomical ultrasound information is used
to indicate blood vessels within the person, wherein this information is used to correct
the modeled temperature distribution.
[0007] Since different blood vessels, which appear to have the same anatomy as indicated
by the anatomical ultrasound information, may influence the temperature distribution
differently, the finally determined corrected temperature distribution may be inaccurate.
SUMMARY OF THE INVENTION
[0008] The scope of the invention is defined by the appended independent claims. Further
preferred embodiments are disclosed in the dependent claims. It is an object of the
present invention to provide a heat sink parameter determination apparatus, and computer
program for determining a parameter of a heat sink within an object, which allows
for an improved determination of a temperature distribution. It is a further object
of the present invention to provide a temperature distribution determination apparatus,
and computer program for determining a temperature distribution within the object,
which use the determined heat sink parameter and therefore allow for a more accurate
determination of the temperature distribution. Moreover, it is an object of the present
invention to provide a heating system for heating the object, which comprises the
temperature distribution determination apparatus.
[0009] In a first aspect of the present invention a heat sink parameter determination apparatus
for determining a parameter of a heat sink within an object is presented, the heat
sink parameter determination apparatus comprises:
- a heat source parameter providing unit for providing a parameter of a heat source,
which defines the heating of the object, the heat source being an ablation element
for ablating a tumor within the object, and the heat sink being a blood vessel within
the object which defines a cooling of the object,
- a model providing unit for providing a model describing a temperature distribution
in the object, while the object is heated by the heat source, wherein the model depends
on the heat sink parameter and the heat source parameter,
- a heat sink location providing unit for providing a heat sink location and an initial
heat sink parameter for initializing the provided model,
- a temperature distribution measuring unit for measuring a temperature distribution
in the object close to the provided heat sink location, the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane,
- a parameter determination unit for determining the heat sink parameter by starting
the determination of the heat sink parameter with the initialized model, and minimizing
a deviation between only parts of the temperature distribution, which have been measured
within the one or several scan planes and having a shortest distance to a border of
the intersection area that are smaller than a predefined distance from the provided
heat sink location, and a modeled temperature distribution, wherein the modeled temperature
distribution is modeled based on the provided heat source parameter and the heat sink
parameter to be determined by using the provided model.
[0010] Since the heat sink parameter is determined by minimizing a deviation between the
measured temperature distribution and the modeled temperature distribution, wherein
the modeled temperature distribution is modeled based on the provided heat source
parameter and the respective heat sink parameter to be determined by using the provided
model, the determination of the heat sink parameter considers the real temperature
distribution within the object and is thus based on its real influence on the temperature
distribution. This can lead to an improved determination of the heat sink parameter
and hence to a more accurate temperature distribution, which may be determined based
on the determined heat sink parameter.
[0011] The temperature distribution measuring unit can be adapted to determine a spatial
and/or temporal temperature distribution. The temperature distribution measuring unit
comprises preferentially an ultrasound probe for acquiring ultrasound data of the
object and an ultrasound thermometry unit for determining the temperature distribution
based on the acquired ultrasound data. The ultrasound probe is adapted to send ultrasound
waves into the object and receive backscattered ultrasound waves from the object,
in order to acquire the ultrasound data. In an embodiment the heat sink is a blood
vessel and the ultrasound probe can be able to acquire the ultrasound data in three
dimensions, wherein for determining the temperature distribution ultrasound data of
one or several planes, which may traverse the blood vessel, may be used. By using
the ultrasound probe, which can acquire the ultrasound data in three dimensions, the
placement of the one or several planes traversing the vessel can be very flexible.
[0012] The heat sink location providing unit can be adapted to determine the heat sink location
based on ultrasound data of the object acquired by the ultrasound probe. In particular,
the ultrasound probe can be adapted to measure ultrasound data in three dimensions,
in order to provide a volume data set in which the heat sink location providing unit
can detect the heat sink, especially a blood vessel, by using, for instance, known
segmentation techniques. The ultrasound probe can therefore be used for at least two
purposes, measuring the temperature distribution and determining the heat sink location,
which can lead to a more compact heat sink determination apparatus.
[0013] In an embodiment, for determining the heat sink parameter the temperature distribution
measuring unit may measure the temperature distribution only close to the provided
heat sink location, wherein the measured temperature distribution may be used for
determining the heat sink parameter. In an embodiment for determining the heat sink
parameter the temperature distribution may be measured close to the provided heat
sink location and also at a larger distance from the heat sink location, wherein for
determining the heat sink parameter only the temperature distribution may be used,
which has been measured close to the heat sink location. The heat source parameter
providing unit may be adapted to determine one or several heat source parameters,
which define the heating of the object, and to provide the determined one or several
heat source parameters, or the heat source parameter providing unit may just be a
storing unit, in which previously determined and/or received one or several heat source
parameters are stored and from which the one or several heat source parameters can
be retrieved for providing the same. The one or several heat source parameters include,
for instance, the location of the heat source, the amount of heat delivered to the
object, in particular, at which times which amount of heat is delivered to the object,
et cetera.
[0014] In an embodiment the object is a living being like a person or an animal, the heat
sink is a blood vessel and the heat source is an ablation element for ablating a tissue
region of the living being. In this case the location of the ablation element may
be determined by a localization system like an x-ray localization system, an electromagnetic
(EM) localization system, an optical shape sensing (OSS) localization system, et cetera,
wherein the location determined by the localization system can be sent to the heat
source parameter providing unit, in order to allow the heat source parameter providing
unit to provide the determined location as a heat source parameter. Moreover, a heating
control unit for controlling the heating by the ablation element can be adapted to
send heating information regarding heating times and heating amounts to the heat source
parameter providing unit, in order to allow the heat source parameter providing unit
to provide this information as heat source parameters. The ablation element can be
an RF ablation element, a microwave ablation element, a high intensity focused ultrasound
(HIFU) ablation element or any other ablation element being adapted to heat the object.
[0015] The parameter determination unit can be adapted to determine one or several heat
sink parameters. Moreover, the model providing unit can be adapted to provide a model,
which depends on one or several heat sink parameters and/or which depends on one or
several heat source parameters. In an embodiment the heat sink parameter determination
apparatus is adapted to determine one or several heat sink parameters which include
a location of the heat sink within the object and/or, if the heat sink comprises a
tubular structure with a flowing fluid, a flow direction and/or a diameter of the
tubular structure and/or a flow rate. These heat sink parameters are very well suited
for determining an improved temperature distribution. However, the heat sink parameter
determination apparatus can also be adapted to alternatively or additionally determine
other heat sink parameters.
[0016] In another embodiment the heat sink parameter determination apparatus comprises an
initial heat sink parameter providing unit for providing an initial heat sink parameter
for initializing the provided model, wherein the parameter determination unit is adapted
to start the determination of the heat sink parameter with the initialized model.
The initial heat sink parameter providing unit can be adapted to provide, for instance,
a location of the heat sink and/or, if the heat sink is a tubular structure with a
flowing fluid, a diameter of the tubular structure and/or a flow rate and/or a flow
direction. These initial heat sink parameters can provide a good start for the process
of determining the one or several heat sink parameters by minimizing the deviation
between the measured temperature distribution and the modeled temperature distribution,
which can lead to a faster determination process and to a reduced likelihood that
the determination process finds a local minimum being larger than an overall minimum.
[0017] The heat sink parameter determination apparatus comprises a heat sink location providing
unit for providing a heat sink location, wherein the temperature distribution measuring
unit is adapted to measure the temperature distribution close to the provided heat
sink location. Thus, the temperature distribution measuring unit is preferentially
adapted to measure the temperature distribution adjacent to or in the vicinity of
the provided heat sink location. In particular, a distance to the heat sink location
may be predefined and for determining the heat sink parameter a temperature distribution
can be measured in parts of the object having a distance to the heat sink location,
which is equal to or smaller than the predefined distance, wherein these parts of
the object may be regarded as being close to the heat sink location. If the temperature
distribution is also measured in other parts of the object, preferentially only the
temperature distribution measured for the close parts of the object is used for determining
the heat sink parameter. The predefined distance is preferentially 3 cm, further preferred
2 cm and even further preferred 1 cm. In an embodiment the temperature distribution
may be measured by using ultrasound data acquired by an ultrasound probe in one or
several scan planes traversing the heat sink, especially a blood vessel, such that
the one or several scan planes each comprise an intersection area where the heat sink
intersects the respective scan plane. For determining the heat sink parameter only
parts of the temperature distribution within the one or several scan planes may be
used, which have a shortest distance to the border of the intersection area being
smaller than the predefined distance, wherein these parts of the temperature distribution
may be regarded as being close to the heat sink. Thus, the parts of the temperature
distribution may be used only, which correspond to, for instance, a margin, especially
an annulus, around the respective intersection area, wherein the thickness of the
margin is defined by the predefined distance. In addition, in an embodiment from the
close parts of the temperature distribution only parts may be used for determining
the heat sink parameter, which have been determined based on ultrasound data which
correspond to backscattered ultrasound waves which have not passed through the heat
sink. Ultrasound data, which correspond to backscattered ultrasound waves which have
not passed through the heat sink, are not affected by, for instance, reflections at
the heat sink and can therefore have a better quality, which can lead to a higher
quality temperature distribution. Since, for instance in the case of ultrasound thermometry
but also when using other temperature measurement techniques, the temperature may
be measurable over a limited temperature range only, during the heating process the
temperature distribution measuring unit may be able to measure the temperature close
to the heat sink over a relatively long time in comparison to a temperature measurement
close to the heat source. More temperature data may therefore be measurable close
to the heat sink and may thus be available for determining the heat sink parameter,
which can lead to a further improved determination of the heat sink parameter. In
an embodiment the heat sink location providing unit and the initial heat sink parameter
providing unit are an integrated unit, wherein in this case the heat sink location
is provided as an initial heat sink parameter.
[0018] It is preferred that the model providing unit is adapted to provide the model such
that it depends also on an object parameter, wherein the parameter determination unit
is adapted to determine the object parameter and the heat sink parameter by minimizing
the deviation between the measured temperature distribution and the modeled temperature
distribution, wherein the modeled temperature distribution is modeled based on the
provided heat source parameter, the heat sink parameter to be determined and the object
parameter to be determined by using the provided model. The model providing unit can
be adapted to provide the model such that it depends on one or several object parameters.
The object parameters may include, for instance, thermal parameters and/or electrical
parameters, i.e. parameters defining thermal and/or electrical properties of the object.
The thermal and electrical parameters can be, for instance, thermal and electrical
conductivities especially of tissue, if the object is a person. Determining also one
or several object parameters, which can be used, in addition to the determined one
or several heat sink parameters, for determining the temperature distribution, can
further improve the accuracy of the temperature distribution, which may be determined
based on the one or several heat sink parameters and on the one or several object
parameters.
[0019] In a further aspect of the present invention a temperature distribution determination
apparatus for determining a temperature distribution within an object comprising a
heat sink is presented, wherein the temperature distribution determination apparatus
comprises:
- a heat sink parameter determination apparatus for determining a parameter of the heat
sink as defined in claim 1,
- a temperature distribution determination unit for determining the temperature distribution
within the object based on the model provided by the model providing unit, a heat
source parameter provided by the heat source parameter providing unit and the determined
heat sink parameter.
[0020] By using the determined heat sink parameter and the model the temperature distribution
within the object can be determined, even if the temperature distribution is not measurable
by the temperature distribution measuring unit, for instance, because the temperatures
are outside of the temperature range measurable by the temperature distribution measuring
unit.
[0021] In an embodiment the heat source is adapted to heat the object in a first heating
period to a lower temperature defined by a provided first heat source parameter and
in a second heating period to a higher temperature defined by a provided second heat
source parameter, wherein a) the temperature distribution measuring unit is adapted
to measure the temperature distribution, which is used for determining the heat sink
parameter, during the first heating period, wherein the parameter determination unit
is adapted to determine the heat sink parameter based on the provided model, the provided
first heat source parameter and the measured temperature distribution, and b) the
temperature distribution determination unit is adapted to determine, during the second
heating period, the temperature distribution within the object based on the provided
model, the provided second heat source parameter and the determined heat sink parameter.
Moreover, it is preferred that the temperature distribution determination apparatus
further comprises a temperature sensing element for sensing, in the second heating
period, a temperature of the object within a temperature range in which a temperature
is not measurable by the temperature distribution measuring unit, wherein the temperature
distribution determination unit is adapted to determine, in the second heating period,
the temperature distribution within the object based on the provided model, the provided
second heat source parameter, the determined heat sink parameter and the temperature
sensed by the temperature sensing element. Using the temperature measured in the second
heating period by the temperature sensing element for determining the temperature
distribution in the second heating period can further improve the quality of the determined
temperature distribution.
[0022] Furthermore, the temperature distribution measuring unit can be adapted to measure,
in the second heating period, a temperature distribution in a first spatial region,
in which a temperature distribution is measurable by the temperature distribution
measuring unit, wherein the temperature distribution determination unit can be adapted
to determine, in the second heating period, the temperature distribution within a
second spatial region, in which a temperature distribution is not measurable by the
temperature distribution measuring unit, based on the provided model, the provided
second heat source parameter, the determined heat sink parameter and the temperature
distribution measured in the second heating period in the first spatial region. In
addition to determining the heat sink parameter also an object parameter can be determined
based on the provided model, the provided first heat source parameter and the temperature
distribution measured during the first heating period, wherein in the second heating
period the temperature distribution within the second spatial region can then be determined
based on the provided model, the provided second heat source parameter, the determined
heat sink parameter and the determined object parameter.
[0023] The object preferentially comprises tissue to be ablated by the heat generated by
the heat source, wherein the heat sink parameter determination apparatus may further
comprise an ablated region determination unit for determining a region in which tissue
has been ablated.
[0024] In a further aspect of the present invention a heating system for heating an object
comprising a heat sink is presented, wherein the heating system comprises:
- a heat source for heating the object, and
- a temperature distribution determination apparatus for determining a temperature distribution
within the object caused by heating the object as defined in claim 6.
[0025] In a preferred embodiment the heating system further comprises a heating control
unit for controlling the heat source depending on the determined temperature distribution.
In particular, a region or interest, which should be ablated by applying heat to the
region of interest, like a tumor region can be provided and compared with a determined
ablated region, wherein the ablated region can be determined depending on the determined
temperature distribution. The heating control unit can then be adapted to control
the heating such that the determined ablated region completely covers the region of
interest, in order to ensure that the region of interest has been ablated completely.
In this way the heating of the object can be improved.
[0026] In another aspect of the present invention a computer program for determining a parameter
of a heat sink within an object is presented, the computer program comprising program
code means for causing a heat sink parameter determination apparatus as defined in
claim 1 to carry out the following steps, when the computer program is run on a computer
controlling the heat sink parameter determination apparatus:
- providing a parameter of a heat source (2), which defines a heating of the object,
by a heat source parameter providing unit (14), the heat source being an ablation
element for ablating a tumor within the object, and the heat sink being a blood vessel
within the object which defines a cooling of the object,
- providing a model describing a temperature distribution in the object, while the object
is heated by the heat source (2) by a model providing unit (15), wherein the model
depends on the heat sink parameter and the heat source parameter,
- providing a heat sink location and an initial heat sink parameter for initializing
the provided model by a heat sink location providing unit (17),
- measuring a temperature distribution in the object close to the provided heat sink
location by a temperature distribution measuring unit (13, 71), the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane,
- determining the heat sink parameter by starting the determination of the heat sink
parameter with the initialized model, and minimizing a deviation between only parts
of the temperature distribution, which have been measured within the one or several
scan planes and having a shortest distance to a border of the intersection area that
are smaller than a predefined distance from the provided heat sink location, and a
modeled temperature distribution by a parameter determination unit (16), wherein the
modeled temperature distribution is modeled based on the provided heat source parameter
and the heat sink parameter to be determined by using the provided model.
[0027] In a further aspect of the present invention a computer program for determining a
temperature distribution within an object comprising a heat sink is presented, the
computer program comprising program code means for causing a temperature distribution
determination apparatus as defined in claim 6 to carry out the following steps, when
the computer program is run on a computer controlling the temperature distribution
determination apparatus:
- determining a parameter of the heat sink by:
- providing a parameter of a heat source (2), which defines a heating of the object,
by a heat source parameter providing unit (14), the heat source being an ablation
element for ablating a tumor within the object, and the heat sink being a blood vessel
within the object which defines a cooling of the object,
- providing a model describing a temperature distribution in the object, while the object
is heated by the heat source (2) by a model providing unit (15), wherein the model
depends on the heat sink parameter and the heat source parameter,
- providing a heat sink location and an initial heat sink parameter for initializing
the provided model by a heat sink location providing unit (17),
- measuring a temperature distribution in the object close to the provided heat sink
location by a temperature distribution measuring unit (13, 71), the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane,
- determining the heat sink parameter by starting the determination of the heat sink
parameter with the initialized model, and minimizing a deviation between only parts
of the temperature distribution, which have been measured within the one or several
scan planes and having a shortest distance to a border of the intersection area that
are smaller than a predefined distance from the provided heat sink location, and a
modeled temperature distribution by a parameter determination unit (16), wherein the
modeled temperature distribution is modeled based on the provided heat source parameter
and the heat sink parameter to be determined by using the provided model, and
- determining the temperature distribution within the object based on the model provided
by the model providing unit (15), a heat source parameter provided by the heat source
parameter providing unit (14) and the determined heat sink parameter by a temperature
distribution determination unit (22).
[0028] It shall be understood that the heat sink parameter determination apparatus of claim
1, the temperature distribution determination apparatus of claim 6, the heating system
of claim 9, the computer program for determining a parameter of a heat sink of claim
11 and the computer program for determining a temperature distribution of claim 12
have similar and/or identical preferred embodiments, in particular, as defined in
the dependent claims.
[0029] These and other aspects of the invention will be apparent from and elucidated with
reference to the embodiments described hereinafter.
BRIEF DESCRIPTION OF THE DRAWINGS
[0030] In the following drawings:
Fig. 1 shows schematically and exemplarily a heating system for heating an object
comprising a heat sink,
Fig. 2 shows schematically and exemplarily a tip of an ablation needle,
Figs. 3 to 6 schematically and exemplarily show mean square difference plots depending
on different parameters,
Fig. 7 illustrates measured temperatures and modeled temperatures,
Fig. 8 shows a flowchart exemplarily illustrating an embodiment of a temperature distribution
method for determining a temperature distribution within the object.
DETAILED DESCRIPTION OF EMBODIMENTS
[0031] Fig. 1 shows schematically and exemplarily an embodiment of a heating system for
heating an object comprising a heat sink. In this embodiment the system 1 is an ablation
system for ablating a tumor within a person 3 lying on a support means 4 like a patient
table. The system 1 comprises a heating source 2 being an energy application element
for applying energy to the person 3, in particular, to a tumor within the person 3.
In this embodiment the energy application element 2 is an ablation needle comprising
an ablation electrode 40 and a temperature sensing element 41 at the tip 5 of the
ablation needle 2 as schematically and exemplarily illustrated in Fig. 2. The temperature
sensing element 41 at the tip 5 of the ablation needle 2 is a thermocouple in this
embodiment, which is electrically connected to a tip temperature measurement determining
unit 18 via an electrical conductor 42 like a wire for determining the temperature
at the tip 5 of the ablation catheter 2 depending on electrical signals received from
the thermocouple.
[0032] The energy applied to the person 3 by the ablation electrode 40 is preferentially
RF energy, wherein the ablation electrode 40 is electrically connected to an ablation
energy control unit 12, which may also be regarded as being a heating control unit
12 and which is adapted to control the application of the RF energy, via an electrical
conductor 43 which may also be a wire. The electrical conductors 42, 43 may be arranged
within the ablation needle 2 and within a cable 72 connecting the ablation needle
2 with the ablation energy control unit 12. In this embodiment the ablation energy
control unit 12 comprises an RF source for providing the RF energy.
[0033] The system 1 further comprises a temperature distribution measuring unit for measuring
a spatially and temporally dependent temperature distribution within the person 3,
while the RF energy is applied to the person 3. The temperature distribution measuring
unit comprises an ultrasound probe 71 for acquiring ultrasound data of the person
3 and an ultrasound thermometry unit 13 for determining the temperature distribution
based on the acquired ultrasound data.
[0034] The system 1 further comprises a model providing unit 15 for providing a model describing
a temperature distribution within the person 3, while the person 3 is heated by using
the ablation needle 2, wherein the model depends on modifiable model parameters including
one or several heat sink parameters, which define a cooling of the person 3 while
the person 3 is heated by the ablation needle 2, and one or several heat source parameters
defining the heating of the person 3 by the ablation needle 2. The modifiable parameters
may further include object parameters like thermal parameters or electrical parameters
of the object. In this embodiment a tumor in the liver 20 of the person 3 should be
ablated such that the modifiable model parameters may include the thermal conductivity
and the electrical conductivity of the liver 20.
[0035] The heat sinks are blood vessels within the liver 20. The heat sink parameters include,
for example, the location of the respective blood vessel, the diameter of the respective
blood vessel, the flow rate through the respective blood vessel, the direction of
flow through the respective blood vessel, et cetera. The heat source parameters include,
for instance, the location of the ablation needle within the person 3, the times at
which energy is applied to the person 3, the amount of energy applied to the person
3, et cetera. The heat source parameters can be provided by a heat source parameter
providing unit 14, which may have received information about the times, at which the
energy is applied, and about the amount of the applied energy from the ablation energy
control unit 12. In particular, there can be a continuous data stream from the ablation
energy control unit 12 to the heat source parameter providing unit 14, in order to
send this information from the ablation control unit 12 to the heat source parameter
providing unit 14. The information about the location of the ablation needle 2 within
the person 3 may be provided to the heat source parameter providing unit 14 by a localization
system 6, which will be described in more detail further below.
[0036] The system 1 further comprises a parameter determination unit 16 for determining
the heat sink parameters by minimizing a deviation between the measured temperature
distribution and a modeled temperature distribution, wherein the modeled temperature
distribution is modeled based on the provided heat source parameters and the heat
sink parameters to be determined by using the provided model. The parameter determination
unit 16 may further consider object parameters like the thermal and electrical conductivities
of the liver tissue during the modeling step, thereby also determining these object
parameters. The parameter determination unit 16 can be adapted to additionally use
the temperature measured by the thermocouple 41 at the tip 5 of the ablation needle
2, wherein in this case the modeled temperature distribution is modeled such that
a deviation between a) the measured temperature distribution, which is measured by
the temperature distribution measuring unit 13, 71, and the temperature measured by
the thermocouple 41 at the tip 5 of the ablation needle 2 and b) the modeled temperature
distribution is minimized.
[0037] The localization system 6 is adapted to detect the position of the tip 5 of the ablation
needle 2 within the person 3. In this embodiment the localization system 6 is an x-ray
fluoroscopy system, in particular, an x-ray C-arm system. The x-ray fluoroscopy system
comprises an x-ray source 7 for generating x-rays 8 which traverse the person 3 on
the support means 4, wherein the x-rays 8, which have traversed the person 3, are
detected by an x-ray detector 9. The x-ray fluoroscopy system 6 further comprises
a fluoroscopy control unit 10 for controlling the x-ray source 7 and the x-ray detector
9. The x-ray detector 9 generates x-ray images of the person 3, which can be shown
on a display 19 and which can be used to determine the position of the tip 5 of the
ablation needle 2.
[0038] The x-ray fluoroscopy system is preferentially an x-ray C-arm system, which allows
acquiring x-ray projection images in different acquisition directions, wherein the
fluoroscopy control unit 10 may be adapted to reconstruct a computed tomography (CT)
image from the x-ray projection images acquired in different acquisition directions.
In the CT image the tip 5 of the ablation needle 2 may be segmented, in order to determine
its position.
[0039] The system 1 preferentially further comprises a heat sink location providing unit
17 for providing a heat sink location based on the reconstructed CT image, wherein
the temperature distribution measuring unit 13, 71 can be adapted to measure the temperature
distribution, which is used for determining the heat sink parameters, close to the
provided heat sink location. For instance, the heat sink location providing unit 17
can be adapted to segment blood vessels within the liver 20 in the reconstructed CT
image, in order to provide the locations of the blood vessels. These locations of
the blood vessels can also be used as initial heat sink parameters, which can be used
to initialize the model before the parameter determination unit 16 determines the
heat sink parameters such that the process of determining the heat sink parameters
is started with the initialized model. The heat sink location providing unit 17 can
therefore also be regarded as being an initial heat sink parameter providing unit.
In other embodiments the heat sink location providing unit 17 can be adapted to determine
the heat sink locations in another way, for instance, based on other images like a
CT image generated by a CT system, a magnetic resonance (MR) image generated by an
MR imaging system, an ultrasound image generated by an ultrasound probe, especially
by the ultrasound probe already used for the ultrasound thermometry, et cetera.
[0040] The ultrasound probe 71 may be adapted to measure the ultrasound data for one or
several scan planes, wherein the heat source parameter providing unit 14 may be adapted
to provide the location of the ablation needle 2 relative to the location of the one
or several scan planes. For instance, the CT image can be used to determine the location
of the ablation needle relative to the one or several scan planes by segmenting the
ablation needle 2 and the ultrasound probe 71 in the CT image and by using a known
spatial relationship between the ultrasound probe 71 and the one or several scan planes.
This known spatial relationship can be provided by the ultrasound probe 71. Also other
images can be used for determining the location of the ablation needle 2 relative
to the location of the one or several scan planes like ultrasound images which may
be generated by using the ultrasound probe 71 or by using another ultrasound probe.
[0041] In other embodiments other localization techniques for detecting the position of
the needle tip within the person can be used like localization techniques which are
based on EM sensors, OSS sensors, et cetera.
[0042] The temperature distribution measuring unit 13, 71, the heat source parameter providing
unit 14, the model providing unit 15, the parameter determination unit 16 and the
heat sink parameter providing unit 17 are used for determining one or several heat
sink parameters, i.e. in this embodiment parameters of blood vessels, within the person
3. These units can therefore be regarded as forming a heat sink parameter determination
apparatus for determining parameters of heat sinks within the person 3.
[0043] The system 1 further comprises a temperature distribution determination unit 22 for
determining a temperature distribution within the person 3 based on the provided model,
a heat source parameter provided by the heat source parameter providing unit 14 and
the determined heat sink parameters. In particular, the ablation needle 2 is preferentially
used to heat the person 3 in a first heating period to a lower temperature defined
by first heat source parameters provided by the heat source parameter providing unit
14 and in a second heating period to a higher temperature defined by second heat source
parameters provided by the heat source parameter providing unit 14. In the second
heating period the temperature is high enough to ablate the tumor region within the
liver 20. The temperature distribution measuring unit 13, 71 is preferentially used
to measure the temperature distribution, which is used for determining the heat sink
parameters, during the first heating period, wherein the parameter determination unit
16 determines the heat sink parameters based on the provided model, the provided first
heat source parameters and the measured temperature distribution. Moreover, the temperature
distribution determination unit 22 preferentially determines, during the second heating
period, the temperature distribution within the person 3 based on the provided model,
the provided second heat source parameters and the determined heat sink parameters.
Moreover, the temperature distribution measuring unit 13, 71 may be adapted to measure,
in the second heating period, a temperature distribution in a first spatial region,
in which a temperature distribution is measurable by the temperature distribution
measuring unit 13, 71, wherein the temperature distribution determination unit 22
may be adapted to determine, in the second heating period, the temperature distribution
within a second spatial region, in which the temperature distribution is not measurable
by the temperature distribution measuring unit 13, 71, based on the provided model,
the provided second heat source parameters, the determined heat sink parameters and
the temperature distribution measured in the second heating period in the first spatial
region. For instance, the second spatial region may be a region close to the location
of the ablation needle 2. Since in this example the second spatial region is close
to the location of the ablation needle 2, the temperatures within the second spatial
region will be relatively high, especially higher than 50 degrees Celsius, and therefore
not be measurable by ultrasound thermometry. The first spatial region has preferentially
a distance to the ablation needle 2, which is large enough to have temperatures within
a temperature range, which is measurable by ultrasound thermometry. The temperature
distribution in the second spatial region can then be determined by using the provided
model, the provided second heat source parameters, the determined heat sink parameters
and the temperature distribution measured in the first spatial region by adapting
the model with the determined heat sink parameters such that a deviation between the
temperature distribution measured in the first spatial region and the modeled temperature
distribution in the first spatial region is minimized. The temperature distribution
provided by the adapted model in the second spatial region is then regarded as being
the temperature distribution in this second spatial region. Thus, a temperature distribution
can be provided in the second spatial region, although the temperature distribution
cannot directly be measured in the second spatial region. The modeling of the temperature
distribution can further consider object parameters like thermal and electrical tissue
conductivities, which may have been determined based on the temperature measurements
performed during the first heating period. Moreover, the modeling can consider temperatures
measured during the second heating period by the thermocouple 41 at the tip 5 of the
ablation needle 2, wherein in this case the modeling can be performed such that a
deviation between a) the temperature distribution measured in the first spatial region
and the temperatures measured by the thermocouple 41 and b) the modeled temperature
distribution in the first spatial region and at the position of the thermocouple 41
is minimized.
[0044] The first heating period is preferentially a heating period, in which the tissue
of the person 3 is heated to a temperature, which is measurable by ultrasound thermometry,
wherein in this low temperature phase heat sink parameters and optionally also object
parameters like thermal and electrical conductivities are determined. During the second
heating period the tissue of the person 3 is preferentially heated to temperatures,
which can be used to ablate regions of the tissue, in particular, a tumor region,
wherein in this case the temperature distribution close to the ablation needle can
be modeled as described above.
[0045] The first spatial region has preferentially a distance to the tip 5 of the ablation
needle 2 such that the temperature in the first spatial region will be smaller than
about 50 degrees Celsius, if the ablation energy is applied to the person 3 in the
second heating period. This ensures that the temperature distribution measuring unit
13, 71 can measure the temperature distribution also during the ablation procedure
in the first spatial region. The second spatial region is closer to the tip 5 of the
ablation needle 2. Preferentially, the second spatial region is adjacent to the tip
5 of the ablation needle 2 and covers a region of interest to be ablated and a surrounding
region surrounding the region of interest. In this example the region of interest
is a tumor region within the liver 20. Thus, the heat sink parameter determination
apparatus 23, 71 and the temperature distribution determination unit 22 are preferentially
adapted such that a temperature distribution can be estimated in the tumor region
and in a surrounding region around the tumor region, wherein these regions may define
the second spatial region. The second spatial region can also be a larger region,
which also covers further regions, which are further away from the tip 5 of the ablation
needle 2 and in which the tissue is also heated to a temperature higher than 50 degrees
Celsius.
[0046] The heat sink parameter determination unit 23, 71 and the temperature distribution
determination unit 22 are part of a temperature distribution determination apparatus
21, 71 for determining a temperature distribution within the person 3, especially
for determining an overall temperature distribution within the first and second spatial
regions.
[0047] The temperature distribution determination apparatus further comprises an ablated
region determination unit 24 for determining an ablated region defining a region within
the person 3 that has been ablated, wherein the ablated region determination unit
24 is adapted to determine the ablated region by determining a part of the person
3 for which the estimated temperature distribution in the second spatial region comprises
or has comprised a temperature being higher than a predefined temperature threshold.
The temperature distribution determination apparatus 21, 71 also comprises a region
of interest providing unit 25 for providing the region of interest being, in this
embodiment, a tumor region, which should be ablated, wherein the determined ablated
region and the tumor region can be shown on the display 19. For instance, an overlay
of the determined ablated region and the tumor region can be shown on the display
19. The predefined temperature threshold is, for instance, 60, 65 or 70 degrees Celsius.
[0048] The ablation energy control unit 12 can be adapted to control the ablation needle
2, i.e. the power of the ablation, depending on the determined temperature distribution.
In particular, the ablation energy control unit 12 can be adapted to control the ablation
power such that the tumor region is completely ablated.
[0049] In another embodiment the modeling may consider temperatures measured during the
second heating period by the thermocouple 41 at the tip 5 of the ablation needle 2,
without considering a thermometry measurement, wherein in this case the modeling can
be performed such that a deviation between a) the temperatures measured by the thermocouple
41 and b) the modeled temperature distribution at the position of the thermocouple
41 is minimized. In particular, in an embodiment the only temperature measurement
being input in the modeling during the second heating period may be the temperature
measurement of the thermocouple 41. Also in this case the ablated region determination
unit 24 can be used to determine an ablated region based on the determined temperature
distribution.
[0050] In this embodiment the ablation needle 2 is navigated directly by hand. In another
embodiment the system can further comprise a navigation unit for navigating the ablation
needle, in particular the needle tip, to a desired location within the person. The
navigation unit can be adapted to allow a user to navigate the ablation needle completely
by hand or semi-automatically. The ablation needle may comprise built-in guiding means,
which can be controlled by the navigation unit. The ablation needle can, for example,
be steered and navigated by the use of steering wires, in order to guide the needle
tip to a desired location within the person.
[0051] Thermal ablation techniques are excellent alternatives to major surgery, which can
pose a risk even with the most experienced surgeon. These techniques are minimally
invasive requiring only needles, which may be adapted to perform an RF therapy, a
cryotherapy or a microwave ablation therapy, or they are non-invasive, wherein, for
instance, a non-invasive heat source such as an ultrasound heating source like a high
intensity focused ultrasound (HIFU) source is used. In most of the procedures, cancerous
tissue is heated to temperatures above 60 degrees Celsius and coagulated.
[0052] For performing an RF ablation (RFA) procedure the system described above with reference
to Fig. 1 comprises a probe with an active electrode tip, i.e. the ablation needle,
through which preferentially a 460 to 500 kHz alternating current is conducted. The
current propagates through the body of the person 3 to grounding pads (not shown in
Fig. 1 for clarity reasons) placed either on the back or the thigh of the person 3.
The current causes ionic agitation and frictional heating. Heat is then dissipated
through thermal conduction to ablate the tumor region. In this embodiment RFA is used
to treat liver cancer.
[0053] In the embodiment described above with reference to Fig. 1 RFA is performed under
x-ray guidance by using an x-ray C-arm system. However, the RFA can also be performed
by using another guidance system, which may be based on ultrasound imaging, CT imaging
or MR imaging guidance. A follow-up examination is preferentially done by using a
CT scan or MRI scan within, for example, a month to assess effectiveness of ablation
and again at three months intervals along with tumor markers to detect residual disease
or recurrence. After state of the art ablation procedures have been performed, relatively
high recurrence rates are often observed because of the often present inability to
monitor and control ablation size sufficiently to adequately kill the tumor cells.
The system described above with reference to Fig. 1 provides therefore real-time feedback
to the clinician by providing a temperature map of the ablated zone. This could also
be achieved with reasonable accuracy with MR based temperature imaging. However, MR
imaging is expensive and may not be readily available. Ultrasound is another modality
that may be used for image guidance during placement of the needle. Due to its ease
of use and availability it may be a preferred method for monitoring the lesions. However,
in the prior art ultrasound is used generally for monitoring the treatment by visualizing
the hyperechoic lesions on a B-mode image. Such visualization is only approximate
and not a good indicator of the treatment efficacy.
[0055] The underlying principle of ultrasound thermometry is that the speed of sound in
tissue changes as a function of temperature which manifests as apparent shifts, i.e.
displacements, in ultrasound echoes. The resulting "temperature-induced strain", which
is mathematically derived by differentiating the displacement along the ultrasound
beam direction, is nominally proportional to the temperature rise in the range up
to 50 degrees Celsius. However, the problem lies in the variation in trend of the
temperature dependence of speed of sound for various tissues. For example, for liver
tissues the speed of sound increases approximately linearly with temperature up to
a temperature range of 50 degrees Celsius, after which the trend plateaus. Hence,
there is no sensitivity to ultrasound echo shifts with temperatures beyond this temperature
range. Also, with the onset of tissue necrosis and the resulting changes in tissue
structure, the signature of the ultrasound echoes changes significantly and makes
the comparison of ultrasound echoes to determine the displacement difficult. Therefore,
for temperatures above 50 degrees Celsius the ultrasound thermometry, which is based
on tracking changes in speed of sound, is not a reliable indicator of temperature
in the tissue.
[0056] In the first heating period the tissue is therefore preferentially heated up to 50
degrees Celsius only, at least in the spatial regions where the ultrasound data for
determining the temperature distribution are measured, wherein these ultrasound data
are preferentially measured in a plane close to a blood vessel. The close plane may
be parallel to the blood vessel or traverse the blood vessel. Moreover, during the
second heating period the first spatial region, in which ultrasound data may be measured,
is preferentially chosen such that the temperature is not higher than 50 degrees Celsius
within the first spatial region.
[0057] Organs such as the liver, to which cancer-related ablative therapies are often applied,
are highly vascularized, wherein blood vessels, especially blood vessels having a
diameter being larger than 3 mm, impact the efficacy of ablative treatments such as
RF ablation. Ablation planning and ablation monitoring can therefore be improved,
if the effects of these blood vessels are considered. Thus, the system described above
with reference to Fig. 1 characterizes their heat sink effect and determines how it
will affect the surrounding tissue and may result in incomplete ablations. This information
may also be used to adapt therapy delivery parameters, before the entire ablation
procedure has been completed. For instance, before the ablation is performed in the
second heating period, in the first heating period the heat sink parameters being
indicative of the heat sink effect can be determined, wherein these heat sink parameters
can be used to determine the temperature distribution also close to the ablation needle
during the second heating period, in which the tumor region should be ablated, wherein
during the ablation procedure heat source parameters like the thermal dose applied
to the tissue can be modified in real-time based on a currently estimated temperature
distribution close to the ablation needle.
[0058] The system described above with reference to Fig. 1 is preferentially adapted to
determine heat sink parameters defining, for instance, the morphology, i.e. geometric
parameters such as position and size, of heat sinks in the vicinity of an ablative
heat source, in order to determine its potential detrimental effect on the therapy
efficacy. The system is preferentially adapted to acquire ultrasound data in the vicinity
of the respective heat sink during a low power heating exercise, i.e. during the first
heating period, where after the ultrasound data are coupled to the thermal model that
includes the ability to model heat sinks to estimate model parameters including at
least heat sink parameters. The thermal model is then preferentially used to infer
temperatures over the entire volume including the core of the heating zone. Finally,
an estimated temperature map, i.e. a temporal and spatial temperature distribution,
is preferentially provided to a physician and the extent of the ablation zone is preferentially
determined based on the estimated temperature map.
[0059] The system described above with reference to Fig. 1 is preferentially configured
to adapt the thermal model, which preferentially simulates ablation heating in the
presence of blood vessels, with real-time temperature measurements obtained by using
ultrasound thermometry. Preferentially, the temperature distribution over the entire
ablation zone is estimated and an undertreatment caused by the presence of blood vessels
is extracted. The system preferentially solves following problems of the prior art.
[0060] The extent of an ablated region may be determined more accurately. Moreover, in the
prior art ultrasound B-mode inspection guided by hyperechoic visualization of the
ablated region is often not accurate, which may render difficult to assess the effectiveness
of therapy. The hyperechoes visualized on B-mode images are caused by gas and vapor
bubbles. In order to generate these bubbles and visualize the treatment region on
ultrasound, an ablation treatment protocol involves heating to temperatures in the
order of 100 degrees Celsius which is overkill for achieving necrosis that only requires
temperatures up to 70 to 80 degrees Celsius. Hence, if ultrasound B-mode imaging is
used for visually monitoring the ablation procedure, the treatment time is longer
than it needs to be. Furthermore, known non-invasive methods do not perform a functional
evaluation of the impact of cooling by blood vessels. While, for instance, CT or ultrasound
Doppler images can provide anatomical information in terms of the location and size
of the blood vessels and the velocities of the blood flowing through the blood vessels,
the thermal impact of these structures is not determined and considered by these known
non-invasive methods. In contrast, the system described above with reference to Fig.
1 evaluates the functional impact of cooling by blood vessels and considers this functional
impact while determining the temperature distribution during the ablation process,
i.e. during the second heating period. The effectiveness of readily available ultrasound
data can therefore be improved during RF ablation treatments and an instant feedback
on the treatment can be provided.
[0061] The system described above with reference to Fig. 1 is adapted to utilize ultrasound
thermometry and estimate the temperature distribution over the entire volume in the
presence of heat sinks such as blood vessels. The system comprises an ultrasound acquisition
system, i.e. the temperature distribution measuring unit, that provides thermometry
estimates preferentially during a low power heating pulse, i.e. during the first heating
period, immediately before the start of the ablation procedure, i.e. immediately before
the second heating period. The ultrasound data are preferentially acquired in one
or several scan planes, wherein two-dimensional ultrasound transducer matrix arrays
may be used, around the heat source being, in this embodiment, the ablation needle
2. Preferably, the one or several scan planes are located close to the blood vessels
so that their cooling effect can be optimally captured. They may be parallel to the
blood vessel or traverse the blood vessel. A thermal model is provided, which can
be a multiphysics model and which is configured to simulate the energy deposition
in tissue based on the used modality which in this embodiment is an RF ablation needle,
but which may also be another heat source like a microwave heat source or a HIFU heat
source. The model is further configured to simulate thermal diffusion resulting from
the heat generation by the heat source and the heat transfer in fluids like blood.
Tissue specific properties such as thermal and electrical parameters and parameters
of the heat sinks like their positions, diameter, et cetera, are preferentially assumed
as unknown parameters and are preferentially estimated
in situ using the ultrasound thermometry data measured by the temperature distribution measuring
unit. Optionally, if anatomical information is available from pre-acquired imaging
data, this information may be used to initialize the locations of the blood vessels
in the model and, if available from, for instance, ultrasound Doppler data, the flow
direction may be initialized in the model. The anatomical information may also be
used to position the ultrasound scan planes for performing the thermometry close to
the blood vessels, wherein the scan planes may be parallel to the blood vessel or
traverse the blood vessel. The parameter determination unit, which may also be regarded
as being an analysis module, is adapted to estimate then the relevant blood vessel
parameters like the position, diameter, flow rate, et cetera, and preferentially to
also estimate tissue parameters like electrical and thermal parameters. Finally, these
parameters maybe used to determine a volumetric temperature distribution optionally
along with the applied thermal dose and/or a lesion size by using the model.
[0062] The system described above with reference to Fig. 1 is adapted to estimate the relevant
parameters in the thermal model that characterize the heat sinks, i.e. in this embodiment
the heat sink parameters characterizing the blood vessels, in the presence of an ablation
heat source, which is the ablation needle in this embodiment. Once the heat sink parameters
and optionally further parameters like tissue parameters have been estimated in the
thermal model, the volumetric temperature distribution can be estimated by a forward
simulation.
[0064] In these equations
V denotes the electrical potential distribution,
T denotes the temperature distribution,
σ denotes the electrical conductivity,
ρ denotes the density,
C denotes the specific heat,
t denotes the time,
k denotes the thermal conductivity and
v denotes the flow velocity in a blood vessel. The variables without subscripts represent
variables of tissue and the variables with the subscript
f represent fluid variables. The density
ρ and the specific heat
C are preferentially assumed to be known and may be based on typical values known from
literature. The thermal conductivity
k and the electrical conductivity
σ are preferentially considered as unknown. Moreover, the heat sink parameters like
the locations and diameters of the blood vessels and the flow rates of the blood flowing
through the blood vessels are also considered as being unknown. These unknown parameters
are determined by the parameter determination unit from a comparison between the model
and measured ultrasound thermometry data. In an embodiment it is assumed that the
electrical conductivity
σ is independent of the temperature such that equation (1) becomes:

[0065] The simulation may consider two different kinds of spatial regions like blood regions
and tissue regions, wherein parameters like the thermal conductivity, the density,
the specific heat and the electrical conductivity may be different in the two different
kinds of spatial regions. However, in an embodiment the simulation may also consider
more different kinds of spatial regions. Moreover, the simulation may consider a tissue
property distribution and a corresponding distribution of tissue parameters as a function
of space.
[0066] The parameter determination unit can be adapted to run the thermal model for different
combinations of the unknown parameters, i.e. of the heat sink parameters and optionally
also of object parameters like electrical and thermal tissue parameters. For the electrical
conductivity
σ a unity value may be used during the simulation with the different combinations,
because the electrical conductivity
σ just linearly affects the temperature distribution as can be seen in above equation
(2) and it can be estimated by a simple scaling operation. For each combination of
parameters a temperature distribution can be modeled and stored in a memory for the
respective ultrasound scan plane. The sweep range for the unknown parameters is preferentially
chosen to be physically meaningful. In order to determine the best combination of
parameters, the parameter determination unit can compare the experimental data, i.e.
the measured temperature distributions, with the stored modeled temperature distribution
and perform a linear fit to estimate the unknown parameter
σ. The comparison can be based on a similarity measure like a mean square difference
(MSD). The parameter determination unit can be adapted to determine the combination
of parameters providing the largest similarity between the respective modeled temperature
distribution and the measured temperature distribution.
[0067] Figs. 3 to 6 schematically and exemplarily show one-dimensional MSD plots for a respective
parameter, wherein the respective other parameters are constant and have values leading
to a respective smallest MSD. In particular, Fig. 3 shows schematically and exemplarily
the MSD depending on the thermal conductivity
k, Fig. 4 shows schematically and exemplarily the MSD depending on a dimension
d of the location of a blood vessel, Fig. 5 schematically and exemplarily shows the
MSD depending on the radius
r of the blood vessel and Fig. 6 schematically and exemplarily shows the MSD depending
on a flow rate
f of a flow of blood through the blood vessel, wherein in each of these plots the parameters
are shown in arbitrary units. Fig. 7 shows schematically and exemplarily a temperature
rise Δ
T along a spatial line along which temperatures have been measured by ultrasound thermometry.
The different positions along the spatial line are indicated by
m , the circles 50 indicate the measured temperatures and the curve 51 represents a
temperature rise as defined by the thermal model with the determined parameters. As
can be seen in Fig. 7, the thermal model with the determined parameters very well
fits the measured temperatures. After the parameter determination unit has determined
these parameters, the temperature distribution determination unit can use these parameters
together with the thermal model, in order to determine a volumetric temperature distribution.
[0068] In the following an embodiment of a temperature distribution determination method
for determining a temperature distribution within an object comprising a heat sink
will exemplarily be described with reference to a flowchart shown in Fig. 8.
[0069] After the heat source has been placed within the object, in step 101 in a first heating
period the object is heated to a relatively low temperature by the heat source, while
a temperature distribution is measured by the temperature distribution measuring unit.
In particular, the heat source is an ablation element placed within a liver tumor
to be ablated, wherein the ablation element is used to heat the liver tumor to a relatively
low temperature, which is preferentially smaller than 50 degrees Celsius, while the
temperature distribution is preferentially measured by ultrasound thermometry.
[0070] In step 102 parameters of the heat source, which define the heating of the object,
are provided by the heat source parameter providing unit, and a thermal model describing
a temperature distribution in the object, while the object is heated by the heat source,
is provided by the model providing unit, wherein the model depends on heat sink parameter
and the heat source parameters. Preferentially, the heat sink parameters include the
location of a blood vessel, the diameter of a blood vessel, the blood flow direction
and the blood flow rate. The heat source parameters include preferentially the times
at which the liver tumor is heated, the amount of heat applied to the liver tumor
and the location of the heat source.
[0071] In step 103 the heat sink parameters are determined by minimizing a deviation between
the measured temperature distribution and a modeled temperature distribution by the
parameter determination unit, wherein the modeled temperature distribution is modeled
based on the provided heat source parameters and the heat sink parameters to be determined
by using the provided model. Steps 101 to 103 can be regarded as being steps of a
heat sink parameter determination method for determining a parameter of a heat sink
within an object.
[0072] In step 104 in a second heating period the object is heated to a higher temperature
and the temperature distribution within the object is determined based on the model
provided by the model providing unit, a heat source parameter defining the heating
during the second heating period provided by the heat source parameter providing unit
and the determined heat sink parameters by the temperature distribution determination
unit. In particular, in step 104 an ablation procedure is performed, wherein during
the ablation procedure the temperature distribution is determined by using the model
with the determined heat sink parameters, in order to monitor the ablation process.
[0073] Although in above described embodiments an ultrasound probe measures ultrasound data
in one or several scan planes, in order measure a temperature distribution by thermometry,
in another embodiment the ultrasound probe may be adapted to measure the ultrasound
data in three dimensions, in order to measure a three-dimensional temperature distribution
by thermometry.
[0074] Although in above described embodiments the ablation element comprises a single ablation
electrode, in other embodiments the ablation element can also comprise more ablation
electrodes. Moreover, although in above described embodiments the ablation element
comprises a single temperature sensing element, i.e. in the described embodiments
a single thermocouple, in other embodiments the ablation element can also comprise
several temperature sensing elements or no temperature sensing element.
[0075] Although in above described embodiments the ablation element is an RF ablation element,
in other embodiments other ablation elements may be used for ablating tissue like
HIFU ablation elements, microwave ablation elements or laser ablation elements.
[0076] Other variations to the disclosed embodiments can be understood and effected by those
skilled in the art in practicing the claimed invention, from a study of the drawings,
the disclosure, and the appended claims.
[0077] In the claims, the word "comprising" does not exclude other elements or steps, and
the indefinite article "a" or "an" does not exclude a plurality.
[0078] A single unit or device may fulfill the functions of several items recited in the
claims. The mere fact that certain measures are recited in mutually different dependent
claims does not indicate that a combination of these measures cannot be used to advantage.
[0079] Procedures like the determination of the heat sink parameters, the determination
of the temperature distributions, the provision of the model, et cetera performed
by one or several units or devices can be performed by any other number of units or
devices. These procedures and/or the control of the temperature distribution determination
apparatus in accordance with the temperature distribution determination method and/or
the control of the heat sink parameter determination apparatus in accordance with
the heat sink parameter determination method can be implemented as program code means
of a computer program and/or as dedicated hardware.
[0080] A computer program may be stored/distributed on a suitable medium, such as an optical
storage medium or a solid-state medium, supplied together with or as part of other
hardware, but may also be distributed in other forms, such as via the Internet or
other wired or wireless telecommunication systems.
[0081] Any reference signs in the claims should not be construed as limiting the scope.
[0082] The invention, as defined in the appended claims, relates to a heat sink parameter
determination apparatus for determining a parameter of a heat sink like a blood vessel
within an object such as a person by minimizing a deviation between a measured temperature
distribution, which has preferentially been measured by ultrasound thermometry, and
a modeled temperature distribution, wherein the modeled temperature distribution is
modeled based on a provided heat source parameter like the location of an ablation
needle and the heat sink parameter to be determined by using a given thermal model.
This determination of heat sink parameters, which may be geometric and/or flow parameters,
considers the real temperature distribution and is thus based on real heat sink influences
on the temperature distribution. This can lead to an improved determination of heat
sink parameters and hence to a more accurate temperature distribution which may be
determined based on the determined heat sink parameters.
1. A heat sink parameter determination apparatus for determining a parameter of a heat
sink within an object, the heat sink parameter determination apparatus (23, 71) comprising:
- a heat source parameter providing unit (14) for providing a parameter of a heat
source (2), which defines a heating of the object, the heat source being an ablation
element for ablating a tumor within the object, and the heat sink being a blood vessel
within the object which defines a cooling of the object,
- a model providing unit (15) for providing a model describing a temperature distribution
in the object, while the object is heated by the heat source (2), wherein the model
depends on the heat sink parameter and the heat source parameter,
- a heat sink location providing unit (17) for providing a heat sink location and
an initial heat sink parameter for initializing the provided model,
- a temperature distribution measuring unit (13, 71) for measuring a temperature distribution
in the object close to the provided heat sink location, the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane, and
- a parameter determination unit (16) for determining the heat sink parameter by starting
the determination of the heat sink parameter with the initialized model, and minimizing
a deviation between only parts of the temperature distribution, which have been measured
within the one or several scan planes and having a shortest distance to a border of
the intersection area that are smaller than a predefined distance from the provided
heat sink location, and a modeled temperature distribution, wherein the modeled temperature
distribution is modeled based on the provided heat source parameter and the heat sink
parameter to be determined by using the provided model.
2. The heat sink parameter determination apparatus as defined in claim 1,
wherein the heat sink parameter determination apparatus (23, 71) is adapted to determine
one or several heat sink parameters which include a location of the heat sink within
the object and/or, if the heat sink includes a tubular structure with a flowing fluid,
a flow direction and/or a radius of the tubular structure and/or a flow rate.
3. The heat sink parameter determination apparatus as defined in claim 1,
wherein the model providing unit (15) is adapted to provide the model such that it
depends also on an object parameter and wherein the parameter determination unit (16)
is adapted to determine the object parameter and the heat sink parameter by minimizing
the deviation between the measured temperature distribution and the modeled temperature
distribution, wherein the modeled temperature distribution is modeled based on the
provided heat source parameter, the heat sink parameter to be determined and the object
parameter to be determined by using the provided model.
4. The heat sink parameter determination apparatus as defined in claim 1,
wherein the temperature distribution measuring unit (13, 71) comprises an ultrasound
probe for acquiring ultrasound data of the object and an ultrasound thermometry unit
for determining the temperature distribution based on the acquired ultrasound data.
5. The heat sink parameter determination apparatus as defined in claim 4,
wherein the heat sink location providing unit (17) is adapted to determine the heat
sink location based on ultrasound data of the object acquired by the ultrasound probe.
6. A temperature distribution determination apparatus for determining a temperature distribution
within an object comprising a heat sink, wherein the temperature distribution determination
apparatus (21, 71) comprises:
- a heat sink parameter determination apparatus (23, 71) for determining a parameter
of the heat sink as defined in claim 1,
- a temperature distribution determination unit (22) for determining the temperature
distribution within the object based on the model provided by the model providing
unit (15), a heat source parameter provided by the heat source parameter providing
unit (14) and the determined heat sink parameter.
7. The temperature distribution determination apparatus as defined in claim 6, wherein
the heat source (2) is adapted to heat the object in a first heating period to a lower
temperature defined by a provided first heat source parameter and in a second heating
period to a higher temperature defined by a provided second heat source parameter,
wherein
- the temperature distribution measuring unit (13, 71) is adapted to measure the temperature
distribution, which is used for determining the heat sink parameter, during the first
heating period, wherein the parameter determination unit (16) is adapted to determine
the heat sink parameter based on the provided model, the provided first heat source
parameter and the measured temperature distribution, and
- the temperature distribution determination unit (22) is adapted to determine, during
the second heating period, the temperature distribution within the object based on
the provided model, the provided second heat source parameter and the determined heat
sink parameter.
8. The temperature distribution determination apparatus as defined in claim 7, wherein
the temperature distribution determination apparatus further comprises a temperature
sensing element (41) for sensing, in the second heating period, a temperature of the
object within a temperature range in which a temperature is not measurable by the
temperature distribution measuring unit (13, 71), wherein the temperature distribution
determination unit (22) is adapted to determine, in the second heating period, the
temperature distribution within the object based on the provided model, the provided
second heat source parameter, the determined heat sink parameter and the temperature
sensed by the temperature sensing element (41).
9. A heating system for heating an object comprising a heat sink, the heating system
(1) comprising:
- a heat source (2) for heating the object, and
- a temperature distribution determination apparatus (21, 71) for determining a temperature
distribution within the object caused by heating the object as defined in claim 6.
10. The heating system as defined in claim 9, wherein the heating system (1) further comprises
a heating control unit (12) for controlling the heat source (2) depending on the determined
temperature distribution.
11. A computer program for determining a parameter of a heat sink within an object, the
computer program comprising program code means for causing a heat sink parameter determination
apparatus (23, 71) as defined in claim 1 to carry out the following steps, when the
computer program is run on a computer controlling the heat sink parameter determination
apparatus (23, 71):
- providing a parameter of a heat source (2), which defines a heating of the object,
by a heat source parameter providing unit (14), the heat source being an ablation
element for ablating a tumor within the object, and the heat sink being a blood vessel
within the object which defines a cooling of the object,
- providing a model describing a temperature distribution in the object, while the
object is heated by the heat source (2) by a model providing unit (15), wherein the
model depends on the heat sink parameter and the heat source parameter,
- providing a heat sink location and an initial heat sink parameter for initializing
the provided model by a heat sink location providing unit (17),
- measuring a temperature distribution in the object close to the provided heat sink
location by a temperature distribution measuring unit (13, 71), the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane,
- determining the heat sink parameter by starting the determination of the heat sink
parameter with the initialized model, and minimizing a deviation between only parts
of the temperature distribution, which have been measured within the one or several
scan planes and having a shortest distance to the border of the intersection area
that are smaller than a predefined distance from the provided heat sink location,
and a modeled temperature distribution by a parameter determination unit (16), wherein
the modeled temperature distribution is modeled based on the provided heat source
parameter and the heat sink parameter to be determined by using the provided model.
12. A computer program for determining a temperature distribution within an object comprising
a heat sink, wherein the heat sink parameter defines a cooling of the object while
the object is heated by a heat source (2), the computer program comprising program
code means for causing a temperature distribution determination apparatus (21, 71)
as defined in claim 6 to carry out the following steps, when the computer program
is run on a computer controlling the temperature distribution determination apparatus
(21, 71):
- determining a parameter of the heat sink by:
- providing a parameter of a heat source (2), which defines a heating of the object,
by a heat source parameter providing unit (14), the heat source being an ablation
element for ablating a tumor within the object, and the heat sink being a blood vessel
within the object which defines a cooling of the object,
- providing a model describing a temperature distribution in the object, while the
object is heated by the heat source (2) by a model providing unit (15), wherein the
model depends on the heat sink parameter and the heat source parameter,
- providing a heat sink location and an initial heat sink parameter for initializing
the provided model by a heat sink location providing unit (17),
- measuring a temperature distribution in the object close to the provided heat sink
location by a temperature distribution measuring unit (13, 71) the temperature distribution
measuring unit (13, 71) comprising an ultrasound probe configured to acquire ultrasound
data in one or more scan planes traversing the heat sink, the one or several scan
planes each comprising an intersection area where the heat sink intersects the respective
scan plane,
- determining the heat sink parameter by starting the determination of the heat sink
parameter with the initialized model, and minimizing a deviation between only parts
of the temperature distribution, which have been measured within the one or several
scan planes and having a shortest distance to the border of the intersection area
that are smaller than a predefined distance from the provided heat sink location and
a modeled temperature distribution by a parameter determination unit (16), wherein
the modeled temperature distribution is modeled based on the provided heat source
parameter and the heat sink parameter to be determined by using the provided model,
and
- determining the temperature distribution within the object based on the model provided
by the model providing unit (15), a heat source parameter provided by the heat source
parameter providing unit (14) and the determined heat sink parameter by a temperature
distribution determination unit (22).
1. Ein Gerät zum Ermitteln eines Kühlkörperparameters in einem Objekt, wobei das Gerät
zum Ermitteln von Kühlkörperparametern (23, 71) Folgendes umfasst:
- eine Einheit zum Bereitstellen von Wärmequellenparametern (14) zum Bereitstellen
eines Parameters einer Wärmequelle (2), der das Erwärmen des Objekts definiert, wobei
es sich bei der Wärmequelle um ein Ablationselement zum Abtragen eines Tumors im Objekt
handelt, und wobei es sich beim Kühlkörper um ein Blutgefäß im Objekt handelt, das
die Kühlung des Objekts definiert,
- eine Modellbereitstellungseinheit (15) zum Bereitstellen eines Modells, das die
Temperaturverteilung im Objekt beschreibt, während das Objekt von der Wärmequelle
(2) erwärmt wird, wobei das Modell vom Kühlkörperparameter und vom Wärmequellenparameter
abhängig ist,
- eine Einheit zum Bereitstellen der Kühlkörperposition (17) zum Bereitstellen der
Kühlkörperposition und eines anfänglichen Kühlkörperparameters zum Initialisieren
des bereitgestellten Modells,
- eine Einheit zum Messen der Temperaturverteilung (13, 71) zum Messen der Temperaturverteilung
im Objekt in der Nähe der bereitgestellten Kühlkörperposition, wobei die Einheit zum
Messen der Temperaturverteilung (13, 71) eine Ultraschallsonde umfasst, die in mindestens
einer Scanfläche Ultraschalldaten erfasst, die den Kühlkörper durchqueren, wobei die
mindestens eine Scanfläche jeweils einen Schnittbereich umfasst, in dem der Kühlkörper
die jeweilige Scanfläche schneidet, und
- eine Parameterermittlungseinheit (16) zum Ermitteln des Kühlkörperparameters durch
Beginnen des Ermittelns des Kühlkörperparameters anhand des initialisierten Modells
sowie durch Minimieren der Abweichung zwischen den Teilen der Temperaturverteilung,
die innerhalb der mindestens einen Scanfläche gemessen wurden, und die den kürzesten
Abstand zur Grenze des Schnittbereichs aufweisen, der kleiner als der vorab festgelegte
Abstand von der bereitgestellten Kühlkörperposition ist, und der modellierten Temperaturverteilung,
wobei die modellierte Temperaturverteilung anhand des bereitgestellten Wärmequellenparameters
und des mithilfe des bereitgestellten Modells zu ermittelnden Kühlkörperparameters
modelliert wird.
2. Das Gerät zum Ermitteln von Kühlkörperparametern gemäß Anspruch 1, wobei das Gerät
zum Ermitteln von Kühlkörperparametern (23, 71) einen oder mehrere Kühlkörperparameter
ermittelt, die die Position des Kühlkörpers im Objekt und/oder, sofern der Kühlkörper
eine röhrenförmige Struktur mit einer strömenden Flüssigkeit umfasst, die Strömungsrichtung
und/oder den Radius der röhrenförmigen Struktur und/oder die Strömungsrate umfassen.
3. Das Gerät zum Ermitteln von Kühlkörperparametern gemäß Anspruch 1,
wobei die Modellbereitstellungseinheit (15) das Modell so bereitstellt, dass dieses
zudem von einem Objektparameter abhängig ist, und wobei die Parameterermittlungseinheit
(16) den Objektparameter und den Kühlkörperparameter durch Minimieren der Abweichung
zwischen der gemessenen Temperaturverteilung und der modellierten Temperaturverteilung
ermittelt, wobei die modellierte Temperaturverteilung anhand des bereitgestellten
Wärmequellenparameters, des zu ermittelnden Kühlkörperparameters und des zu ermittelnden
Objektparameters modelliert wird.
4. Das Gerät zum Ermitteln von Kühlkörperparametern gemäß Anspruch 1,
wobei die Einheit zum Messen der Temperaturverteilung (13, 71) eine Ultraschallsonde
zum Erfassen von Ultraschalldaten des Objekts sowie eine Ultraschall-Thermometrieeinheit
zum Ermitteln der Temperaturverteilung anhand der erfassten Ultraschalldaten umfasst.
5. Das Gerät zum Ermitteln von Kühlkörperparametern gemäß Anspruch 4,
wobei die Einheit zum Ermitteln der Kühlkörperposition (17) die Position des Kühlkörpers
anhand der von der Ultraschallsonde erfassten Ultraschalldaten für das Objekt ermittelt.
6. Ein Gerät zum Ermitteln der Temperaturverteilung zum Ermitteln der Temperaturverteilung
in einem Objekt mit einen Kühlkörper umfasst, wobei das Gerät zum Ermitteln der Temperaturverteilung
(21, 71) Folgendes umfasst:
- ein Gerät (23, 71) zum Ermitteln von Kühlkörperparametern gemäß Anspruch 1,
- eine Einheit zum Ermitteln der Temperaturverteilung (22) zum Ermitteln der Temperaturverteilung
in einem Objekt anhand des von der Modellbereitstellungseinheit (15) bereitgestellten
Modells, des von der Einheit zum Bereitstellen von Wärmequellenparametern (14) bereitgestellten
Wärmequellenparameters und des ermittelten Kühlkörperparameters.
7. Das Gerät zum Ermitteln der Temperaturverteilung gemäß Anspruch 6, wobei die Wärmequelle
(2) das Objekt in einer ersten Heizperiode auf eine niedrigere, anhand des bereitgestellten
ersten Wärmequellenparameters definierte Temperatur und in einer zweiten Heizperiode
auf eine höhere, anhand des bereitgestellten zweiten Wärmequellenparameters definierte
Temperatur erwärmt, wobei
- die Einheit zum Messen der Temperaturverteilung (13, 71) während der ersten Heizperiode
die zum Ermitteln des Wärmequellenparameters verwendete Temperaturverteilung misst,
wobei die Parameterermittlungseinheit (16) den Kühlkörperparameter anhand des bereitgestellten
Modells, des bereitgestellten ersten Wärmequellenparameters und der gemessenen Temperaturverteilung
ermittelt, und
- die Einheit zum Ermitteln der Temperaturverteilung (22) während der zweiten Heizperiode
die Temperaturverteilung im Objekt anhand des bereitgestellten Modells, des bereitgestellten
zweiten Wärmequellenparameters und des ermittelten Kühlkörperparameters ermittelt.
8. Das Gerät zum Ermitteln der Temperaturverteilung gemäß Anspruch 7, wobei das Gerät
zum Ermitteln der Temperaturverteilung zudem ein Temperatursensorelement (41) umfasst,
um in der zweiten Heizperiode die Temperatur des Objekts innerhalb eines Temperaturbereichs
zu erfassen, dessen Temperatur nicht von der Einheit zum Messen der Temperaturverteilung
(13, 71) gemessen werden kann, wobei die Einheit zum Ermitteln der Temperaturverteilung
(22) in der zweiten Heizperiode die Temperaturverteilung im Objekt anhand des bereitgestellten
Modells, des bereitgestellten zweiten Wärmequellenparameters, des ermittelten Kühlkörperparameters
und der vom Temperatursensorelement (41) erfassten Temperatur ermittelt.
9. Ein Heizsystem zum Erwärmen eines Objekts, das einen Kühlkörper umfasst, wobei das
Heizsystem (1) Folgendes umfasst:
- eine Wärmequelle (2) zum Erwärmen des Objekts, und
- ein Gerät zum Ermitteln der durch das Erwärmen des Objekts verursachten Temperaturverteilung
(21, 71) im Objekt gemäß Anspruch 6.
10. Das Heizsystem gemäß Anspruch 9, wobei das Heizsystem (1) zudem eine Heizungssteuerungseinheit
(12) zum Steuern der Wärmequelle (2) anhand der ermittelten Temperaturverteilung umfasst.
11. Ein Computerprogramm zum Ermitteln von Kühlkörperparametern für ein Objekt, wobei
das Computerprogramm Programmcode umfasst, der ein Gerät zum Ermitteln von Kühlkörperparametern
(23, 71) gemäß Anspruch 1 dazu veranlasst, die folgenden Schritte durchzuführen, wenn
das Computerprogramm auf einem Computer ausgeführt wird, der das Gerät zum Ermitteln
von Kühlkörperparametern (23, 71) steuert:
- Bereitstellen eines die Erwärmung des Objekts definierenden Wärmequellenparameters
(2) mit einer Einheit zum Bereitstellen von Wärmequellenparametern (14), wobei es
sich bei der Wärmequelle um ein Ablationselement zum Abtragen eines Tumors im Objekt
handelt, und wobei es sich beim Kühlkörper um ein Blutgefäß im Objekt handelt, das
die Kühlung des Objekts definiert,
- Bereitstellen eines Modells mithilfe einer Modellbereitstellungseinheit (15), das
die Temperaturverteilung im Objekt beschreibt, während das Objekt von der Wärmequelle
(2) erwärmt wird, wobei das Modell vom Kühlkörperparameter und vom Wärmequellenparameter
abhängig ist,
- Bereitstellen der Kühlkörperposition und eines anfänglichen Kühlkörperparameters
zum Initialisieren des bereitgestellten Modells mithilfe einer Einheit zum Bereitstellen
der Kühlkörperposition (17),
- Messen der Temperaturverteilung im Objekt in der Nähe der bereitgestellten Kühlkörperposition
mit einer Einheit zum Messen der Temperaturverteilung (13, 71), wobei die Einheit
zum Messen der Temperaturverteilung (13, 71) eine Ultraschallsonde umfasst, die in
mindestens einer Scanfläche Ultraschalldaten erfasst, die den Kühlkörper durchqueren,
wobei die mindestens eine Scanfläche jeweils einen Schnittbereich umfasst, in dem
der Kühlkörper die jeweilige Scanfläche schneidet,
- Ermitteln des Kühlkörperparameters mithilfe einer einer Parameterermittlungseinheit
(16) durch Beginnen des Ermitteln des Kühlkörperparameters anhand des initialisierten
Modells sowie durch Minimieren der Abweichung zwischen den Teilen der Temperaturverteilung,
die innerhalb der mindestens einen Scanfläche gemessen wurden, und die den kürzesten
Abstand zur Grenze des Schnittbereichs aufweisen, der kleiner als der vorab festgelegte
Abstand von der bereitgestellten Kühlkörperposition ist, und der modellierten Temperaturverteilung,
wobei die modellierte Temperaturverteilung anhand des bereitgestellten Wärmequellenparameters
und des mithilfe des bereitgestellten Modells zu ermittelnden Kühlkörperparameters
modelliert wird.
12. Ein Computerprogramm zum Ermitteln der Temperaturverteilung in einem Objekt mit einen
Kühlkörper umfasst, wobei der Kühlkörperparameter das Kühlen des Objekts definiert,
während das Objekt von einer Wärmequelle (2) erwärmt wird, wobei das Computerprogramm
Programmcode umfasst, der ein Gerät zum Ermitteln der Temperaturverteilung (21, 71)
gemäß Anspruch 6 dazu veranlasst, die folgenden Schritte durchzuführen, wenn das Computerprogramm
auf einem Computer ausgeführt wird, der das Gerät zum Ermitteln der Temperaturverteilung
(21, 71) steuert:
- Ermitteln eines Kühlkörperparameters im Rahmen folgender Schritte:
- Bereitstellen eines die Erwärmung des Objekts definierenden Wärmequellenparameters
(2) mit einer Einheit zum Bereitstellen von Wärmequellenparametern (14), wobei es
sich bei der Wärmequelle um ein Ablationselement zum Abtragen eines Tumors im Objekt
handelt, und wobei es sich beim Kühlkörper um ein Blutgefäß im Objekt handelt, das
die Kühlung des Objekts definiert,
- Bereitstellen eines Modells mithilfe einer Modellbereitstellungseinheit (15), das
die Temperaturverteilung im Objekt beschreibt, während das Objekt von der Wärmequelle
(2) erwärmt wird, wobei das Modell vom Kühlkörperparameter und vom Wärmequellenparameter
abhängig ist,
- Bereitstellen der Kühlkörperposition und eines anfänglichen Kühlkörperparameters
zum Initialisieren des bereitgestellten Modells mithilfe einer Einheit zum Bereitstellen
der Kühlkörperposition (17),
- Messen der Temperaturverteilung im Objekt in der Nähe der bereitgestellten Kühlkörperposition
mit einer Einheit zum Messen der Temperaturverteilung (13, 71), wobei die Einheit
zum Messen der Temperaturverteilung (13, 71) eine Ultraschallsonde umfasst, die in
mindestens einer Scanfläche Ultraschalldaten erfasst, die den Kühlkörper durchqueren,
wobei die mindestens eine Scanfläche jeweils einen Schnittbereich umfasst, in dem
der Kühlkörper die jeweilige Scanfläche schneidet,
- Ermitteln des Kühlkörperparameters mithilfe einer einer Parameterermittlungseinheit
(16) durch Beginnen des Ermitteln des Kühlkörperparameters anhand des initialisierten
Modells sowie durch Minimieren der Abweichung zwischen den Teilen der Temperaturverteilung,
die innerhalb der mindestens einen Scanfläche gemessen wurden, und die den kürzesten
Abstand zur Grenze des Schnittbereichs aufweisen, der kleiner als der vorab festgelegte
Abstand von der bereitgestellten Kühlkörperposition ist, und der modellierten Temperaturverteilung,
wobei die modellierte Temperaturverteilung anhand des bereitgestellten Wärmequellenparameters
und des mithilfe des bereitgestellten Modells zu ermittelnden Kühlkörperparameters
modelliert wird, und
- Ermitteln der Temperaturverteilung in einem Objekt anhand des von der Modellbereitstellungseinheit
(15) bereitgestellten Modells, des von der Einheit zum Bereitstellen von Wärmequellenparametern
(14) bereitgestellten Wärmequellenparameters und des ermittelten Kühlkörperparameters
mithilfe einer Einheit zum Ermitteln der Temperaturverteilung (22).
1. Appareil de détermination de paramètre de dissipateur thermique destiné à la détermination
d'un paramètre d'un dissipateur thermique à l'intérieur d'un objet, ledit appareil
de détermination de paramètre de dissipateur thermique (23, 71) comprenant :
- une unité de fourniture de paramètre de source de chaleur (14) destinée à la fourniture
d'un paramètre d'une source de chaleur (2), lequel définit un chauffage de l'objet,
la source de chaleur étant un élément d'ablation destiné à l'ablation d'une tumeur
à l'intérieur de l'objet, et le dissipateur thermique étant un vaisseau sanguin à
l'intérieur de l'objet, lequel définit un refroidissement de l'objet ;
- une unité de fourniture de modèle (15) destinée à la fourniture d'un modèle décrivant
une distribution de température dans l'objet lorsque l'objet est chauffé par la source
de chaleur (2), dans lequel le modèle dépend du paramètre de dissipateur thermique
et du paramètre de source de chaleur ;
- une unité de fourniture d'emplacement de dissipateur thermique (17) destinée à la
fourniture d'un emplacement de dissipateur thermique et un paramètre initial de dissipateur
thermique pour initialiser le modèle fourni ;
- une unité de mesure de distribution de température (13, 71) destinée à la mesure
d'une distribution de température dans l'objet à proximité de l'emplacement de dissipateur
thermique fourni, l'unité de mesure de distribution de température (13, 71) comprenant
une sonde à ultrasons configurée pour acquérir des données ultrasonores dans au moins
un balayage des plans traversant le dissipateur thermique, l'au moins un plan de balayage
comprenant respectivement une zone d'intersection où le dissipateur thermique intersecte
le plan de balayage respectif ; et
- une unité de détermination de paramètre (16) destinée à la détermination du paramètre
de dissipateur thermique par commencement de la détermination du paramètre de dissipateur
thermique avec le modèle initialisé, et par minimisation d'un écart entre seulement
les parties de la distribution de température, lesquelles ont été mesurées dans l'au
moins un plan de balayage et présentent une distance la plus courte à une bordure
de la zone d'intersection, laquelle est inférieure à une distance prédéfinie de l'emplacement
de dissipateur thermique fourni, et une distribution de température modélisée, dans
lequel la distribution de température modélisée est modélisée en fonction du paramètre
de source de chaleur fourni et du paramètre de dissipateur thermique à déterminer
à l'aide du modèle fourni.
2. Appareil de détermination de paramètre de dissipateur thermique selon la revendication
1,
dans lequel l'appareil de détermination de paramètre de dissipateur thermique (23,
71) est conçu pour déterminer au moins un paramètre de dissipateur thermique, lequel
comprend un emplacement du dissipateur thermique à l'intérieur de l'objet et/ou lorsque
le dissipateur thermique comprend une structure tubulaire comportant un fluide en
circulation, une direction d'écoulement et/ou un rayon de la structure tubulaire et/ou
un débit.
3. Appareil de détermination de paramètre de dissipateur thermique selon la revendication
1,
dans lequel l'unité de fourniture de modèle (15) est conçue pour fournir le modèle
de telle sorte qu'il dépend également d'un paramètre d'objet, et dans lequel l'unité
de détermination de paramètre (16) est conçue pour déterminer le paramètre d'objet
et le paramètre de dissipateur thermique par minimisation de l'écart entre la distribution
de température mesurée et la distribution de température modélisée, dans lequel la
distribution de température modélisée est modélisée en fonction du paramètre de source
de chaleur fourni, du paramètre de dissipateur thermique à déterminer et du paramètre
d'objet à déterminer à l'aide du modèle fourni.
4. Appareil de détermination de paramètre de dissipateur thermique selon la revendication
1,
dans lequel l'unité de mesure de distribution de température (13, 71) comprend une
sonde à ultrasons destinée à l'acquisition des données ultrasonores de l'objet et
une unité de thermométrie par ultrasons pour déterminer la distribution de température
en fonction des données ultrasonores acquises.
5. Appareil de détermination de paramètre de dissipateur thermique selon la revendication
4,
dans lequel l'unité de fourniture d'emplacement de dissipateur thermique (17) est
conçue pour déterminer l'emplacement de dissipateur thermique en fonction des données
ultrasonores de l'objet acquises par la sonde à ultrasons.
6. Appareil de détermination de distribution de température destiné à la détermination
d'une distribution de température à l'intérieur d'un objet comprenant un dissipateur
thermique, dans lequel l'appareil de détermination de distribution de température
(21, 71) comprend :
- un appareil de détermination de paramètre de dissipateur thermique (23, 71) destiné
à la détermination d'un paramètre du dissipateur thermique, tel que défini dans la
revendication 1 ;
- une unité de détermination de distribution de température (22) destinée à la détermination
de distribution de température à l'intérieur de l'objet en fonction du modèle fourni
par l'unité de fourniture de modèle (15), d'un paramètre de source de chaleur fourni
par l'unité de fourniture de paramètre de source de chaleur (14) et du paramètre de
dissipateur thermique déterminé.
7. Appareil de détermination de distribution de température selon la revendication 6,
dans lequel la source de chaleur (2) est conçue pour chauffer l'objet dans une première
période de chauffage à une température inférieure définie par un premier paramètre
de source de chaleur fourni et dans une seconde période de chauffage à une température
supérieure définie par un second paramètre de source de chaleur fourni, dans lequel
- l'unité de mesure de distribution de température (13, 71) est conçue pour mesurer
la distribution de température, laquelle est utilisée pour déterminer le paramètre
de dissipateur thermique, pendant la première période de chauffage, dans lequel l'unité
de détermination de paramètre (16) est conçue pour déterminer le paramètre de dissipateur
thermique en fonction du modèle fourni, du premier paramètre de source de chaleur
fourni et
- de la distribution de température mesurée, et l'unité de détermination de distribution
de température (22) est conçue pour déterminer, pendant la seconde période de chauffage,
la distribution de température à l'intérieur de l'objet en fonction du modèle fourni,
du second paramètre de source de chaleur fourni et du paramètre de dissipateur thermique
déterminé.
8. Appareil de détermination de distribution de température selon la revendication 7,
dans lequel l'appareil de détermination de distribution de température comprend en
outre un élément de détection de température (41) destiné à la détection, dans la
seconde période de chauffage, d'une température de l'objet dans une plage de température
dans laquelle une température n'est pas mesurable par l'unité de mesure de distribution
de température (13, 71), dans lequel l'unité de détermination de distribution de température
(22) est conçue pour déterminer, dans la seconde période de chauffage, la distribution
de température à l'intérieur de l'objet en fonction du modèle fourni, du second paramètre
de source de chaleur fourni, du paramètre de dissipateur thermique déterminé et de
la température détectée par l'élément de détection de température (41).
9. Système de chauffage destiné au chauffage d'un objet comprenant un dissipateur thermique,
ledit système de chauffage (1) comprenant :
- une source de chaleur (2) destinée au chauffage de l'objet ; et
- un appareil de détermination de distribution de température (21, 71) destiné à la
détermination d'une distribution de température à l'intérieur de l'objet provoquée
par le chauffage de l'objet, tel que défini dans la revendication 6.
10. Système de chauffage selon la revendication 9, dans lequel le système de chauffage
(1) comprend en outre une unité de commande de chauffage (12) destinée à la commande
de la source de chaleur (2) en fonction de la distribution de température déterminée.
11. Programme informatique destiné à la détermination d'un paramètre d'un dissipateur
thermique à l'intérieur d'un objet, ledit programme informatique comprenant un moyen
de code de programme destiné à l'amenée d'un appareil de détermination de paramètre
de dissipateur thermique (23, 71), tel que défini dans la revendication 1, à effectuer
les étapes suivantes lorsque ledit programme informatique est exécuté sur un ordinateur
commandant l'appareil de détermination de paramètre de dissipateur thermique (23,
71) :
- la fourniture d'un paramètre d'une source de chaleur (2), lequel définit un chauffage
de l'objet, par une unité de fourniture de paramètre de source de chaleur (14), la
source de chaleur étant un élément d'ablation destiné à l'ablation d'une tumeur à
l'intérieur de l'objet, et le dissipateur thermique étant un vaisseau sanguin à l'intérieur
de l'objet, lequel définit un refroidissement de l'objet ;
- la fourniture d'un modèle décrivant une distribution de température dans l'objet
lorsque l'objet est chauffé par la source de chaleur (2) par une unité de fourniture
de modèle (15), dans lequel le modèle dépend du paramètre de dissipateur thermique
et du paramètre de source de chaleur ;
- la fourniture d'un emplacement de dissipateur thermique et d'un paramètre initial
de dissipateur thermique pour initialiser le modèle fourni par une unité de fourniture
d'emplacement de dissipateur thermique (17) ;
- la mesure d'une distribution de température dans l'objet à proximité de l'emplacement
de dissipateur thermique fourni par une unité de mesure de distribution de température
(13, 71), l'unité de mesure de distribution de température (13, 71) comprenant une
sonde à ultrasons configurée pour acquérir des données d'ultrasons dans au moins un
balayage des plans traversant le dissipateur thermique, l'au moins un plan de balayage
comprenant respectivement une zone d'intersection où le dissipateur thermique intersecte
le plan de balayage respectif ;
- la détermination du paramètre de dissipateur thermique par commencement de la détermination
du paramètre de dissipateur thermique avec le modèle initialisé et par minimisation
d'un écart entre seulement les parties de la distribution de température, lesquelles
ont été mesurées dans l'au moins un plan de balayage et présentent une distance la
plus courte à la bordure de la zone d'intersection, laquelle est inférieure à une
distance prédéfinie de l'emplacement de dissipateur thermique fourni, et une distribution
de température modélisée par une unité de détermination de paramètre (16), dans lequel
la distribution de température modélisée est modélisée en fonction du paramètre de
source de chaleur fourni et du paramètre de dissipateur thermique à déterminer à l'aide
du modèle fourni.
12. Programme informatique destiné à la détermination d'une distribution de température
à l'intérieur d'un objet comprenant un dissipateur thermique, dans lequel le paramètre
de dissipateur thermique définit un refroidissement de l'objet lorsque l'objet est
chauffé par une source de chaleur (2), ledit programme informatique comprenant un
moyen de code de programme destiné à l' amenée d'un appareil de détermination de distribution
de température (21, 71), tel que défini dans la revendication 6, à effectuer les étapes
suivantes lorsque ledit programme informatique est exécuté sur un ordinateur commandant
l'appareil de détermination de distribution de température (21, 71) :
- la détermination d'un paramètre du dissipateur thermique :
- par fourniture d'un paramètre d'une source de chaleur (2), lequel définit un chauffage
de l'objet, par une unité de fourniture de paramètre de source de chaleur (14), la
source de chaleur étant un élément d'ablation destiné à l'ablation d'une tumeur à
l'intérieur de l'objet, et le dissipateur thermique étant un vaisseau sanguin à l'intérieur
de l'objet, lequel définit un refroidissement de l'objet ;
- la fourniture d'un modèle décrivant une distribution de température dans l'objet
lorsque l'objet est chauffé par la source de chaleur (2) par une unité de fourniture
de modèle (15), dans lequel le modèle dépend du paramètre de dissipateur thermique
et du paramètre de source de chaleur ;
- la fourniture d'un emplacement de dissipateur thermique et d'un paramètre initial
de dissipateur thermique pour initialiser le modèle fourni par une unité de fourniture
d'emplacement de dissipateur thermique (17) ;
- la mesure d'une distribution de température dans l'objet à proximité de l'emplacement
de dissipateur thermique fourni par une unité de mesure de distribution de température
(13, 71), l'unité de mesure de distribution de température (13, 71) comprenant une
sonde à ultrasons configurée pour acquérir des données ultrasonores dans au moins
un plan de balayage traversant le dissipateur thermique, l'au moins un plan de balayage
comprenant respectivement une zone d'intersection où le dissipateur thermique intersecte
le plan de balayage respectif ;
- la détermination du paramètre de dissipateur thermique par commencement de la détermination
du paramètre de dissipateur thermique avec le modèle initialisé, et par minimisation
d'un écart entre seulement les parties de la distribution de température, lesquelles
ont été mesurées dans l'au moins un plan de balayage et présentent une distance la
plus courte à la bordure de la zone d'intersection, laquelle est inférieure à une
distance prédéfinie de l'emplacement de dissipateur thermique fourni, et une distribution
de température modélisée par une unité de détermination de paramètre (16), dans lequel
la distribution de température modélisée est modélisée en fonction du paramètre de
source de chaleur fourni et du paramètre de dissipateur thermique à déterminer à l'aide
du modèle fourni ; et
- la détermination de la distribution de température à l'intérieur de l'objet en fonction
du modèle fourni par l'unité de fourniture de modèle (15), d'un paramètre de source
de chaleur fourni par l'unité de fourniture de paramètre de source de chaleur (14)
et du paramètre de dissipateur thermique déterminé par une unité de détermination
de distribution de température (22).