Image processing method, device, apparatus and storage medium

Abstract

The application discloses an image processing method, device and equipment and a storage medium. The method comprises the following steps: acquiring a head contour image, processing the head contour image, and obtaining a head contour edge image; determining a plurality of target feature points on the head contour edge image; wherein the plurality of target feature points comprise a nasal root feature point, an alar point feature point and an occipital protuberance feature point; determining a reference distance based on a first coordinate of the nasal root feature point and a second coordinate of the occipital protuberance feature point; determining a central sulcus position based on the reference distance and the first coordinate; determining a target space position based on the central sulcus position and a third coordinate of the alar point feature point; wherein the target space position comprises a central preoccipital space position and a central postoccipital space position. The technical scheme of the application automatically determines the central preoccipital space position and the central postoccipital position, and the result has high accuracy.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of medical images, and particularly relates to an image processing method and device, equipment and a storage medium. BACKGROUND

[0002] The central sulcus is one of the sulci of the brain that separates the frontal lobe from the parietal lobe in the coronal plane. It runs along the dorsolateral surface of the hemisphere from the precentral gyrus. It originates at the superior margin of the hemisphere and extends inferiorly to the lateral fissure. Above the sulcus is the precentral gyrus in the frontal lobe of the cerebral cortex, and below the sulcus is the postcentral gyrus in the parietal lobe of the cerebral cortex.

[0003] Magnetic resonance (MR) images are relatively clear, and doctors can lock the position of the central sulcus from the images based on years of experience and knowledge of brain anatomy.

[0004] However, the resolution of computed tomography (CT) images is low, and the sulci are not clear, so it is difficult for the human eye to distinguish the spatial positions of the central sulcus and the precentral and postcentral gyri. SUMMARY

[0005] The present application aims to at least partially solve one of the technical problems in the related art. To this end, one object of the present application is to provide an image processing method, device, equipment and storage medium.

[0006] To solve the above technical problems, embodiments of the present application provide the following technical solutions:

[0007] An image processing method comprises:

[0008] Obtaining a head contour image and processing the head contour image to obtain a head contour edge image;

[0009] Determining a plurality of target feature points on the head contour edge image; wherein the plurality of target feature points comprise a nasion feature point, an ala feature point and an inion protuberance feature point;

[0010] Determining a reference distance based on a first coordinate of the nasion feature point and a second coordinate of the inion protuberance feature point;

[0011] Determining a central sulcus position based on the reference distance and the first coordinate;

[0012] Determining a target spatial position based on the central sulcus position and a third coordinate of the ala feature point; wherein the target spatial position comprises a precentral gyrus spatial position and a postcentral gyrus spatial position.

[0013] Optionally, before the obtaining of the head contour edge image, further comprising:

[0014] based on the head contour image, obtaining an initial head contour edge image;

[0015] obtaining a first smoothing parameter and a second smoothing parameter;

[0016] the smoothing processor performs smoothing processing on the initial head contour edge image based on the first smoothing parameter and the second smoothing parameter to obtain the head contour edge image.

[0017] Optionally, the determining of the reference distance based on the first coordinate of the nasal root feature point and the second coordinate of the inion feature point comprises:

[0018] on the head contour edge image, along a preset direction, taking the nasal root feature point as a starting point and the inion feature point as an ending point, all edge pixel points between the nasal root feature point and the inion feature point are sequentially formed into a chain fourth coordinate point set;

[0019] based on the fourth coordinate point set and in combination with the spatial resolution information of the head contour edge image, the reference distance from the nasal root feature point to the inion feature point is obtained.

[0020] Optionally, the obtaining of the reference distance from the nasal root feature point to the inion feature point based on the fourth coordinate point set and in combination with the spatial resolution information of the head contour edge image comprises:

[0021] based on the fourth coordinate point set and the spatial resolution information of the head contour edge image, the distance between every two adjacent fourth coordinate points is accumulated and superimposed to obtain the reference distance from the nasal root feature point to the inion feature point.

[0022] Optionally, the obtaining of the reference distance from the nasal root feature point to the inion feature point by accumulating and superimposing the distance between every two adjacent fourth coordinate points comprises:

[0023] obtaining a first sub-coordinate of every two adjacent fourth coordinate points, and calculating a first term based on each first sub-coordinate point; obtaining a second sub-coordinate of every two adjacent fourth coordinate points, and calculating a second term based on each second sub-coordinate point; obtaining a third sub-coordinate of every two adjacent fourth coordinate points, and calculating a third term based on each third sub-coordinate point;

[0024] based on the first term, the second term and the third term corresponding to every two adjacent fourth coordinate points, a cumulative superimposed term corresponding to every two adjacent fourth coordinate points is respectively calculated.

[0025] calculating a cumulative superposition item corresponding to each two adjacent fourth coordinate points to obtain the reference distance between the nasion feature point and the inion feature point.

[0026] Optionally, the determining the central sulcus position based on the reference distance and the first coordinate includes:

[0027] determining a plurality of central sulcus points on the head contour edge image based on the reference distance and the first coordinate;

[0028] obtaining a fifth coordinate of each of the central sulcus points;

[0029] determining the central sulcus position based on the fifth coordinate of each of the central sulcus points.

[0030] Optionally, the determining the central sulcus position based on the reference distance and the first coordinate includes:

[0031] obtaining a preset distance;

[0032] determining a curve distance based on the preset distance and the reference distance;

[0033] determining a point on the head contour edge image as the central sulcus point, where a distance between the point and the nasion feature point is equal to the curve distance.

[0034] Optionally, the determining the target space position based on the central sulcus position and the third coordinate of the alar point feature point includes:

[0035] selecting a first plane point on an i-th slice of the head contour edge image and a second plane point on a j-th slice of the head contour edge image based on the central sulcus position, where the head contour edge image includes a plurality of slices, i and j are positive integers, and i≠j;

[0036] obtaining a first plane coordinate of the first plane point and a second plane coordinate of the second plane point;

[0037] determining a reference plane equation and a reference plane based on the first plane coordinate, the second plane coordinate, and the third coordinate;

[0038] determining a target space position as a region on the head contour edge image, where a distance between the region and the reference plane is less than or equal to a preset boundary distance, and the central preoccipital space position and the central postoccipital space position are symmetrical about the central sulcus position.

[0039] Embodiments of the present application also provide an image processing apparatus, comprising:

[0040] a processing module configured to acquire a head contour image and process the head contour image to obtain a head contour edge image;

[0041] a first determining module configured to determine a plurality of target feature points on the head contour edge image; wherein the plurality of target feature points comprise a nasion feature point, an ala feature point, and an inion feature point;

[0042] a second determining module configured to determine a reference distance based on a first coordinate of the nasion feature point and a second coordinate of the inion feature point;

[0043] a third determining module configured to determine a central sulcus position based on the reference distance and the first coordinate;

[0044] a fourth determining module configured to determine a target spatial position based on the central sulcus position and a third coordinate of the ala feature point; wherein the target spatial position comprises a precentral sulcus spatial position and a postcentral sulcus spatial position.

[0045] Embodiments of the present application also provide an electronic device, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor implements the method as described above when executing the computer program.

[0046] Embodiments of the present application also provide a computer readable storage medium, comprising a stored computer program, wherein the computer readable storage medium controls a device where the computer readable storage medium is located to execute the method as described above when the computer program runs.

[0047] Embodiments of the present application also provide an electronic device, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor implements the method as described above when executing the computer program.

[0048] Embodiments of the present application also provide a computer readable storage medium, comprising a stored computer program, wherein the computer readable storage medium controls a device where the computer readable storage medium is located to execute the method as described above when the computer program runs.

[0049] Embodiments of the present application have the following technical effects:

[0050] The above technical solutions of the present application automatically identify and determine the location of the central sulcus based on the coordinates of three feature points of the nasion feature point, the inion feature point and the ala feature point on the CT image sequence, and determine the location of the precuneus space and the location of the postcentral gyrus according to the location of the central sulcus, and the result has high accuracy, strong compatibility and fast search speed.

[0051] Additional aspects and advantages of the present application will be in part apparent and in part pointed out hereinafter. BRIEF DESCRIPTION OF DRAWINGS

[0052] Figure 1 is a schematic diagram of an existing CT image structure provided by an embodiment of the present application;

[0053] Figure 2 is a schematic diagram of a flow of an image processing method provided by an embodiment of the present application;

[0054] Figure 3 is a schematic diagram of a head contour image structure provided by an embodiment of the present application;

[0055] Figure 4 is a schematic diagram of a solid head contour image structure provided by an embodiment of the present application;

[0056] Figure 5 is a schematic diagram of a three-dimensional disc-shaped filter structure provided by an embodiment of the present application;

[0057] Figure 6 is a schematic diagram of a working interface of a three-dimensional disc-shaped filter provided by an embodiment of the present application;

[0058] Figure 7 is a schematic diagram of an initial head contour edge image structure provided by an embodiment of the present application;

[0059] Figure 8 is a schematic diagram of an initial head contour edge image structure before smoothing provided by an embodiment of the present application;

[0060] Figure 9 is a schematic diagram of a head contour edge image structure after smoothing provided by an embodiment of the present application;

[0061] Figure 10 is a schematic diagram of a fourth coordinate point determination provided by an embodiment of the present application;

[0062] Figure 11 is a schematic diagram of a central sulcus location determination provided by an embodiment of the present application;

[0063] Figure 12 is a schematic diagram of a target space location determination provided by an embodiment of the present application;

[0064] Figure 13 This is a schematic diagram of a target spatial location provided in an embodiment of this application;

[0065] Figure 14 This is a schematic diagram of another target spatial location provided in an embodiment of this application;

[0066] Figure 15 This is a schematic diagram of the structure of an image processing device provided in an embodiment of this application. Detailed Implementation

[0067] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0068] To facilitate understanding of the embodiments by those skilled in the art, some terms are explained below:

[0069] (1) ITK: Insight Segmentation and Registration Toolkit, which uses advanced multimodal data segmentation and registration algorithms to handle image registration and segmentation problems.

[0070] (2) The precentral gyrus contains a large number of giant pyramidal cells and is the location of the motor center. It controls the movement of the contralateral half of the body and is the origin of the corticospinal tract and the corticobulbar tract.

[0071] The postcentral gyrus, located between the central sulcus and the postcentral sulcus, is the somatosensory area of ​​the cerebral cortex, responsible for sensory functions of the contralateral half of the body (excluding the lower leg and foot).

[0072] The precentral and anterior central gyri are the most important functional areas in the brain besides the vital centers (brainstem). They are crucial areas for radiologists and clinicians to analyze patient conditions, make functional diagnoses, and assess prognosis. MRI images are relatively clear, and based on years of experience and knowledge of brain anatomy, doctors can identify features such as the Ω sign and inverted T sign from the images to pinpoint the location of the central sulcus.

[0073] However, as Figure 1 As shown, CT images have low resolution and unclear sulci and gyri, making it difficult for the human eye to distinguish the spatial location of the central sulcus and anterior and posterior gyri, resulting in low accuracy in determining the anterior and posterior central gyri.

[0074] To address the aforementioned technical problems, embodiments of this application provide the following technical solutions:

[0075] likeFigure 2 As shown, the embodiment of the present application provides an image processing method, comprising:

[0076] Step S21: obtaining a head contour image, and processing the head contour image to obtain a head contour edge image;

[0077] In an optional embodiment of the present application, in order to improve the accuracy of image processing, after obtaining the CT image, the CT image is first processed for three-dimensional correction; for example, using the ITK three-dimensional rigid registration tool to perform three-dimensional correction on the patient's brain.

[0078] In an optional embodiment of the present application, after the three-dimensional correction of the CT image is completed, the head contour of the corrected CT image is extracted to obtain the head contour as shown in Figure 3

[0079] Specifically, the head contour refers to the skin of the head and all the head contours wrapped by the skin V.

[0080] Since the CT value (HU) of air is about -1000, the algorithm of three-dimensional correction itself will normalize the pixel of the head accessory part region to the HU value of 0;

[0081] Therefore, V iSlice = 1; Vimage iSlice > 0;

[0082] In the formula, Vimage is the plain CT image; iSlice is the i-th tomographic image, for example, if the scanning layer thickness of the CT image is 1 mm, there are 165 tomographic images, and the value range of iSlice is {1, 2, …, 165}.

[0083] In an optional embodiment of the present application, as shown in Figure 4 A solid head contour image is obtained based on a morphological hole filling process.

[0084] In an optional embodiment of the present application, as shown in Figure 5 In order to maximize the guarantee of the curvature of the scalp surface without distortion, a three-dimensional disc filter B with a scale of 3*3*3 and a radius of 1 is used on V, and the epidermis edge is extracted on the xoz plane as shown in Figure 6 For example, taking the middle layer of the sagittal position as an example, the initial head contour edge image VC can be obtained as shown in Figure 7

[0085] VC = V-VΘB; wherein Θ is the related working algorithm of B, and the embodiments of the present application do not make specific limitations on this.

[0086] In an optional embodiment of the present application, before the head contour edge image is obtained, the method further comprises: ​​

[0087] Based on the head contour image, an initial head contour edge image is obtained;

[0088] Obtain the first smoothing parameter and the second smoothing parameter;

[0089] The smoothing processor performs smoothing processing on the initial head contour edge image based on the first smoothing parameter and the second smoothing parameter to obtain the head contour edge image.

[0090] In one optional embodiment of this application, since the central groove is located on the curved surface of the line connecting the nasal root feature point and the occipital protuberance feature point, the location of the central groove can be determined by calculating the absolute distance between the nasal root feature point and the occipital protuberance feature point.

[0091] like Figure 8 As shown, the realism of the curved surface directly determines the accuracy of the central groove location. The sagittal projection processing of the data will introduce a certain degree of distortion, which will reduce the accuracy of the central groove location. In order to solve the above problems, the embodiments of this application perform smoothing processing on the initial head contour edge image obtained above based on a smoothing processor.

[0092] In one optional embodiment of this application, a Gaussian filter is selected as the smoothing processor to smooth the initial head contour edge image;

[0093] Specifically, obtain the first smoothing parameter, for example: scale 3*3, and at the same time, obtain the second smoothing parameter, for example: δ=0.5;

[0094] A Gaussian filter smooths the initial head contour edge image based on the obtained 3*3 scale and δ=0, thereby... Figure 8 Smooth the uneven edges in the image to obtain a result like... Figure 9 The image shown shows the outline of the head.

[0095] Step S42: On the head contour edge image, determine multiple target feature points; wherein, the multiple target feature points include nasal root feature points, pterional feature points, and occipital protuberance feature points;

[0096] The embodiments of this application, based on related technologies, determine the nasal root feature points, pterional feature points, and occipital protuberance feature points respectively; however, the specific methods for determining the nasal root feature points, pterional feature points, and occipital protuberance feature points are not within the scope of protection of this application, and therefore will not be described in detail.

[0097] After obtaining the nasion feature point, the ala feature point and the inion feature point, the nasion feature point, the ala feature point and the inion feature point are respectively mapped to corresponding pixel points on the head contour edge image based on the method of minimum Euclidean distance.

[0098] Step S43: determining a reference distance based on the first coordinate of the nasion feature point and the second coordinate of the inion feature point.

[0099] In an optional embodiment of the present application, the determination of the reference distance based on the first coordinate of the nasion feature point and the second coordinate of the inion feature point comprises:

[0100] On the head contour edge image, all edge pixel points between the nasion feature point and the inion feature point are sequentially formed into a chain fourth coordinate point set with the nasion feature point as the starting point and the inion feature point as the ending point along a preset direction.

[0101] The reference distance between the nasion feature point and the inion feature point is obtained based on the fourth coordinate point set and in combination with the spatial resolution information of the head contour edge image.

[0102] In an optional embodiment of the present application, as shown in Figure 10 , the edge portion between the nasion feature point P1(x1, y1, z1) and the inion feature point P M (x M ,y M ,z M ) is intercepted, and a plurality of fourth coordinate points P2, P3, …, P M-1 are sequentially determined between the nasion feature point P1(x1, y1, z1) and the inion feature point P M (x M ,y M ,z M ) on the edge portion along a preset direction.

[0103] The fourth coordinate corresponding to each fourth coordinate point is:

[0104] P2(x2, y2, z2), P3(x3, y3, z3), …, P M-1 (x M-1 , y M-1 , z M-1 ).

[0105] In an optional embodiment of the present application, the reference distance between the nasion feature point and the inion feature point is obtained based on the fourth set of coordinate points and spatial resolution information of the head contour edge image.

[0106] The reference distance between the nasion feature point and the inion feature point is obtained by cumulatively adding distances between every two adjacent fourth coordinate points based on the fourth set of coordinate points and spatial resolution information of the head contour edge image.

[0107] In an optional embodiment of the present application, the spatial resolution information of the head contour edge image comprises:

[0108] The first sub-resolution in the first spatial direction of the head contour image is obtained; the second sub-resolution in the second spatial direction of the head contour image is obtained; and the third sub-resolution in the third spatial direction of the head contour image is obtained.

[0109] In an optional embodiment of the present application, the first spatial direction can correspond to the x-axis direction; the second spatial direction can correspond to the y-axis direction; and the third spatial direction can correspond to the z-axis direction.

[0110] The first sub-resolution can be x space ; the second sub-resolution can be y space ; and the third sub-resolution can be z space .

[0111] In an optional embodiment of the present application, the reference distance between the nasion feature point and the inion feature point is obtained by cumulatively adding distances between every two adjacent fourth coordinate points based on the fourth set of coordinate points and spatial resolution information of the head contour edge image.

[0112] The first sub-coordinate of every two adjacent fourth coordinate points is obtained, and a first term is calculated based on each first sub-coordinate point; the second sub-coordinate of every two adjacent fourth coordinate points is obtained, and a second term is calculated based on each second sub-coordinate point; and the third sub-coordinate of every two adjacent fourth coordinate points is obtained, and a third term is calculated based on each third sub-coordinate point.

[0113] The cumulatively added term corresponding to every two adjacent fourth coordinate points is calculated based on the first term, the second term and the third term corresponding to the two adjacent fourth coordinate points.

[0114] The cumulatively added term corresponding to every two adjacent fourth coordinate points is calculated to obtain the reference distance between the nasion feature point and the inion feature point.

[0115] In an optional embodiment of the present application, the first term, the second term and the third term can be obtained based on the following formulas, respectively:

[0116] First term = ((x t+1 -x t )*x space ) 2 .

[0117] Second term = ((y t+1 -y t )*y space ) 2 .

[0118] Third term = ((z t+1 -z t )*x space ) 2 .

[0119] In an optional embodiment of the present application, the cumulative superposition term can be obtained based on the following formula:

[0120]

[0121] In an optional embodiment of the present application, the reference distance d can be obtained based on the following formula:

[0122]

[0123] Step S44: determining the location of the central sulcus based on the reference distance and the first coordinate;

[0124] In an optional embodiment of the present application, the step of determining the location of the central sulcus based on the reference distance and the first coordinate comprises:

[0125] determining a plurality of central sulcus points on the head profile edge image based on the reference distance and the first coordinate;

[0126] obtaining a fifth coordinate of each of the central sulcus points;

[0127] determining the location of the central sulcus based on the fifth coordinate of each of the central sulcus points.

[0128] In an optional embodiment of the present application, the step of determining a plurality of central sulcus points on the head profile edge image based on the reference distance and the first coordinate comprises:

[0129] obtaining a preset distance;

[0130] determining a curve distance based on the preset distance and the reference distance;

[0131] determining a point on the head profile edge image as the central sulcus point if the distance between the point and the nasal root feature point is equal to the curve distance.

[0132] In an optional embodiment of the present application, after the reference distance is determined, the curve distance is determined based on the reference distance and the preset distance.

[0133] The curve distance = d / 2 + the preset distance.

[0134] For example, the preset distance is 2 cm.

[0135] In an optional embodiment of the present application, as shown in Figure 11 , the first coordinate is taken as a starting point, all points on the curved surface between the nasion feature point and the inion feature point with a distance equal to the curve distance from the nasion feature point are found, and then a plurality of central sulcus points T k are obtained.

[0136] As shown in Figure 12 , the fifth coordinate of each central sulcus point T k is obtained, and the central sulcus position is determined based on each fifth coordinate.

[0137] It should be noted that based on the above method, a plurality of fifth coordinates and central sulcus positions can be determined on each tomographic image.

[0138] Step S45: determining a target space position based on the central sulcus position and the third coordinate of the ala point feature point; wherein the target space position includes a central anterior gyrus space position and a central posterior gyrus space position.

[0139] In an optional embodiment of the present application, as shown in Figure 12 , the determination of the target space position based on the central sulcus position and the third coordinate of the ala point feature point includes:

[0140] selecting a first plane point on the i-th tomographic image of the head contour edge image and a second plane point on the j-th tomographic image of the head contour edge image based on the central sulcus position; wherein the head contour edge image includes a plurality of tomographic images, i and j are both positive integers, and i≠j;

[0141] obtaining a first plane coordinate of the first plane point and a second plane coordinate of the second plane point;

[0142] determining a reference plane equation and a reference plane based on the first plane coordinate, the second plane coordinate, and the third coordinate;

[0143] a preset boundary distance, and determining a region on the head contour edge image with a distance from the reference plane less than or equal to the boundary distance as the target space position; wherein the central anterior gyrus space position and the central posterior gyrus space position are symmetrical about the central sulcus position.

[0144] In an optional embodiment of the present application, a central sulcus point T is selected at the central sulcus position on the i-th tomographic image k ; and another central sulcus point H is selected at the central sulcus position on the j-th tomographic image k ;

[0145] The first plane coordinate corresponding to the central sulcus point T k , the second plane coordinate corresponding to the central sulcus point H k , and the third coordinate corresponding to the wing point feature point are obtained respectively.

[0146] Based on the first plane coordinate, the second plane coordinate, and the third coordinate, a reference plane equation Ax+By+Cz=0 is determined; wherein A, B, and C are coefficients corresponding to the reference plane equation. Because the selected first plane coordinate and the second plane coordinate are different, the coefficients corresponding to the reference plane equation are also different, and thus embodiments of the present application do not specifically limit the coefficients of the reference plane equation.

[0147] Further, based on the reference plane equation, a reference plane can be determined.

[0148] In an optional embodiment of the present application, j=i+5.

[0149] In addition, the size relationship between j and i can also be adjusted according to actual needs.

[0150] In an optional embodiment of the present application, it is assumed that the coordinates of a point on the head contour edge image are Q(x Q ,y Q ,z Q );

[0151] The distance of the point to the above-mentioned reference plane can be calculated based on the following formula:

[0152]

[0153] When D≤boundary distance, Q is a target point.

[0154] Conversely, when D>boundary distance, Q is not a target point.

[0155] By analogy, all target points can be found on the head contour edge image, and the target spatial position can be determined according to all target points.

[0156] For example, the boundary distance can be 2 cm.

[0157] In an optional embodiment of the present application, as Figure 13 and 14As shown, since the spatial positions of the anterior central gyrus and the posterior central gyrus are symmetrical about the position of the central groove, in the target spatial position, the first region on one side of the central groove position is the spatial position of the anterior central gyrus, and the second region on the other side of the central groove position is the spatial position of the posterior central gyrus.

[0158] The embodiments of this application automatically identify and determine the location of the central sulcus based on the coordinates of three feature points on the CT image sequence: the nasal root feature point, the occipital protuberance feature point, and the pterional feature point. Based on the location of the central sulcus, the spatial location of the anterior central gyrus and the location of the posterior central gyrus are determined. The results are highly accurate, highly compatible, and fast in finding the location.

[0159] An optional embodiment of this application, running in an environment with the following hardware configuration (Intel(R) Core(TM) i5-8500 CPU @ 3.00 GHz, memory: 8.0 GB, system type: 64-bit, based on x64 processor), ran 200 sets of data, with an average of 165 tomographic images per set, an average running time of 7 seconds, and an accuracy rate of 95%, achieving automated batch processing.

[0160] like Figure 15 As shown, embodiments of this application also provide an image processing apparatus 150, comprising:

[0161] Processing module 151 is used to acquire a head contour image and process the head contour image to obtain a head contour edge image.

[0162] The first determining module 152 is used to determine multiple target feature points on the head contour edge image; wherein the multiple target feature points include nasal root feature points, pterional feature points and occipital protuberance feature points.

[0163] The second determining module 153 is used to determine a reference distance based on the first coordinates of the nasal root feature point and the second coordinates of the occipital protuberance feature point.

[0164] The third determining module 154 is used to determine the position of the central ditch based on the reference distance and the first coordinates;

[0165] The fourth determining module 155 is used to determine the target spatial position based on the position of the central groove and the third coordinates of the wing feature points; wherein the target spatial position includes the central anterior gyral spatial position and the central posterior gyral spatial position.

[0166] Optionally, before obtaining the head contour edge image, the method further includes:

[0167] Based on the head contour image, an initial head contour edge image is obtained;

[0168] obtaining a first smoothing parameter and a second smoothing parameter;

[0169] The smoothing processor smoothes the initial head contour edge image based on the first smoothing parameter and the second smoothing parameter to obtain the head contour edge image.

[0170] Optionally, the reference distance is determined based on the first coordinate of the nasion feature point and the second coordinate of the inion feature point, and includes:

[0171] On the head contour edge image, all edge pixel points between the nasion feature point and the inion feature point are sequentially formed into a chain fourth coordinate point set along a preset direction with the nasion feature point as a starting point and the inion feature point as an ending point.

[0172] The reference distance between the nasion feature point and the inion feature point is obtained based on the fourth coordinate point set and in combination with spatial resolution information of the head contour edge image.

[0173] Optionally, the reference distance between the nasion feature point and the inion feature point is obtained based on the fourth coordinate point set and in combination with spatial resolution information of the head contour edge image, and includes:

[0174] The distance between every two adjacent fourth coordinate points is accumulated and superimposed based on the fourth coordinate point set and the spatial resolution information of the head contour edge image to obtain the reference distance between the nasion feature point and the inion feature point.

[0175] Optionally, the distance between every two adjacent fourth coordinate points is accumulated and superimposed to obtain the reference distance between the nasion feature point and the inion feature point, and includes:

[0176] A first sub-coordinate of every two adjacent fourth coordinate points is obtained, and a first term is calculated based on each first sub-coordinate point; a second sub-coordinate of every two adjacent fourth coordinate points is obtained, and a second term is calculated based on each second sub-coordinate point; and a third sub-coordinate of every two adjacent fourth coordinate points is obtained, and a third term is calculated based on each third sub-coordinate point.

[0177] A cumulative superimposition term corresponding to every two adjacent fourth coordinate points is calculated based on the first term, the second term and the third term corresponding to the two adjacent fourth coordinate points.

[0178] The cumulative superimposition term corresponding to every two adjacent fourth coordinate points is calculated to obtain the reference distance between the nasion feature point and the inion feature point.

[0179] Optionally, the determining the location of the central sulcus based on the reference distance and the first coordinate comprises:

[0180] determining a plurality of central sulcus points on the head contour edge image based on the reference distance and the first coordinate;

[0181] obtaining a fifth coordinate of each of the central sulcus points;

[0182] determining the location of the central sulcus based on the fifth coordinate of each of the central sulcus points.

[0183] Optionally, the determining the plurality of central sulcus points on the head contour edge image based on the reference distance and the first coordinate comprises:

[0184] obtaining a preset distance;

[0185] determining a curve distance based on the preset distance and the reference distance;

[0186] determining a point on the head contour edge image as the central sulcus point based on the distance between the point and the nasal root feature point being equal to the curve distance.

[0187] Optionally, the determining the target spatial location based on the location of the central sulcus and the third coordinate of the alar point feature point comprises:

[0188] selecting a first plane point on an i-th slice of the head contour edge image and a second plane point on a j-th slice of the head contour edge image based on the location of the central sulcus, wherein the head contour edge image comprises a plurality of slices, i and j are positive integers, and i≠j;

[0189] obtaining a first plane coordinate of the first plane point and a second plane coordinate of the second plane point;

[0190] determining a reference plane equation and a reference plane based on the first plane coordinate, the second plane coordinate, and the third coordinate;

[0191] determining a target spatial location based on a preset boundary distance and a region on the head contour edge image having a distance from the reference plane less than or equal to the boundary distance, wherein the precentral sulcus spatial location and the postcentral sulcus spatial location are symmetrical about the central sulcus location.

[0192] Embodiments of the present application also provide an electronic device comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the method as described above when executing the computer program.

[0193] Embodiments of the present application also provide a computer readable storage medium comprising a stored computer program, wherein the computer readable storage medium, when the computer program runs, controls the device where the computer readable storage medium is located to perform the method as described above.

[0194] In addition, other configurations and functions of the device of the embodiments of the present application are known to those skilled in the art, and to reduce redundancy, they are not described here.

[0195] It should be noted that the logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a list of executable instructions for implementing logic functions, and can be embodied in any computer-readable medium for use by an instruction execution system, apparatus, or device, such as a computer-based system, a system including a processor, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions, or in conjunction with these instructions execution systems, apparatus, or devices. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport the program for use by an instruction execution system, apparatus, or device, or in conjunction with these instruction execution systems, apparatus, or devices. More specific examples (a non-exhaustive list) of the computer-readable medium include the following: an electrical connection having one or more wires (electronic devices), a portable computer diskette (magnetic devices), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber device, and a portable compact disc read-only memory (CD-ROM). In addition, the computer-readable medium can even be paper or other suitable medium on which the program can be printed, as the program can be electronically obtained, for example, by optical scanning of the paper or other medium, followed by electronic conversion into a useable form, and then stored in computer memory.

[0196] It should be understood that parts of the present application can be realized in hardware, software, firmware, or a combination thereof. In the above-described embodiments, a plurality of steps or methods can be realized by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if realized in hardware, and as in another embodiment, it can be realized by any one or a combination of the following technologies known in the art: discrete logic circuit having logic gates for implementing logic functions on data signals, application specific integrated circuit having suitable combination logic gates, programmable gate array (PGA), field programmable gate array (FPGA), etc.

[0197] In the description of the application, the description of the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the application. In the description of the application, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0198] In the description of the application, it should be understood that the orientation or positional relationship indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential" and the like is based on the orientation or positional relationship shown in the drawings, and is only for the purpose of facilitating the description of the application and simplifying the description, and does not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation of the application.

[0199] In addition, the terms "first", "second" are only for the purpose of description, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Therefore, the features defined with "first", "second" can explicitly or implicitly include at least one of the features. In the description of the application, the meaning of "a plurality of" is at least two, for example, two, three, etc., unless otherwise explicitly specified and limited.

[0200] In the present application, unless otherwise explicitly specified and limited, the terms "mounting", "connection", "connection", "fixing" and the like should be understood in a broad sense, for example, it can be fixed connection, or detachable connection, or integral; it can be mechanical connection, or electrical connection; it can be directly connected, or indirectly connected through an intermediate medium, it can be the internal communication of two elements or the interaction relationship between two elements, unless otherwise explicitly limited. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0201] In the present application, unless specifically stated and limited otherwise, a first feature "on" or "under" a second feature can be directly in contact with the second feature, or indirectly in contact with the second feature through an intermediate medium. Also, a first feature "over", "above" and "on top of" a second feature can be directly above or obliquely above the second feature, or simply means that the first feature is horizontally higher than the second feature. A first feature "under", "below" and "underneath" a second feature can be directly below or obliquely below the second feature, or simply means that the first feature is horizontally lower than the second feature.

[0202] Although the embodiments of the present application have been shown and described above, it should be understood by those skilled in the art that the above-described embodiments are exemplary, and cannot be construed as limiting the present application, and those skilled in the art can make changes, modifications, replacements and variations to the above-described embodiments within the scope of the present application.

Claims

1. An image processing method, characterized in that, include: A head contour image is acquired, and the head contour image is processed to obtain a head contour edge image. On the head contour edge image, multiple target feature points are determined; wherein, the multiple target feature points include nasal root feature points, pterional feature points, and occipital protuberance feature points; A reference distance is determined based on the first coordinates of the nasal root feature point and the second coordinates of the occipital protuberance feature point. Based on the reference distance and the first coordinates, the location of the central ditch is determined; The target spatial position is determined based on the location of the central groove and the third coordinates of the wing feature points; wherein, the target spatial position includes the central anterior gyrus spatial position and the central posterior gyrus spatial position; Determining the target spatial location based on the location of the central trench and the third coordinates of the wing feature points includes: Based on the location of the central groove, a first plane point is selected on the i-th tomographic image of the head contour edge image, and a second plane point is selected on the j-th tomographic image of the head contour edge image; wherein, the head contour edge image includes multiple tomographic images, i and j are both positive integers, and i≠j; Obtain the first plane coordinates of the first plane point and the second plane coordinates of the second plane point; Based on the first plane coordinates, the second plane coordinates, and the third coordinates, determine the reference plane equation and the reference plane; A preset boundary distance is set, and the region on the head contour edge image whose distance from the reference plane is less than or equal to the boundary distance is determined as the target spatial position; wherein, the central anterior gyrus spatial position and the central posterior gyrus spatial position are symmetrical about the central groove position.

2. The method according to claim 1, characterized in that, Before obtaining the head contour edge image, the process also includes: Based on the head contour image, an initial head contour edge image is obtained; Obtain the first smoothing parameter and the second smoothing parameter; The smoothing processor performs smoothing processing on the initial head contour edge image based on the first smoothing parameter and the second smoothing parameter to obtain the head contour edge image.

3. The method according to claim 1, characterized in that, The determination of the reference distance based on the first coordinates of the nasal root feature point and the second coordinates of the occipital protuberance feature point includes: On the head contour edge image, along a preset direction, starting from the nasal root feature point and ending at the occipital protuberance feature point, all edge pixels between the nasal root feature point and the occipital protuberance feature point are sequentially chained to form a fourth coordinate point set. Based on the fourth set of coordinate points and combined with the spatial resolution information of the head contour edge image, the reference distance between the nasal root feature point and the occipital protuberance feature point is obtained.

4. The method according to claim 3, characterized in that, The step of obtaining the reference distance from the nasal root feature point to the occipital protuberance feature point based on the fourth coordinate point set and combined with the spatial resolution information of the head contour edge image includes: Based on the spatial resolution information of the fourth coordinate point set and the head contour edge image, the distance between every two adjacent fourth coordinate points is accumulated and superimposed to obtain the reference distance from the nasal root feature point to the occipital protuberance feature point.

5. The method according to claim 4, characterized in that, The cumulative summation of the distances between every two adjacent fourth coordinate points to obtain the reference distance from the nasal root feature point to the occipital protuberance feature point includes: Obtain the first sub-coordinate of every two adjacent fourth coordinate points, and calculate the first item based on each first sub-coordinate point; obtain the second sub-coordinate of every two adjacent fourth coordinate points, and calculate the second item based on each second sub-coordinate point; obtain the third sub-coordinate of every two adjacent fourth coordinate points, and calculate the third item based on each third sub-coordinate point; Based on the first, second, and third terms corresponding to each pair of adjacent fourth coordinate points, the cumulative superposition term corresponding to each pair of adjacent fourth coordinate points is calculated respectively. The cumulative superposition terms corresponding to every two adjacent fourth coordinate points are calculated to obtain the reference distance between the nasal root feature point and the occipital protuberance feature point.

6. The method according to claim 1, characterized in that, Determining the location of the central ditch based on the reference distance and the first coordinates includes: Based on the reference distance and the first coordinates, multiple central groove points are determined on the head contour edge image; Obtain the fifth coordinate of each of the central trench points; The location of the central ditch is determined based on the fifth coordinate of each of the central ditch points.

7. The method according to claim 6, characterized in that, Based on the reference distance and the first coordinates, multiple central groove points are determined on the head contour edge image, including: Get the preset distance; The curve distance is determined based on the preset distance and the reference distance; On the head contour edge image, the point whose distance from the nasal root feature point is equal to the curve distance is determined as the central groove point.

8. An image processing apparatus, characterized in that, include: The processing module is used to acquire a head contour image and process the head contour image to obtain a head contour edge image. The first determining module is used to determine multiple target feature points on the head contour edge image; wherein the multiple target feature points include nasal root feature points, pterional feature points and occipital protuberance feature points. The second determining module is used to determine a reference distance based on the first coordinates of the nasal root feature point and the second coordinates of the occipital protuberance feature point. The third determining module is used to determine the location of the central ditch based on the reference distance and the first coordinates; The fourth determining module is used to determine the target spatial position based on the position of the central groove and the third coordinates of the wing feature points; wherein, the target spatial position includes the central anterior gyral spatial position and the central posterior gyral spatial position; Determining the target spatial location based on the location of the central trench and the third coordinates of the wing feature points includes: Based on the location of the central groove, a first plane point is selected on the i-th tomographic image of the head contour edge image, and a second plane point is selected on the j-th tomographic image of the head contour edge image; wherein, the head contour edge image includes multiple tomographic images, i and j are both positive integers, and i≠j; Obtain the first plane coordinates of the first plane point and the second plane coordinates of the second plane point; Based on the first plane coordinates, the second plane coordinates, and the third coordinates, determine the reference plane equation and the reference plane; A preset boundary distance is set, and the region on the head contour edge image whose distance from the reference plane is less than or equal to the boundary distance is determined as the target spatial position; wherein, the central anterior gyrus spatial position and the central posterior gyrus spatial position are symmetrical about the central groove position.

9. An electronic device, characterized in that, The method includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor, when executing the computer program, implements the method as claimed in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium includes a stored computer program, wherein, when the computer program is executed, it controls the device on which the computer-readable storage medium is located to perform the method as described in any one of claims 1 to 7.

Images