Patent ID: 8583211

Claim:
A method for temperature control in MR-guided administration of ultrasound comprising: administering ultrasound therapy to a patient at an in vivo site by emitting focused ultrasound into the site at multiple foci with a multi-focus ultrasound therapy device and by switching the actual focus position to perform a cycled sonication along a time-invariant volumetric pattern; monitoring temperature in a localized region of an examination subject in which said site is located during said therapy in real-time by MR thermometry; from the MR thermometry, automatically identifying, in a processor, characteristics of a temperature distribution in the monitored region of the examination subject; from the computerized processor, implementing temperature control by regulating the ultrasound therapy device by modulating energy deposition at each individual focal point to homogenize a temperature history of all sonicated foci; in said processor, implementing a calculation calculating a ratio between energy variation δE k (n) at a time point n and an initial value of energy E K,1 , in every sonicated focus k among the multiple locations heated simultaneously, without using any a-priori information on tissue properties, according to the following double-branch master equation that takes into account two regimens, respectively before and after reaching a steady state: δ ⁢ ⁢ E k ⁡ ( n ) E k , 1 = { w · min ⁡ ( n , n max ) · ( 1 - T k , n Θ n ) + w · Ω k , n if ⁢ ⁢ n ≤ n switch - 1 + tanh ⁡ ( α · n switch ) + w · min ⁢ ( n , n max ) · ( 1 - T k , n Θ top ) + w · Ω k , n if ⁢ ⁢ n > n switch or δ ⁢ ⁢ E k ⁡ ( n ) E k , 1 = { w · n max · tanh ⁡ ( n n max ) · ( 1 - T k , n Θ n ) + w · Ω k , n if ⁢ ⁢ n ≤ n switch - 1 + tanh ⁡ ( α · n switch ) + w · n max · tanh ⁡ ( n n max ) · ( 1 - T k , n Θ top ) + w · Ω k , n if ⁢ ⁢ n > n switch wherein tanh is the analytical function hyperbolic tangent, n switch denotes a sampling point where the controller switches irreversibly from the rapid temperature elevation regimen to the steady state regimen, α is a positive number smaller than 0.1, n max is an a-dimentional parameter that limits a maximum strength of the negative reaction, T k,n is the actual temperature at each sonicated focus k for the measurement n, Θ n is an average value of measured temperatures T k,n over the sonicated foci for measurement n, Θ top is a flattop target value of steady state temperature elevation, w is a negative reaction parameter that is a-dimensional and scaled to a square of a temporal resolution, and Ω k,n is a cumulative relative offset of the measured temperature until the measurement n: Ω k , n = { ∑ n ′ = 1 n ⁢ ⁢ ( 1 - T k , n ′ Θ n ′ ) if ⁢ ⁢ n ≤ n switch ∑ n ′ = 1 n switch ⁢ ⁢ ( 1 - T k , n ′ Θ n ′ ) + ∑ n ′ = n switch + 1 n ⁢ ⁢ ( 1 - T k , n ′ Θ top ) if ⁢ ⁢ n > n switch ,