CN103829967B - Estimate device and method, the method determining interest district and the CT equipment in temperature field - Google Patents
Estimate device and method, the method determining interest district and the CT equipment in temperature field Download PDFInfo
- Publication number
- CN103829967B CN103829967B CN201210492598.1A CN201210492598A CN103829967B CN 103829967 B CN103829967 B CN 103829967B CN 201210492598 A CN201210492598 A CN 201210492598A CN 103829967 B CN103829967 B CN 103829967B
- Authority
- CN
- China
- Prior art keywords
- point
- energy
- reference point
- increment
- temperature increment
- 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.)
- Active
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention discloses a kind of device and method, the method determining interest district and the CT equipment of estimating temperature field. The device in this estimation temperature field comprises: one specific some energy estimation unit, and it is for estimating the energy of a specific point; One reference point energy estimation unit, it is for estimating the energy of a reference point, the contiguous described thermal source of wherein said reference point and at described reference point place CT image without artifact; One specific some temperature increment estimation unit, it is for the temperature increment of the increment described specific point of estimation of the CT value according to described specific point; One reference point temperature increment estimation unit, it estimates the temperature increment of described reference point for the increment of the CT value according to described reference point; One specific some temperature increment amending unit, it is for the temperature increment of specific point described in the energy of the energy according to described specific point, described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point.
Description
Technical field
The present invention relates to medical image field, particularly computer tomography.
Background technology
In the many interventional therapy adopting noncontact energy, such as laser, microwave and radio frequency, the too high heat that noncontact energy is introduced is most important. Typically, tumor tissues can be heated to more than 60 DEG C. According to recent research, such as G.D.Pandeya, JH.G.M.Klaessens, M.JW.Grenter, B.Schmidt, T.Flohr, etc., Feasibilityofcomputedtomographybasedthermometryduringint erstitiallaserheatinginbovineliver [J] .EurRadiol (2011) 21:1733-1738, it is feasible for carrying out " observation " temperature distribution by computer tomography in interventional therapy.
But, utilize the precision of the temperature distribution that computer tomography obtains not enough by image artifacts (such as metal artifacts etc.) restriction at present.
Summary of the invention
In view of this, the present invention proposes a kind of device and method, the method determining interest district and the CT equipment of estimating temperature field, improve the precision of Temperature estimate further.
According to the first aspect of the invention, a kind of device estimating temperature field is provided, described temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field, and described device comprises: one specific some energy estimation unit, and it is for estimating the energy of a specific point; One reference point energy estimation unit, it is for estimating the energy of a reference point, the contiguous described thermal source of wherein said reference point and at described reference point place CT image without artifact; One specific some temperature increment estimation unit, it is for the temperature increment of the increment described specific point of estimation of the CT value according to described specific point; One reference point temperature increment estimation unit, it estimates the temperature increment of described reference point for the increment of the CT value according to described reference point; One specific some temperature increment amending unit, it is for the temperature increment of specific point described in the energy of the energy according to described specific point, described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point.
Preferably, in said apparatus, described thermal source is positioned at intervalProbe, L is that pin is long, and described specific some energy estimation unit utilizes formulaEstimate the energy of described specific point, wherein QsIt is the energy of described specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient of energy, ��sL () is the equivalent distances of described specific point to the heat point source at l place; Described reference point energy estimation unit utilizes formulaEstimate the energy of described reference point, wherein QrIt is the energy of described reference point, ��rL () is the equivalent distances of described reference point to the heat point source at l place.
Preferably, in said apparatus, described specific some temperature increment estimation unit utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of described specific point, wherein �� TsIt is the temperature increment of described specific point, �� HUsIt is the increment of the CT value at described specific some place, c1And c0It it is the temperature variant relation coefficient of CT value; Described reference point temperature increment estimation unit utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of described reference point, wherein �� TrIt is the temperature increment of described reference point, �� HUrIt it is the increment of the CT value at described reference point place.
Preferably, in said apparatus, described specific some temperature increment amending unit comprises: a judging unit, it is for utilizing formulaJudge that the temperature increment of described specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of described specific point, the temperature increment of described reference point, the described energy of specific point and the energy of described reference point respectively, �� is thresholding; One amending unit, it is for utilizing formulaRevising the temperature increment of the specific point needing correction, wherein �� is modifying factor.
According to the second aspect of the invention, a kind of CT equipment is provided, comprising the device in an interest district determining unit and estimation temperature field as described above, described interest district determining unit is used for the set of the specific point using the temperature increment thresholding predetermined more than after revising as interest district.
According to the third aspect of the invention we, it is provided that a kind of method estimating temperature field, described temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field, and described method comprises: the energy estimating a specific point; Estimate the energy of a reference point, the contiguous described thermal source of wherein said reference point and at described reference point place CT image without artifact; The increment of the CT value according to described specific point estimates the temperature increment of described specific point; The increment of the CT value according to described reference point estimates the temperature increment of described reference point; The temperature increment of specific point described in energy according to described specific point, the energy of described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point.
Preferably, in the above-mentioned methods, described thermal source is positioned at intervalProbe, L is that pin is long, and the step of the energy of the specific point of described estimation one comprises and utilizes formulaEstimate the energy of described specific point, wherein QsIt is the energy of described specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient of energy, ��sL () is the equivalent distances of described specific point to the heat point source at l place; The step of the energy of described estimation one reference point comprises and utilizes formulaEstimate the energy of described reference point, wherein QrIt is the energy of described reference point, ��rL () is the equivalent distances of described reference point to the heat point source at l place.
Preferably, in the above-mentioned methods, the increment of the described CT value according to described specific point is estimated that the step of temperature increment of described specific point comprises and is utilized formula �� Ts=c1����HUs+c0Estimate the temperature increment of described specific point, wherein �� TsIt is the temperature increment of described specific point, �� UsIt is the increment of the CT value at described specific some place, c1And c0It it is the temperature variant relation coefficient of CT value; The described increment according to the CT value of described reference point is estimated that the step of the temperature increment of described reference point comprises and is utilized formula �� Tr=c1����HUr+c0Estimate the temperature increment of described reference point, wherein �� TrIt is the temperature increment of described reference point, �� HUrIt it is the increment of the CT value at described reference point place.
Preferably, in the above-mentioned methods, described in the energy of the described energy according to described specific point, described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point, the step of the temperature increment of specific point comprises: utilize formulaJudge that the temperature increment of described specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of described specific point, the temperature increment of described reference point, the described energy of specific point and the energy of described reference point respectively, �� is thresholding; If desired revise, then utilize formulaRevising the temperature increment of described specific point, wherein �� is modifying factor.
According to the fourth aspect of the invention, it is provided that a kind of method determining interest district, comprising: the method in an estimation temperature field mentioned above; Using the set of the specific point of the temperature increment thresholding predetermined more than after correction as interest district.
The device in the estimation temperature field of the present invention and method, the method determining interest district and CT equipment utilization CT value estimate temperature increment, and utilize energy to be revised by temperature increment, it is to increase the precision of Temperature estimate. Estimate that the temperature field that obtains can be used for determining interest district, for further clinical application, as taken area-of-interest as the ultralow dosage imaging of target, the Minimally Interventional Therapy etc. that take area-of-interest as monitored object.
Accompanying drawing explanation
The preferred embodiments of the present invention will be described in detail by referring to accompanying drawing below, the person of ordinary skill in the art is more clear that the above-mentioned of the present invention and other feature and advantage, in accompanying drawing:
Fig. 1 is the schematic block diagram of the device in the estimation temperature field according to the first embodiment of the present invention.
Fig. 2 is the schematic block diagram of CT equipment according to the second embodiment of the present invention.
Fig. 3 is the indicative flowchart of the method in estimation temperature field according to the third embodiment of the invention.
Fig. 4 is the indicative flowchart of the method in Dingxing interest district really according to the fourth embodiment of the invention.
In above-mentioned accompanying drawing, the Reference numeral adopted is as follows:
100, device 112,212 judging unit in 201 estimation temperature fields
102,202 specific somes energy estimation unit 114,214 amending units
104,204 reference point energy estimation unit 200CT equipment interest district determining unit
106,206 specific somes temperature increment estimation units 216
108,208 reference point temperature increment estimation unit 300,400,401 method
S302��S304��
S306��S308��
S310��S312��
110,210 specific somes temperature increment amending unit S314; S402, step
S404��S406��
S408��S410��
S412��S414��
S416
Embodiment
For making the object, technical solutions and advantages of the present invention clearly, the present invention is described in more detail by the following examples.
Fig. 1 is the schematic block diagram of the device 100 in the estimation temperature field according to the first embodiment of the present invention. Estimate that the temperature field estimated by device 100 in temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field. As shown in Figure 1, device 100 comprises one specific some energy estimation unit 102, reference point energy estimation unit 104, one specific some temperature increment estimation unit 106, reference point temperature increment estimation unit 108 and one specific some temperature increment amending unit 110.
Specific some energy estimation unit 102 is for estimating the energy of a specific point. At this specific some place, may there is artifact in CT image. In the present embodiment, thermal source is positioned at intervalProbe, L is that pin is long, and specific some energy estimation unit 102 utilizes formulaEstimate the energy of specific point, wherein QsIt is the energy of specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient (relevant with organization type) of energy, ��sL () is the equivalent distances of specific point to the heat point source at l place. Can adopt application number be 201110169757.X Chinese invention application disclosed in method prediction probe position. It is noted that likely there are other shapes due to thermal source, therefore, other formulas above-mentioned formula can be substituted in other embodiments.
Reference point energy estimation unit 104 for estimating the energy of a reference point, the contiguous thermal source of this reference point and at reference point place CT image without artifact. In the present embodiment, reference point energy estimation unit 104 utilizes formulaEstimate the energy of reference point, wherein QrIt is the energy of reference point, ��rL () is the equivalent distances that reference point arrives the heat point source at l place.
Specific some temperature increment estimation unit 106 estimates the temperature increment of specific point for the increment of the CT value according to specific point. In the present embodiment, specific some temperature increment estimation unit 106 utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of specific point, wherein �� TsIt is the temperature increment of specific point, �� HUsIt is the increment of the CT value at specific some place, c1And c0It it is the temperature variant relation coefficient of CT value. �� HUsIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at specific some place.
Reference point temperature increment estimation unit 108 is for the temperature increment of the increment estimation reference point of the CT value according to reference point. In the present embodiment, reference point temperature increment estimation unit 108 utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of reference point, wherein �� TrIt is the temperature increment of reference point, �� HUrIt it is the increment of the CT value at reference point place. �� HUrIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at reference point place.
Specific some temperature increment amending unit 110 is for the temperature increment of the specific point of temperature increment correction of the energy of the energy according to specific point, reference point, the temperature increment of specific point and reference point. Generally speaking the temperature increment at same position place and energy are directly proportional, and available formula is expressed as follows:Wherein, T����T��-����Being the temperature of ��, ��-�� �� time a certain position respectively, Q is the energy of this position, and �� is the tissue density of this position, and C is specific heat. In the present embodiment, specific some temperature increment amending unit 110 comprises judging unit 112 and an amending unit 114. Judging unit 112 utilizes formulaJudge that the temperature increment of specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of specific point, the temperature increment of reference point, the specific energy of point and the energy of reference point respectively, �� is thresholding. If above-mentioned absolute value is greater than ��, so amending unit 114 utilizes formulaRevising the temperature increment of this specific point, wherein �� is modifying factor. �� and �� and therapeutic modality and/or organize relevant, obtain by the mode of experiment.
Fig. 2 is the schematic block diagram of CT equipment 200 according to the second embodiment of the present invention. CT equipment 200 comprises the device 201 in an interest district determining unit 216 and an estimation temperature field.
Estimate that the temperature field estimated by device 201 in temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field. As shown in Figure 2, device 201 comprises one specific some energy estimation unit 202, reference point energy estimation unit 204, one specific some temperature increment estimation unit 206, reference point temperature increment estimation unit 208 and one specific some temperature increment amending unit 210.
Specific some energy estimation unit 202 is for estimating the energy of a specific point. At this specific some place, may there is artifact in CT image. In the present embodiment, thermal source is positioned at intervalProbe, L is that pin is long, and specific some energy estimation unit 202 utilizes formulaEstimate the energy of specific point, wherein QsIt is the energy of specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient (relevant with organization type) of energy, ��sL () is the equivalent distances of specific point to the heat point source at l place. Can adopt application number be 201110169757.X Chinese invention application disclosed in method prediction probe position. It is noted that likely there are other shapes due to thermal source, therefore, other formulas above-mentioned formula can be substituted in other embodiments.
Reference point energy estimation unit 204 for estimating the energy of a reference point, the contiguous thermal source of this reference point and at reference point place CT image without artifact. In the present embodiment, reference point energy estimation unit 204 utilizes formulaEstimate the energy of reference point, wherein QrIt is the energy of reference point, ��rL () is the equivalent distances that reference point arrives the heat point source at l place.
Specific some temperature increment estimation unit 206 estimates the temperature increment of specific point for the increment of the CT value according to specific point. In the present embodiment, specific some temperature increment estimation unit 206 utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of specific point, wherein �� TsIt is the temperature increment of specific point, �� HUsIt is the increment of the CT value at specific some place, c1And c0It it is the temperature variant relation coefficient of CT value. �� HUsIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at specific some place.
Reference point temperature increment estimation unit 208 is for the temperature increment of the increment estimation reference point of the CT value according to reference point. In the present embodiment, reference point temperature increment estimation unit 208 utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of reference point, wherein �� TrIt is the temperature increment of reference point, �� HUrIt it is the increment of the CT value at reference point place. �� HUrIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at reference point place.
Specific some temperature increment amending unit 210 is for the temperature increment of the specific point of temperature increment correction of the energy of the energy according to specific point, reference point, the temperature increment of specific point and reference point. Generally speaking the temperature increment at same position place and energy are directly proportional, and available formula is expressed as follows:Wherein, T����T��-����Being the temperature of ��, ��-�� �� time a certain position respectively, Q is the energy of this position, and �� is the tissue density of this position, and C is specific heat. In the present embodiment, specific some temperature increment amending unit 210 comprises judging unit 212 and an amending unit 214. Judging unit 212 utilizes formulaJudge that the temperature increment of specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of specific point, the temperature increment of reference point, the specific energy of point and the energy of reference point respectively, �� is thresholding. If above-mentioned absolute value is greater than ��, so amending unit 214 utilizes formulaRevising the temperature increment of this specific point, wherein �� is modifying factor. �� and �� and therapeutic modality and/or organize relevant, obtain by the mode of experiment.
Interest district determining unit 216 for using the set of specific point of the temperature increment thresholding predetermined more than after revising as interest district. Such as, this predetermined thresholding can be set to 18 DEG C, 20 DEG C, 23 DEG C. After determining interest district, method disclosed in Chinese invention application that application number is 201110169757.X can be utilized to carry out scanning and image is rebuild.
Fig. 3 is the indicative flowchart of the method 300 in estimation temperature field according to the third embodiment of the invention. This temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field. Method 300 comprises step S302, step S304, step S306, step S308 and step S310.
In step s 302, the energy of a specific point is estimated. At this specific some place, may there is artifact in CT image. In the present embodiment, thermal source is positioned at intervalProbe, L is that pin is long, and step S302 comprises and utilizes formulaEstimate the energy of specific point, wherein QsIt is the energy of specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient (relevant with organization type) of energy, ��sL () is the equivalent distances of specific point to the heat point source at l place. Can adopt application number be 201110169757.X Chinese invention application disclosed in method prediction probe position. It is noted that likely there are other shapes due to thermal source, therefore, other formulas above-mentioned formula can be substituted in other embodiments.
In step s 304, estimate the energy of a reference point, the contiguous thermal source of this reference point and at reference point place CT image without artifact. In the present embodiment, step S304 comprises and utilizes formulaEstimate the energy of reference point, wherein QrIt is the energy of reference point, ��rL () is the equivalent distances that reference point arrives the heat point source at l place.
In step S306, estimate the temperature increment of specific point according to the increment of the CT value of specific point. In the present embodiment, step S306 comprises and utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of specific point, wherein �� TsIt is the temperature increment of specific point, �� HUsIt is the increment of the CT value at specific some place, c1And c0It it is the temperature variant relation coefficient of CT value. �� HUsIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at specific some place.
In step S308, estimate the temperature increment of reference point according to the increment of the CT value of reference point. In the present embodiment, step S308 utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of reference point, wherein �� TrIt is the temperature increment of reference point, �� HUrIt it is the increment of the CT value at reference point place. �� HUrIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at reference point place.
In step S310, according to the temperature increment of the temperature increment of the energy of the energy of specific point, reference point, specific point and the specific point of temperature increment correction of reference point. Generally speaking the temperature increment at same position place and energy are directly proportional, and available formula is expressed as follows:Wherein, T����T��-����Being the temperature of ��, ��-�� �� time a certain position respectively, Q is the energy of this position, and �� is the tissue density of this position, and C is specific heat. In the present embodiment, step S310 comprises step S312 and step S314. In step S312, utilize formulaJudge that the temperature increment of specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of specific point, the temperature increment of reference point, the specific energy of point and the energy of reference point respectively, �� is thresholding. If above-mentioned absolute value is greater than ��, so in step S314, utilize formulaRevising the temperature increment of this specific point, wherein �� is modifying factor. �� and �� and therapeutic modality and/or organize relevant, obtain by the mode of experiment.
Fig. 4 is the indicative flowchart of the method 400 in Dingxing interest district really according to the fourth embodiment of the invention. Determine that the method 400 in interest district comprises the method 401 and step S416 estimating temperature field. This temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field. Method 401 comprises step S402, step S404, step S406, step S408 and step S410.
In step S402, estimate the energy of a specific point. At this specific some place, may there is artifact in CT image. In the present embodiment, thermal source is positioned at intervalProbe, L is that pin is long, and step S402 comprises and utilizes formulaEstimate the energy of specific point, wherein QsIt is the energy of specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient (relevant with organization type) of energy, ��sL () is the equivalent distances of specific point to the heat point source at l place. Can adopt application number be 201110169757.X Chinese invention application disclosed in method prediction probe position. It is noted that likely there are other shapes due to thermal source, therefore, other formulas above-mentioned formula can be substituted in other embodiments.
In step s 404, estimate the energy of a reference point, the contiguous thermal source of this reference point and at reference point place CT image without artifact. In the present embodiment, step S404 comprises and utilizes formulaEstimate the energy of reference point, wherein QrIt is the energy of reference point, ��rL () is the equivalent distances that reference point arrives the heat point source at l place.
In step S406, estimate the temperature increment of specific point according to the increment of the CT value of specific point. In the present embodiment, step S406 comprises and utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of specific point, wherein �� TsIt is the temperature increment of specific point, �� HUsIt is the increment of the CT value at specific some place, c1And c0It it is the temperature variant relation coefficient of CT value. �� HUsIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at specific some place.
In step S408, estimate the temperature increment of reference point according to the increment of the CT value of reference point. In the present embodiment, step S408 utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of reference point, wherein �� TrIt is the temperature increment of reference point, �� HUrIt it is the increment of the CT value at reference point place. �� HUrIt can be such as the increment of the CT value of the CT image of the continuous twice sweep acquisition at reference point place.
In step S410, according to the temperature increment of the temperature increment of the energy of the energy of specific point, reference point, specific point and the specific point of temperature increment correction of reference point. Generally speaking the temperature increment at same position place and energy are directly proportional, and available formula is expressed as follows:Wherein, T����T��-����Being the temperature of ��, ��-�� �� time a certain position respectively, Q is the energy of this position, and �� is the tissue density of this position, and C is specific heat. In the present embodiment, step S410 comprises step S412 and step S414. In step S412, utilize formulaJudge that the temperature increment of specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of specific point, the temperature increment of reference point, the specific energy of point and the energy of reference point respectively, �� is thresholding. If above-mentioned absolute value is greater than ��, so in step S414, utilize formulaRevising the temperature increment of this specific point, wherein �� is modifying factor. �� and �� and therapeutic modality and/or organize relevant, obtain by the mode of experiment.
In step S416, using the set of the specific point of the temperature increment thresholding predetermined more than after correction as interest district. Such as, this predetermined thresholding can be set to 18 DEG C, 20 DEG C, 23 DEG C. After determining interest district, method disclosed in Chinese invention application that application number is 201110169757.X can be utilized to carry out scanning and image is rebuild.
It is noted that so-called energy is the appellation that those skilled in the art are usual in the present invention, strict, it should to be the heating power of unit volume.
The device in the estimation temperature field of the present invention and method, the method determining interest district and CT equipment utilization CT value estimate temperature increment, and utilize energy to be revised by temperature increment, it is to increase the precision of Temperature estimate. Estimate that the temperature field that obtains can be used for determining interest district, for further clinical application, as taken area-of-interest as the ultralow dosage imaging of target, the Minimally Interventional Therapy etc. that take area-of-interest as monitored object.
The foregoing is only the better embodiment of the present invention, not in order to limit the present invention, within the spirit and principles in the present invention all, any amendment of doing, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (8)
1. estimating the device in temperature field, described temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field, and described device comprises:
One specific some energy estimation unit, it is for estimating the energy of a specific point;
One reference point energy estimation unit, it is for estimating the energy of a reference point, the contiguous described thermal source of wherein said reference point and at described reference point place CT image without artifact;
One specific some temperature increment estimation unit, it is for the temperature increment of the increment described specific point of estimation of the CT value according to described specific point;
One reference point temperature increment estimation unit, it estimates the temperature increment of described reference point for the increment of the CT value according to described reference point;
One specific some temperature increment amending unit, it is for the temperature increment of specific point described in the energy of the energy according to described specific point, described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point;
Described specific some temperature increment amending unit comprises:
One judging unit, it is for utilizing formulaJudge that the temperature increment of described specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of described specific point, the temperature increment of described reference point, the described energy of specific point and the energy of described reference point respectively, �� is thresholding;
One amending unit, it is for utilizing formulaRevising the temperature increment of the specific point needing correction, wherein �� is modifying factor.
2. device as claimed in claim 1, is characterized in that, described thermal source is positioned at intervalProbe, L is that pin is long, and described specific some energy estimation unit utilizes formulaEstimate the energy of described specific point, wherein QsIt is the energy of described specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient of energy, ��sL () is the equivalent distances of described specific point to the heat point source at l place;
Described reference point energy estimation unit utilizes formulaEstimate the energy of described reference point, wherein QrIt is the energy of described reference point, ��rL () is the equivalent distances of described reference point to the heat point source at l place.
3. device as claimed in claim 1, is characterized in that, described specific some temperature increment estimation unit utilizes formula �� Ts=c1����HUs+c0Estimate the temperature increment of described specific point, wherein �� TsIt is the temperature increment of described specific point, �� HUsIt is the increment of the CT value at described specific some place, c1And c0It it is the temperature variant relation coefficient of CT value;
Described reference point temperature increment estimation unit utilizes formula �� Tr=c1����HUr+c0Estimate the temperature increment of described reference point, wherein �� TrIt is the temperature increment of described reference point, �� HUrIt it is the increment of the CT value at described reference point place.
4. a CT equipment, comprising an interest district determining unit and the device just like the estimation temperature field described in the arbitrary item of claims 1 to 3, described interest district determining unit is used for the set of the specific point of the temperature increment thresholding predetermined more than after using correction as interest district.
5. estimating the method in temperature field, described temperature field has a thermal source and is positioned at the scanned object in the CT device scan visual field, and described method comprises:
Estimate the energy of a specific point;
Estimate the energy of a reference point, the contiguous described thermal source of wherein said reference point and at described reference point place CT image without artifact;
The increment of the CT value according to described specific point estimates the temperature increment of described specific point;
The increment of the CT value according to described reference point estimates the temperature increment of described reference point;
The temperature increment of specific point described in energy according to described specific point, the energy of described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point;
Described in the energy of the described energy according to described specific point, described reference point, the described temperature increment of specific point and the temperature increment correction of described reference point, the step of the temperature increment of specific point comprises:
Utilize formulaJudge that the temperature increment of described specific point is the need of correction, wherein �� Ts����Tr��QsAnd QrBeing the temperature increment of described specific point, the temperature increment of described reference point, the described energy of specific point and the energy of described reference point respectively, �� is thresholding;
If desired revise, then utilize formulaRevising the temperature increment of described specific point, wherein �� is modifying factor.
6. method as claimed in claim 5, is characterized in that, described thermal source is positioned at intervalProbe, L is that pin is long, and the step of the energy of the specific point of described estimation one comprises and utilizes formulaEstimate the energy of described specific point, wherein QsIt is the energy of described specific point, Q0L () is the energy of the heat point source at l place, �� is the real absorption coefficient of energy, ��sL () is the equivalent distances of described specific point to the heat point source at l place;
The step of the energy of described estimation one reference point comprises and utilizes formulaEstimate the energy of described reference point, wherein QrIt is the energy of described reference point, ��rL () is the equivalent distances of described reference point to the heat point source at l place.
7. method as claimed in claim 5, is characterized in that, the increment of the described CT value according to described specific point is estimated that the step of the temperature increment of described specific point comprises and utilized formula �� Ts=c1����HUs+c0Estimate the temperature increment of described specific point, wherein �� TsIt is the temperature increment of described specific point, �� HUsIt is the increment of the CT value at described specific some place, c1And c0It it is the temperature variant relation coefficient of CT value;
The described increment according to the CT value of described reference point is estimated that the step of the temperature increment of described reference point comprises and is utilized formula �� Tr=c1����HUr+c0Estimate the temperature increment of described reference point, wherein �� TrIt is the temperature increment of described reference point, �� HUrIt it is the increment of the CT value at described reference point place.
8. determine the method in interest district, comprising:
Just like the method in the estimation temperature field according to any one of claim 5 to 7;
Using the set of the specific point of the temperature increment thresholding predetermined more than after correction as interest district.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210492598.1A CN103829967B (en) | 2012-11-28 | 2012-11-28 | Estimate device and method, the method determining interest district and the CT equipment in temperature field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210492598.1A CN103829967B (en) | 2012-11-28 | 2012-11-28 | Estimate device and method, the method determining interest district and the CT equipment in temperature field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103829967A CN103829967A (en) | 2014-06-04 |
| CN103829967B true CN103829967B (en) | 2016-06-01 |
Family
ID=50794087
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210492598.1A Active CN103829967B (en) | 2012-11-28 | 2012-11-28 | Estimate device and method, the method determining interest district and the CT equipment in temperature field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103829967B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101006939A (en) * | 2007-01-18 | 2007-08-01 | 上海交通大学 | Focal prostate cancer conformal radiofrequency ablation electrode |
| CN101859341A (en) * | 2009-04-13 | 2010-10-13 | 盛林 | Image-guided ablation surgery planning device |
| CN102440800A (en) * | 2010-10-13 | 2012-05-09 | 西门子公司 | Method for generating tomographic temperature map, control and computing unit and ct system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012120498A1 (en) * | 2011-03-04 | 2012-09-13 | Technion Research & Development | Non-invasive thermal treatment monitoring |
-
2012
- 2012-11-28 CN CN201210492598.1A patent/CN103829967B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101006939A (en) * | 2007-01-18 | 2007-08-01 | 上海交通大学 | Focal prostate cancer conformal radiofrequency ablation electrode |
| CN101859341A (en) * | 2009-04-13 | 2010-10-13 | 盛林 | Image-guided ablation surgery planning device |
| CN102440800A (en) * | 2010-10-13 | 2012-05-09 | 西门子公司 | Method for generating tomographic temperature map, control and computing unit and ct system |
Non-Patent Citations (2)
| Title |
|---|
| Feasibility of computed tomography based thermometry during interstitial laser heating in bovine liver;G.D.Pandeya,et al;《Eur Radiol》;20110324;第21卷(第8期);1733-1738 * |
| Noninvasive thermometry with a clinical x-ray CT scanner;B.G.Fallone,et al;《Medical Physics》;19820930;第9卷(第5期);715-721 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103829967A (en) | 2014-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kok et al. | Online adaptive hyperthermia treatment planning during locoregional heating to suppress treatment-limiting hot spots | |
| Liu et al. | CT imaging during microwave ablation: analysis of spatial and temporal tissue contraction | |
| CN101859341A (en) | Image-guided ablation surgery planning device | |
| Gavazzi et al. | Advanced patient-specific hyperthermia treatment planning | |
| WO2009118721A1 (en) | Method apparatus and system for analyzing thermal images | |
| Verhaart et al. | The relevance of MRI for patient modeling in head and neck hyperthermia treatment planning: a comparison of CT and CT‐MRI based tissue segmentation on simulated temperature | |
| Brock et al. | Improving image-guided target localization through deformable registration | |
| Li et al. | Improved hyperthermia treatment control using SAR/temperature simulation and PRFS magnetic resonance thermal imaging | |
| Oberacker et al. | Radiofrequency applicator concepts for thermal magnetic resonance of brain tumors at 297 MHz (7.0 Tesla) | |
| Bottauscio et al. | Assessment of computational tools for MRI RF dosimetry by comparison with measurements on a laboratory phantom | |
| Zou et al. | Referenceless MR thermometry—a comparison of five methods | |
| Salahi et al. | Miniature microwave applicator for murine bladder hyperthermia studies | |
| Neufeld et al. | Rapid method for thermal dose‐based safety supervision during MR scans | |
| Lopresto et al. | CT-based investigation of the contraction of ex vivo tissue undergoing microwave thermal ablation | |
| Aklan et al. | Regional deep hyperthermia: Quantitative evaluation of predicted and direct measured temperature distributions in patients with high-risk extremity soft-tissue sarcoma | |
| Numan et al. | Exploration of MR-guided head and neck hyperthermia by phantom testing of a modified prototype applicator for use with proton resonance frequency shift thermometry | |
| Zhang et al. | Simulation of temperature field for temperature-controlled radio frequency ablation using a hyperbolic bioheat equation and temperature-varied voltage calibration: a liver-mimicking phantom study | |
| Zhang et al. | 2D ultrasound and 3D MR image registration of the prostate for brachytherapy surgical navigation | |
| Krywonos et al. | MRI image-based FE modelling of the pelvis system and bladder filling | |
| Kok et al. | Prospective treatment planning to improve locoregional hyperthermia for oesophageal cancer | |
| Tanaka et al. | Cryosurgery planning using bubble packing in 3D | |
| Noetscher et al. | Visible Human Project® female surface based computational phantom (Nelly) for radio-frequency safety evaluation in MRI coils | |
| CN103829967B (en) | Estimate device and method, the method determining interest district and the CT equipment in temperature field | |
| Lin et al. | Detecting early breast tumour by finite element thermal analysis | |
| Roos et al. | A comprehensive motion analysis–consequences for high precision image-guided radiotherapy of esophageal cancer patients |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |