Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP3476871B2 - Fault surface temperature distribution display method and apparatus - Google Patents
[go: Go Back, main page]

JP3476871B2 - Fault surface temperature distribution display method and apparatus - Google Patents

Fault surface temperature distribution display method and apparatus

Info

Publication number
JP3476871B2
JP3476871B2 JP23052893A JP23052893A JP3476871B2 JP 3476871 B2 JP3476871 B2 JP 3476871B2 JP 23052893 A JP23052893 A JP 23052893A JP 23052893 A JP23052893 A JP 23052893A JP 3476871 B2 JP3476871 B2 JP 3476871B2
Authority
JP
Japan
Prior art keywords
temperature
temperature distribution
fault
high frequency
display device
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP23052893A
Other languages
Japanese (ja)
Other versions
JPH0767969A (en
Inventor
五郎 山本
Original Assignee
山本ビニター株式会社
中央電子計測株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 山本ビニター株式会社, 中央電子計測株式会社 filed Critical 山本ビニター株式会社
Priority to JP23052893A priority Critical patent/JP3476871B2/en
Publication of JPH0767969A publication Critical patent/JPH0767969A/en
Application granted granted Critical
Publication of JP3476871B2 publication Critical patent/JP3476871B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Electrotherapy Devices (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、高周波加温治療におけ
る人体等の加温対象物の断層面の温度分布を表示する方
法及び装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and apparatus for displaying a temperature distribution on a tomographic plane of a heated object such as a human body in high frequency heating treatment.

【0002】[0002]

【従来の技術】従来、例えば生体を挾んだ1対の電極か
ら高周波エネルギーを放射供給して、この生体の患部を
加温治療する高周波加温治療装置が知られている(特開
昭55−130675号、特開昭56−80265号公
報)。この高周波加温治療装置は、例えば癌その他の腫
瘍を構成する異常な細胞組織と周囲の正常な細胞組織と
を共に40℃以上の温度範囲で加温すると、前者の異常
な細胞組織が正常な細胞組織に比して2〜2.5℃だけ
高温になるという点に着目し、正常な細胞組織を壊死さ
せない43℃以下に保持する一方、異常な細胞組織を4
5℃前後にまで上昇させて壊死崩壊させんとしたもので
ある。
2. Description of the Related Art Conventionally, there is known a high-frequency heating treatment apparatus for radiating and supplying high-frequency energy from a pair of electrodes sandwiching a living body to heat-treat an affected part of the living body (Japanese Patent Laid-Open No. 55-55). -130675, JP-A-56-80265). This high-frequency hyperthermia treatment apparatus heats both the abnormal cell tissue that constitutes cancer and other tumors and the surrounding normal cell tissue within a temperature range of 40 ° C. or higher to normalize the former abnormal cell tissue. Focusing on the fact that the temperature is higher by 2 to 2.5 ° C. than that of cellular tissue, while maintaining normal cellular tissue at 43 ° C. or below that does not necrosis,
The temperature was raised up to around 5 ° C to prevent necrosis and collapse.

【0003】従って、かかる装置では生体内温度の監視
が極めて重要であり、そのため加温治療の際には、生体
内の温度を監視するために高周波エネルギーが放射され
る体内に所要数の針状温度計を差し込んで、温度計測を
行っていた。
Therefore, in such a device, it is extremely important to monitor the temperature in the living body. Therefore, during the heating treatment, a required number of needle-shaped needles are radiated into the body to radiate high frequency energy in order to monitor the temperature in the living body. The temperature was measured by inserting a thermometer.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、加温治
療装置の電極から放射される高周波エネルギーは場所に
よって均一ではなく、しかも電極構造や生体の患部の位
置等により左右されるため、体温計測用の温度計が少な
いと温度計測が好適に行われない虞れがあり、一方、多
数の体温計を用いると温度監視は好適に行い得るが、患
者の苦痛は勢い増大することとなる。また、実際には数
箇所以上に針を刺すのは非現実的であり、従って、体内
温度の計測には一定の精度的限界があった。
However, since the high-frequency energy radiated from the electrodes of the heating treatment apparatus is not uniform depending on the location and depends on the electrode structure and the position of the affected part of the living body, it is necessary to measure the body temperature. If the number of thermometers is small, there is a possibility that temperature measurement may not be performed appropriately, while if a large number of thermometers are used, temperature monitoring may be performed favorably, but the pain of the patient will increase. In addition, actually, it is unrealistic to puncture the needle at several points or more, and therefore there is a certain accuracy limit in measuring the body temperature.

【0005】本発明は、上記に鑑みてなされたもので、
例えば電極間の高周波エネルギー分布、生体の特性等を
用いて実際に近い生体内の温度分布を求め、これを表示
するようにして加温治療の際の温度監視の重要情報とし
て利用し得るようにした断層面の温度分布表示方法及び
その装置を提供するものである。
The present invention has been made in view of the above,
For example, the in-vivo temperature distribution is obtained using the high-frequency energy distribution between electrodes, the characteristics of the living body, etc., and this can be displayed so that it can be used as important information for temperature monitoring during heating treatment. The present invention provides a method and apparatus for displaying a temperature distribution on a fault plane.

【0006】[0006]

【課題を解決するための手段】本発明は、断層面各部の
温度T(x,y)i+1を、前回の計算値T(x,y)iを元に、造影
剤不使用時の密度による温度上昇率D0(x,y)、高周波エ
ネルギーによる温度上昇率β・e(x,y)及び温度上昇に
伴う冷却効率U(x,y)等を加味して算出し、算出された
各部の温度T(x,y)i+1を表示するようにしたものである
(請求項)。
According to the present invention, the temperature T (x, y) i + 1 at each part of the tomographic plane is calculated based on the previous calculated value T (x, y) i when a contrast agent is not used. The temperature rise rate D0 (x, y) due to the density, the temperature rise rate β · e (x, y) due to the high frequency energy, and the cooling efficiency U (x, y) due to the temperature rise were calculated and calculated. The temperature T (x, y) i + 1 of each part is displayed (claim 1 ).

【0007】なお、上記冷却効率U(x,y)は、造影剤不
使用時の密度による温度上昇寄与率D0(x,y)と造影剤使
用時の密度による温度上昇寄与率D1(x,y)とから実際の
血流による放熱効果に基づく補正係数K(x,y)(=[D0
(x,y)−D1(x,y)]・U(T)/n)を求め、この補正係数K
(x,y)を乗じて、より実際的な冷却効率U(x,y)′=K
(x,y)・U(x,y)に補正するとよい(請求項)。
The cooling efficiency U (x, y) is the temperature increase contribution rate D0 (x, y) due to the density when the contrast agent is not used and the temperature increase contribution rate D1 (x, y due to the density when the contrast agent is used. y) and the correction coefficient K (x, y) (= [D0
(x, y) −D1 (x, y)] · U (T) / n), and the correction coefficient K
Multiplying by (x, y), more practical cooling efficiency U (x, y) ′ = K
It may be corrected to (x, y) · U (x, y) (Claim 2 ).

【0008】また、上記高周波エネルギー量の検出手段
は、被加温体の近傍位置における電磁界の強度を検出す
る電磁界検出手段と、検出された電磁界強度から被加温
体に供給される高周波エネルギー量を算出する演算手段
とから構成するとよい(請求項)。
The high-frequency energy amount detecting means is supplied to the object to be heated from the electromagnetic field detecting means for detecting the strength of the electromagnetic field in the vicinity of the object to be heated and the detected electromagnetic field strength. It may be composed of an arithmetic means for calculating the amount of high frequency energy (claim 3 ).

【0009】また、表示を色信号で行うようにしてもよ
いし(請求項)、濃度信号で行うようにしてもよい
(請求項)。更に、断層面の適所に温度計を設置し、
この温度計からの実際の温度で、演算温度に補正を加え
ることが好ましい(請求項)。
Further, the display may be performed by a color signal (claim 4 ) or by a density signal (claim 5 ). Furthermore, we installed a thermometer in the proper place on the fault plane,
It is preferable to correct the calculated temperature with the actual temperature from this thermometer (claim 6 ).

【0010】[0010]

【作用】本発明によれば、断層面各部の温度T(x,y)i+1
が前回の計算値T(x,y)iを元に、造影剤不使用時の密度
による温度上昇寄与率D0(x,y)、高周波エネルギー量に
基づく温度上昇率β、高周波エネルギー分布による温度
上昇寄与率e(x,y)及び温度上昇に伴う冷却効率U(x,y)
等を加味した所定の演算式により繰り返し算出される。
そして、算出される各部の温度T(x,y)i+1が断層面に重
畳されて逐次、表示手段に表示される(請求項)。
According to the present invention, the temperature T (x, y) i + 1 at each part of the fault plane is
Is based on the previous calculated value T (x, y) i, the temperature rise contribution rate D0 (x, y) due to the density when the contrast agent is not used, the temperature rise rate β based on the high frequency energy amount, and the temperature according to the high frequency energy distribution. Increase contribution rate e (x, y) and cooling efficiency U (x, y) accompanying temperature rise
It is repeatedly calculated by a predetermined arithmetic expression in consideration of the above.
Then, the calculated temperature T (x, y) i + 1 of each part is superimposed on the tomographic plane and sequentially displayed on the display means (claim 1 ).

【0011】また、本発明によれば、造影剤不使時の密
度による温度上昇寄与率D0(x,y)と造影剤使用時の密度
による温度上昇寄与率D1(x,y)とから実際の血流に基づ
く補正係数K(x,y)が算出され、冷却効率U(x,y)にこの
補正係数K(x,y)を乗じて実際の血流に応じた冷却効率
U(x,y)′が算出される。そして、断層面各部の温度T
(x,y)i+1が前回の計算値T(x,y)iを元に、実際の血流に
応じた冷却効率U(x,y)′を加味した所定の演算式によ
り繰り返し算出される(請求項)。
Further, according to the present invention, the temperature increase contribution rate D0 (x, y) due to the density when the contrast agent is not used and the temperature increase contribution rate D1 (x, y) due to the density when the contrast agent is used are actually used. A correction coefficient K (x, y) based on the blood flow of is calculated, and the cooling efficiency U (x, y) is multiplied by the correction coefficient K (x, y) to obtain the cooling efficiency U (x corresponding to the actual blood flow. , y) ′ is calculated. And the temperature T of each part of the fault plane
(x, y) i + 1 is repeatedly calculated based on the previously calculated value T (x, y) i by a predetermined calculation formula that takes into consideration the cooling efficiency U (x, y) ′ according to the actual blood flow. (Claim 2 ).

【0012】また、本発明によれば、被加温体の近傍位
置における電磁界強度を検出し、該検出結果に基づき高
周波エネルギー量が算出され、更にこの算出結果に応じ
て上記温度上昇係数βが演算される(請求項)。
Further, according to the present invention, the electromagnetic field intensity in the vicinity of the object to be heated is detected, the high frequency energy amount is calculated based on the detection result, and the temperature increase coefficient β is calculated according to the calculation result. Is calculated (claim 3 ).

【0013】また、本発明によれば、表示手段への温度
表示は色別で行われる(請求項)。また、表示手段へ
の温度表示は濃度別に行われる(請求項)。更に断層
面の適所に設置された温度計からの実際の温度で、演算
手段による演算温度に補正が加わえられる(請求項
)。
Further, according to the present invention, the temperature display on the display means is performed for each color (claim 4 ). Further, the temperature display on the display means is performed for each concentration (claim 5 ). Further, the actual temperature from a thermometer installed at a proper position on the fault plane can be used to correct the calculated temperature by the calculating means (claim 3
6 ).

【0014】[0014]

【実施例】図1は、本発明に係る断層面温度分布表示装
置が適用される高周波加温治療装置のブロック構成図を
示すものである。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS FIG. 1 shows a block diagram of a high-frequency heating treatment apparatus to which a tomographic plane temperature distribution display apparatus according to the present invention is applied.

【0015】高周波加温治療装置は体内温度を演算する
コンピュータ1、断層面の密度分布情報を取り込むカメ
ラ5及びアンプ6、高周波エネルギーを患部に照射する
ための加温部7、得られた体内温度を表示する表示部8
及び操作部9,91等から構成されている。
The high-frequency heating treatment apparatus is a computer 1 for calculating the temperature inside the body, a camera 5 and an amplifier 6 for taking in the density distribution information of the tomographic plane, a heating section 7 for irradiating the affected area with high-frequency energy, and the obtained body temperature. Display section 8
And the operation units 9 and 91.

【0016】コンピュータ1は、高周波加温治療装置の
全体を統括的に制御するもので、制御プログラムの他、
後述する演算式、生体に関する各種情報及び高周波エネ
ルギーの分布情報等が記憶されたROM2と処理内容を
一時的に保存記憶するためのRAM3とが接続されてい
る。
The computer 1 centrally controls the whole of the high-frequency heating treatment apparatus. In addition to the control program,
A ROM 2 in which arithmetic expressions to be described later, various kinds of information regarding a living body, distribution information of high-frequency energy, and the like are stored, and a RAM 3 for temporarily storing processing contents are connected.

【0017】4は断層面の密度(デンシティ)分布情報
を有する断層写真で、この断層写真4はCT(Computer
Tomography)スキャナーによるX線写真やMRI(Mag
netic Resonance Imaging)による組織の分布写真であ
る。なお、この密度情報は重量密度の他、水素密度に関
するものでもよい。
Reference numeral 4 is a tomographic photograph having density (density) distribution information of the tomographic plane. This tomographic photograph 4 is a CT (Computer).
Tomography) X-ray photograph and MRI (Mag
It is a distribution photograph of tissue by netic Resonance Imaging. The density information may be not only the weight density but also the hydrogen density.

【0018】カメラ5は撮像部に多数配列のCCD(Ch
arge-Coupled Device;固体撮像素子)を有し、上記断
層写真4を密度分布情報として読み取るもので、読み取
られた情報は電気信号に変換されて時系列信号でアンプ
6に入力される。断層写真4の密度分布はカメラ5で、
例えば256階調で読み取られ、アンプ6を介してコン
ピュータ1に導かれて、写真濃度に左右されないように
正規化すべく所要のレベル数、例えば16レベルに変換
されてRAM3に記憶される。
The camera 5 has a large number of CCDs (Ch
It has an arge-Coupled Device (solid-state imaging device) and reads the tomographic photograph 4 as density distribution information. The read information is converted into an electric signal and input to the amplifier 6 as a time series signal. The density distribution of the tomographic photograph 4 is obtained by the camera 5.
For example, it is read in 256 gradations, guided to the computer 1 via the amplifier 6, converted into a required number of levels, for example, 16 levels so as to be normalized so as not to be influenced by the photographic density, and stored in the RAM 3.

【0019】この16レベルの信号はD1〜D16のよ
うに属性が割り付けられ、係数への変更処理が容易に行
えるようにしている。コンピュータ1は取り込んだ撮像
信号を一旦RAM3あるいは別の表示メモリに記憶した
後、これを静止画で表示部8に表示する。また、この1
6レベルの撮像信号は、不図示の外部メモリに記憶さ
れ、後日の利用の便宜が図れるようにしている。
The 16-level signal is assigned an attribute like D1 to D16 so that the process of changing the coefficient can be easily performed. The computer 1 temporarily stores the captured image pickup signal in the RAM 3 or another display memory, and then displays it on the display unit 8 as a still image. Also this 1
The 6-level image pickup signal is stored in an external memory (not shown) for convenience of later use.

【0020】なお、断層の密度分布情報が、例えば通信
手段を介して伝送されるくる場合、あるいは外部メモリ
から取り込まれる場合には、カメラ5に代えてインター
フェースを用いてデータ受信すればよい。
When the density distribution information of the tomographic image is transmitted through, for example, communication means or is fetched from an external memory, the interface may be used instead of the camera 5 to receive the data.

【0021】操作部9はテンキーや操作指示のためのキ
ーを有するキーボード操作部で、取り込んだ断層写真4
の画像を必要に応じて適正な濃淡状態に補正して、後述
の生体特性との対応を取るものである。この濃淡補正
は、例えば断層写真4の最も密度の高い箇所を骨と見做
し、この骨の密度を基準にそれ以下の密度に対して順次
補正を加えるようにするものである。また、キーボード
操作部9は、後述するように各種情報の選択、例えば電
極のタイプを実際に使用される電極の種類に合致させる
べく選択指示するものである。
The operation unit 9 is a keyboard operation unit having a ten-key pad and keys for operating instructions.
The image is corrected to an appropriate shading state as necessary to correspond to the biological characteristics described later. In this grayscale correction, for example, the highest density portion of the tomographic photograph 4 is regarded as a bone, and the density of this bone is used as a reference to sequentially correct the density below that. Further, the keyboard operation section 9 is used to select various information as will be described later, for example, to select and instruct the electrode type to match the type of electrode actually used.

【0022】操作部91は表示部8に表示された画像に
対して操作可能なトラックボール(あるいはマウス)
で、表示部8に静止画表示された断層写真4に対して、
後述の温度演算を施す領域を指定するものである。な
お、演算領域の指定のための生体輪郭の抽出は画像処理
技術を利用してコンピュータ1が行うようにすることも
できる。
The operation unit 91 is a trackball (or mouse) operable on the image displayed on the display unit 8.
Then, for the tomographic photograph 4 displayed as a still image on the display unit 8,
The area is used to specify a temperature calculation to be described later. It should be noted that the extraction of the living body contour for designating the calculation area may be performed by the computer 1 using an image processing technique.

【0023】加温部7は公知の高周波エネルギー発生構
造を有するもので、1対の電極71,72、数MHzの
高周波発生部73、インピーダンス整合部74及び高周
波電力計75を備えてなるものである(図2、参照)。
上記1対の電極71,72は断層写真4に一致する生体
の横断面上で、患部121を挾むようにして密着配設さ
れ、この状態で、高周波発生部からの高周波電界あるい
は磁界等の高周波エネルギーを患部121を中心に照射
するものである。
The heating unit 7 has a known high-frequency energy generating structure, and is provided with a pair of electrodes 71, 72, a high-frequency generating unit 73 of several MHz, an impedance matching unit 74, and a high-frequency power meter 75. Yes (see FIG. 2).
The pair of electrodes 71, 72 are closely arranged so as to sandwich the affected part 121 on the cross section of the living body corresponding to the tomographic photograph 4, and in this state, high frequency energy such as high frequency electric field or magnetic field from the high frequency generating part is applied. The irradiation is performed mainly on the affected area 121.

【0024】上記1対の電極71,72から発せられる
高周波エネルギーの、電極間における分布状態(高周波
電界分布)は予めROM2に取り込まれている。ここ
で、高周波エネルギーの分布状態の予備測定について説
明する。
The distribution state (high-frequency electric field distribution) of the high-frequency energy emitted from the pair of electrodes 71, 72 between the electrodes is stored in the ROM 2 in advance. Here, the preliminary measurement of the distribution state of the high frequency energy will be described.

【0025】ダミーとして密度均一、かつ人体密度に近
似した、例えば寒天等を用い、その両側に電極71,7
2を密着配置した状態で、定格パワーの高周波エネルギ
ーを一定時間だけ放射する。一定時間放射後、高周波エ
ネルギーの照射断面を切断して、この断面をサーモグラ
フィ等の遠隔温度測定装置を用いて温度分布を計測す
る。供給される高周波エネルギーの強弱に応じて被加温
体の上昇温度は一義的に決まるから、密度均一の条件下
では断面温度分布は高周波エネルギーの強度分布に対応
している。そこで、この温度分布を高周波エネルギー分
布に置き換えてROM2に記憶する。
As the dummy, for example, agar, which has a uniform density and approximates the human body density, is used.
In the state where the two are closely arranged, high frequency energy of rated power is radiated for a certain period of time. After radiating for a certain period of time, the irradiation cross section of the high frequency energy is cut, and the temperature distribution is measured on this cross section using a remote temperature measuring device such as thermography. Since the rising temperature of the object to be heated is uniquely determined according to the strength of the supplied high-frequency energy, the cross-sectional temperature distribution corresponds to the strength distribution of the high-frequency energy under the condition of uniform density. Therefore, this temperature distribution is replaced with the high frequency energy distribution and stored in the ROM 2.

【0026】なお、患部121の位置等により大小、形
状等の異なる電極71,72が使用され、更に電極間寸
法が異なる場合には、それぞれの電極の種類、電極間寸
法に対して同様に分布情報を測定し、これらを電極の種
類に対応させてROM2に記憶する。
When electrodes 71, 72 having different sizes, shapes, etc. are used depending on the position of the affected part 121, etc., and the dimension between electrodes is different, the distribution is the same for each type of electrode and dimension between electrodes. Information is measured and stored in the ROM 2 in correspondence with the type of electrode.

【0027】次に、ROM2に予め取り込まれる生体の
各種特性について説明する。ROM2には下記表1に示
す各種特性値が予め取り込まれている。
Next, various characteristics of the living body which are loaded into the ROM 2 in advance will be described. Various characteristic values shown in Table 1 below are previously stored in the ROM 2.

【0028】[0028]

【表1】 [Table 1]

【0029】なお、(x,y)は断層面各部の位置を示すも
ので、例えばRAM3のアドレスに対応するものであ
る。
It should be noted that (x, y) indicates the position of each part of the tomographic plane, and corresponds to the address of the RAM 3, for example.

【0030】表1において、温度上昇寄与率e(x,y)は
高周波エネルギーの分布情報から一義的に決定される値
である。なお、この値e(x,y)は前述したサーモグラフ
ィの測定温度を正規化して得られた数値を採用したもの
でもよい。また、温度上昇寄与率D0(x,y)は密度に応じ
て決定される値で、前記カメラ5から取り込まれた造影
剤不使用時の密度情報(例えば、ρ(x,y))が、後述の
演算に際して変換設定される。
In Table 1, the temperature rise contribution rate e (x, y) is a value uniquely determined from the high frequency energy distribution information. The value e (x, y) may be a value obtained by normalizing the measurement temperature of the thermography described above. The temperature rise contribution rate D0 (x, y) is a value determined according to the density, and the density information (for example, ρ (x, y)) when the contrast agent is not used, which is captured from the camera 5, is Conversion is set in the calculation described later.

【0031】冷却効率U(T)は生体温度に関係する値
で、生体が加温される程、熱発散が大きくなることに着
目して設定されるものである。冷却効率U(T)は、各臓
器について経験的に求められており、この一般的な値を
用いることができる。しかし、上記放熱は主として血流
により行なわれ、しかも実際の血流状態は各臓器により
それぞれ異なるから、望ましくは血流による放熱効果を
加味した、より実際的な冷却効率U(T)′を用いるほう
がよい。なお、この冷却効率U(T)′については後述す
る。
The cooling efficiency U (T) is a value related to the temperature of the living body, and is set by paying attention to the fact that the more the living body is heated, the larger the heat dissipation becomes. The cooling efficiency U (T) is empirically determined for each organ, and this general value can be used. However, since the heat radiation is mainly performed by the blood flow, and the actual blood flow state is different for each organ, it is desirable to use a more practical cooling efficiency U (T) ′ in consideration of the heat radiation effect by the blood flow. Better. The cooling efficiency U (T) 'will be described later.

【0032】また、β0は高周波エネルギーの定格パワ
ーP0に対する温度上昇係数を示すもので、下記数3に
より算出される。
Further, β0 represents a temperature increase coefficient of the high frequency energy with respect to the rated power P0, and is calculated by the following expression 3.

【0033】[0033]

【数3】 [Equation 3]

【0034】因みに、この温度上昇係数β0はe(x,y)=
D0(x,y)=1のとき、1となる。
Incidentally, this temperature rise coefficient β0 is e (x, y) =
It becomes 1 when D0 (x, y) = 1.

【0035】なお、上記加温部7の高周波発生部73の
出力が常に定格パワーP0を出力するものでは、被加温
体の加温部の体積Vに応じて温度上昇係数β0を予め算
出し、この算出結果をROM2に記憶しておけばよい
が、出力が任意に可変できるものでは、実際に出力され
たパワーPに対する温度上昇係数βに変換する必要があ
る。この場合は、上記定格パワーP0に代えて出力パワ
ーPを用いて上記数3により、或いは予め算出された温
度上昇係数β0に定格パワーP0に対する実際の出力パワ
ーPの比率R(=P/P0)を乗じて温度上昇係数β
(=R・β0)を算出することができる。
If the output of the high frequency generator 73 of the heating unit 7 always outputs the rated power P0, the temperature rise coefficient β0 is calculated in advance according to the volume V of the heating unit of the object to be heated. The calculation result may be stored in the ROM 2, but if the output can be arbitrarily changed, it is necessary to convert the temperature rise coefficient β to the actually output power P. In this case, the output power P is used in place of the rated power P0, and the ratio R (= P / P0) of the actual output power P to the rated power P0 is calculated by the above Equation 3 or the temperature rise coefficient β0 calculated in advance. The temperature rise coefficient β
(= R · β0) can be calculated.

【0036】続いて、断層面温度の算出について説明す
る。各部の温度T(x,y)i+1は、下記数4を繰り返し演算
して予め算出される。
Next, the calculation of the cross section surface temperature will be described. The temperature T (x, y) i + 1 of each part is calculated in advance by repeatedly calculating the following Expression 4.

【0037】[0037]

【数4】 [Equation 4]

【0038】なお、上記初期温度T0は通常の体温計等
を使用して測定することができる。また、上記温度上昇
係数βは、上記のように加温体7の高周波発生部73の
出力パワーPに基づき上記数3乃至定格パワーP0に対
する温度上昇係数β0を補正することにより算出され
る。加温体7の高周波発生部73の出力が定格パーP0
の場合は、上記温度上昇係数βとして定格パワーP0に
対する温度上昇係数β0を用いればよい。特に、e(x,y)
=D0(x,y)=1のときは、温度上昇係数β0=1となる
ので、β0・e(x,y)・D0(x,y)=1となり、上記数4の第
2項を考慮することなく温度T(x,y)i+1を演算すること
ができる。
The initial temperature T0 can be measured using a normal thermometer or the like. Further, the temperature rise coefficient β is calculated by correcting the temperature rise coefficient β0 with respect to the equation 3 to the rated power P0 based on the output power P of the high frequency generator 73 of the heating element 7 as described above. The output of the high frequency generator 73 of the heating element 7 is the rated par P0.
In this case, the temperature increase coefficient β0 for the rated power P0 may be used as the temperature increase coefficient β. In particular, e (x, y)
= D0 (x, y) = 1, the temperature rise coefficient β0 = 1, so β0 · e (x, y) · D0 (x, y) = 1, and the second term of the above equation 4 The temperature T (x, y) i + 1 can be calculated without consideration.

【0039】ところで、上記加温部7はインピーダンス
整合部74を備えているが、電極71,72は生体(被
加温体)に圧接され、生体が動くことによりマッチング
状態が変化することから、電極71,72と生体間を完
全にマッチングさせることは困難である。このため、加
温部7の高周波発生部73から上記電極71,72を介
して生体に供給される高周波電力には伝送ロスが生じ易
くなっている。一方、生体の誘電加熱による温度上昇
は、高周波電界の誘電体損による発熱に起因するもの
で、上記電極71,72間の高周波エネルギー分布、す
なわち生体に実際に供給された高周波エネルギー量に直
接、関係するものである。
By the way, the heating section 7 is provided with the impedance matching section 74, but the electrodes 71 and 72 are pressed against the living body (the body to be heated) and the matching state changes due to the movement of the living body. It is difficult to perfectly match the electrodes 71, 72 and the living body. Therefore, the high-frequency power supplied from the high-frequency generator 73 of the heating unit 7 to the living body via the electrodes 71 and 72 is likely to cause a transmission loss. On the other hand, the temperature rise due to the dielectric heating of the living body is caused by the heat generation due to the dielectric loss of the high frequency electric field, and is directly related to the high frequency energy distribution between the electrodes 71 and 72, that is, the amount of the high frequency energy actually supplied to the living body. It is related.

【0040】従って、上記温度上昇係数βは、生体内に
実際に供給された高周波に基づき算出するのが好まし
い。すなわち、生体内に実際に供給された高周波電力
P′により数3を用いて算出された温度上昇係数β′を
数4に適用するのが好ましい。
Therefore, it is preferable to calculate the temperature rise coefficient β based on the high frequency actually supplied to the living body. That is, it is preferable to apply the temperature rise coefficient β ′ calculated by using the equation 3 by the high frequency power P ′ actually supplied to the living body to the equation 4.

【0041】電極71,72に実際に供給された高周波
電力は、一部が放射ロス(電極71,72から空間に電
磁波となって放射することによるロス)となるから、上
記高周波電力の内、生体内に実際に供給される高周波電
力P′、すなわち、誘電加熱に寄与する高周波電力P′
を直接、計測することは、一般に困難である。しかし、
出願人によれば、生体の近傍位置における電磁界の強度
が上記高周波電力P′と良好な相関性をし、該電磁界強
度から算出した温度上昇係数β′を用いて上記温度T
(x,y)i+1を演算すると、その演算結果と実際の温度との
誤差がより小さくなることが確かめられた。
A part of the high frequency power actually supplied to the electrodes 71 and 72 causes a radiation loss (a loss due to radiation from the electrodes 71 and 72 as electromagnetic waves in space). High frequency power P'actually supplied to the living body, that is, high frequency power P'that contributes to dielectric heating.
It is generally difficult to measure directly. But,
According to the applicant, the strength of the electromagnetic field in the vicinity of the living body has a good correlation with the high frequency power P ′, and the temperature T using the temperature rise coefficient β ′ calculated from the strength of the electromagnetic field.
It was confirmed that when (x, y) i + 1 is calculated, the error between the calculation result and the actual temperature becomes smaller.

【0042】従って、生体の近傍位置における電磁界の
強度Wは等価的に上記高周波電力P′として取り扱うこ
とが可能で、該電磁界強度Wと上記温度上昇係数β′と
の相関関数F(=k・W,k;比例係数)を実験等によ
り予め算出しておけば、上記数3に代えて該相関関数F
を用いることにより電磁界強度Wから上記温度上昇係数
β′を算出することができる。
Therefore, the intensity W of the electromagnetic field in the vicinity of the living body can be equivalently treated as the high frequency power P ', and the correlation function F (= F (=) of the electromagnetic field intensity W and the temperature rise coefficient β'. If k · W, k; a proportional coefficient) is calculated in advance by experiments or the like, the correlation function F can be replaced with the above-mentioned equation 3.
The temperature rise coefficient β ′ can be calculated from the electromagnetic field strength W by using

【0043】図2は、上記電磁界強度の計測方法を示す
概略図である。同図において、電磁界強度計13は、微
小コイルからなるセンサ部131と電磁界を受信するレ
シーバ132とから構成され、センサ部131は、生体
14の近傍位置に高周波磁界Hがコイル面を略直角に鎖
交するように配置されている。
FIG. 2 is a schematic diagram showing a method for measuring the electromagnetic field strength. In the figure, the electromagnetic field intensity meter 13 is composed of a sensor unit 131 composed of a minute coil and a receiver 132 that receives an electromagnetic field. The sensor unit 131 has a high-frequency magnetic field H on the coil surface in the vicinity of the living body 14. They are arranged so that they intersect at right angles.

【0044】加温部7の高周波発生部73からインピー
ダンス整合部74及び高周波電力計75を介して生体1
4を挟んで取り付けられた電極71,72に高周波電力
が供給されると、該電極71,72間に高周波電界E
(同図、実線矢印)が生じるとともに、該高周波電界E
の回りに高周波磁界H(同図、点線矢印)が生じる。電
磁界強度計13のセンサ部131は、上記電磁界の磁界
成分をピックアップし、レシーバ132で該高周波磁界
Hの強度Wが検出される。レシーバ132は電磁界強度
Wをエネルギー量(W/m2)、電界強度(V/m)又は磁界
強度(A/m)に換算して出力可能になされ、検出した高
周波磁界Hの強度Wを設定された単位で出力する。な
お、上記電磁界の電界成分を検出して電磁界の強度Wを
算出するようにしてもよい。
From the high frequency generator 73 of the heating unit 7 through the impedance matching unit 74 and the high frequency power meter 75, the living body 1
When high frequency power is supplied to the electrodes 71, 72 attached with the electrode 4 interposed therebetween, a high frequency electric field E is generated between the electrodes 71, 72.
(The same figure, solid line arrow) occurs, the high frequency electric field E
A high-frequency magnetic field H (indicated by a dotted arrow in the figure) is generated around. The sensor unit 131 of the electromagnetic field intensity meter 13 picks up the magnetic field component of the electromagnetic field, and the receiver 132 detects the intensity W of the high frequency magnetic field H. The receiver 132 converts the electromagnetic field strength W into an energy amount (W / m 2 ), an electric field strength (V / m) or a magnetic field strength (A / m), and is capable of outputting the detected strength W of the high frequency magnetic field H. Output in the set unit. The intensity W of the electromagnetic field may be calculated by detecting the electric field component of the electromagnetic field.

【0045】そして、実験等により予め算出された上記
相関関数Fを上記ROM2に記憶しておき、上記レシー
バ132の検出結果Wをコンピュータ1にフィードバッ
クさせるようにすれば、上記相関関数Fから該検出結果
Wを用いて上記温度上昇係数β′が算出され、更に数4
により各部の温度T(x,y)i+1が自動的に演算される。
If the correlation function F calculated in advance by experiments or the like is stored in the ROM 2 and the detection result W of the receiver 132 is fed back to the computer 1, the detection is performed from the correlation function F. The temperature rise coefficient β ′ is calculated using the result W, and
Thus, the temperature T (x, y) i + 1 of each part is automatically calculated.

【0046】なお、上記ROM2に上記相関関数Fに代
えて電磁界強度Wと上記温度上昇係数β′との変換テー
ブルを記憶しておき、該変換テーブルにより電磁界強度
Wから直接、温度上昇係数β′を求めるようにしてもよ
い。
A conversion table for the electromagnetic field strength W and the temperature rise coefficient β'is stored in the ROM 2 instead of the correlation function F, and the temperature rise coefficient is directly calculated from the electromagnetic field strength W by the conversion table. β ′ may be obtained.

【0047】次に、上記冷却効率U(T)′の算出方法に
ついて説明する。上記冷却効率U(T)′は、例えば各部
の実際の血流状態を調べ、実際の血流量に基づく放熱効
果を加味した補正係数K(x,y)を算出し、この補正係数
K(x,y)を上記冷却効率U(T)に乗じて算出することがで
きる。そして、上記補正係数K(x,y)は、例えば造影剤
を注入せずに撮影した断層写真4と造影剤を注入して血
流状態を一層顕在化した断層写真4とをそれぞれカメラ
5で読み取り、読み取った画像を比較して患部付近の臓
器の実際の血流状態を調べることにより算出することが
できる。すなわち、造影剤を使用して撮影した断層写真
4から得られる密度による温度上昇率をD1(x,y)、密
度の分解能をnとすると、補正係数K(x,y)は、下記数
5で表される。
Next, a method of calculating the cooling efficiency U (T) 'will be described. For the cooling efficiency U (T) ′, for example, the actual blood flow state of each part is examined, and the correction coefficient K (x, y) in consideration of the heat radiation effect based on the actual blood flow is calculated, and the correction coefficient K (x , y) can be calculated by multiplying the cooling efficiency U (T). The correction coefficient K (x, y) is calculated by, for example, the tomographic photograph 4 taken without injecting the contrast agent and the tomographic photograph 4 in which the blood flow state is further revealed by injecting the contrast agent by the camera 5. It can be calculated by reading and comparing the read images and examining the actual blood flow state of the organ near the affected area. That is, assuming that the temperature rise rate due to the density obtained from the tomographic photograph 4 taken using the contrast agent is D1 (x, y) and the resolution of the density is n, the correction coefficient K (x, y) is given by It is represented by.

【0048】[0048]

【数5】 [Equation 5]

【0049】従って、上記数4に代えて下記数6により
各部の温度T(x,y)i+1を算出することにより、より実際
に近い温度分布を算出することができる。
Therefore, by calculating the temperature T (x, y) i + 1 of each part by the following expression 6 instead of the above expression 4, a more realistic temperature distribution can be calculated.

【0050】[0050]

【数6】 [Equation 6]

【0051】上記のように、コンピュータ1は、取り込
んだ密度分布情報、すなわちD0(x,y)と、予め記憶され
た高周波エネルギー分布、すなわちe(x,y)に、加温部
7から供給される高周波エネルギー量に相当する値、す
なわち温度上昇係数β又はβ′等に基づいて断層面各部
の温度T(x,y)i+1を直前の温度T(x,y)iを元に算出す
る。
As described above, the computer 1 supplies the imported density distribution information, that is, D0 (x, y) and the prestored high frequency energy distribution, that is, e (x, y) from the heating unit 7. Based on the temperature T (x, y) i + 1 of each part of the fault plane based on the value corresponding to the amount of high frequency energy, that is, the temperature rise coefficient β or β ', etc. calculate.

【0052】算出された温度T(x,y)i+1は表示部8の対
応する位置(x,y)に表示される。この表示は、例えば温
度に対応する白黒のトーンで階調表示するものでもよい
し、あるいは温度に対応させて、高温側から所定温度幅
毎に、白、赤、ピンク、黄、緑、青、水色となるように
色別表示するようにしてもよい。また、特定位置、ある
いは患部とその周辺の指定位置を数値データのまま、例
えば表形式で表示するようにしてもよい。
The calculated temperature T (x, y) i + 1 is displayed at the corresponding position (x, y) on the display unit 8. This display may be, for example, gradation display in black and white tones corresponding to temperature, or in accordance with temperature, white, red, pink, yellow, green, blue You may make it display by color so that it may become light blue. Further, the specific position, or the designated position around the affected area and its periphery may be displayed as it is in numerical data, for example, in a tabular form.

【0053】図3は、表示部8の表示状態を示す図であ
る。図において、71,72は電極で、10は生体断層
面である。この生体断層面10内には、脊髄11、各臓
器12,12、…があり、その一部に患部121があ
る。なお、図中、トーンあるいは色別表示は省略してあ
るが、断層面の表示は前述した輪郭表示のみならず、密
度がある程度異なる境界に対しても輪郭表示するように
してもよい。これにより、各臓器間の密度の相違、また
臓器と該臓器内の患部とを輪郭線で明確に表示すること
ができ、臓器の位置や形状等の把握性も向上する。
FIG. 3 is a diagram showing a display state of the display section 8. In the figure, 71 and 72 are electrodes, and 10 is a biological tomographic plane. In the living body tomographic plane 10, there are a spinal cord 11, respective organs 12, 12, ... It should be noted that although the display by tone or color is omitted in the drawing, the display of the tomographic plane may be performed not only by the contour display described above but also by the contour display for a boundary having a certain density difference. As a result, it is possible to clearly display the difference in density between the organs, the organs and the affected part in the organs by the contour lines, and the graspability of the position, shape, etc. of the organs is improved.

【0054】図4は、高周波加温治療装置による温度分
布表示の動作を説明するフローチャートである。
FIG. 4 is a flow chart for explaining the operation of displaying the temperature distribution by the high frequency heating treatment apparatus.

【0055】先ず、電源(パワー)が投入され(S1)、
この状態で、例えば断層写真4であるCTフィルムがカ
メラ5の前面適所にセットされる(S3)。そして、カ
メラ5は内蔵する測光手段や自動焦点調節手段により、
ピント、ズーム及び絞りが自動調整される(S5)。こ
の作業は操作者がマニュアルで設定するようにしてもよ
い。そして、キーボード操作部9からの操作によって、
CTフィルムの断層写真4が画像データとして取り込ま
れ(S7)、ここで必要に応じて濃淡補正を施した後、
コンピュータ1により、256階調が16階調に変換さ
れる(S9)。
First, the power source is turned on (S1),
In this state, for example, a CT film, which is a tomographic photograph 4, is set at an appropriate position on the front surface of the camera 5 (S3). Then, the camera 5 has a built-in photometric unit and an automatic focus adjusting unit.
The focus, zoom and aperture are automatically adjusted (S5). This operation may be manually set by the operator. Then, by the operation from the keyboard operation unit 9,
The tomographic photograph 4 of the CT film is taken in as image data (S7), and is subjected to shading correction if necessary,
The computer 1 converts 256 gradations into 16 gradations (S9).

【0056】次いで、取り込まれた画像からトラックボ
ール91による指示により不要な画像データが削除さ
れ、前記数4の演算領域を指定するための生体の輪郭が
抽出される(S11)。この画像輪郭抽出処理は画像処
理技術を用いて自動的に行うようにしてもよい。
Then, unnecessary image data is deleted from the captured image according to an instruction from the trackball 91, and the contour of the living body for designating the calculation region of the equation 4 is extracted (S11). This image contour extraction processing may be automatically performed using an image processing technique.

【0057】続いて、実際に使用される電極71,72
の種類をキーボード操作部9で選択設定し、更に電極位
置、傾きを指示する(S15)。電極位置は、電極間寸
法であり、傾きはROM2に記憶されている高周波エネ
ルギー分布に電極傾斜に合った回転を与えるためであ
る。
Next, the electrodes 71, 72 actually used
The type of is selected and set by the keyboard operating unit 9, and the electrode position and inclination are further instructed (S15). The electrode position is the dimension between the electrodes, and the inclination is for giving the high-frequency energy distribution stored in the ROM 2 rotation suitable for the electrode inclination.

【0058】この後、前記数4に基づいて温度算出のた
めの計算が開始される(S17)。すなわち、先ず、加
温部7を駆動させ、パワー調整が行われた後、この高周
波エネルギーP、すなわち投入パワーPのデータを受信
し、それを元に入力されている、例えば電界分布、生体
の密度分布により、S11で指定された領域内に対して
各画素毎に温度を計算する(S19)。なお、高周波エ
ネルギーが常に定格出力される装置である場合には、前
記係数βを一定値、例えば“1”として、すなわち考慮
することなく、演算処理するようにしてもよい。
Thereafter, the calculation for calculating the temperature is started based on the equation 4 (S17). That is, first, after the heating unit 7 is driven and power adjustment is performed, the high-frequency energy P, that is, the input power P data is received, and is input based on the received data, such as the electric field distribution and the biological body. Using the density distribution, the temperature is calculated for each pixel in the area specified in S11 (S19). In the case of a device in which high-frequency energy is constantly output at a rated value, the coefficient β may be set to a constant value, for example, “1”, that is, may be calculated without consideration.

【0059】次いで、上記数4により得られた温度分布
は対応する色別信号に変換されてモニターである表示部
8に表示される(S21)。このとき、画像輪郭と温度
分布との位置関係を対応させておくと、重畳表示ができ
る。そして、指定範囲に対して温度計算が終了したら、
最初の画素位置に戻って計算を繰り返し実行する(S2
3)。このように計算を繰り返すことで、前回の温度を
元に今回の温度を継続的に算出することが可能となる。
Next, the temperature distribution obtained by the above equation 4 is converted into corresponding color signals and displayed on the display unit 8 which is a monitor (S21). At this time, if the positional relationship between the image contour and the temperature distribution is made to correspond to each other, the superimposed display can be performed. And when the temperature calculation is completed for the specified range,
The calculation is repeated by returning to the first pixel position (S2).
3). By repeating the calculation in this way, it is possible to continuously calculate the current temperature based on the previous temperature.

【0060】なお、上記数6に基づいて各部の温度を演
算する場合は、上記ステップS3〜S9において、造影
剤を使用していないCTフィルムの画像データを取り込
り、ステップS9とステップS11間で造影剤を使用し
ているCTフィルムの画像データを取り込むとともに、
両画像データから補正係数K(x,y)を算出し、ステップ
S17及びステップS19で数4に代えて数6に基づい
て温度算出のための計算を行なうようにするとよい。
When calculating the temperature of each part based on the above equation 6, in steps S3 to S9, the image data of the CT film which does not use the contrast agent is taken in, and between steps S9 and S11. While capturing the image data of the CT film that uses the contrast agent in
The correction coefficient K (x, y) may be calculated from both image data, and the calculation for temperature calculation may be performed based on the equation 6 instead of the equation 4 in steps S17 and S19.

【0061】また、ステップS19〜S23の処理を所
定時間、例えば1分置きに行うようにすれば、温度上昇
変化を定量的に把握することが容易となる。あるいは、
ステップS23の終了後直ちに最初の画素に戻って温度
計算を開始するようにしてもよいし、所要の計算周期を
マニュアル設定し得るようにしてもよい。更に、温度の
上昇勾配に応じて計算周期を適宜可変にして、温度変化
に効果的に追尾し得るようにしてもよい。このようにす
ることで、よりリアルタイムに近い温度表示が可能とな
る。この場合、前記β,β′、e(x,y)、D(x,y)、D0
(x,y)、D1(x,y)、U(T)の各値は計算周期に対応した複
数の数値テーブルが記憶されていることが必要となる。
Further, if the processes of steps S19 to S23 are performed every predetermined time, for example, every one minute, it becomes easy to quantitatively grasp the temperature rise change. Alternatively,
Immediately after the end of step S23, the temperature calculation may be started by returning to the first pixel, or the required calculation cycle may be manually set. Further, the calculation cycle may be appropriately changed according to the rising gradient of the temperature so that the change in the temperature can be effectively tracked. By doing so, it is possible to display the temperature in more real time. In this case, β, β ', e (x, y), D (x, y), D0
For each value of (x, y), D1 (x, y), and U (T), it is necessary to store a plurality of numerical tables corresponding to the calculation cycle.

【0062】なお、本実施例では、今回の温度を前回の
温度を元に計算したが、本発明は、かかる計算方法に限
定されず、例えばその都度、初期温度を元に現在の温度
を算出する方法を採用してもよい。
In the present embodiment, the current temperature is calculated based on the previous temperature, but the present invention is not limited to this calculation method, and the current temperature is calculated based on the initial temperature each time. The method of doing may be adopted.

【0063】また、断層面内の適所(x,y)に1本あるい
は少数本だけ、例えば針状温度計を刺し込むようにし、
この温度計からの実際の温度と演算による求まる温度と
のずれを比較し、ずれがあるときは、そのずれを補正す
るようにすれば、より正確な温度表示が可能となる。
In addition, one or a small number, for example, a needle thermometer, is inserted at an appropriate position (x, y) in the fault plane,
By comparing the deviation between the actual temperature from this thermometer and the temperature obtained by calculation, and if there is a deviation, it is possible to correct the deviation so that a more accurate temperature display is possible.

【0064】[0064]

【発明の効果】以上説明したように、本発明によれば、
高周波により誘電加熱される被加温体の温度分布を該被
加温体の断層面に重畳して表示する断層面温度分布表示
装置であって、被加温体に供給される高周波エネルギー
量に応じた温度上昇係数β、冷却効率U(T)を加味し
て、前回の温度T(x,y)iを元に今回の温度T(x,y)i+1を
繰り返し求めるようにしたので、より実際の温度に近い
ものにすることができる。
As described above, according to the present invention,
A tomographic plane temperature distribution display device for superimposing and displaying a temperature distribution of a heated body, which is dielectrically heated by a high frequency, on a tomographic plane of the heated body. Since the temperature increase coefficient β and the cooling efficiency U (T) are taken into consideration, the current temperature T (x, y) i + 1 is repeatedly obtained based on the previous temperature T (x, y) i. , Can be closer to the actual temperature.

【0065】また、上記冷却効率U(T)を、造影剤不使
用時の密度による温度上昇寄与率D0(x,y)と造影剤使用
時の密度による温度上昇寄与率D1(x,y)とから得られる
補正係数K(x,y)を乗じて、実際の血流による放熱効果
を加味したより実際的な冷却効率U(T)′に補正するよ
うにしたので、加温体の各部の温度の算出精度がより向
上する。
Further, the cooling efficiency U (T) is calculated by using the density increase contribution ratio D0 (x, y) due to the density when the contrast agent is not used and the temperature increase contribution ratio D1 (x, y) due to the density when the contrast agent is used. Since the correction coefficient K (x, y) obtained from is multiplied to correct the cooling efficiency U (T) 'in consideration of the heat radiation effect of the actual blood flow, each part of the heating element is corrected. The calculation accuracy of the temperature is further improved.

【0066】また、被加温体の近傍位置における電磁界
の強度を検出し、該検出結果に基づき被加温体に供給さ
れる高周波エネルギー量を算出するようにしたので、こ
の高周波エネルギー量に対する温度上昇係数βはより実
際のものに近くなり、被加温体各部の温度の算出精度が
更に向上する。
Further, the intensity of the electromagnetic field in the vicinity of the object to be heated is detected, and the amount of high frequency energy supplied to the object to be heated is calculated based on the detection result. The temperature rise coefficient β becomes closer to the actual one, and the calculation accuracy of the temperature of each part to be heated is further improved.

【0067】また、求めた温度を濃度別や色別で段階表
示するので、温度把握が容易となる。更に、実際の温度
との補正ができるようにしているので、より正確な温度
表示が可能となる。
Further, since the obtained temperature is displayed in stages by density and color, it is easy to grasp the temperature. Further, since it is possible to correct the actual temperature, more accurate temperature display can be performed.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明に係る断層面温度分布表示装置が適用さ
れる高周波加温治療装置のブロック構成図を示すもので
ある。
FIG. 1 is a block configuration diagram of a high-frequency heating treatment device to which a tomographic plane temperature distribution display device according to the present invention is applied.

【図2】生体の近傍位置における電磁界の強度の計測方
法を示す概略図である。
FIG. 2 is a schematic diagram showing a method for measuring the strength of an electromagnetic field at a position near a living body.

【図3】表示部の表示状態を示す図である。FIG. 3 is a diagram showing a display state of a display unit.

【図4】本発明が適用される高周波加温治療装置による
温度分布表示の動作を説明するフローチャートである。
FIG. 4 is a flowchart illustrating an operation of displaying a temperature distribution by the high-frequency heating treatment device to which the present invention is applied.

【符号の説明】[Explanation of symbols]

1 コンピュータ 2 ROM 3 RAM 4 断層写真 5 カメラ 6 アンプ 7 加温部 71,72 電極 73 高周波発生部 74 インピーダンス整合部 75 電力計 8 表示部 9,91 操作部 13 電磁界強度計 131 センサ部 132 レシーバ部 14 生体 1 computer 2 ROM 3 RAM 4 tomography 5 camera 6 amplifier 7 heating section 71,72 electrodes 73 High frequency generator 74 Impedance matching section 75 electricity meter 8 Display 9,91 Operation part 13 Electromagnetic field strength meter 131 Sensor part 132 Receiver section 14 living body

フロントページの続き (56)参考文献 特開 平4−49977(JP,A) 特開 平5−245217(JP,A) 特開 平5−31080(JP,A) 特開 平5−95909(JP,A) 特開 昭63−252172(JP,A) 実開 平1−107307(JP,U) (58)調査した分野(Int.Cl.7,DB名) A61N 1/40 A61B 5/04 Continuation of front page (56) Reference JP-A-4-49977 (JP, A) JP-A-5-245217 (JP, A) JP-A-5-31080 (JP, A) JP-A-5-95909 (JP , A) JP-A-63-252172 (JP, A) Actual Kaihei 1-107307 (JP, U) (58) Fields investigated (Int.Cl. 7 , DB name) A61N 1/40 A61B 5/04

Claims (6)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 高周波により誘電加熱される被加温体の
温度分布を該被加温体の断層面に重畳して表示する断層
面温度分布表示装置であって、上記被加温体の断層面の
密度分布情報を取り込む情報取込手段と、この情報取込
手段で取り込んだ密度分布情報を温度上昇寄与率D0(x,
y)に変換して、記憶する第1の記憶手段と、高周波を印
加する電極間の高周波エネルギー分布におけるエネルギ
ー強度に起因する温度上昇寄与率e(x,y)を記憶する第
2の記憶手段と、上記被加温体に供給される高周波エネ
ルギー量を検出する検出手段と、検出された高周波エネ
ルギー量に応じた温度上昇係数βを算出する第1の演算
手段と、断層面各部の温度T(x,y)i+1を下記数1の演算
式を繰り返し実行して算出する第2の演算手段と、算出
された各部の温度T(x,y)i+1を表示する表示手段とを備
えてなる断層面温度分布表示装置。 【数1】
1. A heating target object that is dielectrically heated by high frequency
A fault whose temperature distribution is superimposed and displayed on the fault plane of the heated body
A surface temperature distribution display device, comprising:
Information capturing means for capturing density distribution information, and this information capturing
The density increase information captured by the means is used to calculate the temperature rise contribution ratio D0 (x,
The first storage means for converting to y) and storing the high frequency
Energy in high-frequency energy distribution between applied electrodes
-The temperature rise contribution rate e (x, y) due to the intensity is stored
2 storage means and high-frequency energy supplied to the heated body
The detection means for detecting the amount of energy and the detected high frequency energy
First calculation for calculating the temperature rise coefficient β according to the amount of rugies
Means and the temperature T (x, y) i + 1 of each part of the fault plane is calculated by the following equation 1.
Second calculation means for repeatedly executing the formula to calculate, and
And a display means for displaying the temperature T (x, y) i + 1 of each part
A fault surface temperature distribution display device. [Equation 1]
【請求項2】 請求項1記載の断層面温度分布表示装置
において、上記情報取込手段で取り込んだ造影剤使用時
の密度分布情報を温度上昇寄与率D1(x,y)に変換して、
記憶する第3の記憶手段を備え、上記第2の演算手段
は、断層面の各部の温度T(x,y)i+1を上記数1に代えて
下記数2の演算式を繰り返し実行して算 出することを特
徴とする断層面温度分布表示装置。 【数2】
2. A fault surface temperature distribution display device according to claim 1.
At the time of using the contrast agent captured by the above information capture means,
The density distribution information of is converted into the temperature rise contribution rate D1 (x, y),
A second storage means for storing the third storage means;
Replaces the temperature T (x, y) i + 1 of each part of the fault plane with the above equation 1.
Japanese to repeatedly output calculated by executing the following Equation 2 arithmetic expression
Display device for temperature distribution of fault surface. [Equation 2]
【請求項3】 請求項1又は2記載の断層面温度分布表
示装置において、上記高周波エネルギーの検出手段は、
被加温体の近傍位置における電磁界の強度を検出する電
磁界検出手段と、検出された電磁界強度から上記高周波
エネルギー量を算出する演算手段とからなることを特徴
とする断層面温度分布表示装置。
3. A fault surface temperature distribution table according to claim 1 or 2.
In the apparatus shown, the high-frequency energy detecting means,
An electromagnetic wave that detects the strength of the electromagnetic field in the vicinity of the object to be heated.
From the magnetic field detection means and the detected electromagnetic field strength, the high frequency
It is characterized by comprising an arithmetic means for calculating the amount of energy
A device for displaying temperature distribution on a fault plane.
【請求項4】 前記表示手段への表示は前記演算結果を
色信号で行うものであることを特徴とする請求項1〜3
のいずれかに記載の断層面温度分布表示装置。
4. The calculation result is displayed on the display means.
4. The method according to claim 1, wherein the color signal is used.
The fault surface temperature distribution display device according to any one of 1.
【請求項5】 前記表示手段への表示は前記演算結果を
濃度信号で行うものであることを特徴とする請求項1〜
3のいずれかに記載の断層面温度分布表示装置。
5. The calculation result is displayed on the display means.
It is performed by using a density signal.
3. The fault surface temperature distribution display device according to any one of 3 above.
【請求項6】 請求項1〜3のいずれかに記載の断層面
温度分布表示装置において、前記断層面の適所(x,y)に
1または複数の温度計を設置し、この測定温度値で前記
演算手段の演算結果に補正を加えることを特徴とする断
層面温度分布表示装置。
6. A tomographic plane according to claim 1.
In the temperature distribution display device, at the appropriate place (x, y) on the fault plane
Install one or more thermometers and use the measured temperature value
A feature characterized by adding a correction to the calculation result of the calculation means.
Layer surface temperature distribution display device.
JP23052893A 1992-09-16 1993-09-16 Fault surface temperature distribution display method and apparatus Expired - Lifetime JP3476871B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP23052893A JP3476871B2 (en) 1992-09-16 1993-09-16 Fault surface temperature distribution display method and apparatus

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
JP24644792 1992-09-16
JP4-246447 1992-09-16
JP15937393 1993-06-29
JP5-159373 1993-06-29
JP23052893A JP3476871B2 (en) 1992-09-16 1993-09-16 Fault surface temperature distribution display method and apparatus

Publications (2)

Publication Number Publication Date
JPH0767969A JPH0767969A (en) 1995-03-14
JP3476871B2 true JP3476871B2 (en) 2003-12-10

Family

ID=27321532

Family Applications (1)

Application Number Title Priority Date Filing Date
JP23052893A Expired - Lifetime JP3476871B2 (en) 1992-09-16 1993-09-16 Fault surface temperature distribution display method and apparatus

Country Status (1)

Country Link
JP (1) JP3476871B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4543219B2 (en) * 2004-02-19 2010-09-15 学校法人明治大学 Method and apparatus for specifying heating mode for local heating
JP4534053B2 (en) * 2004-02-19 2010-09-01 学校法人明治大学 Local heating device
CN114618083A (en) * 2020-12-14 2022-06-14 澳美力科技(成都)有限公司 Novel adaptive radio frequency variable subcutaneous tissue radiation volume control method

Also Published As

Publication number Publication date
JPH0767969A (en) 1995-03-14

Similar Documents

Publication Publication Date Title
JP4545291B2 (en) Method and apparatus for controlling the dynamic range of digital diagnostic images
US7996066B2 (en) Topographic optical infrared tomography system for biophysical imaging with infrared diagnostic exploratory algorithm sequencing (IDEAS) scripting language
EP0904732B1 (en) X-Ray computed tomography method and apparatus
US20030149380A1 (en) Ultrasound treatment apparatus
US8030922B2 (en) Magnetic resonance tomography method and apparatus for correcting drift of the B0 field in a temperature map
JP2000107159A (en) X-ray imaging system, X-ray imaging image display method, and storage medium
JP2003515825A (en) Method of generating a gray scale transfer function for use in displaying digital x-ray images
JP7209552B2 (en) X-RAY DIAGNOSTIC APPARATUS, CONTROL METHOD OF X-RAY DIAGNOSTIC APPARATUS, AND IMAGE PROCESSING APPARATUS
JP3476871B2 (en) Fault surface temperature distribution display method and apparatus
CN108024774A (en) Spectral Imaging Phantoms and Methods
JPH10328318A (en) Radiotherapy system
JP2004135976A (en) Image acquisition device, image acquisition method, recording medium, and program
Blad et al. Quantitative assessment of impedance tomography for temperature measurements in hyperthermia
Cetas et al. Methods of hyperthermia control
JPH08146540A (en) Image processing device and image data originating device
CN119655867A (en) Real-time temperature and image feedback based data ablation monitoring system and method
JPH11319123A (en) Radiotherapy system
JP2714807B2 (en) Radiation image processing device
JP3827851B2 (en) Image information network
JP2006311922A (en) X-ray equipment
JPH06209906A (en) Optical imaging device
JPH02289229A (en) Temperature distribution measuring device
Soceanu et al. [P047] Ultrasound imaging: Towards smaller portable systems using signal and image processing tools
JPH08103440A (en) Photographing control method for x-ray ct system
JP2000157534A (en) X-ray CT system

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20030520

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20030826

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070926

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080926

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080926

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090926

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100926

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110926

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120926

Year of fee payment: 9