A method for functionally specific target localization of individual rTMS without navigation

Abstract

The present invention relates to a method for positioning individual rTMS functional specific targets without navigation, and belongs to the field of target positioning. The method comprises S1, calculating the position of the cortical target; S2, segmenting the brain tissue of the standard brain template to obtain a standard and complete skull and scalp image; S3, outlining the edge of the skull and scalp image to obtain an edge image; S4, outlining the outermost edge on the edge image to obtain a standard scalp mask; S5, transferring the standard scalp mask to the individual space, and then registering with the individual image to obtain an individual scalp mask; S6, mapping the cortical target obtained in step S1 to the individual scalp mask to obtain a scalp target; S7, creating a coordinate system on the individual scalp mask, reading the coordinates of the scalp target, the highest points of the two auricles, the root of the nose and the occipital protuberance; S8, determining the position of the scalp target according to the positional relationship between the scalp target and the highest points of the two auricles, the root of the nose and the occipital protuberance. The present invention does not need to adopt a professional navigation and positioning system, and has low implementation cost.

Description

Technical Field

[0001] The invention belongs to the field of target location, and in particular to a method for non-navigation individual rTMS function-specific target location. Background Art

[0002] In the scientific research and clinical application of TMS (transcranial magnetic stimulation) and rTMS (repetitive transcranial magnetic stimulation), the selection and positioning of stimulation targets are crucial and closely related to the efficacy of TMS and rTMS. Targets calculated based on fMRI (functional magnetic resonance imaging) are targeted at specific functions, which can accelerate the effectiveness and improve the efficacy. In institutions with conditions, after the cortical targets are calculated, the navigation and positioning system can be used to accurately locate them on the subject's scalp. However, not all clinical and medical institutions are equipped with navigation systems. In the absence of navigation, how to locate the cortical targets defined by fMRI, that is, functionally specific targets, directly on the scalp is a problem that needs to be solved. Summary of the invention

[0003] The technical problem to be solved by the present invention is to provide a method for positioning individual rTMS function-specific targets without navigation, which can also perform rTMS function-specific target positioning on the scalp without a navigation positioning system.

[0004] In order to solve the above problems, the technical solution adopted by the present invention is: a method for positioning individual rTMS functional specific targets without navigation, comprising the following steps:

[0005] S1. Calculate the location of cortical targets using fMRI.

[0006] S2. Segment the brain tissue of the standard brain template to obtain a standard and complete skull and scalp image;

[0007] S3, outlining the edges of the skull and scalp images to obtain edge images;

[0008] S4, outlining the outermost edge on the edge image to obtain a standard scalp mask;

[0009] S5, converting the standard scalp mask to the individual space, and then registering it with the individual image to obtain the individual scalp mask;

[0010] S6, mapping the cortical target point obtained in step S1 onto the individual scalp mask to obtain the scalp target point;

[0011] S7. Create a coordinate system on the individual scalp mask, with the line connecting the occipital protuberance and the root of the nose as the x-axis, the line connecting the highest points of the two auricles as the y-axis, the intersection of the x-axis and the y-axis as the origin of the coordinates, the axis passing through the origin of the coordinates and perpendicular to both the x-axis and the y-axis as the z-axis, and the intersection of the z-axis and the scalp as the scalp landmark point; then read out the coordinates of the scalp target point, the highest points of the two auricles, the root of the nose, and the occipital protuberance;

[0012] S8. Determine the position of the scalp target point on the individual's head based on the positional relationship between the scalp target point and the highest points of the two auricles, the root of the nose, and the occipital protuberance.

[0013] Further, step S1 includes:

[0014] S11. Identify functionally specific targets in the primary motor cortex of the brain;

[0015] S12. rTMS intervention on functionally specific targets:

[0016] S121, preprocessing the task state data:

[0017] Correction of inter-layer acquisition time difference;

[0018] Rearrange the rigid bodies between layers to correct for head motion;

[0019] Align the functional image to the T1 structural image and register the T1 weighted image with the functional image;

[0020] The T1 structural image is segmented, and an inverse matrix is ​​generated during the segmentation process;

[0021] The functional images were spatially smoothed with a 6 mm half-width Gaussian smoothing kernel; covariates that affected the analysis results were removed;

[0022] S122. Calculate individual task activation points using a dynamic causal model: Set Units for design, Interscan interval, Microtime resolution, Microtime onset, Multiple regressors, and High-pass filter in the dynamic causal model to generate a design matrix;

[0023] The design matrix is ​​estimated using the estimation module in the dynamic causal model;

[0024] The estimated design matrix was statistically tested to obtain the individual task activation areas;

[0025] Find the individual task activation point in the individual task activation area, and the individual task activation point serves as the cortical target.

[0026] Further, step S3 includes:

[0027] S31, calculating the gradient and gradient change rate of the skull and scalp 3D images;

[0028] S32, calculating the intensity of the gradient;

[0029] S33, extracting the parts with gradient intensity in the top 5% as possible edge points;

[0030] S34. Among the possible edge points, the points whose edge strength exceeds the threshold are marked as 1, and the points whose edge strength does not exceed the threshold are marked as 0, so as to generate a binary image.

[0031] Further, step S4 includes:

[0032] S41, performing connected region detection on the binary image obtained in step S34, finding all connected regions and storing the information in a structure;

[0033] S42, calculating the number of all pixels in each connected area;

[0034] S43, finding the connected region containing the most pixels;

[0035] S44. Create an all-zero matrix according to the index of the connected region containing the most pixels. The connected region is stored in the all-zero matrix to obtain a standard scalp mask.

[0036] Furthermore, in step S5, the standard scalp mask is transferred to the individual space by using the inverse matrix generated by segmenting the T1 structural image in S121.

[0037] Furthermore, in step S2, the standard brain template adopts the MNI standard brain template, and step S6 includes:

[0038] S61, reading all non-zero voxels in the individual scalp mask, and converting the non-zero voxels into three-dimensional coordinates and storing them;

[0039] S62, converting the MNI coordinates into image matrix coordinates;

[0040] S63, calculating the distance between the cortical target point and all edge points of the individual scalp mask;

[0041] S64, find the nearest edge point, which is the mapped scalp target point, and convert the image matrix coordinates of the edge point into MNI coordinates.

[0042] Further, in step S8, a measuring tool is used to determine the position of the scalp target point on the individual's head, the measuring tool comprising a handle, the lower end of the handle is provided with an angle measuring disk, two arc-shaped fixed arc rulers and an arc-shaped movable arc ruler, the upper end of the fixed arc ruler is fixed to the angle measuring disk, and the two arc-shaped fixed arc rulers are symmetrically arranged relative to the angle measuring disk; the upper end of the movable arc ruler is slidably matched with the angle measuring disk, so that the movable arc ruler can rotate around the center line of the angle measuring disk; the fixed arc ruler and the movable arc ruler are provided with scale lines and readings, and the readings are arc distances from the corresponding scale lines to the center of the angle measuring disk;

[0043] The process of determining the location of the scalp target point on an individual's head is:

[0044] S81, selecting two reference points, the reference points being the highest points of the two auricles, or the reference points being the nasion and the occipital protuberance, and respectively calculating arc distances from the scalp target point and the two reference points to the scalp landmark point, wherein the arc is the shortest connecting line on the individual scalp mask;

[0045] S82, connecting the scalp landmark point and the scalp target point with a straight line, projecting the straight line onto the xy plane to obtain a projection line, and calculating an angle between the projection line and the reference axis;

[0046] S83. Use two fixed arc rulers of the measuring tool to clamp the individual's head. The two fixed arc rulers pass through the reference point respectively. According to the arc distance from the reference point to the scalp landmark point, adjust the position of the fixed arc ruler so that the reading of the fixed arc ruler at the reference point is the same as the arc distance from the reference point to the scalp landmark point. At this time, the center of the angle measurement disk is at the scalp landmark point; rotate the movable arc ruler so that the angle between the movable arc ruler and the reference axis is the angle calculated in step S82; then, according to the arc distance from the scalp target point to the scalp landmark point, find a point on the movable arc ruler where the reading is equal to the arc distance. This point is the individual scalp target point.

[0047] Furthermore, in S81, the calculation process of the arc distance is: given the coordinates of any two points A and B on the individual scalp mask, assuming that the path from point A to point B generates multiple path points through linear interpolation, the Euclidean distance between two adjacent path points is calculated and accumulated to obtain the actual arc distance.

[0048] The beneficial effect of the present invention is that since the cortical target point calculated by fMRI is located in the skull and cannot be directly observed or measured, the present invention adopts a vertical mapping method to map the cortical target point to the scalp surface. After being mapped to the scalp, it can only be seen on the magnetic resonance image, and it is still impossible to specifically locate this point on the scalp. Therefore, a coordinate system is established on the scalp, and the position mapped to the scalp in the image is implemented on the subject's scalp through measurement.

[0049] The present invention does not require the use of a professional navigation and positioning system and has low implementation cost. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 is a flow chart of the present invention;

[0051] Figure 2 is a schematic diagram of a standard brain template;

[0052] Figure 3 It is a standard and complete skull and scalp image;

[0053] Figure 4 is a schematic diagram of an edge image;

[0054] Figure 5 is a schematic diagram of a standard scalp mask;

[0055] Figure 6 is a schematic diagram of the standard scalp mask before registration with the individual image;

[0056] Figure 7 is a schematic diagram of the standard scalp mask after registration with the individual image;

[0057] Figure 8 is a side view schematic diagram of creating a coordinate system on an individual scalp mask;

[0058] Fig. 9 It is a top-down schematic diagram of creating a coordinate system on an individual scalp mask;

[0059] Fig.10 It is a schematic diagram of the measuring tool;

[0060] Fig.11 is a top view schematic diagram of the angle measuring disc;

[0061] Fig.12 is a schematic diagram of calculating the angle between the projection line and the reference axis in step S82;

[0062] Figure numerals: 1—handle; 2—angle measuring disc; 3—fixed arc ruler; 4—movable arc ruler. DETAILED DESCRIPTION

[0063] The present invention is further described below in conjunction with the accompanying drawings and embodiments.

[0064] The method of the present invention for the non-navigated individual rTMS function-specific target localization is as follows: Figure 1 As shown, the following steps are included:

[0065] S1. Calculate the location of cortical targets using fMRI.

[0066] Specifically include:

[0067] S11. Identify functionally specific targets in the primary motor cortex of the brain.

[0068] The cortical targets are defined by functional-specific tasks (i.e., functional-specific targets). The functional-specific tasks can be manual tasks of the right thumb, or other fMRI indicators such as ALFF (low-frequency amplitude), FC (functional connectivity), or other existing computational methods to define functional-specific targets on the cortex.

[0069] S12. After calculating the location of the functional specific target, perform rTMS intervention on the functional specific target:

[0070] S121, preprocessing the task state data, where the task state data is the obtained magnetic resonance imaging machine image, and the specific contents include:

[0071] Correct the acquisition time difference between layers. Since the scanning is done by layer-by-layer scanning, there is a gap between the scanning times of different layers in the same brain region. Therefore, the scanning times of different layers in the same brain region need to be corrected to the same time point.

[0072] Rearrange the rigid bodies between layers to correct head movement. During the scanning process, it is difficult for the individual (the person being scanned) to keep his head completely fixed and may shake slightly. Correction is needed to reduce the error caused by head movement. For the sake of correction, head movement is considered to be less than 2.5 degrees of rotation in any direction and less than 2.5 mm of translation.

[0073] The functional image is aligned to the T1 structural image, and the T1 weighted image is registered with the functional image. The scanned images include the functional image, T1 structural image, and T1 weighted image, etc. After registration, it is convenient for joint analysis.

[0074] The T1 structural image is segmented, and an inverse matrix is ​​generated during the segmentation process, which can be used to convert the standard scalp mask into the individual space.

[0075] The functional images were spatially smoothed with a Gaussian smoothing kernel of 6 mm at half maximum width to reduce noise.

[0076] Remove covariates that affect the analysis results.

[0077] Step S121 can be performed in the DPABI software using the software's built-in functions.

[0078] S122. Calculate individual task activation points using a dynamic causal model. The dynamic causal model is SPM12 software, and existing technology can be used to calculate individual task activation areas using the Specify 1st-level module in the SPM12 software to obtain function-specific targets.

[0079] Specifically, set Units for design, Interscan interval, Microtime resolution, Microtime onset, Multiple regressors, and High-pass filter in the dynamic causal model to generate a design matrix. Units for design, Interscan interval, Microtime resolution, Microtime onset, Multiple regressors, and High-pass filter are all submodules in the Specify 1st-level module. When operating, open each submodule to set parameters, and the software automatically generates a design matrix, which is the SPM.mat file.

[0080] The design matrix is ​​estimated using the estimation module in the dynamic causal model. The estimation module (Estimate) is also a submodule in the Specify 1st-level module. The SPM.mat file generated in the previous step is used for estimation. The software automatically completes the estimation.

[0081] Perform statistical tests on the estimated design matrix to obtain the individual task activation area. Statistical test (Result—contrast) is also a submodule in the Specify 1st-level module. Import the estimated SPM.mat file, select t-contrasts, and click done to generate the spmT_000x.nii file, which is the individual task activation area.

[0082] The individual task activation point is found in the individual task activation area using existing technology, and the individual task activation point is used as the cortical target.

[0083] In the case of navigation, the coordinates of the individual task activation points can be directly input into the navigation system and relevant operations can be performed to achieve positioning. In the case of no navigation, the data needs to be further processed.

[0084] S2. Segment the brain tissue of the standard brain template to obtain a standard and complete skull and scalp image.

[0085] The standard brain template uses the MNI standard brain template, such as Figure 2 Specifically, the segment function module in the SPM12 software is used to segment the standard brain template, and the obtained standard and complete skull and scalp images are shown in FIG. Figure 3 shown.

[0086] S3, outlining the edges of the skull and scalp images to obtain edge images. The edge images are as follows: Figure 4 shown.

[0087] Specifically, step S3 includes:

[0088] S31. Calculate the gradient and gradient change rate of the skull and scalp 3D images, and use the gradient function to calculate the gradient change rate.

[0089] S32. Calculate the strength of the gradient. Specifically, use the gradient_magnitude function to calculate the strength of the gradient.

[0090] S33. Extract the parts with gradient intensity in the top 5% as possible edge points.

[0091] S34. Among the possible edge points, the points whose edge strength exceeds the threshold are marked as 1, and the points whose edge strength does not exceed the threshold are marked as 0, so as to generate a binary image.

[0092] S4. Since the present invention only wants to obtain the outermost edge information, the outermost edge is outlined on the edge image to obtain a standard scalp mask.

[0093] Step S4 includes:

[0094] S41, performing connected region detection on the binary image obtained in step S34, finding all connected regions and storing the information in a structure.

[0095] S42, calculating the number of all pixels in each connected area;

[0096] S43, finding the connected region containing the most pixels;

[0097] S44. Create an all-zero matrix according to the index of the connected region containing the most pixels. The connected region is stored in the all-zero matrix to obtain a standard scalp mask. The standard scalp mask is as follows: Figure 5 shown.

[0098] S5. Convert the standard scalp mask to the individual space, and then register it with the individual image to obtain the individual scalp mask. Figure 6 As shown, after registration Figure 7 shown.

[0099] Specifically, in S121 , the T1 structural image is segmented to generate an inverse matrix to transfer the standard scalp mask to the individual space.

[0100] S6. Map the cortical target point obtained in step S1 onto the individual scalp mask to obtain the scalp target point.

[0101] Step S6 includes:

[0102] S61. Read all non-zero voxels in the individual scalp mask, where the non-zero voxels are scalp edge points, and convert the non-zero voxels into three-dimensional coordinates and store them.

[0103] S62, converting the MNI coordinates into image matrix coordinates. Since the standard brain template uses the MNI standard brain template, the coordinates of the non-zero voxels obtained on the basis of the MNI standard brain template are MNI coordinates, which usually represent anatomical positions and are coordinates of the standard space. The voxel values ​​of the image matrix of the image processed in step S1 correspond to the image matrix coordinates. Therefore, the MNI coordinates of the non-zero voxels are converted into image matrix coordinates to unify the coordinates for subsequent calculations.

[0104] S63. Calculate the distance between the cortical target point and all edge points of the individual scalp mask.

[0105] S64, find the edge point closest to the cortical target point, which is the mapped scalp target point, and convert the image matrix coordinates of this edge point into MNI coordinates. After the conversion into image matrix coordinates in step S62, the coordinates of the scalp target point found are also image matrix coordinates, and in subsequent positioning, MNI coordinates need to be used, so the image matrix coordinates are converted into MNI coordinates again.

[0106] Since the cortical target points calculated by fMRI are located inside the skull and cannot be directly observed or measured, the present invention uses a vertical mapping method to map the cortical target points to the scalp surface and obtain the scalp target points on the individual scalp mask. After being mapped to the scalp, it can only be seen on the magnetic resonance image, and it is still impossible to specifically locate this point on the human scalp. Therefore, a coordinate system is established on the scalp, and the points mapped to the scalp in the image are implemented on the subject's scalp through measurement.

[0107] Specifically, S7, create a coordinate system on the individual scalp mask: the line connecting the occipital protuberance and the nasion is the x-axis, the line connecting the highest points of the two auricles is the y-axis, the intersection of the x-axis and the y-axis is the coordinate origin, and the axis passing through the coordinate origin and perpendicular to the x-axis and the y-axis is the z-axis, such as Figure 8 and Fig. 9 As shown. It is necessary to determine the coordinates of the highest points of the two auricles, the nasal root and the occipital protuberance on the T1 structural image. The intersection of the z-axis and the scalp is the scalp landmark point. At this time, the coordinates of the scalp target point, the highest points of the two auricles, the nasal root and the occipital protuberance can be directly read.

[0108] S8. Determine the position of the scalp target point on the individual's head based on the positional relationship between the scalp target point and the highest points of the two auricles, the root of the nose, and the occipital protuberance.

[0109] Specifically, a measuring tool is used to determine the location of the scalp target point on the individual's head. The measuring tool is as follows: Fig.10As shown, it includes a handle 1, and the lower end of the handle 1 is provided with an angle measuring disc 2, two arc-shaped fixed arc rulers 3 and an arc-shaped movable arc ruler 4. The upper end of the fixed arc ruler 3 is fixed to the angle measuring disc 2, and the two arc-shaped fixed arc rulers 3 are symmetrically arranged relative to the angle measuring disc 2; the movable arc ruler 4 can be one or two. The upper end of the movable arc ruler 4 is slidably matched with the angle measuring disc 2, so that the movable arc ruler 4 can rotate around the center line of the angle measuring disc 2; the fixed arc ruler 3 and the movable arc ruler 4 are provided with scale lines and readings, and the readings are the arc distances from the corresponding scale lines to the center of the angle measuring disc 2. The fixed arc ruler 3 and the movable arc ruler 4 can be made of silicone rulers.

[0110] The angle measuring disc 2 can be an existing circular protractor. The angle measuring disc 2 is fixedly connected to the handle 1. The entire measuring tool can be moved by holding the handle 1. The measuring range of the angle measuring disc 2 is 0 to 360 degrees. Fig.11 As shown, when the movable arc ruler 4 rotates around the center line of the angle measuring disc 2, the rotation angle of the movable arc ruler 4 can be detected. The center of the angle measuring disc 2 simulates the scalp landmark point, and two straight lines passing through the center of the angle measuring disc 2 in the radial direction and perpendicular to each other can simulate the x-axis and the y-axis. For example, the projection of two fixed arc rulers 3 on the angle measuring disc 2 can simulate the x-axis or the y-axis. On the angle measuring disc 2, the angle between the projection of the movable arc ruler 4 on the angle measuring disc 2 and the x-axis or the y-axis can be read.

[0111] The process of determining the location of the scalp target point on an individual's head is:

[0112] S81. Select two reference points. The reference points can be the highest points of the two auricles, or the root of the nose and the occipital protuberance. Calculate the arc distances from the scalp target point and the two reference points to the scalp landmark point respectively. The arc is the shortest connecting line on the individual scalp mask. The shapes of the fixed arc ruler 3 and the movable arc ruler 4 are the same as the shape of the arc. In actual positioning, the fixed arc ruler 3 and the movable arc ruler 4 can be used to simulate the arc.

[0113] Specifically, the arc distance is calculated by linear interpolation: given the coordinates of any two points A and B on the individual scalp mask, assume that the path from point A to point B generates multiple path points through linear interpolation, calculate the Euclidean distance between two adjacent path points and accumulate them to obtain the actual arc distance. The path points are inserted with equal length, and the specific number of path points can be 100. The number of path points is determined according to the actual distance. The more path points there are, the smaller the calculation error.

[0114] S82, connect the scalp landmark point and the scalp target point with a straight line, project the straight line onto the xy plane to obtain a projection line, and calculate the angle between the projection line and the reference axis, such as Fig.12As shown. The reference axis can be the positive direction of the x-axis, the negative direction of the x-axis, the positive direction of the y-axis or the negative direction of the y-axis. According to the quadrant in which the projection line is located, the closer axis is selected as the reference axis. The projection line rotates the calculated angle toward the reference axis with the z-axis as the rotation center, that is, coincides with the reference axis. In actual positioning, the center line of the angle measurement disk 2 simulates the z-axis.

[0115] S83, using two fixed arc rulers 3 of the measuring tool to clamp the individual's head, the two fixed arc rulers 3 pass through reference points respectively, the two reference points are the highest points of the two auricles of the human head, or the two reference points are the root of the nose and the occipital protuberance of the human body. According to the arc distance from the reference point to the scalp landmark point calculated in step S81, the position of the fixed arc ruler 3 is adjusted so that the reading of the fixed arc ruler 3 at the reference point is the same as the arc distance from the reference point to the scalp landmark point. At this time, the center of the angle measuring disk 2 is at the scalp landmark point, the angle measuring disk 2 is used as the xy plane, the center line of the angle measuring disk 2 is the z axis, the projection of the fixed arc ruler 3 on the angle measuring disk 2 is the x axis or the y axis, and the positive or negative direction of the x axis or the y axis is used as the reference axis.

[0116] Then, the movable arc ruler 4 is rotated around the center line (z axis) of the angle measuring disk 2, and the angle between the movable arc ruler 4 and the reference axis is measured by the angle measuring disk 2, so that the angle between the projection of the movable arc ruler 4 on the angle measuring disk 2 and the reference axis is the angle calculated in step S82. Finally, according to the arc distance from the scalp target point to the scalp landmark point calculated in step S81, a point on the movable arc ruler 4 where the reading is equal to the arc distance is found, and this point is the individual scalp target point.

[0117] The present invention does not need to adopt a professional navigation and positioning system and has low implementation cost.

[0118] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for non-navigated individual rTMS functional specific target localization, characterized in that: The following steps are involved: S1. Calculate the location of cortical targets using fMRI; including: S11. Identify functionally specific targets in the primary motor cortex of the brain; S12. rTMS intervention on functionally specific targets: S121, preprocessing the task state data: Correction of inter-layer acquisition time difference; Rearrange the rigid bodies between layers to correct for head motion; Align the functional image to the T1 structural image and register the T1 weighted image with the functional image; The T1 structural image is segmented, and an inverse matrix is ​​generated during the segmentation process; The functional images were spatially smoothed with a 6 mm half-width Gaussian smoothing kernel; covariates that affected the analysis results were removed; S122. Calculate individual task activation points using a dynamic causal model: Set Units for design, Interscan interval, Microtime resolution, Microtime onset, Multiple regressors, and High-pass filter in the dynamic causal model to generate a design matrix; The design matrix is ​​estimated using the estimation module in the dynamic causal model; The estimated design matrix was statistically tested to obtain the individual task activation areas; Find the individual task activation point in the individual task activation area, and the individual task activation point will be used as the cortical target point; S2. Segment the brain tissue of the standard brain template to obtain a standard and complete skull and scalp image; S3, outlining the edges of the skull and scalp images to obtain edge images; S4, outlining the outermost edge on the edge image to obtain a standard scalp mask; S5, converting the standard scalp mask to the individual space, and then registering it with the individual image to obtain the individual scalp mask; S6, mapping the cortical target point obtained in step S1 onto the individual scalp mask to obtain the scalp target point; S7. Create a coordinate system on the individual scalp mask, with the line connecting the occipital protuberance and the root of the nose as the x-axis, the line connecting the highest points of the two auricles as the y-axis, the intersection of the x-axis and the y-axis as the origin of the coordinates, the axis passing through the origin of the coordinates and perpendicular to both the x-axis and the y-axis as the z-axis, and the intersection of the z-axis and the scalp as the scalp landmark point; then read out the coordinates of the scalp target point, the highest points of the two auricles, the root of the nose, and the occipital protuberance; S8. Determine the position of the scalp target point on the individual's head based on the positional relationship between the scalp target point and the highest points of the two auricles, the root of the nose, and the occipital protuberance.

2. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 1, characterized in that: Step S3 includes: S31, calculating the gradient and gradient change rate of the skull and scalp 3D images; S32, calculating the intensity of the gradient; S33, extracting the parts with gradient intensity in the top 5% as possible edge points; S34. Among the possible edge points, the points whose edge strength exceeds the threshold are marked as 1, and the points whose edge strength does not exceed the threshold are marked as 0, so as to generate a binary image.

3. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 2, characterized in that: Step S4 includes: S41, performing connected region detection on the binary image obtained in step S34, finding all connected regions and storing the information in a structure; S42, calculating the number of all pixels in each connected area; S43, finding the connected region containing the most pixels; S44. Create an all-zero matrix according to the index of the connected region containing the most pixels. The connected region is stored in the all-zero matrix to obtain a standard scalp mask.

4. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 1, characterized in that: In step S5, the standard scalp mask is transferred to the individual space by using the inverse matrix generated by segmenting the T1 structural image in step S121.

5. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 3, characterized in that: In step S2, the standard brain template adopts the MNI standard brain template, and step S6 includes: S61, reading all non-zero voxels in the individual scalp mask, and converting the non-zero voxels into three-dimensional coordinates and storing them; S62, converting the MNI coordinates into image matrix coordinates; S63, calculating the distance between the cortical target point and all edge points of the individual scalp mask; S64, find the nearest edge point, which is the mapped scalp target point, and convert the image matrix coordinates of the edge point into MNI coordinates.

6. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 1, characterized in that: In step S8, a measuring tool is used to determine the position of a scalp target point on the head of an individual, the measuring tool comprising a handle (1), the lower end of the handle (1) being provided with an angle measuring disc (2), two arc-shaped fixed arc rulers (3) and an arc-shaped movable arc ruler (4), the upper end of the fixed arc ruler (3) being fixed to the angle measuring disc (2), and the two arc-shaped fixed arc rulers (3) being symmetrically arranged relative to the angle measuring disc (2); the upper end of the movable arc ruler (4) being slidably matched with the angle measuring disc (2), so that the movable arc ruler (4) can rotate around the center line of the angle measuring disc (2); the fixed arc ruler (3) and the movable arc ruler (4) are provided with scale lines and readings, the readings being the arc distances from the corresponding scale lines to the center of the angle measuring disc (2); The process of determining the location of the scalp target point on an individual's head is: S81, selecting two reference points, the reference points being the highest points of the two auricles, or the reference points being the nasion and the occipital protuberance, and respectively calculating arc distances from the scalp target point and the two reference points to the scalp landmark point, wherein the arc is the shortest connecting line on the individual scalp mask; S82, connecting the scalp landmark point and the scalp target point with a straight line, projecting the straight line onto the xy plane to obtain a projection line, and calculating an angle between the projection line and the reference axis; S83, using two fixed arc rulers (3) of the measuring tool to clamp the head of the individual, the two fixed arc rulers (3) respectively pass through the reference point, and according to the arc distance from the reference point to the scalp mark point, adjust the position of the fixed arc ruler (3) so that the reading of the fixed arc ruler (3) at the reference point is the same as the arc distance from the reference point to the scalp mark point, at this time, the center of the angle measurement disk (2) is at the scalp mark point; rotate the movable arc ruler (4) so ​​that the angle between the movable arc ruler (4) and the reference axis is the angle calculated in step S82; then, according to the arc distance from the scalp target point to the scalp mark point, find a point on the movable arc ruler (4) where the reading is equal to the arc distance, and this point is the individual scalp target point.

7. The method for non-navigated individual rTMS functional specific target localization as claimed in claim 6, characterized in that: In S81, the calculation process of the arc distance is: given the coordinates of any two points A and B on the individual scalp mask, assuming that the path from point A to point B generates multiple path points through linear interpolation, the Euclidean distance between two adjacent path points is calculated and accumulated to obtain the actual arc distance.

Images