Patent Application
20230030917
The present Patent application for industrial invention relates to a method for the treatment of a tumor area by means of high intensity focused ultrasound (HIFU—High Intensity Focused Ultrasound) combined with monitoring the temperature distribution and the actually treated area.
There are many documents describing tumor mass treatments by means of high intensity focused ultrasound (HIFU).
During a HIFU treatment, one of the main problems is to determine the actual temperature spatial distribution in real time during the treatment.
An example of temperature monitoring by means of the analysis of the ultrasound signal reflected by tissues is described in US2015005634, where a method is described which uses the variation of acoustic attenuation in tissues while increasing temperature, and, on the basis of this effect, measures temperature by detecting temperature variation of the harmonics of the ultrasound signal reflected by tissues as a response to the high temperature pulses.
However, the method described requires a calibration on fictitious tissue portions (phantom) in order to quantify the variations in harmonic parameters as a function of temperature variations.
It is clear that the phantoms used to study the problem (i) do not allow to consider the heat diffusive phenomena linked to body fluids circulation and, however (ii) approach only partially the real tissue features, since tissues of different body parts have different behavior and, in case of equal body part interested, the tissues of different subjects have a different behavior, in particular with reference to the actual acoustic heating.
Other embodiments of temperature non-invasive measuring methods are described in documents US2020237234, US2019142513, US2019262074, WO2019173138, but none of these ones solves the technical problem to measure a temperature spatial distribution by simply considering the specific features of tissues undergoing a HIFU treatment, all by means of a unique ultrasound probe and in real time.
More generally, to the best knowledge of the current inventors, some problems linked to the usage of ultrasound signals for diagnosis and treatment of tumor masses in various organs of the human body remain unsolved.
In fact, embodiments which, by means of a unique ultrasound probe, allow: (i) to individuate the tumor mass; (ii) to verify the correct positioning of the probe with respect to the tumor mass; (iii) to treat the tumor mass by means of HIFU; (iv) to verify the extension of the area actually treated by temperature monitoring in real time are not known.
Therefore, aim of the present invention is to provide a method for tumor mass treatment, which, by means of a unique ultrasound probe, allows: (i) to individuate the tumor mass; (ii) to verify the correct positioning of the probe with respect to the tumor mass; (iii) to treat the tumor mass by means of HIFU; (iv) to verify the extension of the area actually treated by temperature monitoring in real time.
More in particular, aim of the present invention is to provide an ultrasound device and a method which allow to measure a temperature spatial distribution by simply considering the specific features of the tissues undergoing a HIFU treatment, and in particular the actual acoustic heating rate, all by using a unique ultrasound probe and in real time.
These and other advantages will be clear from the detailed description of the invention, which will refer to the appended drawings 1 to 6. In
Firstly, it is to be said that in the context of the present Patent application, with “raw ultrasound signal” or “radiofrequency ultrasound signal” it is intended the ultrasound signal emitted by the probe and reflected by the human body towards the same probe, before the same is processed in order to obtain the ultrasound image; yet, if not otherwise specified, it is to be specified that with “ultrasound image” it is intended an ultrasound image of B-mode type, obtained along the propagation plane of the ultrasound beam emitted by the probe.
Yet, it is to be specified that with “raw ultrasound signal corresponding to a determined ROI” it is intended the portion of the raw ultrasound signal which, suitably treated, originated a segment of ultrasound image relating to the area of interest (ROI).
In fact, as it is known, the correlation between the position of each pixel in the ultrasound image and the ultrasound signal is carried out as a function of the time interval between the ultrasound pulse emission and the relative echo reception (reflected signal), since the signal reflected by the tissues positioned at greater depths uses more time to reach the probe after being reflected.
Therefore, regardless of the nature of the possible following processing, the segmentation of the “raw ultrasound signal corresponding to a determined ROI” occurs in the time domain, in order to isolate the raw ultrasound signal portion that, in the ultrasound image, originated a determined segment.
Yet, if not otherwise specified, it is intended that, according to what commonly known at the state of the art, the ultrasound probe used is a probe containing an array of piezoelectric transducers or CMUT configured to emit a plurality of ultrasound signals arranged side by side, so that to each signal a vertical line (“line of sight”) of the ultrasound image corresponds and the assembly of the lines arranged side by side allows to make up the ultrasound image again. It is specified that, for brevity and clarity of description, the processing carried out on ultrasound signals will be described with reference to a single ultrasound signal. Obviously, even where it is not explicitly specified, it is clear that all the processing can be conveniently applied to a plurality of raw ultrasound signals, each one received by one of the piezoelectric transducers comprised in the ultrasound probe.
Yet, it is specified that with “point” of raw signal (radiofrequency ultrasound signal) it is intended the value of the raw signal in a single sampling point: in a 20 MHz sampling, as a way of purely example, 20000 points per millisecond of acquired signal are obtained.
It is also specified that all the processing described in the following is carried out by means of an echographic device provided with at least an ultrasound probe comprising at least an array of piezoelectric transducers or CMUT and suitable piloting means of said probe, with computing means for processing the signal configured to generate the signal to be sent by means of said probe and to analyze the signal received by said probe to obtain an ultrasound image, with user interaction means comprising a graphical interface and controlling means such as for example keyboards and/or pointing devices.
The device used comprises a HIFU ultrasound probe of the type known at the state of the art. Preferably, said probe comprises an array of piezoelectric transducers or CMUT configured to emit and receive ultrasound signals of frequency and intensity useful for the ultrasound imaging and an array of piezoelectric transducers configured to emit high intensity focused ultrasound signals.
These devices, also as a whole, are known at the state of the art and commonly used in the ultrasound technique. However, it is specified that the ultrasound device used for carrying out the methods described in the following is configured to process the raw ultrasound signal (radiofrequency ultrasound signal) for determining an ultrasound image, but also to store the raw ultrasound signal in order to carry out following processing: it is a function not available in the ultrasound imaging devices available on the market.
As yet stated, a problem solved by the present invention is to provide a method for tumor mass treatment which, by means of a unique ultrasound probe, allows
The point (i) of tumor mass individuation can be carried out by means of an automatic tissue classification method, or by means of manual individuation by the doctor.
The definition of the area to be treated in point (ii) is proposed to the doctor as area identified as tumor, in addition to a buffer of possible further extension determinable by the doctor. It is clear that the definition of the area can be also defined by the doctor and highlighted on the ultrasound image by means of suitable commands provided on the echographic device.
As it is known, the effects of a HIFU treatment (positive if the treatment is applied to a tumor mass, negative if the treatment is applied to a healthy tissue) are correlated to the temperature reached by the tissues after the treatment; therefore, it is necessary to have a mean for measuring the temperature reached by the tissues treated during the treatment.
In fact, the key goal of all the ablative therapies is to reach and keep a temperature between 60 and 100° C. in the whole volume of tissue to be treated to induce coagulative necrosis.
It has been said that in the treatment systems known at the state of the art, the temperature reached is not directly assessable in a reliable manner as a function of the treatment time, since the response to ultrasound insonification depends on the specific tissue considered and is variable from site to site of the human body, and in case of same site, from patient to patient.
In other words, it is not perfectly known, a priori, the entity of acoustic heating of each tissue as a response to a HIFU treatment.
The method according to the present invention allows instead to monitor the actual temperature distribution in ROI, during the whole duration of the HIFU treatment, in a non-invasive way and by means of exclusive ultrasound signal analysis.
In
In
It has to be specified that in the following, for simplicity, we will refer to the focal point (11) of the HIFU beam. With this term it is meant the centre of the focusing volume. The HIFU beam focusing volume can be schematized in a bidimensional ultrasound image by means of a focusing area. The focal point (11) indicated in the following is the centre of such focusing area.
As it is possible to observe, the probe (1) is configured:
Conveniently, the probe can be configured to use a unique array of piezoelectric or CMUT transducers both to obtain the ultrasound image and to detect—clearly in discrete moments—the ultrasound signals reflected and/or emitted by the tissues of the patient after the HIFU treatment. As an alternative, the probe can comprise two arrays of transducers, one dedicated to the HIFU treatment and the other one to the acquisition of signals for the ultrasound imaging.
In other words, the device is configured to carry out cycles of:
For each reflected raw ultrasound signal, the device is also configured to calculate the frequency spectrum (200s).
In
It is to be noted how during each HIFU cycle the temperature increases, and so the following low intensity pulses are applied to a tissue with gradually increasing temperatures.
In
Therefore, in the ultrasound image plane it is possible to define a second coordinate system (X, Y) originated in the repere point P. This second coordinate system (X, Y) is fixed to the body of the patient, while the first two systems are integral to the probe (1). It is clear that each point of the ROI (15) can be identified in each of the three mentioned coordinate systems, and that by means of coordinate transformations with formulae already known at the state of the art, each point can be identified in each of the three mentioned coordinate systems.
After defining the coordinate systems, it is possible to determine firstly a method for monitoring the temperature distribution after the HIFU treatment, assuming that the ultrasound probe (1) is still in a static position.
Therefore, it is assumed that being a starting temperature distribution in ROI assigned, the time temperature variation in each point of the ROI after HIFU treatment can be expressed by the heating diffusion equation:
Wherein:
While I0(r,z) is a feature of the ultrasound probe, and so it is known once the used ultrasound probe is characterized, Q depends on a plurality of specific parameters, such for example the tissue absorbing rate and the tissue attenuation rate, which are variable from patient to patient and, for each patient, inside each site, and so it is an unknown value.
It is specified that the possibility to determine correctly and in real time the value of the actual rate of acoustic heating Q for the specific patient can allow to calculate by means of equation (1) the temperature spatial and time distribution in tissues interested by the treatment, to view it in real time by means of the graphical interface of the sonographer and to allow the operator to carry out the treatment to avoid exceeding temperature thresholds that are critical for specific tissues close to the ones interested by the treatment.
The present invention allows to determine in real time and by using exclusively the just described ultrasound probe the acoustic heating rate value Q for the specific tissue object of treatment of each patient, and, so, to determine the temperature spatial and time distribution.
In order to obtain the temperature spatial distribution, the device according to the present invention is configured to carry out the following method:
In this way, it was obtained both the determination of an actual heating rate value Q specific for tissues under treatment and a starting temperature distribution for the next treatment step.
It is specified that in a first embodiment said peaks (2021) can be individuated by carrying out the next steps for each iteration:
Moreover, preferably, in order to classify a frequency (f1) as peak, it is to be verified that the difference or the ratio between said lower (f1′) and higher (f1″) frequencies is lower than a predetermined amount, thus introducing a test condition of the peak amplitude.
It is to be specified that the temperature of 100° C. at a significant change of the spectrum shape is justified by the fact that this value is a good assessment of the temperature at which localized boiling/cavitation phenomena can occur.
However, the whole operation is carried out with the device still and with the HIFU beam constantly focused on a unique point inside the ROI.
This means that only a very little region (ideally a point) of tissue inside the ROI, preferably chosen at the geometric center of gravity of the same, reaches the temperature of 100° C., while all the volume surrounding it reaches lower temperatures, thus guaranteeing the safety of treatment.
The device is also configured to calculate a plurality of parameters of shape of the frequency spectrum relating to the raw ultrasound signal reflected in response to the broadband ultrasound signal emitted during each break between emission of two focused ultrasound beams of HIFU treatment, and in particular of the portion of raw signal corresponding to a propagation path of the ultrasound signal comprising the focal point.
After determining the temperature time distribution, it is therefore possible to associate the respective focal point temperature to each calculated frequency spectrum (200s, 202s).
Therefore, the HIFU treatment can continue in other points localized inside the region to be treated. In this case, the device is configured to measure the temperature at the focal point (1) relative to the new position of the probe:
As temperature value at the focus, it is assumed the temperature value for whose spectrum said correlation coefficient is higher.
In another embodiment said correlation coefficients are provided in decreasing order, and as temperature value it is assumed the average of temperature values associated to the highest correlation coefficients, for example to the 5 highest correlation coefficients.
It is specified that the condition verified in point f) needs to verify if boiling is observed at the focus (11). Since the treated tissues are aqueous media, in the present method it is in fact assumed that boiling starts at about 100° C. The presence of boiling causes not only that the tissue reflects the HIFU ultrasound pulses according to its own features, but that it has also an ultrasound emission of its own due to pulses emitted by water vapour bubbles forming, which can be detected by means of the frequency spectrum analysis. It is clear that the focus has to be inside the area to be treated, to avoid damaging healthy tissues.
Moreover, it is specified that for solving equation (1), in steps i) and j) a literature value for the K parameter can be conveniently used, or that the same can be determined directly by suitable tissue-mimicking phantoms.
Without this being limiting for the aims of the invention, the K value used is preferably between 1,1E-07 m2/s and 1,6E-07 m2/s.
Yet, it is specified that the iterative solution of equation (1) at point i) is carried out assuming that boiling occurs in primis at the focus.
Equation (1) is solved with an attempt value for Q, by calculating the temperature distribution for all the time intervals between moment t=0 of treatment start, for which the temperature starting distribution was assigned, up to moment t=tboil in which boiling was detected by means of the frequency spectrum analysis. The Q value is then modified depending on the fact that the temperature value calculated at the focus for t=tboil is higher or lower than 100° C., and the procedure is repeated up to when the calculated value differs from 100° C. by an amount lower than a predetermined threshold, for example 1° C.
As it is shown in
Moreover, the method proposed allows to verify also the actual temperature distribution inside the treated area, also in case of probe displacement.
As it is shown in
However, the identification on the new ultrasound image of the new position of the repere point (P) allows to correctly position the temperature spatial distribution values calculated in the previous probe position, since their coordinates (X, Y) in the coordinates system originated in the repere point are not changed.
Therefore, in the method according to the present invention, after each updating of temperature spatial distribution, conveniently carried out in the system (x, y) integral to the ultrasound probe, said temperature spatial distribution is transformed in the system (X, Y) originated in the repere point and integral thereto.
After probe displacement, the inverted operation will be carried out by transforming the temperature spatial distribution from the system integral to the repere point to the system integral to the ultrasound probe again.
| Number | Date | Country | Kind |
|---|---|---|---|
| 102019000025306 | Dec 2019 | IT | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/IB2020/062408 | 12/23/2020 | WO |
