CN115659767A - Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium - Google Patents

Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium Download PDF

Info

Publication number
CN115659767A
CN115659767A CN202211702466.7A CN202211702466A CN115659767A CN 115659767 A CN115659767 A CN 115659767A CN 202211702466 A CN202211702466 A CN 202211702466A CN 115659767 A CN115659767 A CN 115659767A
Authority
CN
China
Prior art keywords
target
particle
structural parameters
parameters
coil
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
CN202211702466.7A
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.)
Institute of Advanced Technology University of Science and Technology of China
Original Assignee
Institute of Advanced Technology University of Science and Technology of China
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 Institute of Advanced Technology University of Science and Technology of China filed Critical Institute of Advanced Technology University of Science and Technology of China
Priority to CN202211702466.7A priority Critical patent/CN115659767A/en
Publication of CN115659767A publication Critical patent/CN115659767A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

The invention discloses a magnetic resonance radio frequency coil design method, a device, equipment and a readable storage medium, and belongs to the technical field of magnetic resonance imaging. The method comprises the following steps: constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target, and acquiring the structural parameters of the target coil based on the initial geometric model; simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil; and acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation and analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil. The coil is optimally designed through the joint simulation of COMSOL software and MATLAB software, and the purpose of improving the optimization efficiency of the coil is achieved.

Description

Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium
Technical Field
The present invention relates to the field of magnetic resonance imaging, and in particular, to a method, an apparatus, a device, and a readable storage medium for designing a magnetic resonance radio frequency coil.
Background
Magnetic Resonance Imaging (MRI) has the advantages of high resolution, no ionizing radiation, multi-parameter Imaging and the like, and is widely applied to medical clinical examination. Radio frequency coils (RF coils) are a key component of magnetic resonance imaging scanners (MRI scanners) and are an important part of obtaining images with high signal strength and good signal-to-noise ratio. Good RF coils can receive MRI signals uniformly in the imaging region, i.e. B of the coil 1 The magnetic field should be spatially uniform within the ROI. Meanwhile, in order to improve the signal-to-noise ratio of a magnetic resonance image and improve the image quality of thin-layer scanning, high-resolution scanning and a low-field machine, a multi-channel radio frequency coil is a main basic magnetic resonance coil model.
To obtain the optimal coil structure, the three-dimensional model of the coil is usually B-modeled using electromagnetic simulation software such as HFSS, CST, FEKO, COMSOL, etc 1 Simulating the radio frequency magnetic field, calculating the electromagnetic field distribution of coils in the region of interest by using a finite element numerical analysis method, and optimizing the structure to obtain spatially uniform B 1 A magnetic field.
However, although the electromagnetic field simulation software can simulate the radio frequency magnetic field distribution in the region of interest, the optimization capability is weak, and the application of the electromagnetic field simulation software in engineering optimization design is limited, specifically, the electromagnetic field simulation software has the characteristics of low speed, long time consumption, large workload, poor optimization effect and the like.
The above is only for the purpose of assisting understanding of the technical aspects of the present invention, and does not represent an admission that the above is prior art.
Disclosure of Invention
The invention mainly aims to provide a magnetic resonance radio frequency coil design method, a magnetic resonance radio frequency coil design device, magnetic resonance radio frequency coil design equipment and a readable storage medium, and aims to solve the technical problem that the existing coil optimization method is low in efficiency.
In order to achieve the above object, the present invention provides a method for designing a magnetic resonance radio frequency coil, comprising the steps of:
constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target, and acquiring the structural parameters of the target coil based on the initial geometric model;
simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
and acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation and analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil.
Optionally, the step of optimizing the structural parameter based on MATLAB scientific computational analysis software, the structural parameter, the performance parameter, and the constraint condition to obtain the target structural parameter of the target coil includes:
setting a population scale and a spatial dimension according to the structural parameters, and respectively corresponding the structural parameters to each particle in the spatial dimension;
optimizing the structural parameters by a particle swarm algorithm based on the MATLAB scientific calculation analysis software, the structural parameters, the performance parameters, the population scale, the space dimensions and the constraint conditions to obtain target structural parameters of the target coil.
Optionally, the constraint condition includes a structural constraint condition and a performance constraint condition, and the step of optimizing the structural parameter by a particle swarm algorithm based on the MATLAB scientific computational analysis software, the structural parameter, the performance parameter, the population scale, the spatial dimension, and the constraint condition to obtain the target structural parameter of the target coil includes:
determining the particle position of each particle based on the MATLAB scientific calculation analysis software, the structure parameters and the performance parameters, and acquiring a preset position range of each particle;
calculating the fitness value of each particle through a fitness function according to the particle position of each particle;
for each particle, traversing to obtain a historical optimal position corresponding to each particle, comparing the fitness value of the particle position with the fitness value of the historical optimal position, and determining the particle position with a high fitness value or the historical optimal position as the optimal position of the current particle;
selecting a target optimal position with the largest fitness value from the particle positions aiming at all the particles, acquiring historical population optimal positions corresponding to all the particles, comparing the fitness value of the target optimal position with the fitness value of the historical population optimal position, and determining the target optimal position with the largest fitness value or the historical population optimal position as the current population optimal position;
determining the optimal position of each particle according to the optimal position of the current population;
judging whether the optimized position meets the structural constraint condition or not based on the preset position range;
if so, judging whether the performance parameter of each particle at the optimized position meets the performance constraint condition; if so, taking the optimized position as a target structure parameter;
and if not, updating the optimized position according to the preset position range, and executing the step of judging whether the performance parameter of each particle on the optimized position meets the performance constraint condition.
Optionally, the step of determining the optimal position of each particle according to the optimal position of the current population includes:
acquiring the particle speed of each particle and a preset speed range corresponding to each particle, and constraining the particle speed in the preset speed range;
and updating the particle position as an optimized position based on the particle speed and the current population optimal position.
Optionally, the step of calculating the fitness value of each particle through a fitness function according to the particle position of each particle includes:
performing electromagnetic simulation through a fitness function based on the COMSOL software and the particle position of each particle to obtain a target simulation result;
and determining the fitness value of each particle based on the MATLAB scientific calculation analysis software and the target simulation result.
Optionally, the step of obtaining a constraint condition corresponding to the imaging target, and optimizing the structural parameter based on MATLAB scientific computational analysis software, the structural parameter, the performance parameter, and the constraint condition to obtain a target structural parameter of the target coil includes:
acquiring constraint conditions corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions;
performing electromagnetic simulation based on the COMSOL software and the optimized structural parameters to obtain optimized performance parameters;
and if the optimized performance parameters meet the performance parameter targets corresponding to the imaging targets, taking the optimized structure parameters as the target structure parameters of the target coil.
Optionally, after the step of performing electromagnetic simulation based on the COMSOL software and the optimized structure parameters to obtain optimized performance parameters, the method includes:
if the optimized structural parameters do not meet the performance parameter target corresponding to the imaging target, the COMSOL multi-physical field simulation-based software and the structural parameters are repeatedly executed, the target coil is simulated to obtain the performance parameters of the target coil until the constraint conditions corresponding to the imaging target are obtained, and the structural parameters are optimized to obtain the target structural parameters of the target coil based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions.
In addition, to achieve the above object, the present invention also provides a magnetic resonance radio frequency coil designing apparatus, including:
the acquisition module is used for constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target and acquiring the structural parameters of the target coil based on the initial geometric model;
the simulation module is used for simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
and the optimization module is used for acquiring a constraint condition corresponding to the imaging target and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil.
In addition, to achieve the above object, the present invention also provides a magnetic resonance radio frequency coil designing apparatus, including: a memory, a processor and a computer program stored on the memory and executable on the processor, the computer program being configured to implement the steps of the magnetic resonance radio frequency coil design method.
Furthermore, to achieve the above object, the present invention further provides a readable storage medium having stored thereon a computer program, which when executed by a processor, implements the steps of the magnetic resonance radio frequency coil design method.
In one technical solution provided in the embodiment of the present invention, an initial geometric model of a target coil is constructed based on a main magnetic field direction of a magnetic resonance system and design requirements of an imaging target, and structural parameters of the target coil are obtained based on the initial geometric model; then, simulating the target coil based on COMSOL software and structural parameters to obtain performance parameters of the target coil; and acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil. Compared with the method that the simulation and optimization are completed in only one simulation software, the technical scheme provided by the embodiment of the invention fully considers the advantages of the COMSOL software and the MATLAB software, determines the performance parameters by using the COMSOL software after determining the initial geometric model and the structural parameters of the target coil, and then realizes the optimization design of the structural parameters of the coil by using the MATLAB software to realize algorithm programming, so that not only can the highly accurate electromagnetic simulation be realized, but also the automatic optimization of the algorithm can be realized, and the multi-parameter debugging is not required to be performed step by depending on the design experience of workers, so that the workload is greatly reduced, the optimization speed is accelerated, and the optimization precision is improved.
Drawings
FIG. 1 is a schematic diagram of a magnetic resonance radio frequency coil design apparatus in a hardware operating environment according to an embodiment of the present invention;
FIG. 2 is a schematic flow chart of a first embodiment of a method for designing a magnetic resonance RF coil according to the present invention;
FIG. 3 is a schematic diagram of a finite element geometric model according to a first embodiment of the design method of an MR RF coil of the present invention;
FIG. 4 is a flow chart illustrating a seventh embodiment of a method of designing a magnetic resonance RF coil in accordance with the present invention;
FIG. 5 is a functional block diagram of an embodiment of an apparatus for designing an MR RF coil according to the present invention.
The implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Referring to fig. 1, fig. 1 is a schematic structural diagram of a magnetic resonance radio frequency coil design device in a hardware operating environment according to an embodiment of the present invention.
As shown in fig. 1, the magnetic resonance radio frequency coil designing apparatus may include: a processor 1001, such as a Central Processing Unit (CPU), a communication bus 1002, a user interface 1003, a network interface 1004, and a memory 1005. Wherein a communication bus 1002 is used to enable connective communication between these components. The user interface 1003 may include a Display (Display), an input unit such as a Keyboard (Keyboard), and the optional user interface 1003 may also include a standard wired interface, a wireless interface. The network interface 1004 may optionally include a standard wired interface, a WIreless interface (e.g., a WIreless-FIdelity (WI-FI) interface). The Memory 1005 may be a Random Access Memory (RAM) Memory, or may be a Non-Volatile Memory (NVM), such as a disk Memory. The memory 1005 may alternatively be a storage device separate from the processor 1001.
Those skilled in the art will appreciate that the configuration shown in figure 1 does not constitute a limitation of the magnetic resonance radio frequency coil design apparatus and may include more or fewer components than shown, or some components may be combined, or a different arrangement of components.
As shown in fig. 1, a memory 1005, which is a kind of storage medium, may include therein an operating system, a data storage module, a network communication module, a user interface module, and a computer program.
In the magnetic resonance radio frequency coil design apparatus shown in fig. 1, the network interface 1004 is mainly used for data communication with other apparatuses; the user interface 1003 is mainly used for data interaction with a user; the processor 1001 and the memory 1005 of the magnetic resonance radio frequency coil design apparatus of the present invention may be provided in the magnetic resonance radio frequency coil design apparatus, which calls the computer program stored in the memory 1005 through the processor 1001 and executes the magnetic resonance radio frequency coil design method provided by the embodiment of the present invention.
An embodiment of the present invention provides a method for designing a magnetic resonance radio frequency coil, and referring to fig. 2, fig. 2 is a schematic flow chart of a first embodiment of the method for designing a magnetic resonance radio frequency coil according to the present invention.
In this embodiment, the method for designing a magnetic resonance radio frequency coil includes:
step S11: constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target, and acquiring the structural parameters of the target coil based on the initial geometric model;
it can be understood that COMSOL is a large-scale high-level numerical simulation software, and based on Finite Element Method (FEM), the simulation of real physical phenomena is realized by solving partial differential equations, and highly accurate numerical simulation is realized with high-efficiency computational performance and outstanding multi-field bidirectional direct coupling analysis capability. MATLAB as a professional mathematic software has strong data analysis and numerical calculation capacity, and provides a comprehensive solution for scientific research, engineering design and a plurality of scientific fields which need to carry out effective numerical calculation. Moreover, the COMSOL software is actually a Toolbox originated from MATLAB software and originally named as Toolbox1.0, so that the two types of software are compatible with each other, and the respective advantages of the two types of software can be fully utilized for joint simulation.
Specifically, an initial geometric model of a target coil is constructed based on a main magnetic field direction of a magnetic resonance system and design requirements of an imaging target, where the main magnetic field direction of the magnetic resonance system and the design requirements of the imaging target are specifically set by technicians according to factors such as detection purposes, device precision, and patient groups, and the present embodiment is not limited to this, and the target coil includes, but is not limited to, a breast intervention coil, an ankle joint coil, a head and neck joint coil, a body coil, a spine coil, and the like. It should be noted that the step of specifically constructing the initial geometric model of the target coil may be directly performed in the COMSOL software, or may be implemented in other software having a modeling function.
Further, structural parameters of the target coil, including but not limited to width, arc, height, spacing between multi-channel coils, are obtained based on the initial geometric model. It should be noted that the structural parameters of the coil are used to characterize the structural condition of the coil and the related background condition, and since the structural parameters of the radio frequency coil are closely related to the signal-to-noise ratio and the uniformity of the image reconstructed by the coil receiving the magnetic resonance signal, the structural parameters of the radio frequency coil are a core problem in the coil design process.
Exemplarily, fig. 3 is a finite element geometric model of a magnetic resonance breast interventional coil built in COMSOL software, wherein 1-a single-channel magnetic resonance breast interventional coil unit; 2-magnetic resonance mammary gland interventional coil radian; 3-two channel interval of the magnetic resonance mammary gland interventional coil; 4-set breast interventional coil target radio frequency magnetic field interested region.
Step S12: simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
it can be understood that after the structural parameters of the target coil are obtained based on the initial geometric model, the performance evaluation needs to be performed on the structural parameters in the COMSOL software, and specifically, the COMSOL software performs electromagnetic simulation on the structural parameters of the target coil to obtain the performance parameters of the target coil, including but not limited to scattering parameters, magnetic field distribution, and minimum magnetic field strength. It should be noted that the performance parameter is used to characterize the radio frequency magnetic field generated by the target coil, and may not only be used to determine whether the current target coil is up to standard, but also play a role in constraint in the subsequent optimization process.
Step S13: acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil;
it can be understood that, in order to ensure that the target coil does not deform or distort excessively in the specific optimization process, certain constraint conditions need to be set for the structural parameters and/or the performance parameters of the target coil, and the constraint conditions can be further determined according to the imaging target, for example, the radian range of the target coil is a-B, the interval between two channels is C-D, and the magnetic field uniformity is E.
It can be understood that COMSOL software itself is less capable of optimization, and therefore, the step of coil optimization needs to be completed by MATLAB software. Specifically, the compatibility between the COMSOL software and the MATLAB software is utilized, and the data in the COMSOL software is directly saved as MATLAB code and imported into the MATLAB software.
Further, based on MATLAB scientific calculation analysis software, structural parameters, performance parameters and constraint conditions, the structural parameters are optimized, and the specific process is that the MATLAB software acquires the structural parameters, the performance parameters and the constraint conditions from MATLAB codes, substitutes the parameters into an optimization algorithm, determines the optimal structural parameters through continuous optimization iteration, and finally obtains the target structural parameters of the target coil as the optimization result so as to ensure that the final magnetic resonance coil structure has a good radio frequency magnetic field meeting the optimization target in an interested region (namely an MRI imaging region). The optimization algorithm comprises but is not limited to a particle swarm algorithm, a genetic algorithm, a simulated annealing algorithm, taboo search and a neural network, the particle swarm algorithm simulates the foraging behavior of bird groups, the simulated annealing thought is derived from the annealing process of solid matters in physics, the genetic algorithm refers to the evolutionary thought of superior and inferior natural world, the taboo search simulates the intelligence process of human beings with a memory process, the neural network directly simulates the human brain, the embodiment is not specifically limited as to the selection of the specific algorithm, technicians can select the algorithm according to actual conditions, and a plurality of optimization algorithms can be organically integrated together to make up for deficiencies of the advantages, so that the performance is more excellent.
It should be noted that after the structural parameters are optimized by the MATLAB software to obtain the target structural parameters, the target structural parameters can be further verified by using the MATLAB software or the COMSOL software to ensure the rationality of the target structural parameters.
In a technical solution provided in this embodiment, an initial geometric model of a target coil is constructed based on a main magnetic field direction of a magnetic resonance system and design requirements of an imaging target, and structural parameters of the target coil are obtained based on the initial geometric model; then, simulating the target coil based on COMSOL software and structural parameters to obtain performance parameters of the target coil; and then acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil. Compared with the method that simulation and optimization are completed in only one simulation software, the technical scheme provided by the embodiment fully considers respective advantages of the COMSOL software and the MATLAB software, determines the performance parameters by using the COMSOL software after determining the initial geometric model and the structural parameters of the target coil, and then performs optimization design on the structural parameters of the coil by using the MATLAB software to realize algorithm programming, so that not only can high-precision electromagnetic simulation be realized, but also automatic optimization of the algorithm can be realized, multi-parameter debugging is not required to be performed step by depending on design experience of workers, the workload is greatly reduced, the optimization speed is increased, and the optimization accuracy is improved.
Further, a second embodiment of the method of designing a magnetic resonance radio frequency coil of the present invention is presented. Based on the embodiment shown in fig. 2, the step of optimizing the structure parameter based on MATLAB scientific computational analysis software, the structure parameter, the performance parameter, and the constraint condition to obtain the target structure parameter of the target coil includes:
step S21: setting a population scale and a spatial dimension according to the structural parameters, and respectively corresponding the structural parameters to each particle in the spatial dimension;
it can be understood that MATLAB software can optimize structural parameters using a particle swarm algorithm, which is a random search algorithm based on swarm cooperation developed by simulating the foraging behavior of a bird swarm, and the basic idea is to find an optimal solution through cooperation and information sharing among individuals in the swarm. Specifically, the technician sets basic parameters such as population size N (solution range, i.e., total number of particles) and spatial dimension D (three-direction coordinate optimization) according to actual conditions, and sets the structural parameters corresponding to each particle in the spatial dimension as independent variables to be optimized.
Step S22: optimizing the structural parameters by a particle swarm algorithm based on the MATLAB scientific calculation analysis software, the structural parameters, the performance parameters, the population scale, the space dimensions and the constraint conditions to obtain target structural parameters of the target coil.
Specifically, the structural parameters, the performance parameters, the population scale, the spatial dimensions and the constraint conditions are used as input parameters, the structural parameters are optimized through a particle swarm algorithm to obtain target structural parameters of the target coil, the specific process can be understood as solving the fitness function to obtain the optimal value of the fitness function, and the current structural parameters in the optimal value are used as the target structural parameters of the target coil.
It should be noted that, not only the structural parameters, the performance parameters, the population size, the spatial dimensions and the constraint conditions may be used as input parameters, but also the weight coefficients, the correction coefficients and the like may be added to ensure the accuracy of the output result.
In a technical solution provided in this embodiment, a population scale and a spatial dimension are set according to a structural parameter, the structural parameter is respectively corresponding to each particle in the spatial dimension, and then the structural parameter is optimized by a particle swarm optimization based on MATLAB software, the structural parameter, a performance parameter, the population scale, the spatial dimension, and a constraint condition, so as to obtain a target structural parameter of a target coil. Compared with other optimization algorithms, the technical scheme provided by the embodiment optimizes the structural parameters by adopting the particle swarm optimization, has the advantages of simple principle of the particle swarm optimization, easiness in realization and less parameters needing to be adjusted, and has better effects on higher convergence speed and avoidance of premature trapping in a local optimal solution by setting the self-cognition item and the swarm cognition item.
Further, a third embodiment of the method of designing a magnetic resonance radio frequency coil of the present invention is presented. Based on the second embodiment, the constraint conditions include structural constraint conditions and performance constraint conditions, and the step of optimizing the structural parameters by using a particle swarm algorithm based on the MATLAB scientific computational analysis software, the structural parameters, the performance parameters, the population scale, the spatial dimensions, and the constraint conditions to obtain the target structural parameters of the target coil includes:
step S31: determining the particle position of each particle based on the MATLAB scientific calculation analysis software, the structure parameters and the performance parameters, and acquiring a preset position range of each particle;
specifically, the particle position of each particle is determined based on MATLAB software, structural parameters, and performance parameters, where each particle corresponds to one structural parameter, for example, particle 1 is a coil length, particle 2 is a coil width, and particle 3 is a coil radian, which is not specifically limited in this embodiment. The position of a particle is used to indicate the numerical value corresponding to the particle.
Further, in order to ensure the rationality of the subsequent optimization process and the optimization result, the preset position range of each particle is obtained and used for limiting the variation range of the position. It should be noted that the preset position range is specifically set by a technician according to factors such as a detection purpose, equipment precision, a patient group, and the like, and the embodiment is not particularly limited.
Step S32: calculating the fitness value of each particle through a fitness function according to the particle position of each particle;
specifically, the fitness value corresponding to the particle position of each particle is calculated through a fitness function according to the particle position of each particle, it should be noted that the step of calculating the fitness value may be implemented through numerical calculation in MATLAB software, or may be implemented through electromagnetic simulation in COMSOL software, and this embodiment is not particularly limited.
Step S33: for each particle, traversing and acquiring a historical optimal position corresponding to each particle, comparing the fitness value of the particle position with the fitness value of the historical optimal position, and determining the particle position with a large fitness value or the historical optimal position as the optimal position of the current particle;
specifically, for each particle, the historical optimal position corresponding to each particle is obtained in a traversing manner, the fitness value of the particle position is compared with the fitness value of the historical optimal position, if the fitness value of the particle position is larger, the particle position is used as the optimal position of the current particle, and if the fitness value of the historical optimal position is larger, the historical optimal position is used as the optimal position of the current particle. It should be noted that, if the historical optimal position corresponding to each particle cannot be obtained initially, the comparison is not required, and the particle position is directly used as the optimal position of the current particle.
Step S34: selecting a target optimal position with the largest fitness value from the particle positions aiming at all the particles, acquiring historical population optimal positions corresponding to all the particles, comparing the fitness value of the target optimal position with the fitness value of the historical population optimal position, and determining the target optimal position with the largest fitness value or the historical population optimal position as the current population optimal position;
specifically, for all particles, the optimal target position with the largest fitness value is selected from the particle positions, the optimal historical population positions corresponding to all the particles are obtained, the fitness value of the optimal target position is compared with the fitness value of the optimal historical population position, if the fitness value of the optimal target position is larger, the optimal target position is determined as the optimal current population position, and if the fitness value of the optimal historical population position is larger, the optimal historical population position is determined as the optimal current population position. It should be noted that, if the historical optimal population positions corresponding to all the particles cannot be obtained initially, comparison is not needed, and the target optimal position is directly determined as the current optimal population position.
Step S35: determining the optimal position of each particle according to the optimal position of the current population;
specifically, the optimal position of each particle is determined according to the current population optimal position, where a specific optimization process may directly use the current population optimal position as the optimal position of the particle, or may additionally set an optimization speed, and gradually adjust the optimal position to the current population optimal position according to the optimization speed, which is not specifically limited in this embodiment.
Step S36: judging whether the optimized position meets the structural constraint condition or not based on the preset position range;
specifically, based on a preset position range, whether the optimized position meets a structural constraint condition is judged, so as to analyze whether the optimized position of each particle fluctuates within a reasonable range.
Step S37: if yes, judging whether the performance parameter of each particle at the optimized position meets the performance constraint condition; if so, taking the optimized position as a target structure parameter;
specifically, if the optimized position meets the structural constraint condition, it indicates that the optimized position of the particle is still within a reasonable range, and further, calculates the performance parameter of each particle at the optimized position, if the performance parameter meets the performance constraint condition, it indicates that the performance parameter of the particle at the optimized position has reached the target, and at this time, the optimized position may be directly used as the target structural parameter.
Step S38: if not, updating the optimized position according to the preset position range, and executing a step of judging whether the performance parameter of each particle on the optimized position meets a performance constraint condition.
Specifically, if the optimized position does not satisfy the structural constraint condition, it indicates that the optimized position of the particle is not within the reasonable range, and at this time, the optimized position needs to be updated directly according to the preset position range.
In one technical solution provided in this embodiment, the constraint condition is refined into a structural constraint condition and a performance constraint condition, and a specific step of performing optimization by using a particle swarm algorithm is provided, so that it can be ensured that the accurate target structural parameter of the target coil is obtained after the structural parameter, the performance parameter, and the like are input into the particle swarm algorithm.
Further, a fourth embodiment of the method of designing a magnetic resonance radio frequency coil of the present invention is presented. Based on the third embodiment, the step of determining the optimal position of each particle according to the optimal position of the current population includes:
step S41: acquiring the particle speed of each particle and a preset speed range corresponding to each particle, and constraining the particle speed in the preset speed range;
specifically, the particle velocity of each particle and a preset velocity range corresponding to each particle are obtained, where the particle velocity and the preset velocity range are preset by a technician according to actual conditions. Further, the particle velocity is constrained according to a predetermined velocity range, which is specifically that if the particle velocity is greater than the upper limit of the velocity range, the original value of the particle velocity is replaced by an upper limit, and if the particle velocity is less than the lower limit of the velocity range, the original value of the particle velocity is replaced by a lower limit, so as to ensure that the particle velocity is kept within a reasonable range.
Step S42: and updating the particle position as an optimized position based on the particle speed and the current population optimal position.
Specifically, on the premise that the particle speed is within a reasonable range, the particle position is updated based on the particle speed and the current optimal position of the population, and the updated particle position is used as the optimal position.
In a technical solution provided in this embodiment, a particle speed and a preset speed range of each particle are obtained, the particle speed is constrained according to the preset speed range, and then a particle position is updated as an optimized position according to the particle speed and a current optimal position of a population. According to the technical scheme provided by the embodiment, the optimal position of each particle is adjusted by increasing the particle speed on the basis of the optimal position of the current population, so that each particle can be gradually adjusted to the optimal position of the current population, large-amplitude error adjustment caused by error of the optimal position of the current population is reduced as much as possible, and the optimal design efficiency of the coil is improved.
Further, a fifth embodiment of the method of designing a magnetic resonance radio frequency coil of the present invention is presented. Based on the third embodiment, the step of calculating the fitness value of each particle through a fitness function according to the position of the particle of each particle includes:
step S51: performing electromagnetic simulation through a fitness function based on the COMSOL software and the particle position of each particle to obtain a target simulation result;
step S52: and determining the fitness value of each particle based on the MATLAB scientific calculation analysis software and the target simulation result.
Specifically, the COMSOL software has good simulation capability, so that a target simulation result is obtained by performing electromagnetic simulation through a fitness function based on the COMSOL software and the position of each particle. And because the MATLAB software has good data analysis capability, the adaptability value of each particle is determined based on the MATLAB software and the target simulation result.
In one technical solution provided in this embodiment, electromagnetic simulation is performed through a fitness function based on COMSOL software and a particle position of each particle to obtain a target simulation result, and then a fitness value of each particle is determined based on MATLAB software and the target simulation result. Compared with the method that only MATLAB software is used for obtaining a pure numerical value or COMSOL software is used for obtaining a pure simulation display, the technical scheme provided by the embodiment fully utilizes the simulation capability of the COMSOL software and the data analysis capability of the MATLAB software, combines the two software to complete the calculation of the fitness value, not only can the accuracy of the result be ensured, but also the result can be presented in a three-dimensional model mode.
Furthermore, a sixth embodiment of the design method of the magnetic resonance radio frequency coil is provided. Based on the first embodiment shown in fig. 2, the step of obtaining the constraint condition corresponding to the imaging target, and optimizing the structural parameter based on MATLAB scientific computational analysis software, the structural parameter, the performance parameter, and the constraint condition to obtain the target structural parameter of the target coil includes:
step S61: acquiring constraint conditions corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions;
specifically, the constraint conditions corresponding to the imaging target are obtained, and the constraint conditions are further set by a technician according to the imaging target, such as the magnetic field uniformity is 90%, and the minimum magnetic field strength is greater than 1A/m.
Further, continuously optimizing and iterating the structural parameters based on MATLAB software, the structural parameters, the performance parameters and the constraint conditions, and determining the optimal structural parameters as the optimization result.
Step S62: performing electromagnetic simulation based on the COMSOL software and the optimized structural parameters to obtain optimized performance parameters;
specifically, electromagnetic simulation is performed on the optimized structural parameters in COMSOL software, and corresponding optimized performance parameters are determined according to simulation results.
Step S63: and if the optimized performance parameters meet the performance parameter targets corresponding to the imaging targets, taking the optimized structure parameters as the target structure parameters of the target coil.
Further, if the optimized performance parameters meet the performance parameter targets corresponding to the imaging targets, if the performance parameter targets are that the magnetic field uniformity is 90% and the minimum magnetic field strength is greater than 1A/m, and the magnetic field uniformity of the optimized performance parameters is 95% and the minimum magnetic field strength is 1.2A/m, the performance parameter targets are met, which indicates that the optimized structural parameters reach the imaging targets at the moment, and the optimized structural parameters can also play a better role in subsequent practical application, so the optimized structural parameters can be directly used as the target structural parameters of the target coil.
In a technical solution provided in this embodiment, after the structural parameters are optimized, electromagnetic simulation is performed based on the COMSOL software and the optimized structural parameters to obtain optimized performance parameters, and if the optimized performance parameters meet a performance parameter target corresponding to an imaging target, the optimized structural parameters are used as target structural parameters of a target coil.
Further, with reference to fig. 4, a seventh embodiment of the magnetic resonance radio frequency coil design method of the present invention is presented. Based on the sixth embodiment, after the step of performing electromagnetic simulation based on the COMSOL software and the optimized structural parameters to obtain the optimized performance parameters, the method includes:
step S71: if the optimized structural parameters do not meet the performance parameter target corresponding to the imaging target, the COMSOL multi-physical field simulation-based software and the structural parameters are repeatedly executed, the target coil is simulated to obtain the performance parameters of the target coil until the constraint conditions corresponding to the imaging target are obtained, and the structural parameters are optimized to obtain the target structural parameters of the target coil based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions.
Specifically, after the structural parameters of the coil are optimized, the performance parameters under the optimized structural parameters are further analyzed, and if the optimized structural parameters do not meet the performance parameter target corresponding to the imaging target, it is indicated that the performance of the coil does not reach the standard after the optimization of this time, and the coil cannot be directly put into practical application, but the structure of the coil should be optimized again based on the data optimized last time. Further, repeatedly executing the COMSOL-based multi-physical field simulation software and the structural parameters, simulating the target coil to obtain the performance parameters of the target coil until obtaining the constraint conditions corresponding to the imaging target, optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions to obtain the target structural parameters of the target coil, and determining the target structural parameters of the target coil until the optimized structural parameters meet the performance parameter target.
In a technical solution provided in this embodiment, if the optimized structural parameter does not meet a performance parameter target corresponding to an imaging target, the COMSOL-based multi-physical field simulation software and the structural parameter are repeatedly executed to simulate the target coil to obtain the performance parameter of the target coil until the constraint condition corresponding to the imaging target is obtained, and the structural parameter is optimized based on MATLAB scientific calculation analysis software, the structural parameter, the performance parameter, and the constraint condition to obtain the target structural parameter of the target coil. According to the technical scheme, the processing mode under the condition that the optimized structural parameters do not meet the performance parameter target is provided, and the magnetic induction coil is optimized in a closed loop mode through continuous circulation of coil structure optimization, electromagnetic simulation and condition judgment until the appropriate structural parameters are output and serve as the target structural parameters of the target coil, so that the design efficiency and accuracy of the coil are effectively improved.
An embodiment of the present invention provides a magnetic resonance radio frequency coil design device, and referring to fig. 5, a functional module diagram of an embodiment of the magnetic resonance radio frequency coil design device of the present invention is shown.
The acquisition module is used for constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target and acquiring the structural parameters of the target coil based on the initial geometric model;
the simulation module is used for simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
and the optimization module is used for acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil.
Since the embodiments of the apparatus portion and the method portion correspond to each other, please refer to the description of the embodiments of the method portion for the embodiments of the apparatus portion, which is not repeated herein.
The embodiment of the invention provides a magnetic resonance radio frequency coil design device, which comprises: a memory, a processor and a computer program stored on the memory and executable on the processor, the computer program being configured to implement the steps of any of the embodiments of the magnetic resonance radio frequency coil design method.
Since the embodiment of the magnetic resonance radio frequency coil design apparatus portion and the embodiment of the method portion correspond to each other, please refer to the description of the embodiment of the method portion for the embodiment of the magnetic resonance radio frequency coil design apparatus portion, which is not repeated herein.
Embodiments of the present invention provide a readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps in any of the above-mentioned magnetic resonance radio frequency coil design methods.
Since the embodiment of the readable storage medium portion corresponds to the embodiment of the method portion, please refer to the description of the embodiment of the method portion for the embodiment of the readable storage medium portion, and details are not repeated here again.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or system. Without further limitation, an element defined by the phrases "comprising one of 8230; \8230;" 8230; "does not exclude the presence of additional like elements in a process, method, article, or system that comprises the element.
The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments.
Through the description of the foregoing embodiments, it is clear to those skilled in the art that the method of the foregoing embodiments may be implemented by software plus a necessary general hardware platform, and certainly may also be implemented by hardware, but in many cases, the former is a better implementation. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium (e.g., ROM/RAM, magnetic disk, optical disk) as described above and includes instructions for enabling a terminal device (e.g., a mobile phone, a computer, a server, or a network device) to execute the method according to the embodiments of the present invention.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.

Claims (10)

1. A magnetic resonance radio frequency coil design method is characterized by comprising the following steps:
constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target, and acquiring the structural parameters of the target coil based on the initial geometric model;
simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
and acquiring constraint conditions corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions to obtain target structural parameters of the target coil.
2. The magnetic resonance radio frequency coil design method as set forth in claim 1, wherein the step of optimizing the structural parameters based on MATLAB scientific computational analysis software, the structural parameters, the performance parameters, and the constraints to obtain target structural parameters of the target coil comprises:
setting a population scale and a spatial dimension according to the structural parameters, and respectively corresponding the structural parameters to each particle in the spatial dimension;
optimizing the structural parameters through a particle swarm algorithm based on the MATLAB scientific calculation analysis software, the structural parameters, the performance parameters, the population scale, the space dimensions and the constraint conditions to obtain target structural parameters of the target coil.
3. The magnetic resonance radio frequency coil design method according to claim 2, wherein the constraints include structural constraints and performance constraints, and the step of optimizing the structural parameters by a particle swarm algorithm based on the MATLAB scientific computational analysis software, the structural parameters, the performance parameters, the population size, the spatial dimensions, and the constraints to obtain target structural parameters of the target coil comprises:
determining the particle position of each particle based on the MATLAB scientific calculation analysis software, the structure parameters and the performance parameters, and acquiring a preset position range of each particle;
calculating the fitness value of each particle through a fitness function according to the particle position of each particle;
for each particle, traversing and acquiring a historical optimal position corresponding to each particle, comparing the fitness value of the particle position with the fitness value of the historical optimal position, and determining the particle position with a large fitness value or the historical optimal position as the optimal position of the current particle;
selecting a target optimal position with the largest fitness value from the particle positions aiming at all the particles, acquiring historical population optimal positions corresponding to all the particles, comparing the fitness value of the target optimal position with the fitness value of the historical population optimal position, and determining the target optimal position with the largest fitness value or the historical population optimal position as the current population optimal position;
determining the optimal position of each particle according to the optimal position of the current population;
judging whether the optimized position meets the structural constraint condition or not based on the preset position range;
if yes, judging whether the performance parameter of each particle at the optimized position meets the performance constraint condition; if so, taking the optimized position as a target structure parameter;
if not, updating the optimized position according to the preset position range, and executing a step of judging whether the performance parameter of each particle on the optimized position meets a performance constraint condition.
4. The method of claim 3, wherein the step of determining the optimal location of each particle based on the current population optimal location comprises:
obtaining the particle speed of each particle and a preset speed range corresponding to each particle, and constraining the particle speed in the preset speed range;
and updating the particle position as an optimized position based on the particle speed and the current population optimal position.
5. The magnetic resonance radio frequency coil design method as set forth in claim 3, wherein the step of calculating the fitness value of each particle by a fitness function based on the particle position of each particle comprises:
performing electromagnetic simulation through a fitness function based on the COMSOL software and the particle position of each particle to obtain a target simulation result;
and determining the fitness value of each particle based on the MATLAB scientific calculation analysis software and the target simulation result.
6. The method of claim 1, wherein the step of obtaining constraints corresponding to the imaging target and optimizing the structural parameters based on MATLAB scientific computational analysis software, the structural parameters, the performance parameters, and the constraints to obtain target structural parameters of the target coil comprises:
acquiring constraint conditions corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions;
performing electromagnetic simulation based on the COMSOL software and the optimized structural parameters to obtain optimized performance parameters;
and if the optimized performance parameters meet the performance parameter targets corresponding to the imaging targets, taking the optimized structure parameters as the target structure parameters of the target coils.
7. The method of claim 6, wherein said step of performing electromagnetic simulation based on said COMSOL software and optimized structural parameters to obtain optimized performance parameters comprises, after said step of:
if the optimized structural parameters do not meet the performance parameter targets corresponding to the imaging targets, the COMSOL multi-physical-field simulation-based software and the structural parameters are executed repeatedly, the target coil is simulated to obtain the performance parameters of the target coil until the constraint conditions corresponding to the imaging targets are obtained, and the structural parameters are optimized based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint conditions to obtain the target structural parameters of the target coil.
8. An apparatus for magnetic resonance radio frequency coil design, the apparatus comprising:
the acquisition module is used for constructing an initial geometric model of a target coil based on the main magnetic field direction of a magnetic resonance system and the design requirement of an imaging target and acquiring the structural parameters of the target coil based on the initial geometric model;
the simulation module is used for simulating the target coil based on COMSOL multi-physical field simulation software and the structural parameters to obtain performance parameters of the target coil;
and the optimization module is used for acquiring a constraint condition corresponding to the imaging target, and optimizing the structural parameters based on MATLAB scientific calculation analysis software, the structural parameters, the performance parameters and the constraint condition to obtain target structural parameters of the target coil.
9. A magnetic resonance radio frequency coil design apparatus, the apparatus comprising: memory, a processor and a computer program stored on the memory and executable on the processor, the computer program being configured to implement the steps of the magnetic resonance radio frequency coil design method as claimed in any one of claims 1 to 7.
10. A readable storage medium, having stored thereon a computer program which, when being executed by a processor, carries out the steps of the magnetic resonance radio frequency coil design method as set forth in any one of claims 1 to 7.
CN202211702466.7A 2022-12-29 2022-12-29 Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium Pending CN115659767A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211702466.7A CN115659767A (en) 2022-12-29 2022-12-29 Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211702466.7A CN115659767A (en) 2022-12-29 2022-12-29 Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium

Publications (1)

Publication Number Publication Date
CN115659767A true CN115659767A (en) 2023-01-31

Family

ID=85022427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211702466.7A Pending CN115659767A (en) 2022-12-29 2022-12-29 Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium

Country Status (1)

Country Link
CN (1) CN115659767A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116579284A (en) * 2023-07-12 2023-08-11 深圳飞骧科技股份有限公司 Method, system and related equipment for designing frequency band of surface acoustic wave filter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897487A (en) * 2017-01-13 2017-06-27 华中科技大学 A kind of modeling method of high-temperature superconducting magnet nonlinear analysis
CN111695292A (en) * 2020-05-15 2020-09-22 东南大学 Superconducting magnetic energy storage magnet optimization design method based on MATLAB and COMSOL joint simulation
CN111695275A (en) * 2020-05-11 2020-09-22 东南大学 Method for improving optimal calculation speed of superconducting magnetic energy storage coil
CN115526141A (en) * 2022-08-19 2022-12-27 西安电子科技大学芜湖研究院 Thermal stress coupling optimization design method for coaxial through-silicon via

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897487A (en) * 2017-01-13 2017-06-27 华中科技大学 A kind of modeling method of high-temperature superconducting magnet nonlinear analysis
CN111695275A (en) * 2020-05-11 2020-09-22 东南大学 Method for improving optimal calculation speed of superconducting magnetic energy storage coil
CN111695292A (en) * 2020-05-15 2020-09-22 东南大学 Superconducting magnetic energy storage magnet optimization design method based on MATLAB and COMSOL joint simulation
CN115526141A (en) * 2022-08-19 2022-12-27 西安电子科技大学芜湖研究院 Thermal stress coupling optimization design method for coaxial through-silicon via

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
尹斯星, 北京邮电大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116579284A (en) * 2023-07-12 2023-08-11 深圳飞骧科技股份有限公司 Method, system and related equipment for designing frequency band of surface acoustic wave filter
CN116579284B (en) * 2023-07-12 2024-04-12 深圳飞骧科技股份有限公司 Method, system and related equipment for designing frequency band of surface acoustic wave filter

Similar Documents

Publication Publication Date Title
CN109255830B (en) Three-dimensional face reconstruction method and device
US7657071B2 (en) Systems, devices, and methods for bundle segmentation in diffusion tensor magnetic resonance imaging
Jeng et al. Simulated annealing in compound Gaussian random fields (image processing)
CN111339695B (en) Apparatus and method for electromagnetic field simulation
CN107665736A (en) Method and apparatus for generating information
JP2022519678A (en) Methods and equipment for automated target and tissue segmentation using multimodal imaging and ensemble machine learning models
Chen et al. MRI brain tissue classification using unsupervised optimized extenics-based methods
Subramanian et al. Simulation of glioblastoma growth using a 3D multispecies tumor model with mass effect
CN109410253B (en) For generating method, apparatus, electronic equipment and the computer-readable medium of information
CN115659767A (en) Magnetic resonance radio frequency coil design method, device, equipment and readable storage medium
CN111028884B (en) Filling method, device and server for genotype data deletion
CN111179247A (en) Three-dimensional target detection method, training method of model thereof, and related device and equipment
AU744253B2 (en) A 3-D electromagnetic infinite element
CN101310302A (en) A method, a system and a computer program for volumetric registration
Zhou et al. Learning stochastic object models from medical imaging measurements by use of advanced ambient generative adversarial networks
CN117392119B (en) Tumor lesion area detection method and device based on position priori and feature perception
CN118133403A (en) City planning design drawing generation method, device, equipment, medium and product
CN112101438B (en) Left-right eye classification method, device, server and storage medium
CN109124635B (en) Model generation method, magnetic resonance imaging scanning method and system
CN116403210A (en) Dinner plate information identification method, dinner plate information identification device, electronic equipment and computer readable medium
CN108256451B (en) Method and device for detecting human face
CN114066986B (en) Three-dimensional coordinate determination method and device, electronic equipment and storage medium
CN116167320A (en) Automatic layout method and device for symmetrical constraint of analog integrated circuit
Winter et al. Open‐source magnetic resonance imaging: Improving access, science, and education through global collaboration
CN112037886B (en) Radiotherapy plan making device, method and storage medium

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20230131