CN116227277A - Method, apparatus, device, medium and program product for generating electrode layout scheme - Google Patents

Method, apparatus, device, medium and program product for generating electrode layout scheme Download PDF

Info

Publication number
CN116227277A
CN116227277A CN202310085786.0A CN202310085786A CN116227277A CN 116227277 A CN116227277 A CN 116227277A CN 202310085786 A CN202310085786 A CN 202310085786A CN 116227277 A CN116227277 A CN 116227277A
Authority
CN
China
Prior art keywords
electrode
electrode point
field intensity
pair
point position
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.)
Pending
Application number
CN202310085786.0A
Other languages
Chinese (zh)
Inventor
林汲
杨鑫
胡春华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202310085786.0A priority Critical patent/CN116227277A/en
Publication of CN116227277A publication Critical patent/CN116227277A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00—Electrotherapy; Circuits therefor
    • A61N1/02—Details
    • A61N1/04—Electrodes
    • A61N1/0404—Electrodes for external use
    • A61N1/0408—Use-related aspects
    • A61N1/0456—Specially adapted for transcutaneous electrical nerve stimulation [TENS]
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Epidemiology (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Electrotherapy Devices (AREA)

Abstract

The present application relates to a method, apparatus, device, medium and program product for generating an electrode layout scheme. The method comprises the following steps: obtaining a standard finite element model of a target part; the standard finite element model comprises a region of interest; then, according to a plurality of preset electrode points on a standard finite element model, acquiring the average field intensity of the region of interest when the electric stimulation is applied to different electrode point pairs; and finally, generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point pair. Each electrode point position pair is any two electrode point positions in a plurality of preset electrode point positions. By adopting the method, an electrode layout scheme which can optimize the field intensity distribution of a tumor area can be obtained.

Description

Method, apparatus, device, medium and program product for generating electrode layout scheme
Technical Field
The present invention relates to the field of computer technologies, and in particular, to a method, an apparatus, a device, a medium, and a program product for generating an electrode layout scheme.
Background
Electric field therapy is a therapy performed by portable, non-invasive medical devices.
For example, in tumor electric field therapy, in general, when tumor electric field therapy is performed, a plurality of electrode sheets are required to be attached to the brain skin to continuously provide a low-intensity, medium-frequency alternating electric field to a tumor area. While the electric field is distributed in a non-uniform manner in the brain, different array electrode layouts can result in large field strength differences in the tumor area.
Therefore, how to determine an electrode layout scheme for optimizing the field intensity distribution of a tumor region is a technical problem to be solved.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, apparatus, device, medium, and program product for generating an electrode layout scheme that can optimize the field strength distribution of a tumor region.
In a first aspect, there is provided a method of generating an electrode layout scheme, the method comprising:
obtaining a standard finite element model of a target part; the standard finite element model comprises a region of interest;
acquiring average field intensity of the region of interest when electric stimulation is applied at different electrode point pairs according to a plurality of preset electrode points on the standard finite element model; each electrode point position pair is any two electrode point positions in a plurality of preset electrode point positions;
And generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point pair.
In one embodiment, obtaining a standard finite element model of a target site includes:
acquiring medical image data of a target part in at least one mode;
performing tissue segmentation on each piece of medical image data to obtain tissue segmentation information of a target part;
generating a three-dimensional model of the target part according to the tissue segmentation information of the target part;
and filling the tetrahedral mesh into the three-dimensional model of the target part to obtain a standard finite element model of the target part.
In one embodiment, obtaining the average field strength of the region of interest when the electrical stimulus is applied at different pairs of electrode points from a plurality of preset electrode points on a standard finite element model comprises:
for any electrode point position pair, acquiring a field intensity vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest when electric stimulation is applied to the electrode point position pair;
and determining the average field intensity of the region of interest when the electric stimulus is applied to the electrode point position pair according to the field intensity vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest.
In one embodiment, determining the average field strength of the region of interest when the electrical stimulus is applied at the electrode point pair according to the field strength vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest comprises:
determining the ratio of each grid cell in the total volume of all grid cells as the volume weight of each grid cell;
and determining the average field intensity of the region of interest when the electric stimulation is applied to the electrode point position pair according to the field intensity vector distribution value of each grid cell and the volume weight of each grid cell.
In one embodiment, the generating an optimal electrode layout scheme of the target site according to the average field intensity of each electrode point pair comprises:
determining a first electrode layout scheme according to the average field intensity of each electrode point position pair;
removing electrode points corresponding to the first electrode layout scheme from a plurality of preset electrode points on the standard finite element model to obtain residual preset electrode points on the standard finite element model;
determining a second electrode layout scheme according to the remaining preset electrode points;
and determining the first electrode layout scheme and the second electrode layout scheme as the optimal electrode layout scheme of the target part.
In one embodiment, determining the first electrode layout scheme according to the average field intensity of each electrode point pair comprises:
sequencing the average field intensity of each electrode point position pair, and acquiring a target electrode point position pair with the maximum average field intensity according to a sequencing result;
if the average field intensity of the target electrode point pair is greater than a preset field intensity threshold value, determining the electrode point position and the electric field direction of the target electrode point pair as a first electrode layout scheme;
if the average field intensity of the target electrode point pair is smaller than the field intensity threshold value, a first electrode layout scheme is determined by performing field intensity accumulation operation on the average field intensity of each electrode point pair.
In one embodiment, determining the first electrode layout scheme by performing a field strength summation operation on the average field strengths of the electrode point pairs includes:
acquiring a plurality of candidate electrode point pair sets according to the sequence of the average field intensity of each electrode point pair from large to small; each candidate electrode point position pair set comprises at least two candidate electrode point position pairs;
obtaining the average field intensity of a plurality of candidate electrode point pair sets by performing field intensity accumulation operation on the candidate electrode point pairs in each candidate electrode point pair set;
Determining a target candidate electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate electrode point position pair sets;
a set of target candidate electrode point pairs is determined as a first electrode layout scheme.
In one embodiment, the method for obtaining the plurality of candidate electrode point pair sets according to the order of the average field intensity of each electrode point pair from large to small comprises the following steps:
obtaining an electrode point pair sequence according to the sequence of the average field intensity of each electrode point pair from large to small;
performing the operation of acquiring combined electrode point pairs in parallel on each electrode point pair in the electrode point pair sequence to acquire a target combined electrode point pair of each electrode point pair; the target combined electrode point pair represents the combined electrode point pair with the largest average field intensity;
and correspondingly determining the target combined electrode point pair corresponding to each electrode point pair as a candidate electrode point pair set.
In one embodiment, determining the second electrode layout scheme according to the remaining preset electrode points includes:
according to the residual preset electrode point positions, acquiring the average field intensity of the region of interest when the electric stimulation is applied to different residual electrode point position pairs; each remaining electrode point pair is any two electrode points in the remaining preset electrode points;
And generating a second electrode layout scheme according to the average field intensity of each remaining electrode point position pair.
In one embodiment, the generating a second electrode layout scheme according to the average field intensity of each remaining electrode point pair comprises:
sequencing the average field intensity of each residual electrode point position pair, and acquiring a target residual electrode point position pair with the maximum average field intensity according to a sequencing result;
if the average field intensity of the target residual electrode point position pair is larger than a preset field intensity threshold value, determining the electrode point position and the electric field direction of the target residual electrode point position pair, and determining the electrode point position and the electric field direction as a second electrode layout scheme;
if the average field intensity of the target remaining electrode point pairs is smaller than the field intensity threshold value, determining a second electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each remaining electrode point pair.
In one embodiment, determining the second electrode layout scheme by performing a field strength summation operation on the average field strengths of each remaining electrode point pair comprises:
acquiring a plurality of candidate residual electrode point position pair sets according to the sequence of the average field intensity of each residual electrode point position pair from large to small; each set of candidate remaining electrode point pairs comprises at least two candidate remaining electrode point pairs;
Obtaining the average field intensity of a plurality of candidate residual electrode point position pair sets by performing field intensity accumulation operation on the candidate residual electrode point position pairs in each candidate residual electrode point position pair set;
determining a target candidate residual electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate residual electrode point position pair sets;
and determining the target candidate residual electrode point position pair set as a second electrode layout scheme.
In one embodiment, the obtaining a plurality of candidate remaining electrode point pair sets according to the order of the average field intensity of each remaining electrode point pair from big to small includes:
obtaining a sequence of residual electrode point pairs according to the sequence of the average field intensity of each residual electrode point pair from large to small;
performing the operation of acquiring combined residual electrode point pairs in parallel for each residual electrode point pair in the residual electrode point pair sequence to acquire target combined residual electrode point pairs of each residual electrode point pair; the target combined residual electrode point position pair represents the combined residual electrode point position pair with the largest average field intensity;
and correspondingly determining the target combination residual electrode point position pairs corresponding to each residual electrode point position pair as a candidate residual electrode point position pair set.
In one embodiment, the manner in which each electrode point pair is determined on a standard finite element model includes any of the following:
determining any two different preset electrode points on the surface of the standard finite element model as electrode point pairs;
any two different preset electrode points inside the standard finite element model are determined as electrode point pairs.
In a second aspect, there is provided an apparatus for generating an electrode layout scheme, the apparatus comprising:
the model acquisition module is used for acquiring a standard finite element model of the target part; the standard finite element model comprises a region of interest;
the field intensity acquisition module is used for acquiring the average field intensity of the region of interest when the electric stimulation is applied to different electrode point pairs according to a plurality of preset electrode point positions on the standard finite element model; each electrode point position pair is any two electrode point positions in a plurality of preset electrode point positions;
the scheme generating module is used for generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point position pair.
In a third aspect, the present application also provides a computer device. The computer device comprises a memory storing a computer program and a processor implementing the steps of the method in any of the embodiments of the first aspect described above when the computer program is executed.
In a fourth aspect, the present application also provides a computer-readable storage medium. The computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the method in any of the embodiments of the first aspect described above.
In a fifth aspect, the present application also provides a computer program product. The computer program product comprising a computer program which, when executed by a processor, implements the steps of the method in any of the embodiments of the first aspect described above.
The method, the device, the equipment, the medium and the program product for generating the electrode layout scheme firstly acquire a standard finite element model of a target part; the standard finite element model comprises a region of interest; then, according to a plurality of preset electrode points on a standard finite element model, acquiring the average field intensity of the region of interest when the electric stimulation is applied to different electrode point pairs; wherein each electrode point pair is any two electrode points in a plurality of preset electrode points; and finally, generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point pair. The method is based on a standard finite element model, and obtains an optimal electrode layout scheme according to the average field intensity of different electrode point pairs, which is equivalent to taking all electrode point pairs of the standard finite element model into consideration in the process of generating the electrode layout scheme, and applying electric stimulation to all electrode point pairs, wherein the average field intensity generated by the region of interest at the moment is used as the basis for selecting the optimal electrode layout scheme, so that the determined electrode layout scheme can ensure that the field intensity distribution of the region of interest is optimal.
Drawings
FIG. 1 is a block diagram of a computer device of a method of generating an electrode layout scheme in one embodiment;
FIG. 2 is a flow chart of a method of generating an electrode layout scheme in one embodiment;
FIG. 3 is a flow diagram of a method of obtaining a standard finite element model in one embodiment;
FIG. 4 is a flow chart of a method of obtaining average field strength in one embodiment;
FIG. 5 is a flow chart of a method for obtaining average field strength in another embodiment;
FIG. 6 is a flow chart of a method of generating an electrode layout scheme in another embodiment;
FIG. 7 is a flow chart of a method for generating a first electrode layout scheme according to another embodiment;
FIG. 8 is a flow chart of a method for generating a first electrode layout scheme according to another embodiment;
FIG. 9 is a flow chart of a method of candidate electrode point pair set determination in one embodiment;
FIG. 10 is a flow chart of a method for generating a first electrode layout scheme according to another embodiment;
FIG. 11 is a flow chart of a method of generating a second electrode layout scheme in one embodiment;
FIG. 12 is a flow chart of a method for generating a second electrode layout scheme according to another embodiment;
FIG. 13 is a flow chart of a method for generating a second electrode layout scheme according to another embodiment;
FIG. 14 is a flow chart of a method for determining a set of candidate remaining electrode point pairs in one embodiment;
FIG. 15 is a flow chart of a method for determining electrode point pairs in one embodiment;
FIG. 16 is a flow chart of a method of generating an electrode layout scheme in another embodiment;
fig. 17 is a schematic structural view of an electrode layout scheme generating device in one embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application.
The method for generating the electrode layout scheme can be applied to computer equipment. The computer device may be a terminal, and its internal structure may be as shown in fig. 1. The computer device includes a processor, a memory, an input/output interface, a communication interface, a display unit, and an input means. The processor, the memory and the input/output interface are connected through a system bus, and the communication interface, the display unit and the input device are connected to the system bus through the input/output interface. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The input/output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless mode can be realized through WIFI, a mobile cellular network, NFC (near field communication) or other technologies. The computer program is executed by a processor to implement a method of generating an electrode layout scheme. The display unit of the computer device is used for forming a visual picture, and can be a display screen, a projection device or a virtual reality imaging device. The display screen can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, can also be a key, a track ball or a touch pad arranged on the shell of the computer equipment, and can also be an external keyboard, a touch pad or a mouse and the like. It will be appreciated by those skilled in the art that the structure shown in fig. 1 is merely a block diagram of some of the structures associated with the present application and is not limiting of the computer device to which the present application may be applied, and that a particular computer device may include more or fewer components than shown, or may combine certain components, or have a different arrangement of components.
Electric field therapy is a therapy performed by portable, non-invasive medical devices. Taking tumor treatment electric field (Tumor Treatment Fields, TTF) as an example, in general, when tumor electric field treatment is performed, it is required to attach a plurality of electrode sheets to the brain skin to continuously provide a low-intensity, medium-frequency alternating electric field to a tumor area. While the electric field is distributed in a non-uniform manner in the brain, different array electrode layouts can result in large field strength differences in the tumor area.
In the related art, using a preset head model, the following drawbacks exist in the calculation electrode layout scheme:
(1) And the head is reconstructed by a 3D high-precision conductivity model without using a brain tissue segmentation algorithm, so that the electric field simulation calculation precision is affected.
(2) The model generated by using the deformed elliptical structure has larger difference from the model of the actual target position, and the array electrode layout can not be ensured to be effectively applied to the actual target position in practice.
(3) The cost of manually measuring the tumor area is excessive and coarse errors may exist.
(4) The corresponding array layout schemes are fixed at nine, and the layout is required to be a 3×3 array, so that the flexibility is low.
In summary, the electrode plate array mode in the electrode layout scheme calculated by the related technology has the problems of low singleness and low precision, so that the field intensity distribution of the region of interest cannot be determined to be optimal. Therefore, how to determine an electrode layout scheme for optimizing the field intensity distribution of the region of interest is a technical problem to be solved. Based on this, the present application proposes a method for generating an electrode layout scheme for determining that the field intensity distribution of the region of interest reaches the optimum, and the method for generating an electrode layout scheme will be described below by way of an embodiment.
In one embodiment, as shown in fig. 2, there is provided a method for generating an electrode layout scheme, including the steps of:
s202, acquiring a standard finite element model of a target part; the region of interest is included in the standard finite element model.
The finite element model is a model established by applying a finite element analysis method, and is a group of unit assemblies which are connected only at nodes, are transferred only by the nodes and are constrained only at the nodes. In the embodiment of the present application, the standard finite element model of the target site is a set of units in which the geometric mechanism of the continuous target site is discretized into a finite number of units, and a finite number of nodes are set in each unit, so that the continuum is regarded as an aggregate of a group of units connected only at the nodes. Meanwhile, the node value of the field function is selected as a basic unknown quantity, an approximate interpolation function is assumed in each unit to represent the distribution rule of the field function in the unit, and a finite element equation set for solving the node unknown quantity is built, so that the infinite degree of freedom problem in one continuous domain is converted into the finite degree of freedom problem in a discrete domain. The region of interest refers to a sample, mask, crop area, or other operation in the image that is typically used as an image classification, generated by selecting or using a method such as thresholding or conversion from other files. In the embodiment of the application, the standard finite element model comprises an interested region and a non-interested region, wherein the interested region is a tissue region obtained in an image of a target part according to processing means such as image recognition, image segmentation and the like; the non-region of interest refers to a region of tissue in the target site other than the tissue corresponding to the region of interest.
For example, if a brain tumor is to be treated, the target area is the head, the region of interest is a tumor region and a surrounding portion of the brain parenchyma region within the head, the region of non-interest is each tissue region of the head, such as the skull, spinal fluid, and the tissue regions such as gray matter and white matter that remove the tumor region and the surrounding portion of the brain parenchyma region, and the standard finite element model of the target area is the finite element model corresponding to the head containing the tumor region.
Specifically, the standard finite element model is a three-dimensional model of the target part by establishing a series of approximation functions to perform simulation on the basis of the image corresponding to the target part. It should be noted that, the standard finite element model may select different acquisition modes based on actual requirements. Optionally, if the target part structure is complex and the simulation period is long, the simulated standard finite element model can be obtained in an offline mode through preset target part image information and related function relations. Optionally, if the target site has a simple structure, the information scanning device and the model building module of the target site can be integrated in the same system to acquire the information of the target site in real time in an online manner, and further generate a standard finite element model of the target site.
S204, according to a plurality of preset electrode points on a standard finite element model, acquiring average field intensity of an interested region when electric stimulation is applied to different electrode point pairs; each electrode point pair is any two electrode points in a plurality of preset electrode points.
The preset electrode point positions are positions of a plurality of electrode plates established on a standard finite element model through different standards and are used for placing the electrode plates. The common preset electrode point positions can be determined according to empirical data, such as setting standards of EEG10-10, EEG10-20 and the like, and can also be calculated based on real-time data, such as real-time dividing the target position according to preset proportion and the like. Because the electric stimulation is generated based on a loop consisting of an anode and a cathode and acts on the whole standard finite element model, two electrode points are selected as electrode point pairs from a plurality of preset electrode points on the standard finite element model.
The average field strength is a physical quantity representing the electric field vector distribution obtained by the region of interest in the standard finite element model when the electrode point pairs are electrically stimulated. It should be noted that, although the field intensity distribution generated by the electrical stimulation acts on the whole standard finite element model, the main purpose is to optimize the field intensity distribution of the region of interest, so that a positive-negative electrical stimulation is applied to the electrode point pair, and only the average field intensity of the region of interest in the standard finite element model needs to be considered.
Further, in the standard finite element model, the region of interest is constructed from different cubic grid cells, and the electric field distribution of different grid cells in the region of interest is different for the same electrical stimulus due to the difference in distance between the different grid cells and the electrode point pair to which the electrical stimulus is applied. Thus, for the same electrical stimulus, the average field strength of the region of interest is obtained based on the individual electric field profiles of the individual grid cells.
Exemplary, in the standard finite element model, if the preset electrode point is represented as N 1 To N 8 If the region of interest comprises five cube grid cells, any two electrode points can be selected as electrode point pairs, and N can be selected 1 And N 8 To form electrode point pairs, N can also be 1 And N 14 And forming an electrode point position pair. When N is selected 1 And N 8 As electrode point position pairs, at N 1 And N 8 When positive and negative electric stimulation is applied, the five cube grid units of the region of interest generate different field intensity vector distributions respectively, and the average field intensity of the region of interest is obtained according to the field intensity vector distribution values generated by the five cube grid units of the region of interest.
S206, generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point position pair.
The average field strength of an electrode point pair refers to the average field strength generated by the region of interest when an electrical stimulus is applied to the electrode point pair. The electrode point position pair applying the electric stimulation acts on a plurality of grid cells in the region of interest, field intensity distribution corresponding to the grid cells is generated, and average field intensity is obtained according to the field intensity distribution corresponding to the grid cells. When electrical stimulation is applied to different pairs of electrode points, the average field strengths of a plurality of regions of interest are correspondingly generated.
The optimal electrode layout scheme is an electrode layout scheme for optimizing field intensity distribution of a region of interest of a target portion, and comprises the number of electrode point pairs, the electric field direction, the position and layout information of each electrode plate in the electrode point pairs. It should be noted that, in the electrode layout scheme, a plurality of groups of electrode point pairs are generally required to act on the target site at the same time, so the number of electrode point pairs in the electrode layout scheme may be one or more.
After the average field intensity of each electrode point pair is obtained, firstly, selecting one group or a plurality of groups of electrode point pairs which enable the field intensity distribution of the region of interest to be optimal according to preset conditions, such as a field intensity threshold value, a quantity threshold value or an array position; and integrating the electric field direction, the position and the layout information of one or more groups of electrode point pairs meeting preset conditions to serve as an optimal electrode layout scheme of the target part.
The electrode layout generating method provided by the embodiment of the application firstly acquires a standard finite element model of a target part; the standard finite element model comprises a region of interest; then, according to a plurality of preset electrode points on a standard finite element model, acquiring the average field intensity of the region of interest when the electric stimulation is applied to different electrode point pairs; wherein each electrode point pair is any two electrode points in a plurality of preset electrode points; and finally, generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point pair. The method is characterized in that an optimal electrode layout scheme is obtained according to the average field intensity of different electrode point pairs on the basis of a standard finite element model, and is equivalent to that in the process of generating the electrode layout scheme, all the electrode point pairs of the standard finite element model are considered, electric stimulation is applied to all the electrode point pairs, the average field intensity generated by an interested region at the moment is used as the basis for selecting the optimal electrode layout scheme, and the determined optimal electrode layout scheme can enable the field intensity distribution of the interested region to be optimal.
In the process of generating an optimal electrode layout scheme, the optimal electrode layout scheme is generally calculated on the basis of a standard finite element model so as to accurately acquire the average field intensity of the region of interest and generate the optimal electrode layout scheme. In addition, when the standard finite element model of the target part is constructed, the model can be obtained through simulation in a function construction mode, and can also be directly generated in an image processing mode. Based on this, a method of acquiring the standard finite element model is described below by way of one embodiment.
In one embodiment, as shown in FIG. 3, obtaining a standard finite element model of a target site includes:
s302, acquiring medical image data of at least one mode of a target part.
Wherein the medical image data is a discrete series of medical images, the medical images comprising the region of interest may be one or more. After acquiring the medical image data, first a medical image containing the region of interest is extracted from the medical image data, and then information in the region of interest is acquired in the medical image containing the region of interest.
Medical image data under different modalities refers to different types of medical image data generated by a scanning imaging device due to different scanning modes. For example, medical image data acquired by a nuclear magnetic resonance scanning apparatus is MRI data generated by imaging using the principles of nuclear magnetism and magnetic fields. Further, MRI data may be further divided into 2D MRI data and 3D MRI data according to different scan modes.
For a target site, a standard finite element model can be constructed by acquiring medical image data under one mode. When acquiring medical image data of two or more modes, the medical image data of multiple modes can be fused to generate fused medical image data, and then a standard finite element model is constructed according to the fused medical image data; the initial model can be generated according to the medical image data of one mode, then the initial model is corrected according to the medical image data of other modes, and the corrected initial model is used as a standard finite element model.
It should be noted that, the number of modes of the medical image data and the accuracy of the initial finite element model are positively correlated, the more the acquired medical image data modes are, the higher the accuracy of the standard finite element model correspondingly constructed is, the more the structural information of the target part can be truly restored, and the fewer the medical image data modes are, the lower the accuracy of the standard finite element model is.
S304, performing tissue segmentation on each medical image data to obtain tissue segmentation information of the target part.
Because the medical image data under the single mode is a series of two-dimensional images in nature, the medical image data under the single mode is subjected to tissue segmentation, and the tissue segmentation information of the target part is obtained and comprises image pixel position information and a real coordinate transformation matrix. The real coordinate conversion matrix is used to convert image data of a pixel unit into site area information of an actual unit.
Taking a tumor region of the segmented head medical image data as an example, if the medical image data includes N images, and the (3, 4) pixel point of the 10 th image belongs to the tumor region, acquiring (3, 4, 10) and a corresponding pixel-to-reality coordinate conversion matrix, and converting each pixel into a voxel. Then, taking the coordinate of each voxel as the center of the cube, and expanding a specified length, such as 5mm, in the x, y and z directions at the same time to form a cube region corresponding to each voxel, wherein the specification of each cube is 10mm by 10mm. It should be noted that, the medical image data may include a tumor area formed by a plurality of pixels, and then a plurality of cubes are corresponding to each tumor area, and when a plurality of cube areas overlap, the union of the cubes is regarded as the tumor area.
And in the medical image data of the target part, carrying out tissue segmentation on the medical image data in each mode, and dividing an interested region and a non-interested region in the medical image data to obtain tissue segmentation information corresponding to the medical image data in each mode. Illustratively, if the target site is the head, the region of interest is a tumor region within the head, and the regions of non-interest refer to various tissue regions of the head, such as gray matter, white matter, skull, spinal fluid, and the like.
It should be noted that the target site includes a plurality of tissue sections, and the conductivities of the tissue sections are not the same. That is, after tissue segmentation of the medical image data, the conductivity of each tissue is different. Alternatively, the spm12 device may be used for tissue segmentation of medical images. The manner in which the tissue is segmented is not limited in this application.
S306, generating a three-dimensional model of the target part according to the tissue segmentation information of the target part.
The tissue segmentation information of the target part comprises partitioned tissue partitions, and the tissue segmentation information of each tissue partition is discrete, so that the tissue segmentation information corresponding to medical image data of each mode is required to be integrated by utilizing an image fusion technology, the segmentation information of each tissue partition can be fused by utilizing different conductivities corresponding to each tissue partition, a finite element conductivity model corresponding to the target part is generated, and the finite element conductivity model is determined as a three-dimensional model of the target part. The three-dimensional model at this time has restored the structural information of the target site from the mesh layer.
S308, filling tetrahedral grids into the three-dimensional model of the target part to obtain a standard finite element model of the target part.
And filling the tetrahedral mesh into the three-dimensional model of the target part, and filling each tissue partition of the target part in a tetrahedral form so as to conveniently and intuitively acquire the field intensity distribution condition of the region of interest. Alternatively, a tetrahedral mesh may be filled using a three-dimensional finite element mesh generator Gmsh. The filling mode is not limited in the present application.
In the embodiment of the application, the initial finite element model is obtained by carrying out tissue segmentation on medical image data of different modes, so that each tissue partition information of the target part can be subdivided, and the standard finite element model is closer to the real structure of the target part. The standard finite element model obtained in this way considers the differences among different tissue partitions of the target part, and can be flexibly and pertinently constructed for the target part, so that the actual effectiveness of the electrode layout scheme is ensured. In addition, the three-dimensional model of the target part is filled with grids, so that the field intensity distribution condition of the region of interest can be intuitively acquired.
When the optimal electrode layout scheme is obtained, the average field intensity of interest is usually used as a criterion to ensure that the electrode layout scheme enables the field intensity distribution of the region of interest to be optimal. Based on this, a method of acquiring the average field strength will be described below by way of one embodiment.
In one embodiment, as shown in fig. 4, obtaining the average field strength of the region of interest when the electrical stimulus is applied at different pairs of electrode points from a plurality of preset electrode points on a standard finite element model, includes:
s402, for any electrode point position pair, acquiring a field intensity vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest when the electric stimulation is applied at the electrode point position pair.
When the electrode point position pair applies electric stimulation, the generated electric field is distributed in each area in the standard finite element model, and the electric stimulation intensities received by different areas are different. Since the present application aims to optimize the field strength distribution of the region of interest, it is necessary to obtain the field strength distribution of the region of interest when applying electrical stimulation at the electrode point pair.
In the standard finite element model, the region of interest is made up of a plurality of grid cells characterized in the form of cubes. Obviously, when the electric stimulus is applied to the electrode point position pairs, the stimulation to each grid unit in the region of interest is different, and the corresponding generated field intensity vector distribution is also different. It is therefore necessary to obtain the field strength vector distribution values of each grid cell in the region of interest when an electrical stimulus is applied at the electrode point pairs. In addition, in order to obtain the average field strength of the region of interest, the field strengths of all grid cells, that is, the average field strength of the region of interest, need to be obtained by combining the volume of each grid cell and the field strength vector distribution value of each grid cell.
S404, determining the average field intensity of the region of interest when the electric stimulus is applied at the electrode point position pair according to the field intensity vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest.
After acquiring the field strength vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest, it is necessary to acquire the average field strength of the region of interest. The field intensity vector distribution value of each grid cell can be represented by a 1×n distribution matrix, the number of columns of the distribution matrix depends on the number of grid cells, the volume of each grid cell can be represented by an n×1 volume matrix, the number of rows of the volume matrix corresponds to the number of columns of the distribution matrix, and the distribution matrix also depends on the data of the grid cells. Obviously, a specific value can be obtained according to the distribution matrix of each grid cell and the volume matrix of each grid cell, and the obtained value is taken as the average field intensity of the region of interest.
In the embodiment of the application, the region of interest is further divided into a plurality of grid cells, and different field intensity vectors of the grid cells are combined according to the volumes of the grid cells, so that the obtained average field intensity of the region of interest is high in interpretation and reliability.
When the average field intensity is obtained, the scalar field intensity distribution value of each grid unit and the volume of each grid unit can be directly multiplied to obtain a specific value as the average field intensity, or the volume of each grid unit can be further processed and then multiplied by the scalar field intensity distribution value of each grid unit to obtain the specific value as the average field intensity. Based on this, the manner of obtaining the average field strength will be described below by way of an example.
In one embodiment, as shown in fig. 5, determining the average field strength of the region of interest when the electrical stimulus is applied at the electrode point pair according to the field strength vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest includes:
s502, determining the ratio of each grid cell in the total volume of all grid cells as the volume weight of each grid cell.
After the volume of each grid cell is acquired, the volume weight of each grid cell may be acquired according to the ratio of each grid cell in the total volume of all grid cells.
S504, determining the average field intensity of the region of interest when the electric stimulation is applied to the electrode point position pair according to the field intensity vector distribution value of each grid cell and the volume weight of each grid cell.
The field intensity vector distribution value of each grid unit is in the three-dimensional directions of x, y and z, and obviously, the volume weight is a scalar quantity, and the average field intensity is also a scalar quantity. Therefore, it is necessary to first convert the field intensity vector distribution value of each grid cell into a scalar value of each grid cell, then perform a weighted sum of the scalar value of each grid cell and the volume weight of each grid cell, and finally use the weighted sum result as the average field intensity of the region of interest when the electric stimulation is applied at the electrode point pair.
Exemplary, if the region of interest comprises three grid cells, the corresponding volume values are V 1 、V 2 And V 3 When the electrode point pair is electrically stimulated, the electric field vector distribution generated by the three grid cells is (E x1 ，E y1 ，E z1 )，(E x2 ，E y2 ，E z2 ) And (E) x3 ，E y3 ，E z3 ) Scalar of electric field generated by three grid cells is E 1 、E 2 And E is 3 The representation, wherein,
Figure BDA0004068815150000081
Figure BDA0004068815150000082
then the volume weights of the three grid cells are V 1 /(V 1 +V 2 +V 3 )、V 2 /(V 1 +V 2 +V 3 )、V 3 /(V 1 +V 2 +V 3 ) The method comprises the steps of carrying out a first treatment on the surface of the The average field strength of the region of interest is: (E) 1 ×V 1 +E 2 ×V 2 +E 3 ×V 3 )/(V 1 +V 2 +V 3 )。
According to the embodiment of the application, the average field intensity of the region of interest is obtained according to the volume weight and the field intensity vector distribution value of each grid unit, and the vector result is reasonably subjected to scalar conversion, so that the defect that the field intensity vector distribution values cannot be compared is overcome.
In determining the optimal electrode layout scheme, the average field strength is generally used as a preset condition to determine the electrode layout scheme that optimizes the field strength distribution of the region of interest. Based on this, a method of generating an electrode layout scheme will be described below by way of one embodiment.
In one embodiment, as shown in fig. 6, the generating an optimal electrode layout scheme of the target site according to the average field intensity of each electrode point pair includes:
s602, determining a first electrode layout scheme according to the average field intensity of each electrode point position pair.
The standard finite element model of the target part comprises a plurality of preset electrode points, and two preset electrode points form a pair of electrode point pairs, namely the standard finite element model comprises a plurality of groups of electrode point pairs. The electric stimulation is applied to the electrode point position pairs, so that field intensity is generated in the region of interest of the target part, and the aim of possibly optimizing the field intensity distribution of the region of interest is fulfilled. This means that each electrode point pair may be an electrode layout scheme when applying electrical stimulation or when applying electrical stimulation in combination with a plurality of electrode point pairs.
In the process of determining the optimal electrode layout scheme, the preset conditions for meeting the first electrode layout scheme include conditions such as average field strength of electrode point pairs, the number of the electrode point pairs, an array of the electrode point pairs and the like. In the embodiments of the present application, a method of generating an electrode layout scheme based on an average field strength will be described in the case where other conditions are satisfied. Specifically, when electric stimulation is applied to each electrode point position pair, the average field intensity is calculated according to the field intensity generated by the region of interest, the calculated result is used as the average field intensity of each electrode pair, then the average field intensity is compared with the preset average field intensity, and the position and the direction of the electrode point position pair meeting the preset average field intensity are used as a first electrode layout scheme.
For example, when a group of electrode point pairs apply electric stimulation, and the generated average field intensity meets a preset condition, the position of the electrode point pairs and the direction of an electric field applying the electric stimulation are used as a first electrode layout scheme; when the electric stimulation is applied to the electrode point pairs, and the generated average field intensity meets the preset condition, the positions of the electrode point pairs and the electric field direction of the electric stimulation are used as a first electrode layout scheme.
S604, eliminating the electrode points corresponding to the first electrode layout scheme from a plurality of preset electrode points on the standard finite element model to obtain the remaining preset electrode points on the standard finite element model.
Although the first electrode layout solution can optimize the field intensity distribution of the region of interest, other preset electrode point pairs or combinations of other preset electrode point pairs that can optimize the field intensity distribution of the region of interest are not excluded from the remaining preset electrode points. Therefore, after the first electrode layout scheme is obtained, the electrode point positions corresponding to the first electrode layout scheme are removed from the preset electrode point positions, and the remaining preset electrode point positions on the standard finite element model are obtained and used for obtaining other electrode point position layouts which enable the field intensity distribution of the region of interest to be optimal.
S606, determining a second electrode layout scheme according to the remaining preset electrode points.
And (3) obtaining the average field intensity of the region of interest when the electric stimulation is applied to different preset residual electrode point pairs by applying the electric stimulation to the residual preset electrode point pairs on the standard finite element model, and comparing the average field intensity with preset conditions comprising the preset field intensity to obtain a second electrode layout scheme.
And S608, determining the first electrode layout scheme and the second electrode layout scheme as the optimal electrode layout scheme of the target part.
In consideration of the side effects of applying continuous electrical stimulation to the same electrode point pair in a single layout scheme, the embodiment of the application combines the practical situations, obtains two first electrode layout schemes and two second electrode layout schemes which reach preset conditions, and determines the first electrode layout scheme and the second electrode layout scheme as the optimal electrode layout scheme of the target part.
It should be noted that the first electrode layout scheme and the second electrode layout scheme work alternately with a certain period rule, so that the field intensity distribution of the region of interest in the target part can still be continuously optimized under the electrical stimulation of different electrode point pairs.
In the embodiment of the application, the first electrode layout scheme and the second electrode layout scheme are determined to be the optimal layout scheme of the target part, and the electrode point positions of the layout schemes are different in pairs and can act on the target part at the same time, so that the distribution of the electrode point positions of the target part is more uniform.
When determining the first electrode arrangement, the first electrode arrangement can be obtained by means of the average field strength. Based on this, a determination method of how to determine the electrode layout scheme from the average field strength is explained below by way of an example.
In one embodiment, as shown in fig. 7, determining the first electrode layout scheme according to the average field strength of each electrode point pair includes:
s702, sorting the average field intensity of each electrode point position pair, and acquiring the target electrode point position pair with the maximum average field intensity according to the sorting result.
The average field strength of an electrode point represents the overall field strength distribution of the region of interest when the electrode point pair is electrically stimulated. The average field intensity of the plurality of electrode point pairs corresponds to a plurality of field intensity distribution conditions of the region of interest when the electric stimulus is applied to the different electrode point pairs. The average field intensity of each electrode point pair is a scalar which is obtained based on the field intensity vector distribution value of each grid unit of the region of interest, the electrode point pairs with the largest average field intensity can be sequenced, and the electrode point pair with the largest average field intensity is obtained according to the sequencing result and is taken as the target electrode point pair.
And S704, if the average field intensity of the target electrode point pair is greater than a preset field intensity threshold, determining the electrode point position and the electric field direction of the target electrode point pair as a first electrode layout scheme.
The field intensity threshold is used for judging whether the field intensity distribution condition of the region of interest reaches the optimal condition.
If the average field strength of the target electrode point pair is greater than a preset field strength threshold, the fact that electric stimulation is applied to the target electrode point pair is indicated that the field strength distribution condition of the region of interest can be optimized, and then the electrode point position and the electric field direction of the target electrode point pair are determined to be the first electrode layout scheme.
S706, if the average field intensity of the target electrode point pair is smaller than the field intensity threshold value, determining a first electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each electrode point pair.
If the average field intensity of the target electrode point pair is smaller than a preset field intensity threshold value, which indicates that the electric stimulation is applied to the target electrode point pair, so that the field intensity distribution situation of the region of interest cannot be optimized, then the electric stimulation is required to be applied to the target electrode point pair and other preset electrode point pairs at the same time, the field intensity distribution situation of the region of interest is optimized by performing field intensity accumulation operation on the average field intensity of each electrode point pair, and the position and the electric field direction of each electrode point pair, which enables the field intensity distribution situation of the region of interest to be optimized, are determined as a first electrode layout scheme. In this case, the first electrode layout scheme includes a plurality of electrode dot pairs, and in the first electrode layout scheme, the directions of the plurality of electrode dot pairs are recorded with the electric field direction of the target electrode dot pair as the reference electric field direction.
In the embodiment of the application, the target electrode point position pair is determined according to the sequencing result of the average field intensity, the logic is clear, the position and the electric field direction of the target electrode point position pair are used as the first electrode layout scheme, and the obtained first electrode layout scheme is reasonable.
The above embodiment uses the electrode point position and the electric field direction of the target electrode point pair as the first electrode layout scheme by sorting the average field strengths, if the maximum average field strength reaches the preset field strength threshold, and if the maximum average field strength does not reach the preset field strength threshold, the field strength accumulating operation is performed on the average field strengths of the electrode point pair to obtain the first electrode layout scheme. In view of this, the following describes, by way of an embodiment, the manner in which the first electrode arrangement is obtained in the case where the maximum average field strength does not reach the preset field strength threshold.
In one embodiment, as shown in fig. 8, the determining the first electrode layout scheme by performing a field strength summation operation on the average field strengths of the electrode point pairs includes:
s802, acquiring a plurality of candidate electrode point pair sets according to the sequence of the average field intensity of each electrode point pair from large to small; each set of candidate electrode point pairs comprises at least two candidate electrode point pairs.
Because the average field intensity of one electrode point pair does not reach the preset field intensity threshold value, at least two electrode point pairs need to be combined to obtain a candidate electrode point pair set, and the average field intensity of the candidate electrode point pair set reaches or approaches to the preset field intensity threshold value. The average field intensity of each electrode point pair is ordered from big to small, and each electrode point pair is used as the first electrode point pair of the candidate electrode point pair set. That is, the number of candidate electrode point pair sets corresponds to the number of electrode point pairs.
S804, performing field intensity accumulation operation on the candidate electrode point position pairs in each candidate electrode point position pair set to obtain the average field intensity of the plurality of candidate electrode point position pair sets.
The candidate electrode point position pair set comprises a plurality of candidate electrode point position pairs, and each candidate electrode point position pair corresponds to field intensity vectors of the plurality of sub-grid cells and volume weights of the plurality of sub-grid cells respectively. And (3) describing a candidate electrode point position pair set, sequentially carrying out superposition/superposition on field intensity vectors of each candidate electrode point position pair in the candidate electrode point position pair set to obtain a field intensity vector superposition/superposition result, converting the superposition/superposition result into a set of scalar quantities, carrying out weighted summation on the scalar quantities and a set of volume weights to obtain corresponding average field intensity, comparing the average field intensity corresponding to the superposition/superposition result, and selecting each candidate electrode point position pair corresponding to the largest average field intensity for combination to be used as the candidate electrode point position pair set. It should be noted that the superposition result indicates that the electric field directions of the two candidate electrode point pairs are consistent, and the superposition result indicates that the electric field directions of the two candidate point pairs are inconsistent.
S806, determining a target candidate electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate electrode point position pair sets.
After the plurality of candidate electrode point position pair sets are obtained, determining an electric field treatment requirement according to actual requirements, wherein the electric field treatment requirement can be that the average field intensity of the candidate electrode point position pair sets meets a preset field intensity threshold value, or that the average field intensity of the candidate electrode point position pair sets does not meet the preset field intensity threshold value, but the candidate electrode point position pair numbers of the candidate electrode point position pair sets meet a preset quantity.
And determining a target candidate electrode point pair set meeting the electric field treatment requirement from the average field intensity of the plurality of candidate electrode point pair sets.
S808, determining the target candidate electrode point position pair set as a first electrode layout scheme.
It should be noted that the first electrode layout scheme in the embodiment of the present application is different from and includes different conditions for the determination of the first electrode layout scheme obtained in the above-described embodiment S704. The first electrode layout solution obtained in the above embodiment S704 is generated when the average field strength meets the preset field strength threshold, and includes the electrode point position and the electric field direction of one electrode point pair. The first electrode layout scheme in the embodiment of the application is generated when the average field intensity meets the preset field intensity threshold, or when the average field intensity does not meet the preset field intensity threshold but meets the preset quantity, and comprises the electrode point position and the electric field direction of at least two electrode point pairs.
In the embodiment of the application, the average field intensity of each electrode point pair is sequenced, and a plurality of candidate electrode point pair sets are constructed based on each electrode point pair, so that the candidate electrode point pair sets can cover all the electrode point pair combination modes. Further, after the plurality of candidate electrode point position pair sets are obtained, the mode of determining the target candidate electrode point position pair set considers the determining mode meeting the preset field intensity threshold value and the determining mode not meeting the preset field intensity threshold value, so that the determined target candidate electrode point position pair set, namely the first electrode layout scheme, is more reasonable.
When the candidate electrode point position pair set is obtained, a plurality of electrode point position pairs are combined in sequence according to a preset strategy, so that the candidate electrode point position pair set is ensured to be effective. In view of this, a description will be given below of a method of acquiring a candidate electrode point pair set by way of an embodiment.
In one embodiment, as shown in fig. 9, the obtaining a plurality of candidate electrode point pair sets according to the order of the average field intensity of each electrode point pair from big to small includes:
s902, obtaining an electrode point pair sequence according to the sequence of the average field intensity of each electrode point pair from large to small.
Each electrode point pair corresponds to the field intensity vector of the grid cells and the volume weight of the grid cells, and the average field intensity of the electrode point pair is obtained according to the weighted sum of the field intensity vector and the volume weight. For a plurality of electrode point pairs, a plurality of average field strengths exist, and the average field strengths of the electrode point pairs are sequenced to obtain an electrode point pair sequence.
For example, if there are three electrode point pairs A1, A2 and A3, the average field strengths corresponding to A1 to A3 are E1, E2, E3, and E2> E1> E3, then the electrode point pair sequence is: a2, A1, A3.
S904, executing the operation of acquiring combined electrode point pairs in parallel for each electrode point pair in the electrode point pair sequence to acquire a target combined electrode point pair of each electrode point pair; the target combined electrode point pair represents the combined electrode point pair with the largest average field strength.
The process of combining electrode point pairs in parallel of each electrode point pair of the electrode point pair sequence is to update the average field intensity of the electrode point pairs of the electrode point pair sequence and the corresponding combined electrode point pairs, and one combined electrode point pair is added during each update. Obviously, the average field strength of the combined electrode point pair is increased every time it is updated.
And when the average field intensity of the electrode point pairs of the electrode point pair sequence and the corresponding combined electrode point pairs meet the preset field intensity threshold value, the average field intensity is not updated any more.
When the average field intensity of the electrode point pairs of the electrode point pair sequence and the corresponding combined electrode point pairs do not meet the preset field intensity threshold, adding one combined electrode point pair, calculating the average field intensity of the electrode point pairs of the electrode point pair sequence and all updated combined electrode point pairs, comparing the average field intensity with the preset field intensity threshold again, if yes, not updating any more, otherwise, continuing updating until the number of the combined electrode point pairs reaches the preset number.
S906, determining the target combined electrode point pair corresponding to each electrode point pair as a candidate electrode point pair set.
The candidate electrode point position pair sets are candidates of the first electrode layout scheme.
The electrode point position pair sequence comprises a plurality of electrode point position pairs, and each electrode point position pair corresponds to one target combined electrode point position pair respectively through the operation of the combined electrode point position pair. Thus, the electrode point pair sequence corresponds to a plurality of candidate electrode point pair sets.
When each candidate electrode point position pair set is obtained, the number of electrode point position pairs in the electrode point position pair set is increased in sequence, and the construction process can ensure that the diversity of each candidate electrode point position pair set is allowed while different candidate electrode point position pair sets meet preset conditions. And the acquisition process of the plurality of candidate electrode point position pair sets is operated in parallel, and the acquisition time is shortened under the condition of computer resource permission.
In one embodiment, in the case that the average field strength of a single pair of electrode point pairs does not meet the preset field strength threshold, the obtaining manner of the first electrode layout scheme is as shown in fig. 10, taking n electrode point pairs as an example, and obtaining a sequence of electrode point pairs according to the order of the average field strengths of the electrode point pairs from large to small, and obtaining combined electrode point pairs for the electrode point pairs in parallel. The process of obtaining the combined electrode point pair and generating the candidate electrode point pair set is described for the electrode point pair A1 with the first high average field intensity. Firstly, calculating the electric field vector distribution of the target area corresponding to the A1 and the electric field vector distribution of the target area corresponding to a plurality of coordinate point pairs which do not contain the A1 respectively, and obtaining a plurality of average field strengths; acquiring the maximum average field intensity from the plurality of average field intensities, and an electrode point position pair A1 and a first combined electrode point position pair corresponding to the maximum average field intensity; then comparing the maximum average field strength with a preset field strength threshold, if the maximum average field strength reaches the preset field strength threshold, taking the A1 and the first combined electrode point pair as a candidate electrode point pair set, and if the maximum average field strength does not reach the preset field strength threshold, searching for a second combined electrode point pair in the electrode point pair which does not comprise the first combined electrode point pair and the A1; and respectively calculating the electric field vector distribution of the target area corresponding to the A1 and the first combined electrode point pair with the electric field vector distribution of the target area corresponding to a plurality of electrode point pairs which do not contain the A1 and the first combined electrode point pair, obtaining a plurality of average field strengths, selecting the maximum average field strength from the average field strengths and comparing with a preset field strength threshold value, taking the A1, the first combined electrode point pair and the second combined electrode point pair as a candidate electrode point pair set if the maximum average field strength reaches the preset field strength threshold value, continuously searching a third combined electrode point pair if the maximum average field strength does not reach the preset field strength threshold value, and so on until the number of the combined electrode point pairs reaches the preset number, and generating the candidate electrode point pair set. After all candidate electrode point position pair sets are obtained, the number and average field intensity of candidate electrode point position pairs in each candidate electrode point position pair set are obtained, and a target candidate electrode point position pair set is obtained according to treatment requirements and is used as a first electrode layout scheme.
In the process of obtaining the first electrode layout scheme, the embodiment of the application is used for sequencing the average field intensity of each electrode point pair and constructing a plurality of candidate electrode point pair sets based on each electrode point pair, so that the candidate electrode point pair sets can cover the combination mode of all the electrode point pairs. Further, after the plurality of candidate electrode point position pair sets are obtained, the mode of determining the target candidate electrode point position pair set considers the determining mode meeting the preset field intensity threshold value and the determining mode not meeting the preset field intensity threshold value, so that the determined target candidate electrode point position pair set, namely the first electrode layout scheme, is more reasonable.
When the optimal electrode scheme is obtained, two electrode layout schemes are usually determined, and the field intensity distribution of the target area is ensured to be uniform through the rotation work of the two electrode schemes. Based on this, the acquisition step of the second electrode layout scheme will be described below by way of one embodiment.
In one embodiment, as shown in fig. 11, determining the second electrode layout scheme according to the remaining preset electrode points includes:
s1102, according to the residual preset electrode point positions, acquiring the average field intensity of the region of interest when the electric stimulation is applied to different residual electrode point position pairs; each remaining electrode point pair is any two electrode points in the remaining preset electrode points.
S1104, generating a second electrode layout scheme according to the average field intensity of each remaining electrode point position pair.
Similar to the determination of the first electrode layout scheme, but the preset electrode points of the second electrode layout scheme are in preset point pairs of the first electrode layout scheme, and the electrode point pairs of the first electrode layout scheme are eliminated.
In the embodiment of the application, the second electrode layout scheme is calculated based on the same logic reasoning and field intensity vector distribution operation as the first electrode layout scheme, so that the whole optimal electrode layout scheme is clear in hierarchy and easy to realize on computer equipment.
When determining the second electrode arrangement, the second electrode arrangement can be obtained by averaging the basis of the field strengths. Based on this, a determination method of how to determine the electrode layout scheme from the average field strength is explained below by way of an example.
In one embodiment, as shown in fig. 12, the generating a second electrode layout scheme according to the average field strength of each remaining electrode point pair includes:
s1202, sorting the average field intensity of each residual electrode point position pair, and acquiring a target residual electrode point position pair with the maximum average field intensity according to a sorting result.
The average field strength of an electrode point represents the overall field strength distribution of the region of interest when the electrode point pair is electrically stimulated. The average field intensity of the plurality of remaining electrode point pairs corresponds to a plurality of field intensity distribution conditions of the region of interest when the electrical stimulation is applied to the different remaining electrode point pairs. The average field intensity of each residual electrode point position pair is a scalar which is obtained based on the field intensity vector distribution value of each grid unit of the region of interest, the residual electrode point positions with the maximum average field intensity can be sequenced, the residual electrode point position pair with the maximum average field intensity is obtained according to the sequencing result, and the residual electrode point position pair with the maximum average field intensity is taken as the target residual electrode point position pair.
And S1204, if the average field intensity of the target residual electrode point pair is greater than a preset field intensity threshold, determining the electrode point position and the electric field direction of the target residual electrode point pair, and determining the second electrode layout scheme.
If the average field intensity of the target remaining electrode point pair is larger than a preset field intensity threshold value, which indicates that the electric stimulation is applied to the target remaining electrode point pair, the field intensity distribution condition of the region of interest can be optimized, and then the electrode point position and the electric field direction of the target remaining electrode point pair are determined to be the second electrode layout scheme.
S1206, if the average field intensity of the target remaining electrode point pairs is smaller than the field intensity threshold, determining a second electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each remaining electrode point pair.
If the average field intensity of the target remaining electrode point position pair is smaller than a preset field intensity threshold value, the fact that the electric stimulation is applied to the target remaining electrode point position pair is indicated that the field intensity distribution condition of the region of interest cannot be optimized, then the electric stimulation is required to be applied to the target remaining electrode point position pair and other preset electrode point positions at the same time, the field intensity distribution condition of the region of interest is optimized by performing field intensity accumulation operation on the average field intensity of each remaining electrode point position pair, and the position and the electric field direction of each remaining electrode point position pair, which are optimal, of the field intensity distribution condition of the region of interest are determined to be the second electrode layout scheme. In this case, the second electrode layout scheme includes a plurality of sets of remaining electrode dot pairs, and in the second electrode layout scheme, the directions of the plurality of sets of remaining electrode dot pairs are recorded with the electric field direction of the target remaining electrode dot pair as the reference electric field direction.
In the embodiment of the application, the second electrode layout scheme is calculated based on the same logic reasoning and field intensity vector distribution operation as the first electrode layout scheme, so that the whole optimal electrode layout scheme is clear in hierarchy and easy to realize on computer equipment.
The above embodiment uses the positions of the remaining electrode points and the electric field directions of the target remaining electrode point pairs as the first electrode layout scheme by sorting the average field strengths, if the maximum average field strength reaches the preset field strength threshold, and if the maximum average field strength does not reach the preset field strength threshold, performs the field strength accumulation operation on the average field strengths of the remaining electrode point pairs, thereby obtaining the second electrode layout scheme. Based on this, the following describes, by way of an embodiment, the manner in which the second electrode arrangement is obtained in the case where the maximum average field strength does not reach the preset field strength threshold.
In one embodiment, as shown in fig. 13, the determining the second electrode layout scheme by performing a field strength summation operation on the average field strengths of each remaining electrode point pair includes:
s1302, acquiring a plurality of candidate residual electrode point position pair sets according to the sequence of the average field intensity of each residual electrode point position pair from large to small; each set of candidate remaining electrode point pairs comprises at least two candidate remaining electrode point pairs.
Because the average field intensity of one residual electrode point position pair does not reach the preset field intensity threshold value, at least two residual electrode point position pairs need to be combined to obtain a candidate residual electrode point position pair set, and the average field intensity of the candidate residual electrode point position pair set reaches or approaches to the preset field intensity threshold value. The average field intensity of each residual electrode point pair is ordered according to the order from big to small, and each residual electrode point pair is used as the first residual electrode point pair of the candidate residual electrode point pair set. That is, the number of candidate remaining electrode point pair sets corresponds to the number of remaining electrode point pairs.
S1304, performing field intensity accumulation operation on the candidate residual electrode point pairs in each candidate residual electrode point pair set to acquire the average field intensity of the plurality of candidate residual electrode point pair sets.
The candidate remaining electrode point position pair set comprises a plurality of candidate remaining electrode point position pairs, and each candidate remaining electrode point position pair corresponds to the field intensity vector of the plurality of sub-grid cells and the volume weight of the plurality of sub-grid cells respectively. And (3) describing a set of candidate residual electrode point position pairs, sequentially carrying out superposition/superposition on field intensity vectors of each candidate residual electrode point position pair in the set of candidate residual electrode point position pairs to obtain a field intensity vector superposition/superposition result, converting the superposition/superposition result into a group of scalar quantities, carrying out weighted summation on the group of scalar quantities and a group of volume weights to obtain corresponding average field intensity, comparing the average field intensity corresponding to the superposition/superposition result, and selecting each candidate residual electrode point position pair corresponding to the largest average field intensity for combination to be used as the set of candidate residual electrode point position pairs. It should be noted that the superposition result indicates that the electric field directions of the two candidate residual electrode point pairs are consistent, and the superposition result indicates that the electric field directions of the two candidate residual electrode point pairs are inconsistent.
S1306, determining a target candidate residual electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate residual electrode point position pair sets.
After the plurality of candidate residual electrode point position pair sets are obtained, determining an electric field treatment requirement according to actual requirements, wherein the electric field treatment requirement can be that the average field intensity of the candidate residual electrode point position pair sets meets a preset field intensity threshold, or that the average field intensity of the candidate residual electrode point position pair sets does not meet the preset field intensity threshold, but the candidate residual electrode point position pair numbers of the candidate residual electrode point position pair sets meet a preset quantity, or that the included angle between the candidate residual electrode point position pair sets and the first electrode layout scheme meets a preset angle. The included angle between the candidate remaining electrode point position pair set and the first electrode layout scheme is obtained through cosine theorem, firstly, a group of field intensity vectors of a plurality of grid units corresponding to the remaining electrode point position pair set and a group of field intensity vectors of a plurality of grid units corresponding to the first electrode layout scheme are obtained, then the included angle between each grid unit is obtained through cosine theorem according to two field intensity vectors corresponding to each grid unit, then an absolute value is obtained according to difference between 90 degrees and the included angle between each grid unit, an orthogonal offset angle corresponding to each grid unit is obtained, and the orthogonal offset angle of each grid unit and the volume weight corresponding to each grid unit are subjected to weighted summation operation, wherein the operation result is the included angle between the candidate remaining electrode point position pair set and the first electrode layout scheme.
S1308, the target candidate remaining electrode point pair set is determined as the second electrode layout scheme.
It should be noted that the second electrode layout scheme in the embodiment of the present application is different from and includes different conditions for the determination of the second electrode layout scheme obtained in the above-described embodiment S1204. The second electrode layout scheme obtained in the above embodiment S1204 is generated when the average field strength meets the preset field strength threshold, and includes the candidate electrode point position and the electric field direction of one candidate electrode point pair. The second electrode layout scheme in the embodiment of the application is generated when the average field intensity meets a preset field intensity threshold, or the average field intensity does not meet the preset field intensity threshold but meets the preset quantity, or the average field intensity does not meet the preset field intensity threshold, but the included angle between the candidate residual electrode point pair set and the first electrode layout scheme meets the preset angle, and the second electrode layout scheme comprises residual electrode point positions and electric field directions of at least two residual electrode point pairs.
In the embodiment of the application, the average field intensity of each residual electrode point position pair is sequenced, and a plurality of candidate residual electrode point position pair sets are constructed based on each residual electrode point position pair, so that the candidate residual electrode point position pair sets can cover the combination mode of all the residual electrode point position pairs. Further, after the plurality of candidate remaining electrode point position pair sets are obtained, the mode of determining the target candidate remaining electrode point position pair set considers the determining mode meeting the preset field intensity threshold value and the determining mode not meeting the preset field intensity threshold value, so that the determined target candidate remaining electrode point position pair set, namely the second electrode layout scheme is more reasonable.
When the candidate remaining electrode point pair set is obtained, a plurality of remaining electrode point pairs are combined in sequence according to a preset strategy, so that the candidate remaining electrode point pair set is ensured to be effective. Based on this, the following describes a method of acquiring the candidate remaining electrode point pair set by an embodiment.
In one embodiment, as shown in fig. 14, obtaining a plurality of candidate remaining electrode point pair sets in order of increasing average field strength of each remaining electrode point pair includes:
s1402, obtaining a sequence of the residual electrode point pairs according to the sequence of the average field intensity of each residual electrode point pair from large to small.
And each remaining electrode point position pair corresponds to the field intensity vector of the grid cells and the volume weight of the grid cells, and the average field intensity of the remaining electrode point position pairs is obtained according to the weighted sum of the field intensity vector and the volume weight. For a plurality of residual electrode point pairs, a plurality of average field strengths exist, and the average field strengths of the residual electrode point pairs are sequenced to obtain a residual electrode point pair sequence.
S1404, performing the operation of obtaining combined residual electrode point pairs in parallel for each residual electrode point pair in the residual electrode point pair sequence to obtain target combined residual electrode point pairs of each residual electrode point pair; the target combined remaining electrode point pair represents the combined remaining electrode point pair having the largest average field strength.
The method comprises the steps of combining the residual electrode point pairs in parallel for each residual electrode point pair of the residual electrode point pair sequence, wherein the process of combining the residual electrode point pairs is to update according to the residual electrode point pairs of the residual electrode point pair sequence and the average field intensity of the corresponding combined residual electrode point pairs, and one combined residual electrode point pair is added during each updating. Obviously, the average field intensity of the combined remaining electrode point pairs is increased every time it is updated.
And when the average field intensity of the residual electrode point pairs of the residual electrode point pair sequence and the corresponding combined residual electrode point pairs meet the preset field intensity threshold value, the residual electrode point pairs are not updated.
When the average field intensity of the residual electrode point pairs of the residual electrode point pair sequence and the corresponding combined residual electrode point pairs do not meet the preset field intensity threshold, adding one combined residual electrode point pair, calculating the average field intensity of the residual electrode point pairs of the residual electrode point pair sequence and all updated combined residual electrode point pairs, comparing with the preset field intensity threshold again, if yes, not updating, otherwise, continuing updating until the number of the combined residual electrode point pairs reaches the preset number.
And S1406, correspondingly determining the target combination residual electrode point position pairs corresponding to each residual electrode point position pair as a candidate residual electrode point position pair set.
The candidate remaining electrode point position pair sets are candidates of the second electrode layout scheme.
The residual electrode point position pair sequence comprises a plurality of residual electrode point position pairs, and each residual electrode point position pair corresponds to one target combined residual electrode point position pair through the operation of combining the residual electrode point position pairs. Thus, the sequence of remaining electrode point pairs corresponds to a plurality of candidate sets of remaining electrode point pairs.
It should be appreciated that the first electrode layout scheme is performed sequentially with the second electrode layout scheme, the distinction comprising two parts of the number of electrode site pairs and the treatment requirements. In terms of the number of electrode dot pairs, the first electrode layout scheme is acquired in all of the electrode dot pairs, and the second electrode layout scheme is acquired in the remaining dot pairs excluding the electrode dot corresponding to the first electrode layout scheme. In terms of treatment requirements, the first electrode layout scheme is determined according to the number of electrode point pairs or the magnitude of average field intensity, and the second electrode layout scheme is determined according to the number of electrode point pairs, the magnitude of average field intensity or the included angle between the first electrode layout scheme and the candidate remaining electrode point pair set.
When each candidate residual electrode point position pair set is obtained, the number of residual electrode point position pairs in the residual electrode point position pair set is increased in sequence, and the construction process can ensure that the diversity of each candidate residual electrode point position pair set is allowed while different candidate residual electrode point position pair sets meet preset conditions. And the acquisition process of the plurality of candidate residual electrode point position pair sets is operated in parallel, and the acquisition time is shortened under the condition of computer resource permission. In addition, in the embodiment of the application, the second electrode layout scheme is calculated based on the same logic reasoning and field intensity vector distribution operation as the first electrode layout scheme, so that the whole optimal electrode layout scheme is clear in level and easy to realize on computer equipment.
The foregoing embodiments illustrate how to determine the electrode layout based on the average field strength of the electrode point pairs, and it is apparent that the manner in which the electrode point pairs are selected is also important. In view of this, a manner of determining the electrode point pairs will be described below by way of an example.
In one embodiment, as shown in fig. 15, the manner in which each electrode point pair is determined on the standard finite element model includes any one of the following:
S1502, determining any two different preset electrode points on the surface of the standard finite element model as electrode point pairs.
It should be appreciated that when applying electrical stimulation to a standard finite element model, electrode pads are typically placed at two different locations (i.e., electrode sites) of the standard finite element model and an electrical stimulus of plus-minus is applied. Thus, the electrode site pair is composed of two electrode sites to which electrical stimulation is applied.
On a standard finite element model, firstly, a plurality of preset electrode points are determined, each preset electrode point corresponds to a different position on the surface of the standard finite element model, then, among the plurality of preset electrode points, two preset electrode point pairs with different positions are selected at will to serve as a group of electrode point pairs, and finally, all the preset electrode points are traversed to determine the electrode point pairs.
Taking n preset electrode points on the surface of a standard finite element model as an example (n>2) Then according to n preset electrode points, it can be determined
Figure BDA0004068815150000161
Pairs of electrode sites.
Alternatively, the preset electrode point position may be determined according to the existing standard, and common setting standards include EEG10-10, EEG10-20, and other standards, or may be calculated based on real-time data, for example, dividing the surface of the standard finite element model according to a preset proportion.
S1504, determining any two different preset electrode points in the standard finite element model as electrode point pairs.
On a standard finite element model, firstly, a plurality of preset electrode points are determined, each preset electrode point corresponds to different positions in the standard finite element model, then, among the plurality of preset electrode points, two preset electrode point pairs with different positions are selected at will to serve as a group of electrode point pairs, and finally, all the preset electrode points are traversed to determine the electrode point pairs.
Optionally, after determining the position of the preset electrode point on the surface of the standard finite element model, acquiring the internal position of the electrode point placed in the standard finite element model, preprocessing the inside of the standard finite element model according to a preset rule, and determining the corresponding position of the preset electrode point in the preprocessed standard finite element model.
Taking a standard finite element model as an example of a head model, the head model comprises a series of tissues such as scalp, skull and the like, after the internal positions of electrode points are determined, the tissues corresponding to the interior of the head model can be removed to realize pretreatment of the head model, so that the electrode points can be placed in the head model after pretreatment.
In the embodiment of the application, the mode of determining the electrode point position pairs on the surface of the standard finite element model and in the standard finite element model is provided respectively, so that the determined mode is convenient for a computer to select according to actual conditions, and the determined electrode point position pairs are more closely related to the standard finite element model based on the mode.
In one embodiment, as shown in fig. 16, an embodiment is provided in combination with the method for generating an electrode layout scheme, taking a target site as a head, and a region of interest as a tumor region as an example, and the embodiment includes the following steps:
s1601, acquiring head MRI image data: at least acquiring T1 enhanced 3D MRI data and T2Flair data with resolution not less than 1 x 1mm, additionally, 3D MRI data of T2 and 2D MRI data of T1 can be acquired, so as to further correct errors of a head model and improve tumor recognition accuracy.
S1602, extracting the target tumor region: the tumor region is extracted using a head tumor recognition model based on machine learning, and checked with a doctor, and the doctor ensures that the tumor region contains a tumor and reduces the erroneous recognition region as much as possible as a target region.
S1603, segmenting the head MRI image: brain tissue segmentation was performed using spm12 for head MRI images.
S1604, generating a head finite element model: and generating a 3d model according to the brain tissue segmentation information, filling the tetrahedral mesh by using Gmsh, and generating a head finite element model.
S1605, the target tumor region is mapped to a target volume: and acquiring corresponding registration matrixes according to the corresponding 2D and 3D MRI data, and mapping the target tumor area to the target volume of the head finite element model according to the registration matrixes.
S1606, a personalized electrode layout is calculated.
S1607, outputting a personalized electrode layout.
According to the embodiment of the application, the optimal electrode layout scheme is obtained according to the average field intensity of different electrode point pairs on the basis of the standard finite element model, which is equivalent to that in the process of generating the electrode layout scheme, all the electrode point pairs of the standard finite element model are considered, electric stimulation is applied to all the electrode point pairs, and the average field intensity generated by the region of interest at the moment is taken as the basis for selecting the optimal electrode layout scheme, so that the determined optimal electrode layout scheme can enable the field intensity distribution of the region of interest to be optimal.
In order to verify the average field strength of the region of interest corresponding to the electrode layout scheme in the embodiment of the present application, in one embodiment, as shown in table 1, the average field strengths of the regions of interest corresponding to 18 electrode layouts are counted, where 1 to 17 are the average field strengths of the regions of interest corresponding to the electrode layout obtained by the conventional method, and 18 is the average field strength of the region of interest corresponding to the first electrode layout scheme in the embodiment of the present application.
TABLE 1
Sequence number Average field intensity (V/M)
1 71.46
2 66.69
3 46.78
4 30,57
5 146.59
6 155.55
7 134.67
8 122.08
9 120.46
10 114.59
11 114.00
12 108.51
13 84.36
14 105.53
15 101.72
16 74.24
17 71.25
18 199.32
As can be seen from table 1, the average field intensity of the region of interest corresponding to the first electrode layout scheme in the embodiment of the present application is 199.32V/M, which is significantly better than the average field intensity corresponding to the layout scheme in the conventional method, and further demonstrates the effectiveness of the method for generating the electrode layout scheme in the embodiment of the present application.
It should be understood that, although the steps in the flowcharts related to the above embodiments are sequentially shown as indicated by arrows, these steps are not necessarily sequentially performed in the order indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in the flowcharts described in the above embodiments may include a plurality of steps or a plurality of stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of the steps or stages is not necessarily performed sequentially, but may be performed alternately or alternately with at least some of the other steps or stages.
Based on the same inventive concept, the embodiment of the application also provides a generating device of the electrode layout scheme for realizing the generating method of the electrode layout scheme. The implementation of the solution provided by the device is similar to the implementation described in the above method, so the specific limitation in the embodiment of the generating device of one or more electrode layout schemes provided below may refer to the limitation of the generating method of the electrode layout scheme hereinabove, and will not be repeated herein.
In one embodiment, as shown in fig. 17, there is provided a generating apparatus 1700 of an electrode layout scheme, comprising: a model acquisition module 1702, a field strength acquisition module 1704, and a scenario generation module 1706, wherein:
the model acquisition module 1702 is configured to acquire a standard finite element model of a target portion; the standard finite element model comprises a region of interest;
the field intensity obtaining module 1704 is configured to obtain, according to a plurality of preset electrode points on the standard finite element model, an average field intensity of the region of interest when the electrical stimulation is applied at different electrode point pairs; each electrode point position pair is any two electrode point positions in a plurality of preset electrode point positions;
The scheme generating module 1706 is configured to generate an optimal electrode layout scheme of the target site according to the average field intensity of each electrode point pair.
In one embodiment, the model acquisition module 1702 includes a first acquisition unit, a second acquisition unit, a first generation unit, and a third acquisition unit, where:
the first acquisition unit is used for acquiring medical image data of the target part in at least one mode;
the second acquisition unit is used for carrying out tissue segmentation on each medical image data to obtain tissue segmentation information of the target part;
the first generation unit is used for generating a three-dimensional model of the target part according to the tissue segmentation information of the target part;
and the third acquisition unit is used for filling the tetrahedral mesh into the three-dimensional model of the target part to obtain a standard finite element model of the target part.
In one embodiment, the field strength acquisition module 1704 includes a fourth acquisition unit and a first determination unit, wherein:
a fourth acquisition unit for acquiring, for any one of the electrode point pairs, a field intensity vector distribution value of each grid cell in the region of interest and a volume of each grid cell in the region of interest when the electrical stimulation is applied at the electrode point pair;
And the first determining unit is used for determining the average field intensity of the region of interest when the electric stimulus is applied to the electrode point position pair according to the field intensity vector distribution value of each grid unit in the region of interest and the volume of each grid unit in the region of interest.
In one embodiment, the first determining unit comprises a first determining subunit and a first obtaining subunit, wherein:
a first determining subunit configured to determine, as a volume weight of each grid cell, a ratio of each grid cell to a total volume of all grid cells;
the first acquisition subunit is used for determining the average field intensity of the region of interest when the electric stimulation is applied to the electrode point position pair according to the field intensity vector distribution value of each grid unit and the volume weight of each grid unit.
In one embodiment, the scenario generation module 1706 includes a second determination unit, a fifth acquisition unit, a third determination unit, and a fourth determination unit, wherein:
the second determining unit is used for determining a first electrode layout scheme according to the average field intensity of each electrode point position pair;
a fifth obtaining unit, configured to reject, from a plurality of preset electrode points on the standard finite element model, the electrode points corresponding to the first electrode layout scheme, so as to obtain remaining preset electrode points on the standard finite element model;
The third determining unit is used for determining a second electrode layout scheme according to the remaining preset electrode points;
and a fourth determining unit configured to determine the first electrode layout scheme and the second electrode layout scheme as optimal electrode layout schemes for the target region.
In one embodiment, the second determining unit comprises a second acquisition subunit, a second determining subunit, and a third determining subunit, wherein:
the second acquisition subunit is used for sequencing the average field intensity of each electrode point position pair, and acquiring a target electrode point position pair with the maximum average field intensity according to the sequencing result;
the second determining subunit is configured to determine, as a first electrode layout scheme, an electrode point position and an electric field direction of the target electrode point pair if an average field strength of the target electrode point pair is greater than a preset field strength threshold;
and the third determination subunit is used for determining the first electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each electrode point pair if the average field intensity of the target electrode point pair is smaller than the field intensity threshold value.
In one embodiment, the third determining subunit is further configured to obtain a plurality of candidate electrode point pair sets according to an order from the big to the small of average field strength of each electrode point pair; each candidate electrode point position pair set comprises at least two candidate electrode point position pairs; obtaining the average field intensity of a plurality of candidate electrode point pair sets by performing field intensity accumulation operation on the candidate electrode point pairs in each candidate electrode point pair set; determining a target candidate electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate electrode point position pair sets; a set of target candidate electrode point pairs is determined as a first electrode layout scheme.
In one embodiment, the third determining subunit is further configured to obtain a sequence of electrode point pairs according to the order of from the big to the small of the average field strength of each electrode point pair; performing the operation of acquiring combined electrode point pairs in parallel on each electrode point pair in the electrode point pair sequence to acquire a target combined electrode point pair of each electrode point pair; the target combined electrode point pair represents the combined electrode point pair with the largest average field intensity; and correspondingly determining the target combined electrode point pair corresponding to each electrode point pair as a candidate electrode point pair set.
In one embodiment, the third determining unit comprises a third acquisition subunit and a fourth acquisition subunit, wherein:
the third acquisition subunit is used for acquiring the average field intensity of the region of interest when the electric stimulation is applied to different residual electrode point pairs according to the residual preset electrode point positions; each remaining electrode point pair is any two electrode points in the remaining preset electrode points;
and the fourth acquisition subunit is used for generating a second electrode layout scheme according to the average field intensity of each remaining electrode point position pair.
In one embodiment, the fourth obtaining subunit is further configured to sort the average field strengths of the remaining electrode point pairs, and obtain, according to the sorting result, a target remaining electrode point pair with the maximum average field strength; if the average field intensity of the target residual electrode point position pair is larger than a preset field intensity threshold value, determining the electrode point position and the electric field direction of the target residual electrode point position pair, and determining the electrode point position and the electric field direction as a second electrode layout scheme; if the average field intensity of the target remaining electrode point pairs is smaller than the field intensity threshold value, determining a second electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each remaining electrode point pair.
In one embodiment, the fourth obtaining subunit is further configured to obtain a plurality of candidate remaining electrode point pair sets according to an order from the greater average field strength of each remaining electrode point pair to the lesser average field strength; each set of candidate remaining electrode point pairs comprises at least two candidate remaining electrode point pairs; obtaining the average field intensity of a plurality of candidate residual electrode point position pair sets by performing field intensity accumulation operation on the candidate residual electrode point position pairs in each candidate residual electrode point position pair set; determining a target candidate residual electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate residual electrode point position pair sets; and determining the target candidate residual electrode point position pair set as a second electrode layout scheme.
In one embodiment, the fourth obtaining subunit is further configured to obtain a sequence of remaining electrode point pairs according to a sequence from the greater average field strength of each remaining electrode point pair to the lesser average field strength; performing the operation of acquiring combined residual electrode point pairs in parallel for each residual electrode point pair in the residual electrode point pair sequence to acquire target combined residual electrode point pairs of each residual electrode point pair; the target combined residual electrode point position pair represents the combined residual electrode point position pair with the largest average field intensity; and correspondingly determining the target combination residual electrode point position pairs corresponding to each residual electrode point position pair as a candidate residual electrode point position pair set.
In one embodiment, the field strength acquisition module further comprises a fifth determination unit and a sixth determination unit, wherein:
a fifth determining unit, configured to determine any two different preset electrode points on the surface of the standard finite element model as an electrode point pair;
and a sixth determining unit, configured to determine any two different preset electrode points inside the standard finite element model as an electrode point pair.
The respective modules in the generation means of the above-described electrode layout scheme may be implemented in whole or in part by software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, including a memory and a processor, where the memory stores a computer program, and the processor implements the technical solution of the method for generating an electrode layout solution provided in any of the foregoing embodiments when executing the computer program.
The computer device provided in the foregoing embodiments has similar implementation principles and technical effects to those of the foregoing method embodiments, and will not be described herein in detail.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which when executed by a processor implements the technical solution of the method for generating an electrode layout solution provided in any of the above embodiments.
The foregoing embodiment provides a computer readable storage medium, which has similar principles and technical effects to those of the foregoing method embodiment, and will not be described herein.
In an embodiment, a computer program product is provided, which comprises a computer program, which when executed by a processor implements the technical solution of the method for generating an electrode layout solution provided in any of the above embodiments.
The foregoing embodiment provides a computer program product, which has similar principles and technical effects to those of the foregoing method embodiment, and will not be described herein.
It should be noted that, the user information (including, but not limited to, user equipment information, user personal information, etc.) and the data (including, but not limited to, data for analysis, stored data, presented data, etc.) referred to in the present application are information and data authorized by the user or sufficiently authorized by each party, and the collection, use and processing of the related data are required to comply with the related laws and regulations and standards of the related countries and regions.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, database, or other medium used in the various embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical Memory, high density embedded nonvolatile Memory, resistive random access Memory (ReRAM), magnetic random access Memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric Memory (Ferroelectric Random Access Memory, FRAM), phase change Memory (Phase Change Memory, PCM), graphene Memory, and the like. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, and the like. By way of illustration, and not limitation, RAM can be in the form of a variety of forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), and the like. The databases referred to in the various embodiments provided herein may include at least one of relational databases and non-relational databases. The non-relational database may include, but is not limited to, a blockchain-based distributed database, and the like. The processors referred to in the embodiments provided herein may be general purpose processors, central processing units, graphics processors, digital signal processors, programmable logic units, quantum computing-based data processing logic units, etc., without being limited thereto.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The foregoing examples represent only a few embodiments of the present application, which are described in more detail and are not thereby to be construed as limiting the scope of the present application. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application shall be subject to the appended claims.

Claims (17)

1. A method of generating an electrode layout scheme, the method comprising:
obtaining a standard finite element model of a target part; the standard finite element model comprises a region of interest;
acquiring average field intensity of the region of interest when electric stimulation is applied to different electrode point pairs according to a plurality of preset electrode point positions on the standard finite element model; each electrode point position pair is any two electrode point positions in the plurality of preset electrode point positions;
And generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point position pair.
2. The method of claim 1, wherein the obtaining a standard finite element model of the target site comprises:
acquiring medical image data of a target part in at least one mode;
performing tissue segmentation on each piece of medical image data to obtain tissue segmentation information of the target part;
generating a three-dimensional model of the target part according to the tissue segmentation information of the target part;
and filling the tetrahedral mesh into the three-dimensional model of the target part to obtain a standard finite element model of the target part.
3. The method according to claim 1 or 2, wherein the obtaining the average field strength of the region of interest when applying electrical stimulation at different pairs of electrode points from a plurality of preset electrode points on the standard finite element model comprises:
for any electrode point position pair, acquiring a field intensity vector distribution value of each grid cell in the region of interest and a volume of each grid cell in the region of interest when electric stimulation is applied to the electrode point position pair;
And determining the average field intensity of the region of interest when the electric stimulus is applied to the electrode point position pair according to the field intensity vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest.
4. A method according to claim 3, wherein said determining the average field strength of the region of interest when the electrical stimulus is applied at the pair of electrode sites from the field strength vector distribution value of each grid cell in the region of interest and the volume of each grid cell in the region of interest comprises:
determining the ratio of each grid cell in the total volume of all grid cells as the volume weight of each grid cell;
and determining the average field intensity of the region of interest when the electric stimulus is applied to the electrode point position pair according to the field intensity vector distribution value of each grid cell and the volume weight of each grid cell.
5. The method according to claim 1 or 2, wherein generating an optimal electrode layout scheme for the target site based on the average field strength of each of the electrode point pairs comprises:
determining a first electrode layout scheme according to the average field intensity of each electrode point position pair;
Removing electrode points corresponding to the first electrode layout scheme from a plurality of preset electrode points on the standard finite element model to obtain residual preset electrode points on the standard finite element model;
determining a second electrode layout scheme according to the residual preset electrode points;
and determining the first electrode layout scheme and the second electrode layout scheme as the optimal electrode layout scheme of the target part.
6. The method of claim 5, wherein determining a first electrode layout scheme based on the average field strength of each of the electrode point pairs comprises:
sequencing the average field intensity of each electrode point pair, and acquiring a target electrode point pair with the maximum average field intensity according to a sequencing result;
if the average field intensity of the target electrode point pair is larger than a preset field intensity threshold value, determining the electrode point position and the electric field direction of the target electrode point pair as the first electrode layout scheme;
and if the average field intensity of the target electrode point pair is smaller than the field intensity threshold value, determining the first electrode layout scheme by executing field intensity accumulation operation on the average field intensity of each electrode point pair.
7. The method of claim 6, wherein said determining said first electrode layout scheme by performing a field strength summation operation on the average field strengths of each of said electrode point pairs comprises:
acquiring a plurality of candidate electrode point position pair sets according to the sequence of the average field intensity of each electrode point position pair from large to small; each candidate electrode point position pair set comprises at least two candidate electrode point position pairs;
obtaining average field intensity of a plurality of candidate electrode point pair sets by performing field intensity accumulation operation on the candidate electrode point pairs in each candidate electrode point pair set;
determining a target candidate electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate electrode point position pair sets;
and determining the target candidate electrode point position pair set as the first electrode layout scheme.
8. The method of claim 7, wherein the obtaining a plurality of candidate electrode point pair sets in order of increasing average field strength of each of the electrode point pairs comprises:
obtaining an electrode point pair sequence according to the sequence of the average field intensity of each electrode point pair from large to small;
Executing the operation of acquiring combined electrode point pairs in parallel for each electrode point pair in the electrode point pair sequence to acquire a target combined electrode point pair of each electrode point pair; the target combined electrode point position represents the combined electrode point position pair with the largest average field intensity;
and correspondingly determining the target combined electrode point pair corresponding to each electrode point pair as a candidate electrode point pair set.
9. The method of claim 5, wherein determining a second electrode layout scheme based on the remaining preset electrode points comprises:
acquiring average field intensity of the region of interest when electric stimulation is applied to different residual electrode point pairs according to the residual preset electrode point positions; each remaining electrode point pair is any two electrode points in the remaining preset electrode points;
and generating the second electrode layout scheme according to the average field intensity of each remaining electrode point position pair.
10. The method of claim 9, wherein generating the second electrode layout scheme based on the average field strength of each of the remaining electrode point pairs comprises:
sequencing the average field intensity of each residual electrode point position pair, and acquiring a target residual electrode point position pair with the maximum average field intensity according to a sequencing result;
If the average field intensity of the target residual electrode point pair is larger than a preset field intensity threshold value, determining the electrode point position and the electric field direction of the target residual electrode point pair, and determining the electrode point position and the electric field direction as the second electrode layout scheme;
and if the average field intensity of the target remaining electrode point position pairs is smaller than the field intensity threshold value, determining the second electrode layout scheme by performing field intensity accumulation operation on the average field intensity of each remaining electrode point position pair.
11. The method of claim 10, wherein said determining said second electrode layout scheme by performing a field strength summation operation on the average field strengths of each of said remaining electrode point pairs comprises:
acquiring a plurality of candidate residual electrode point position pair sets according to the sequence of the average field intensity of each residual electrode point position pair from large to small; each candidate residual electrode point position pair set comprises at least two candidate residual electrode point position pairs;
obtaining average field intensity of a plurality of candidate residual electrode point pair sets by performing field intensity accumulation operation on the candidate residual electrode point pairs in each candidate residual electrode point pair set;
determining a target candidate residual electrode point position pair set meeting the electric field treatment requirement from the average field intensity of the candidate residual electrode point position pair sets;
And determining the target candidate residual electrode point position pair set as the second electrode layout scheme.
12. The method of claim 11, wherein the obtaining a plurality of candidate sets of remaining electrode point pairs in order of increasing average field strength for each remaining electrode point pair comprises:
obtaining a sequence of residual electrode point pairs according to the sequence of the average field intensity of each residual electrode point pair from large to small;
performing the operation of obtaining combined residual electrode point pairs in parallel for each residual electrode point pair in the residual electrode point pair sequence to obtain target combined residual electrode point pairs of each residual electrode point pair; the target combined residual electrode point position pair represents a combined residual electrode point position pair with the largest average field intensity;
and correspondingly determining the target combination residual electrode point position pairs corresponding to each residual electrode point position pair as a candidate residual electrode point position pair set.
13. The method of claim 1, wherein determining each of the electrode point pairs on the standard finite element model comprises any one of:
determining any two different preset electrode points on the surface of the standard finite element model as electrode point pairs;
And determining any two different preset electrode points in the standard finite element model as electrode point pairs.
14. An apparatus for generating an electrode layout scheme, the apparatus comprising:
the model acquisition module is used for acquiring a standard finite element model of the target part; the standard finite element model comprises a region of interest;
a field intensity acquisition module; the method comprises the steps of obtaining average field intensity of a region of interest when electric stimulation is applied to different electrode point pairs according to a plurality of preset electrode point positions on the standard finite element model; each electrode point position pair is any two electrode point positions in the plurality of preset electrode point positions;
and the scheme generating module is used for generating an optimal electrode layout scheme of the target part according to the average field intensity of each electrode point position pair.
15. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any one of claims 1 to 13 when the computer program is executed.
16. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 13.
17. A computer program product comprising a computer program, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any one of claims 1 to 13.
CN202310085786.0A 2023-01-19 2023-01-19 Method, apparatus, device, medium and program product for generating electrode layout scheme Pending CN116227277A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310085786.0A CN116227277A (en) 2023-01-19 2023-01-19 Method, apparatus, device, medium and program product for generating electrode layout scheme

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310085786.0A CN116227277A (en) 2023-01-19 2023-01-19 Method, apparatus, device, medium and program product for generating electrode layout scheme

Publications (1)

Publication Number Publication Date
CN116227277A true CN116227277A (en) 2023-06-06

Family

ID=86572484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310085786.0A Pending CN116227277A (en) 2023-01-19 2023-01-19 Method, apparatus, device, medium and program product for generating electrode layout scheme

Country Status (1)

Country Link
CN (1) CN116227277A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117438039A (en) * 2023-12-19 2024-01-23 湖南安泰康成生物科技有限公司 Method and device for determining application distribution of electrode plates
CN117954050A (en) * 2024-03-25 2024-04-30 湖南安泰康成生物科技有限公司 Electrode slice application scheme determining method and device, equipment and storage medium

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117438039A (en) * 2023-12-19 2024-01-23 湖南安泰康成生物科技有限公司 Method and device for determining application distribution of electrode plates
CN117438039B (en) * 2023-12-19 2024-03-22 湖南安泰康成生物科技有限公司 Method and device for determining application distribution of electrode plates
CN117954050A (en) * 2024-03-25 2024-04-30 湖南安泰康成生物科技有限公司 Electrode slice application scheme determining method and device, equipment and storage medium

Similar Documents

Publication Publication Date Title
CN116227277A (en) Method, apparatus, device, medium and program product for generating electrode layout scheme
CN110235175B (en) Online learning enhanced map-based automatic segmentation
CN116127808A (en) Electrode layout scheme generation method, device, computer equipment and storage medium
CN109903268A (en) It determines the method for the Exception Type of spine image collection and calculates equipment
CN109903269A (en) It determines the method for the Exception Type of backbone cross-sectional image and calculates equipment
CN111210465A (en) Image registration method and device, computer equipment and readable storage medium
US11308662B2 (en) System and method for image reconstruction
CN113129418B (en) Target surface reconstruction method, device, equipment and medium based on three-dimensional image
CN116350958A (en) Method, device, equipment, medium and program product for setting radiotherapy plan parameters
US11455755B2 (en) Methods and apparatus for neural network based image reconstruction
CN112950684A (en) Target feature extraction method, device, equipment and medium based on surface registration
CN116630239A (en) Image analysis method, device and computer equipment
CN116071404A (en) Image registration method, device, computer equipment and storage medium
CN112669450B (en) Human body model construction method and personalized human body model construction method
CN116128895A (en) Medical image segmentation method, apparatus and computer readable storage medium
CN114169468A (en) Image classification method, computer device, storage medium, and computer program product
CN115861248A (en) Medical image segmentation method, medical model training method, medical image segmentation device and storage medium
CN114723723A (en) Medical image processing method, computer device and storage medium
CN114299010A (en) Method and device for segmenting brain tumor image, computer equipment and storage medium
CN114972026A (en) Image processing method and storage medium
CN114668498A (en) Sequence recognition method of mark points, surgical robot system and storage medium
CN113780519A (en) Method and device for generating confrontation network training, computer equipment and storage medium
CN114388100B (en) Method and device for establishing muscle mapping model and computer equipment
CN116797726B (en) Organ three-dimensional reconstruction method, device, electronic equipment and storage medium
CN110853043B (en) Image segmentation method and device, readable storage medium and electronic equipment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination