Intelligent control method and system for radio frequency tumor endogenous field thermal therapy system

Abstract

The invention relates to the technical field of medical equipment control, in particular to an intelligent control method and system for a radio frequency tumor endogenous field thermal therapy system, and the method comprises the steps: collecting the power and thermal response data of each transmitting channel, analyzing the trend, jointly judging the thermal state, generating a thermal state parameter set, sorting and suppressing the priority, optimizing a power plan and a channel activation sequence, and achieving the output control. According to the method, the change trend of power input in different time periods is dynamically collected, and linkage analysis of the temperature rise amplitude and the thermal diffusion rate of the target area is combined, so that the thermal response state of the target tissue can be judged in real time, and then thermal state information of different emission channels in different periods is identified and labeled; on the basis, deviation evaluation is carried out on a power request value and a historical output value, multi-factor sorting is carried out through inter-channel power request deviation and heat conduction characteristics, priority is suppressed to determine regulation and control of power distribution and an activation sequence, and power control optimization under thermal response dynamic changes of different organization areas is achieved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of medical equipment control, in particular to a radio frequency tumor endogenous field hyperthermia system intelligent control method and system. BACKGROUND

[0002] The radio frequency tumor endogenous field hyperthermia system intelligent control method belongs to the technical field of medical equipment control, which mainly covers the operation control method, control algorithm implementation form, device state management and monitoring, data acquisition and analysis process, and dynamic adjustment based on feedback mechanism of various medical equipment used in the disease diagnosis and treatment process.

[0003] Among them, the radio frequency tumor endogenous field hyperthermia system intelligent control method refers to setting fixed treatment parameters to control radio frequency power output, hyperthermia time and treatment area positioning in the radio frequency hyperthermia process, and using power preset and channel sequence setting to realize transmission control of radio frequency signal according to the given physical properties of the target tissue such as electrical conductivity, heat capacity and depth.

[0004] The prior art only controls by setting fixed radio frequency power parameters and channel sequence, and relies on preset values to implement power output on the target tissue in the hyperthermia process, which lacks real-time analysis capability of instantaneous power change trend and actual heat expansion dynamics. When the electrical or thermal properties of the target tissue differ or change, the fixed parameters cannot effectively cope with the uneven heat conduction or unstable temperature rise caused by individual differences, especially under the condition of multi-channel combined treatment, the different channels are prone to heat expansion interference and power resource allocation imbalance due to different tissue response characteristics. For example, heat accumulation occurs in a channel due to weak heat conductivity of the target tissue, but the system cannot dynamically adjust the output according to the feedback, which easily causes local tissue overheating or even burns, while other channels are inhibited too much, resulting in insufficient temperature rise during treatment, thereby affecting the overall efficacy and treatment safety. The one-way nature and lack of feedback regulation capability of the fixed control strategy are the main reasons for its response lag and insufficient regulation in complex tissue scenarios. SUMMARY

[0005] The purpose of the present application is to solve the shortcomings in the prior art and propose a radio frequency tumor endogenous field hyperthermia system intelligent control method and system.

[0006] In order to achieve the above purpose, the present application adopts the following technical scheme: a radio frequency tumor endogenous field hyperthermia system intelligent control method, comprising the following steps: S1: Under the endogenous field, collect the instantaneous power input of the specified period of the hyperthermia emission channel corresponding region, judge the power trend change in the adjacent period, and obtain the power trend judgment result; S2: under the endogenous field, the temperature rise amplitude and the heat diffusion rate of the current emission channel targeted treatment tissue area are collected, and the criticality judgment is performed in combination with the power trend judgment result in the corresponding period to obtain the targeted treatment tissue area heat response analysis result; S3: according to the targeted treatment tissue area heat response analysis result, the thermal state information of the radio frequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate is labeled, the power request value of the current period and the output power value of the last period of the target channel are collected, and the channel thermal state parameter set is generated; S4: the power request deviation between the power request value of the current period and the output power value of the last period in the channel thermal state parameter set is calculated, the suppression priority among channels is sorted according to the power request deviation, and the power suppression priority table is generated; S5: according to the power suppression priority table, the next period power emission plan is adjusted, and the radio frequency hyperthermia emission channel activation sequence is set, and the channel output behavior control result is obtained.

[0007] As a further scheme of the application, the power trend judgment result includes power input change rate, adjacent period power difference value, power trend type, the targeted treatment tissue area heat response analysis result includes temperature rise amplitude stability, heat diffusion rate characteristic, joint criticality state, the channel thermal state parameter set includes thermal state information label, current period power request value, last period output power value, the power suppression priority table includes channel order sequence, suppression priority level, and the number of channels participating in sorting, and the channel output behavior control result includes next period power emission value, channel activation sequence, and control adjustment state.

[0008] As a further scheme of the application, the power trend judgment result acquisition step is specifically: S111: the radio frequency hyperthermia emission channel currently in the active state is acquired, the power signal of the targeted treatment tissue area acted on by the channel is collected, the instantaneous power input value in the continuous period is acquired in the specified sampling period, and the continuous period power input value sequence is generated; S112: based on the power input pair of every two adjacent periods in the continuous period power input value sequence, the power change relationship between each period segment is established in order, the progressive trend of the power input value of the next period and the corresponding value of the previous period is judged, and the adjacent period power change state sequence is generated; S113: in combination with whether there is a trend segment meeting the continuous decline condition in the adjacent period power change state sequence, the overall trend structure of the current power change is identified, and the power trend judgment result is obtained.

[0009] As a further scheme of the application, the targeted treatment tissue area heat response analysis result acquisition step is specifically: S211: Based on the power trend judgment result continuously descending radio frequency hyperthermia emission channel, the temperature value in the current period is collected in the target channel corresponding to the targeted treatment tissue area, and the temperature change data of each recording point is extracted, the temperature spatial distribution gradient of the adjacent area is combined, the temperature rise amplitude and the heat diffusion rate in the specified time are calculated, and the heat conduction characteristic data of the targeted treatment tissue area is obtained; S212: Based on the heat conduction characteristic data of the targeted treatment tissue area and the power trend judgment result, the synchronism between the power input trend and the temperature diffusion behavior is identified, and the current period data combination is screened under the consistent trend condition, the period sample meeting the temperature rise stable and diffusion enhancement characteristics is extracted, and the heat expansion dominant trend interval is generated; S213: According to the heat expansion dominant trend interval, it is judged whether the targeted treatment tissue area is in the heat expansion dominant state, and the change amplitude and distribution characteristics of each trend data item in the corresponding state are combined to judge whether the target tissue area shows stable temperature rise and diffusion coordinated conduction characteristics, and the heat response analysis result of the targeted treatment tissue area is obtained.

[0010] As a further scheme of the application, the acquisition step of the channel thermal state parameter set is specifically: S311: According to the heat response analysis result of the targeted treatment tissue area, the radio frequency hyperthermia emission channel with stable state is marked as a to-be-switched state, and the temperature rise amplitude and heat diffusion rate stable channel are obtained; S312: Based on the temperature rise amplitude and heat diffusion rate stable channel, the power request value in the current period and the output power value of the last period are collected for each channel, and the two power data are respectively associated to the corresponding channel, and the radio frequency hyperthermia emission channel power request and output value group is generated. S313: According to the radio frequency hyperthermia emission channel power request and output value group, the power input sequence, the temperature rise amplitude and the heat diffusion rate of each channel corresponding to the three recording data are integrated, and the channel thermal state parameter set is generated.

[0011] As a further scheme of the application, the acquisition step of the power suppression priority table is specifically: S411: Extract all radio frequency hyperthermia emission channels marked as to-be-switched state from the channel thermal state parameter set, calculate the difference value between the current period power request value and the last period output power value corresponding to each channel, and generate the to-be-switched state channel power request deviation sequence; S412: Based on the to-be-switched state channel power request deviation sequence, the power request deviation of each channel and the parameters of the temperature rise amplitude and the heat diffusion rate in the channel thermal state parameter set are combined, and the channel sorting input set is generated. S413: According to the power request deviation, the temperature rise amplitude and the heat diffusion rate corresponding to each channel in the channel ranking input set, a multi-parameter ranking is performed by a TOPSIS optimization algorithm, a position mapping is performed on the ranking result and the original channel number, and a power suppression priority table is generated by sequence arrangement.

[0012] As a further scheme of the present application, the channel output behavior control result obtaining step specifically comprises: S511: Extracting the power suppression priority table specified ranking radio frequency hyperthermia emission channel, configuring the limit parameter, setting the next cycle power emission upper limit and delayed activation time, generating the power limit and activation timing setting table; S512: Based on the power limit and activation timing setting table, combining the power input sequence and the temperature rise amplitude value in the channel thermal state parameter set, performing parameter weighted distribution on the remaining power resources, and setting the target power of each remaining channel according to the distribution weight, generating a next cycle power emission plan structure set; S513: According to the next cycle power emission plan structure set and the control parameters of the corresponding channel in the power limit and activation timing setting table, integrating the power setting value and the activation sequence of each channel, outputting the control instruction set in a predetermined format, and obtaining the channel output behavior control result.

[0013] An intelligent control system of a radio frequency endogenous field hyperthermia system, which is used to execute the intelligent control method of the radio frequency endogenous field hyperthermia system, and comprises: A power trend identification module, which collects the instantaneous power input of a specified period of a hyperthermia emission channel corresponding region under an endogenous field, judges the power trend change in adjacent periods, and obtains a power trend judgment result; A thermal response joint determination module, which collects the temperature rise amplitude and the heat diffusion rate of a target treatment tissue region of a current emission channel under an endogenous field, and performs joint criticality determination with the power trend judgment result in the corresponding period, and obtains a target treatment tissue region thermal response analysis result; A thermal state information labeling module, which labels the thermal state information of a radio frequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate according to the target treatment tissue region thermal response analysis result, collects the power request value of the target channel in the current period and the output power value in the last period, and generates a channel thermal state parameter set; A priority ranking module, which calculates the power request deviation between the power request value of the target channel in the current period and the output power value in the last period in the channel thermal state parameter set, performs suppression priority ranking between channels according to the power request deviation, and generates a power suppression priority table; An output behavior control module adjusts a next cycle power emission plan according to the power suppression priority table, and sets a radio frequency hyperthermia emission channel activation sequence to obtain a channel output behavior control result.

[0014] Compared with the prior art, the application has the advantages and positive effects that: In the application, under the endogenous field, the real-time judgment of the thermal response state of the target tissue can be realized by dynamically collecting the variation trend of the power input in different time periods and combining the linkage analysis of the target region temperature rise amplitude and the heat diffusion rate, and then the thermal state information of different emission channels in different periods is identified and labeled, and the deviation of the power request value and the historical output value is evaluated, the multi-factor sorting of the power request deviation between channels and the heat conduction characteristics is performed, the power allocation and the activation sequence are regulated and controlled by the suppression priority, the power control optimization under the dynamic change of the thermal response of different tissue regions is realized, the power configuration in the overall hyperthermia process is more consistent with the actual thermal effect distribution state, the resource waste and the local overheating situation are avoided, the heat diffusion efficiency and the stability of the tissue temperature rise are improved, the controllability and the targeting of the hyperthermia process are guaranteed, the response accuracy and the self-adaptive control ability of the system under the complex heat conduction scene of multiple channels are enhanced, the thermal action consistency and the conduction coordination of the targeted treatment tissue region are improved, and the risk of heat diffusion deviation and heat accumulation imbalance is effectively avoided. BRIEF DESCRIPTION OF DRAWINGS Figure 1 It is a work flow diagram of the application; Figure 2 It is a flow chart of step S1 of the application; Figure 3 It is a flow chart of step S2 of the application; Figure 4 It is a flow chart of step S3 of the application; Figure 5 It is a flow chart of step S4 of the application; Figure 6 It is a flow chart of step S5 of the application. DETAILED DESCRIPTION

[0015] In order to make the purpose, technical scheme and advantages of the application more clear, the application is further described in detail below in combination with the drawings and examples. It should be understood that the specific examples described herein are only used to explain the application and do not limit the application.

[0016] Please refer to Figure 1 The application provides a technical scheme: an intelligent control method of a radio frequency endogenous field hyperthermia system for tumors, comprising the following steps: S1: Under the endogenous field, collect the instantaneous power input of the specified period of the corresponding region of the hyperthermia emission channel, judge the power trend change in the adjacent period, and obtain the power trend judgment result; S2: Under the endogenous field, collect the temperature rise amplitude and heat diffusion rate of the targeted treatment tissue area of the current emission channel, and jointly determine the criticality with the power trend judgment result in the corresponding period to obtain the thermal response analysis result of the targeted treatment tissue area; S3: According to the thermal response analysis result of the targeted treatment tissue area, label the thermal state information of the radiofrequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate, collect the power request value of the target channel in the current period and the output power value in the last period, and generate a channel thermal parameter set; S4: Calculate the power request deviation between the power request value of the target channel in the current period and the output power value in the last period in the channel thermal parameter set, and perform suppression priority sorting between channels according to the power request deviation to generate a power suppression priority table; S5: According to the power suppression priority table, adjust the next period power emission plan, and set the radiofrequency hyperthermia emission channel activation sequence to obtain the channel output behavior control result; The power trend judgment result includes power input rate of change, adjacent period power difference, and power trend type. The targeted treatment tissue area thermal response analysis result includes temperature rise amplitude stability, heat diffusion rate characteristics, and joint criticality state. The channel thermal parameter set includes thermal state information label, current period power request value, and last period output power value. The power suppression priority table includes channel order, suppression priority level, and participating channel number. The channel output behavior control result includes next period power emission value, channel activation sequence, and control adjustment state.

[0017] Please refer to Figure 2 The acquisition step of the power trend judgment result is specifically: S111: Obtain the radiofrequency hyperthermia emission channel currently in the active state, collect the power signal of the targeted treatment tissue area acted on by the channel, obtain the instantaneous power input value in the continuous period within the specified sampling period, and generate a continuous period power input value sequence; To obtain the radio frequency hyperthermia transmitting channel currently in the active state, first, the real-time state of the multiple radio frequency channels needs to be read at the control end of the treatment device system. The activation state of each channel is identified by a flag bit on the control panel. For example, if the state bit of channel number 3 is "1", it indicates that it is in the active state. Then the system calls the power collection module corresponding to the channel to detect the output signal of the electrode connection area, collects the instantaneous power input value of the targeted treatment tissue area acted on by the channel, sets the sampling period to 1 second, and the system samples every 0.2 seconds within this period to obtain 5 groups of instantaneous power values. For example, the values obtained in a certain period are 13.2, 13.1, 13.4, 13.0 and 13.3 watts respectively. The system stores them as the first group of power input sequence, and generates the second group by re-sampling after the next second. According to this rule, the sampling value groups of multiple periods are continuously recorded to form a complete continuous period power input value sequence. In the application scenario, if the radio frequency hyperthermia system is used to treat liver tumor target area, the radio frequency energy applied through the electrode patch will penetrate the tissue, and the collection system will record the power value received by the penetrated area in real time. After each sampling, the system will automatically label the sampling time and associate it with the channel number and time point to form a continuous period power input value sequence.

[0018] S112: Based on each pair of adjacent periods in the continuous period power input value sequence, the power change relationship between each period segment is established in order, the progressive trend of the power input value of the next period and the corresponding value of the previous period is judged, and a sequence of adjacent period power change states is generated; If the average power of the first period is 13.2 watts and the second period is 12.8 watts, the system calculates the power change difference between the two periods as 0.4 watts, judges whether the difference is positive, negative or zero, and generates a state identifier accordingly. This process is performed in pairs throughout the sequence to record the change state between multiple periods. To ensure the accuracy of state identification, the system defines a judgment standard for "change". If the difference exceeds 0.3 watts, it is considered to be "significant change". This standard is derived from the analysis of the temperature change data of the tissue caused by power fluctuations above 0.3 watts in hyperthermia treatment. It is believed that this value can produce a significant response to the tissue. Therefore, when the difference is 0.4 watts, it is judged to be a significant decrease. The complete change trend sequence is obtained by processing all period pairs in turn. For example, the state sequence appears "decrease", "decrease", "unchanged", "increase" and so on. Each state value reflects the power relationship change between the corresponding two periods. In hyperthermia applications, this processing method is used to continuously monitor whether the received power change of the tissue is stable.

[0019] S113: Combine the trend segment that meets the continuous decrease condition in the adjacent period power change state sequence to identify the overall trend structure of the current power change, and obtain the power trend judgment result. To determine whether the current power is in an overall downward trend, a criterion needs to be set. A criterion is set that a continuous 3 cycles or more in a downward state is considered as a downward trend segment. This criterion is set with reference to the tissue temperature response delay time in medical literature. Generally, the tissue temperature rise needs to be continuously decreased for more than 3 seconds to show thermal reaction attenuation. Therefore, 3 cycles are taken as the minimum length of continuous decrease. The system scans the state sequence from the beginning in a sliding window manner to detect whether there is a sub-sequence of "decrease" for 3 or more times, for example, the state of "decrease" appears 4 times continuously from the 4th cycle. This segment is identified as a trend segment, the start and end positions are recorded and the length is accumulated. Continue to scan backward until the end of the sequence. Finally, the total length of all downward trend segments and the proportion of the total length of the sequence are counted. If the total downward trend segment length exceeds 40% of the total sequence, it is considered that the overall is in a downward trend. The 40% criterion is based on the experience of the total proportion of the actual effective downward phase in the simulation treatment process. If a total of 20 cycles are collected in a sampling, 8 cycles of downward trend segments are accumulated, the proportion is 40%, which meets the determination condition. The final output power trend judgment result is downward.

[0020] Please refer to Figure 3 The acquisition step of the targeted treatment tissue area thermal response analysis result is specifically: S211: Based on the power trend judgment result of continuous decrease, the radio frequency hyperthermia emission channel collects the temperature value in the current cycle in the target channel corresponding to the targeted treatment tissue area, extracts the temperature change data of each record point, combines the temperature spatial distribution gradient of the adjacent area, calculates the temperature rise amplitude and heat diffusion rate in a specified time, and obtains the thermal conduction characteristic data of the targeted treatment tissue area. Based on the power trend judgment result for continuous decline of radio frequency hyperthermia emission channel, the system needs to collect temperature data in the target therapy tissue region corresponding to the channel. First, the thermocouple sensor embedded in the tissue or the surface infrared temperature control probe is used for real-time monitoring in the current period. Under the setting of 1 second sampling period, the temperature value is obtained once per second, and the value is recorded as 39.5 degrees Celsius. The temperature value in the last period is 38.9 degrees Celsius. The calculation method of temperature rise amplitude is current temperature minus last period temperature, i.e. temperature rise amplitude = current period temperature - last period temperature. For example: 39.5 degrees Celsius - 38.9 degrees Celsius = 0.6 degrees Celsius. This value is the current period temperature rise amplitude. Then, the temperature difference between adjacent detection points (for example, arranged at an interval of 3 mm) is analyzed. For example, the center point is 39.5 degrees Celsius, and the outer point is 38.3 degrees Celsius. The calculation method of temperature spatial distribution gradient is center point temperature minus edge point temperature, i.e. temperature gradient = center point temperature - edge point temperature. For example: 39.5 degrees Celsius - 38.3 degrees Celsius = 1.2 degrees Celsius. The calculation method of heat diffusion rate is temperature gradient divided by detection point spacing, i.e. heat diffusion rate = temperature gradient ÷ detection point spacing. For example: 1.2 degrees Celsius ÷ 3 mm = 0.4 degrees Celsius per millimeter. The system records the temperature rise amplitude and heat diffusion rate of the period, establishes a time sequence, compares the temperature rise and diffusion rate in 5 consecutive periods, and records them into the database to construct a time-temperature response characteristic curve and form the heat conduction characteristic data of the target therapy tissue region.

[0021] S212: Based on the heat conduction characteristic data of the target therapy tissue region and the power trend judgment result, the synchronization between the power input trend and the temperature diffusion behavior is identified, and the current period data combination is selected under the consistent trend condition. The period sample meeting the characteristics of stable temperature rise and diffusion enhancement is extracted, and the heat diffusion dominant trend interval is generated. Correlation analysis is performed on the targeted treatment tissue region heat conduction characteristic data and the corresponding power trend judgment results. The goal is to identify whether the power decline trend and the temperature change behavior are synchronized. If the temperature rise amplitude does not decrease or slightly increases and the heat diffusion rate remains relatively stable during continuous power decline, it is determined that the two are synchronized. The synchronization standard is that the temperature rise amplitude fluctuation range is not more than 0.2 degrees Celsius and the diffusion rate variation is not more than 0.1 degrees Celsius per millimeter. For example, in the case of continuous power decline for three periods, the temperature rise amplitude is recorded as 0.5, 0.6, and 0.5 degrees Celsius, and the diffusion rate is 0.35, 0.36, and 0.34 degrees Celsius per millimeter. The fluctuations are all less than the standard judgment range, so it is considered to have a synchronized trend. The system further screens such periodic samples and extracts a sample combination that meets the stable temperature rise and diffusion enhancement conditions in the current period data as a representative, generating a data section called "heat expansion dominant trend interval". In this section, the power decreases but the temperature continues to stabilize or the diffusion speed slightly increases. Such sample segments are marked by the system as typical heat expansion dominant reaction regions.

[0022] S213: According to the heat expansion dominant trend interval, it is determined whether the targeted treatment tissue region is in a heat expansion dominant state, and the change amplitude and distribution characteristics of each trend data item in the corresponding state are combined to determine whether the target tissue region exhibits stable temperature rise and diffusion coordination conduction characteristics, and the targeted treatment tissue region heat response analysis result is obtained; According to the heat expansion dominant trend interval, it is further analyzed and judged whether the targeted treatment tissue region is in a heat expansion dominant state. The judgment condition is set as follows: the temperature change amplitude is stable within 0.5 degrees Celsius for more than 3 consecutive periods in this interval, and the heat diffusion rate remains above 0.3 degrees Celsius per millimeter with a fluctuation range within 0.1 degrees Celsius. If this condition is met, it is determined that the tissue region is in a heat expansion dominant state. This standard is based on the time window of heat expansion behavior in the literature for temperature control response in the tumor edge region. By quantitatively processing the change amplitude of each trend data item and jointly analyzing the spatial distribution characteristics, for example, if the temperature rise is 0.6, 0.5, 0.5, and 0.4 degrees Celsius for 4 consecutive periods, and the heat diffusion rate is 0.32, 0.34, 0.31, and 0.30 degrees Celsius per millimeter, the fluctuation range is within the allowed threshold. If it meets the heat expansion dominant state judgment, the number of power decline periods in this section and the heat expansion rate growth trend are combined to confirm that the target tissue region exhibits stable temperature rise and diffusion coordination heat conduction behavior in this stage, and the targeted treatment tissue region heat response analysis result is output.

[0023] Please refer to Figure 4 The acquisition steps of the channel thermal parameter set are as follows: S311: According to the targeted treatment tissue region heat response analysis result, the radio frequency hyperthermia transmission channel with stable state is marked as the state to be switched, and the temperature rise amplitude and heat diffusion rate stable channel are obtained. According to the results of the targeted tissue area thermal response analysis, the system selects those channels that meet the "state stable" determination condition in the radio frequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate, and marks them as the state to be switched. The determination criteria include that the temperature change amplitude is between 0.4 and 0.6 degrees Celsius for more than three consecutive periods, the diffusion rate is maintained between 0.3 and 0.4 degrees Celsius per millimeter, and the fluctuation amplitude is not more than 0.05 degrees Celsius per millimeter. Such criteria are set with reference to the tumor area heat conduction simulation experiment. When the channel continuously meets the above conditions, the system performs the marking operation, writes the channel ID and its corresponding current state into the "switching state list" in the control module, for example, channel No. 4, whose temperature rise amplitude is 0.5, 0.6, and 0.5 degrees Celsius in the last three periods, the diffusion rate is 0.32, 0.33, and 0.34 degrees Celsius per millimeter, and the fluctuation amplitude is 0.02, all of which meet the threshold requirements. The system marks it as a "stable channel" and enters the state to be switched. Such operations will scan the temperature rise and diffusion data of all active channels one by one, judge whether they meet the temperature stability condition, and perform the marking operation. Finally, the system generates a "temperature rise amplitude and heat diffusion rate stable channel" list based on all the stable state channels marked.

[0024] S312: Based on the temperature rise amplitude and heat diffusion rate stable channel, the current period power request value and the last period output power value of each channel are collected, and the two power data are associated with the corresponding channel to generate a radio frequency hyperthermia emission channel power request and output value group; Based on the temperature rise amplitude and heat diffusion rate stable channel, the current period power request value and the last period output power value of each channel are collected, and the two power data are associated with the corresponding channel to generate a radio frequency hyperthermia emission channel power request and output value group;

[0025] S313: According to the radio frequency hyperthermia emission channel power request and output value group, integrate each channel corresponding power input sequence, temperature rise amplitude and heat diffusion rate three record data, generate channel thermal parameter set; According to the radio frequency hyperthermia emission channel power request and output value group, continue to integrate the other two types of key data in the current channel, that is, the power input sequence, the temperature rise amplitude sequence and the heat diffusion rate sequence. The collection frequency of the three data needs to be consistent, such as recording once every second. The system sets the window length to five cycles. Five groups of continuous data are combined and packed to form a complete thermal parameter sequence. The system first retrieves the channel number where the current cycle is located, and extracts the power input value, corresponding temperature rise amplitude and heat diffusion rate three-dimensional data of the channel in the last 5 seconds from the database in time sequence. Each type of data needs to be confirmed to be consistent with the time tag of the power request and output value group. The system performs binding operation after synchronous check of the data, for example, the power recorded by channel No. 4 in the 10th to 14th second is 13.8, 14.1, 13.9, 14.0 and 14.2 watts, the temperature rise amplitude is 0.5, 0.4, 0.6, 0.5 and 0.5 degrees Celsius, and the heat diffusion rate is 0.34, 0.36, 0.35, 0.33 and 0.32 degrees Celsius per millimeter. The three data correspond in unit time. The system writes them into the thermal parameter structure to form the thermal parameter set of the channel in this period of time.

[0026] Please refer to Figure 5 , the acquisition steps of the power suppression priority table are as follows: S411: Extract all radio frequency hyperthermia emission channels marked as switching state from the channel thermal parameter set, calculate the difference between the current cycle power request value and the last cycle output power value of each channel, and generate a switching state channel power request deviation sequence; Extract all radio frequency hyperthermia emission channels marked as switching state from the channel thermal parameter set. The system first calls the switching channel identification list, matches the channel numbers in the parameter set one by one, selects the corresponding power data field for analysis, extracts the power request value of the current cycle and the output power value of the last cycle for each selected channel, and then calculates the difference between the two values. The calculation formula is described in words as follows: power request deviation = current cycle power request value - last cycle output power value. For example, the request value of channel No. 6 is 15.0 watts, and the output value is 14.2 watts. The power request deviation is 0.8 watts. Similarly, the same calculation operation is performed on other channels. The system records the numbers of all channels and the corresponding power deviations in a structured manner to form a switching state channel power request deviation sequence. In actual application, all channel power deviation values will be used to judge the adjustment pressure degree. The larger the value, the stronger the current power adjustment demand of the channel.

[0027] S412: Based on the channel power request deviation sequence to be switched, combine the power request deviation of each channel with the temperature rise amplitude and heat diffusion rate parameters in the channel thermal state parameter set to generate a channel ranking input set; Based on the channel power request deviation sequence to be switched, the system continues to extract the temperature rise amplitude and heat diffusion rate of the corresponding channel from the channel thermal state parameter set, and combines the three data to form a complete ranking input item. Each channel forms an input group containing three indicators, and the combination method is: power request deviation + temperature rise amplitude + heat diffusion rate, stored in a structure array. During the combination process, the system needs to verify whether the timestamps of the three data correspond to ensure that they come from the same cycle data source. For example, the request deviation of channel number 6 is 0.8 watts, the temperature rise amplitude is 0.6 degrees Celsius, and the diffusion rate is 0.35 degrees Celsius per millimeter. The channel ranking input group is a combination structure composed of three numerical values. After all channels complete data combination according to the same rule, the system forms a channel ranking input set.

[0028] S413: According to the power request deviation, temperature rise amplitude and heat diffusion rate corresponding to each channel in the channel ranking input set, perform multi-parameter sorting through the TOPSIS optimization algorithm, and map the sorting results to the original channel number to generate a power suppression priority table in sequence; According to the power request deviation, temperature rise amplitude and heat diffusion rate corresponding to each channel in the channel ranking input set, the system uses the TOPSIS optimization algorithm with a preference adjustment factor to perform multi-parameter sorting. The essence of the TOPSIS method is to construct a positive ideal solution (i.e. the combination of the best values of all evaluation indicators) and a negative ideal solution (i.e. the combination of the worst values of all evaluation indicators), then calculate the distance between each channel and the two extreme solutions, and determine the ranking score of each channel through the distance ratio.

[0029] The evaluation indicators include: power request deviation (unit: watts): indicates the difference between the current required power of the channel and the actual output of the previous period, the larger the value, the higher the power adjustment pressure of the channel, and therefore the larger the value, the better the indicator; temperature rise amplitude (unit: degrees Celsius): reflects the temperature rise degree of the organization per unit time, the faster the temperature rises, the stronger the heat efficiency, and the larger the value, the better the indicator; heat diffusion rate (unit: degrees Celsius / mm): indicates the heat diffusion speed to the surrounding tissue, the higher the value, the more active the thermal response of the tissue, and the larger the value, the better the indicator. Therefore, the three indicators are all "the larger the better" type, and the positive ideal solution constructed in TOPSIS is the combination of the maximum values of each indicator, while the negative ideal solution ​is the minimum combination of all indicators. Let the total number of channels be n, and the total number of indicators be m = 3. For the jth indicator, the positive ideal solution takes the maximum value of the indicator in all channels, and the negative ideal solution takes the minimum value.

[0030] To deal with the problem of inconsistent dimensions between indicators, the system first normalizes each indicator data, unifies the dimension, and the specific method is: divide the value of each channel on a certain indicator by the square root of the sum of squares of all channels on the indicator. Finally, each column of data is mapped to 0 to 1, while keeping the relative size relationship between channels unchanged.

[0031] Subsequently, the weight term is introduced Importance level is given to each indicator. The weight setting principle is as follows: power request deviation weight : The greater the adjustment demand, the more the channel needs to be suppressed in priority, so the maximum weight is given; temperature rise amplitude weight : Fast temperature rise indicates that the channel responds strongly to heat treatment, which is less than power demand; heat diffusion rate weight : As an auxiliary evaluation item, the weight is the lowest but still meaningful.

[0032] The system also introduces a preference adjustment factor , which adjusts the sorting result to be more biased towards the ideal solution or away from the negative ideal solution. In this scenario, it is set to , indicating a bias towards the positive ideal solution. Its significance lies in: when the channel performs well, the system is more inclined to promote its priority; when the channel performs poorly, the sorting impact is small. Finally, the TOPSIS relative closeness score calculation formula is as follows: ; Where, : The sorting score of the ith channel, the value range is between 0 and 1, the higher the priority, the earlier the priority; : Preference adjustment factor, set to 0.6; : The weight of the jth indicator, , , ; : The value of the ith channel on the jth indicator; , : The maximum and minimum values of the jth indicator of all channels; m: total number of indicators, here it is 3.

[0033] Let the first channel indicator value be: power request deviation: 1.2 watts, temperature rise amplitude: 0.5 degrees Celsius, heat diffusion rate: 0.35 degrees Celsius / mm.

[0034] Assume that in all channels, the maximum and minimum values of the three indicators are: Power request deviation maximum , minimum value , maximum value of temperature rise amplitude , minimum value , maximum value of heat diffusion rate , minimum value . The index weight is set as follows: (power request deviation), (temperature rise amplitude), (heat diffusion rate). The preference adjustment factor is set as: .

[0035] Forward distance numerator item calculation: power request deviation item: , , ; temperature rise amplitude item: , , ; heat diffusion rate item: , , ; forward distance square sum is: , .

[0036] negative distance denominator item calculation: power request deviation item: , , ; negative distance square sum is: , .

[0037] relative closeness calculation: , the ranking score of the channel is , the value will be compared in all channel rankings, and decide its final ranking position in the power suppression priority table.

[0038] Firstly, the three indicators of power request deviation, temperature rise amplitude and heat diffusion rate of each channel are standardized to ensure their comparability in the same dimension. Then, the importance of different indicators is reflected by setting weights, among which the power request deviation is given the highest priority, representing the urgency of energy adjustment, and the temperature rise amplitude and heat diffusion rate reflect the thermal response level of the target tissue. On this basis, the system calculates the closeness (i.e. positive distance) between each channel and the ideal optimal indicator combination and the distance (i.e. negative distance) between each channel and the worst indicator combination. These two distances measure the positive and negative levels of the current channel performance respectively. Finally, the system calculates a relative closeness score by ratio fusion of the two distances through a comprehensive score formula. The score is between zero and one, and the larger the value, the closer the channel is to the best state and the higher the ranking. The core of the whole process is to express multi-dimensional evaluation data in a quantitative way. Through the four steps of standardization, weighting, distance calculation and ratio conversion, a complete, executable and logically clear channel priority ranking model is constructed to support objective selection in power suppression scheduling.

[0039] Please refer to Figure 6 , the acquisition step of the channel output behavior control result is: S511: Extract the power suppression priority table specified sorting radio frequency hyperthermia emission channel, configure the limit parameter, and set the next cycle power emission upper limit and delayed activation time, generate the power limit and activation timing setting table; Extract the radio frequency hyperthermia emission channel determined according to the sorting result in the power suppression priority table. The system first reads the channel number and its corresponding priority level value in the sorting list, selects the channel with high ranking as the priority restriction object, configures the limit parameter for these channels, including the power emission upper limit value and the activation delay time. The setting of limit parameter needs to be based on the current thermal state response characteristics and historical fluctuation of the channel, for example, if the power request deviation of a channel is higher than 1.0 watt in the last five cycles, and the temperature rise speed is above 0.6 degrees Celsius per second, the system limits the power emission upper limit of the channel to the current average output power minus 1.5 watts, and sets the activation delay time to 2 seconds, which is based on the thermal stability recovery time set in the tissue cooling response model. If the current average power of the channel is 14.2 watts, the power emission upper limit is set to 12.7 watts. All parameters and channel numbers are written into the control table item and indexed according to the channel number. Finally, the system generates the power limit and activation timing setting table, which contains the channel number, limit power value, delay time and limit state identifier field.

[0040] S512: Based on the power limit and activation timing setting table, combined with the power input sequence and temperature rise amplitude value in the channel thermal parameter set, the remaining power resources are allocated with parameter weighting, and the target power of each remaining channel is set according to the allocation weight to generate the power transmission plan structure set for the next cycle. Based on the power limiting and activation timing settings table, before entering the next power allocation cycle, the system divides all unrestricted radio frequency hyperthermia transmission channels into two categories: thermal response driven channels (based on temperature rise amplitude) and load-bearing channels (based on historical power input), each using an independent power allocation strategy. The system proportionally splits the remaining allocable total power resource R into two parts, for example, into temperature rise channel resources. Wattage and load-type channel resources Tiles, and then proportionally distribute them respectively.

[0041] Category 1: Power allocation formula for temperature rise driven channels: ; in, Assigned to the first The target power for each temperature-driven channel is expressed in watts. The total power allocated to temperature rise channels, in watts. : No. The temperature rise of each channel in the most recent period, in degrees Celsius. The sum of the temperature rise amplitudes of all temperature rise channels. The total number of temperature rise channels. : No. The temperature rise amplitude of each temperature-driven channel, in degrees Celsius. : The index of the current temperature rise driven channel. : The first general index variable used when traversing the channel.

[0042] Let channels A and B belong to the temperature rise driven type, with temperature rise amplitudes of respectively. , The temperature ranges from 0.9 degrees Celsius to 0.9 degrees Celsius. This is a resource related to temperature rise. Tile, then: Channel A: watt; Channel B: watt.

[0043] Category 2: Power allocation formula for load-bearing channels: ; in, : Assigned to the The target power of each load-bearing channel is expressed in watts. : Total power of all load type channels, unit: watt, : Historical average power input value of the : Historical average power input value of the : Sum of average power value of all load type channels, : Total number of load type channels, : Average power input value of the : Average power input value of the : Index of current load type channel, : Second general index variable used when traversing channels.

[0044] Suppose channels C and D belong to load carrying type, with average power input of , 24 watts, and load type resource 24 watts, then: Channel C: 24 watts; Channel D: 24 watts.

[0045] The system writes the target power value of all channels , into the next cycle power emission plan structure set, and binds it with the channel number to form a data structure, ensuring that the downstream control logic can identify and execute scheduling commands according to channel type. This strategy follows the principle of separation of different indicators and independent allocation of different resources, avoiding logical errors caused by direct combination of different unit data, while also considering the dual regulation goals of organizing temperature control response and device power stability.

[0046] First, the unrestricted channels are divided into two categories according to their response characteristics: one is the temperature rise driven channel, which is used to identify those tissue regions that have a faster temperature rise in the last cycle; the other is the load bearing channel, which is used to identify those output paths that have stable historical power input and strong energy carrying capacity. The system reserves different proportions of power resources for these two types of channels, and uses two completely independent proportional allocation formulas to calculate the target power respectively. The first formula is used for temperature rise channels. It compares the proportion of the temperature rise amplitude of each channel to the total temperature rise amplitude to determine the proportion of the corresponding channel in the temperature rise driven resource. Finally, the target power is calculated. The larger the value, the faster the corresponding tissue region of the channel rises in temperature, and it will receive more energy supply. The second formula is used for load channels. It compares the proportion of the historical average power of each channel to the total input to determine its allocation proportion in the load resource. Finally, the target power is determined. The larger the value, the higher the long-term power input of the channel, and the more capable it is of bearing the power task of the new cycle. The two parts of the results are finally written into the transmission plan data structure for the output control of the channels in the next cycle. The overall process strictly follows the physical quantity classification calculation principle, avoiding mixed units that can cause calculation errors, while ensuring that the power resources are allocated to the most needed locations.

[0047] S513: According to the control parameters of the corresponding channels in the next cycle power transmission plan structure set and the power limit and activation timing setting table, the power setting value and activation sequence of each channel are integrated, and the control instruction set is output in a predetermined format to obtain the channel output behavior control result; According to the control parameters of each channel in the next cycle power transmission plan structure set and the power limit and activation timing setting table, the system indexes by channel number, cross-compares and integrates the target power setting value of each channel with the power upper limit value and activation delay time in the limit parameter. First, it is determined whether the current set target power of the channel exceeds the set upper limit. If it exceeds, the set value is modified to the power upper limit value of the corresponding channel. For example, the target power of channel number 5 is 15.5 watts, but the limit upper limit is 14 watts, then the system automatically modifies the channel power setting value to 14 watts. Then, combined with the activation delay time field configured for the channel, all channels to be activated are arranged in ascending order of activation time to form an activation sequence queue. For example, channel 3 has a delay of 1 second, channel 5 has a delay of 3 seconds, and channel 7 has a delay of 2 seconds. The order is channel 3, 7, 5. The system combines the final channel number, the modified target power value, the delay time value and the activation sequence position information to construct a structured instruction format, generating a series of data instruction lines that conform to the output control protocol. Finally, all channel instructions are combined to form a control instruction set, which is used as the direct output entry of the system controller to drive each radio frequency hyperthermia transmission channel to execute the output behavior according to the specified power and timing in the next cycle, generating a complete channel output behavior control result.

[0048] The application discloses an intelligent control system of a radio frequency intratumoral endogenous field hyperthermia system. The power trend identification module collects the instantaneous power input of a specified period of a hyperthermia emission channel corresponding region under the endogenous field, judges the power trend change in adjacent periods, and obtains a power trend judgment result. The thermal response joint determination module collects the temperature rise amplitude and thermal diffusion rate of a targeted treatment tissue region of a current emission channel under the endogenous field, and jointly determines the criticality with the power trend judgment result in the corresponding period to obtain a thermal response analysis result of the targeted treatment tissue region. The thermal state information labeling module labels the thermal state information of the radio frequency hyperthermia emission channel with stable temperature rise amplitude and thermal diffusion rate according to the thermal response analysis result of the targeted treatment tissue region, collects the power request value of the target channel in the current period and the output power value in the last period, and generates a channel thermal state parameter set. The priority sorting module calculates the power request deviation between the power request value of the target channel in the current period and the output power value in the last period in the channel thermal state parameter set, performs suppression priority sorting between channels according to the power request deviation, generates a power suppression priority table, and adjusts the next period power emission plan according to the power suppression priority table. The output behavior control module adjusts the next period power emission plan according to the power suppression priority table, sets the activation sequence of the radio frequency hyperthermia emission channel, and obtains a channel output behavior control result.

[0049] The above is only a preferred embodiment of the application, and does not limit the application in other forms. Any person skilled in the art can modify or change the above disclosed technical content to equivalent embodiments applied to other fields, but any simple modification, equivalent change and modification made according to the technical essence of the application to the above embodiments still belongs to the protection scope of the application technical scheme.

Claims

1. A method for intelligent control of a radiofrequency intratumoral endogenous field hyperthermia system, characterized in that, The method comprises the following steps: S1: under the endogenous field, collecting the instantaneous power input of the specified period of the corresponding region of the radiofrequency hyperthermia emission channel, judging the power trend change in the adjacent period, and obtaining the power trend judgment result; S2: under the endogenous field, collecting the temperature rise amplitude and heat diffusion rate of the targeted treatment tissue area of the current emission channel, and jointly judging the criticality with the power trend judgment result in the corresponding period, to obtain the thermal response analysis result of the targeted treatment tissue area; S3: according to the thermal response analysis result of the targeted treatment tissue area, labeling the thermal state information of the radiofrequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate, collecting the power request value of the target channel in the current period and the output power value in the last period, and generating a channel thermal state parameter set; S4: calculating the power request deviation between the power request value of the target channel in the current period and the output power value in the last period in the channel thermal state parameter set, performing suppression priority sorting between channels according to the power request deviation, and generating a power suppression priority table; S5: according to the power suppression priority table, adjusting the next period power emission plan, and setting the radiofrequency hyperthermia emission channel activation sequence, to obtain the channel output behavior control result.

2. The intelligent control method of the radiofrequency intratumoral field thermotherapy system according to claim 1, characterized in that, The power trend judgment result includes power input change rate, adjacent period power difference, and power trend type. The thermal response analysis result of the targeted treatment tissue area includes temperature rise amplitude stability, heat diffusion rate characteristics, and joint criticality state. The channel thermal state parameter set includes thermal state information label, current period power request value, and last period output power value. The power suppression priority table includes channel order sequence, suppression priority level, and participating sorting channel number. The channel output behavior control result includes next period power emission value, channel activation sequence, and control adjustment state. 3.The intelligent control method of the radiofrequency intratumoral field thermotherapy system according to claim 1, characterized in that, The acquisition step of the power trend judgment result is specifically: S111: acquiring the radiofrequency hyperthermia emission channel currently in the activated state, collecting the power signal of the targeted treatment tissue area acted on by the channel under the endogenous field, acquiring the instantaneous power input value in the continuous period within the specified sampling period, and generating a continuous period power input value sequence; S112: based on the power input pairs of every two adjacent periods in the continuous period power input value sequence, establishing the power change relationship between each period segment in order, judging the progressive trend of the power input value of the next period and the corresponding value of the previous period, and generating an adjacent period power change state sequence; S113: combining whether there is a trend segment meeting the continuous descending condition in the adjacent period power change state sequence, identifying the overall trend structure of the current power change, and obtaining the power trend judgment result.

4. The intelligent control method of the radiofrequency intratumoral endogenous field hyperthermia system according to claim 3, characterized in that, The acquisition step of the thermal response analysis result of the targeted treatment tissue area is specifically: S211: based on the radiofrequency hyperthermia emission channel with continuous descending power trend judgment result, collecting the temperature value in the current period in the targeted treatment tissue area corresponding to the target channel, and extracting the temperature change data of each recording point, combining the temperature spatial distribution gradient of the adjacent region, calculating the temperature rise amplitude and heat diffusion rate in the specified time, and obtaining the thermal conduction characteristic data of the targeted treatment tissue area; S212: Based on the target therapy tissue area heat conduction characteristic data and the power trend judgment result, the synchronization between the power input trend and the temperature diffusion behavior is identified, and under the trend consistency condition, the current period data combination is screened, the period sample conforming to the temperature rise stable and diffusion enhancement characteristics is extracted, and the heat diffusion dominant trend interval is generated; S213: According to the heat diffusion dominant trend interval, it is judged whether the target therapy tissue area is in the heat diffusion dominant state, and the change amplitude and distribution characteristics of each trend data item in the corresponding state are combined to judge, whether the target tissue area shows stable temperature rise and diffusion coordinated conduction characteristics, and the target therapy tissue area heat response analysis result is obtained.

5. The intelligent control method of radiofrequency intratumoral endogenous field hyperthermia system according to claim 4, characterized in that, The acquisition step of the channel thermal state parameter set is specifically: S311: According to the target therapy tissue area heat response analysis result, the radio frequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate is marked as a state to be switched, and a channel with stable temperature rise amplitude and heat diffusion rate is obtained; S312: Based on the channel with stable temperature rise amplitude and heat diffusion rate, the power request value in the current period and the output power value in the last period are collected for each channel, and the two power data are associated with the corresponding channel respectively, and a radio frequency hyperthermia emission channel power request and output value group is generated; S313: According to the radio frequency hyperthermia emission channel power request and output value group, the power input sequence, temperature rise amplitude and heat diffusion rate of each channel are integrated, and a channel thermal state parameter set is generated.

6. The intelligent control method of radiofrequency intratumoral endogenous field hyperthermia system according to claim 5, characterized in that, The acquisition step of the power suppression priority table is specifically: S411: Extract all radio frequency hyperthermia emission channels marked as a state to be switched from the channel thermal state parameter set, calculate the difference between the current period power request value and the last period output power value corresponding to each channel, and generate a channel power request deviation sequence in a state to be switched; S412: Based on the channel power request deviation sequence in a state to be switched, the power request deviation of each channel is combined with the temperature rise amplitude and heat diffusion rate parameters in the channel thermal state parameter set to generate a channel sorting input set; S413: According to the power request deviation, temperature rise amplitude and heat diffusion rate corresponding to each channel in the channel sorting input set, multi-parameter sorting is performed through TOPSIS optimization algorithm, the sorting result is mapped with the original channel number, and the power suppression priority table is generated in sequence.

7. The intelligent control method of radiofrequency intratumoral endogenous field hyperthermia system according to claim 6, characterized in that, The acquisition step of the channel output behavior control result is specifically: S511: Extract the radio frequency hyperthermia emission channel specified in the power suppression priority table, configure the limit parameter, set the next period power emission upper limit and delayed activation time, and generate a power limit and activation time sequence setting table; S512: Based on the power limit and activation time sequence setting table, combining the power input sequence and temperature rise amplitude value in the channel thermal state parameter set, the remaining power resources are executed parameter weighted distribution, and the target power of each remaining channel is set according to the distribution weight, and a next period power emission plan structure set is generated; S513: According to the next cycle power emission plan structure set and the corresponding channel control parameters in the power limit and activation timing setting table, the power setting value and the activation sequence of each channel are integrated, the control instruction set is output in a predetermined format, and the channel output behavior control result is obtained.

8. An intelligent control system for a radiofrequency tumor endogenous hyperthermia system, characterized in that, The intelligent control method of the radio frequency intratumor endogenous field hyperthermia system according to any one of claims 1-7, wherein the system comprises: A power trend identification module, under the endogenous field, collects the instantaneous power input of the specified period of the region corresponding to the hyperthermia emission channel, judges the power trend change in the adjacent period, and obtains the power trend judgment result; A thermal response joint determination module, under the endogenous field, collects the temperature rise amplitude and heat diffusion rate of the targeted treatment tissue region of the current emission channel, and performs joint criticality determination with the power trend judgment result in the corresponding period to obtain the thermal response analysis result of the targeted treatment tissue region; A thermal state information labeling module labels the thermal state information of the radio frequency hyperthermia emission channel with stable temperature rise amplitude and heat diffusion rate according to the thermal response analysis result of the targeted treatment tissue region, collects the power request value of the target channel in the current period and the output power value in the last period, and generates a channel thermal state parameter set; A priority sorting module calculates the power request deviation between the power request value of the target channel in the current period and the output power value in the last period in the channel thermal state parameter set, performs suppression priority sorting between channels according to the power request deviation, and generates a power suppression priority table; An output behavior control module adjusts the next cycle power emission plan and sets the activation sequence of the radio frequency hyperthermia emission channel according to the power suppression priority table, and obtains the channel output behavior control result.

Images