[go: up one dir, main page]

CN118365632B - Image processing method and device and computer-aided diagnosis method of brain diseases - Google Patents

Image processing method and device and computer-aided diagnosis method of brain diseases Download PDF

Info

Publication number
CN118365632B
CN118365632B CN202410766838.5A CN202410766838A CN118365632B CN 118365632 B CN118365632 B CN 118365632B CN 202410766838 A CN202410766838 A CN 202410766838A CN 118365632 B CN118365632 B CN 118365632B
Authority
CN
China
Prior art keywords
vertex
target
information
brain
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202410766838.5A
Other languages
Chinese (zh)
Other versions
CN118365632A (en
Inventor
赵奋强
唐禹行
张灵
吕乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Alibaba China Co Ltd
Original Assignee
Alibaba China Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Alibaba China Co Ltd filed Critical Alibaba China Co Ltd
Priority to CN202410766838.5A priority Critical patent/CN118365632B/en
Publication of CN118365632A publication Critical patent/CN118365632A/en
Application granted granted Critical
Publication of CN118365632B publication Critical patent/CN118365632B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Pathology (AREA)
  • General Engineering & Computer Science (AREA)
  • Epidemiology (AREA)
  • Computer Hardware Design (AREA)
  • Primary Health Care (AREA)
  • Computer Graphics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本说明书实施例提供图像处理方法及装置、脑部疾病的计算机辅助诊断方法,其中图像处理方法包括:接收图像处理任务,其中,图像处理任务携带三维脑图像;将三维脑图像输入至解剖分割模型,获得解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;根据灰质信息生成目标大脑皮层,并确定目标大脑皮层上顶点的曲率特征信息;根据灰质信息和目标大脑皮层上顶点的曲率特征信息,确定目标大脑皮层对应的皮层表面分割结果;根据皮层表面分割结果、白质信息和皮层下结构信息,确定三维脑图像对应的全脑分割结果。本方法提供了从粗粒度到细粒度的两步分割方法,加快了数据处理速度,又能提升图像处理的准确度。

The embodiments of this specification provide an image processing method and device, and a computer-aided diagnosis method for brain diseases, wherein the image processing method includes: receiving an image processing task, wherein the image processing task carries a three-dimensional brain image; inputting the three-dimensional brain image into an anatomical segmentation model, obtaining gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on the contrast information of the three-dimensional brain image; generating a target cerebral cortex according to the gray matter information, and determining the curvature characteristic information of the vertices on the target cerebral cortex; determining the cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and the curvature characteristic information of the vertices on the target cerebral cortex; determining the whole brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, white matter information, and subcortical structure information. This method provides a two-step segmentation method from coarse-grained to fine-grained, which speeds up data processing and improves the accuracy of image processing.

Description

图像处理方法及装置、脑部疾病的计算机辅助诊断方法Image processing method and device, computer-aided diagnosis method for brain diseases

技术领域Technical Field

本说明书实施例涉及计算机技术领域,特别涉及一种图像处理方法。The embodiments of the present specification relate to the field of computer technology, and in particular to an image processing method.

背景技术Background Art

在医学领域的脑部磁共振成像(Magnetic Resonance Imaging,MRI)分析中,全脑结构的分割是一个基本且关键的步骤。它为每个体素分配一个与特定神经解剖结构相对应的语义标签,从而可以获得诸如体积、厚度、各种感兴趣区域的面积等定量属性值。这些属性值对于许多临床应用和科学研究至关重要。In the analysis of brain magnetic resonance imaging (MRI) in the medical field, the segmentation of whole brain structure is a basic and critical step. It assigns a semantic label corresponding to a specific neuroanatomical structure to each voxel, so that quantitative attribute values such as volume, thickness, and area of various regions of interest can be obtained. These attribute values are crucial for many clinical applications and scientific research.

目前的全脑分割算法通常分为两类,第一类是基于脑图谱的方法,第二类是利用深度学习技术开发端到端的分割网络。第一类方法在将灰质细分为更小区域很难将不同区域对齐,且算法较为复杂和耗时。第二类方法处理速度较快,但是其准确性较差。因此,亟需一种新的全脑分割方法解决上述问题。The current whole-brain segmentation algorithms are generally divided into two categories. The first category is based on brain maps, and the second category is to develop an end-to-end segmentation network using deep learning technology. The first method is difficult to align different areas when subdividing gray matter into smaller areas, and the algorithm is more complex and time-consuming. The second method has a faster processing speed, but its accuracy is poor. Therefore, a new whole-brain segmentation method is urgently needed to solve the above problems.

发明内容Summary of the invention

有鉴于此,本说明书实施例提供了一种图像处理方法。本说明书一个或者多个实施例同时涉及一种脑部疾病的计算机辅助诊断方法、一种图像处理装置,一种计算设备,一种计算机可读存储介质以及一种计算机程序产品,以解决现有技术中存在的技术缺陷。In view of this, an embodiment of this specification provides an image processing method. One or more embodiments of this specification also relate to a computer-aided diagnosis method for brain diseases, an image processing device, a computing device, a computer-readable storage medium, and a computer program product to solve the technical defects existing in the prior art.

根据本说明书实施例的第一方面,提供了一种图像处理方法,包括:According to a first aspect of an embodiment of this specification, there is provided an image processing method, including:

接收图像处理任务,其中,所述图像处理任务携带三维脑图像;receiving an image processing task, wherein the image processing task carries a three-dimensional brain image;

将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image;

根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex;

根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex;

根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。A whole-brain segmentation result corresponding to the three-dimensional brain image is determined according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

根据本说明书实施例的第二方面,提供了一种图像处理方法,应用于云侧设备,包括:According to a second aspect of an embodiment of this specification, there is provided an image processing method, which is applied to a cloud-side device, and includes:

接收端侧发送的图像处理任务,其中,所述图像处理任务携带三维脑图像;An image processing task sent by a receiving end, wherein the image processing task carries a three-dimensional brain image;

将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image;

根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex;

根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex;

根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果;Determining a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information, and the subcortical structure information;

向所述端侧设备发送所述全脑分割结果。The whole-brain segmentation result is sent to the client-side device.

根据本说明书实施例的第三方面,提供了一种脑部疾病的计算机辅助诊断方法,包括:According to a third aspect of the embodiments of this specification, a computer-aided diagnosis method for brain diseases is provided, comprising:

接收脑部疾病检测任务,其中,所述脑部疾病检测任务携带三维脑图像和目标脑结构标识;receiving a brain disease detection task, wherein the brain disease detection task carries a three-dimensional brain image and a target brain structure identifier;

将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image;

根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex;

根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex;

根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果;Determining a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information, and the subcortical structure information;

根据所述目标脑结构标识在所述全脑分割结果中确定目标脑结构,针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果。The target brain structure is determined in the whole-brain segmentation result according to the target brain structure identifier, and the target brain structure is detected to obtain the segmentation and detection results corresponding to the target brain structure.

根据本说明书实施例的第四方面,提供了一种图像处理装置,包括:According to a fourth aspect of the embodiments of this specification, there is provided an image processing device, including:

接收模块,被配置为接收图像处理任务,其中,所述图像处理任务携带三维脑图像;A receiving module, configured to receive an image processing task, wherein the image processing task carries a three-dimensional brain image;

第一分割模块,被配置为将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;A first segmentation module is configured to input the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image;

生成模块,被配置为根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;A generating module, configured to generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex;

第二分割模块,被配置为根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;A second segmentation module is configured to determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex;

确定模块,被配置为根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。The determination module is configured to determine the whole brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

根据本说明书实施例的第五方面,提供了一种计算设备,包括:According to a fifth aspect of an embodiment of this specification, a computing device is provided, including:

存储器和处理器;Memory and processor;

所述存储器用于存储计算机程序/指令,所述处理器用于执行所述计算机程序/指令,该计算机程序/指令被处理器执行时实现上述方法的步骤。The memory is used to store computer programs/instructions, and the processor is used to execute the computer programs/instructions. When the computer programs/instructions are executed by the processor, the steps of the above method are implemented.

根据本说明书实施例的第六方面,提供了一种计算机可读存储介质,其存储有计算机程序/指令,该计算机程序/指令被处理器执行时实现上述方法的步骤。According to a sixth aspect of the embodiments of this specification, a computer-readable storage medium is provided, which stores a computer program/instruction, and the steps of the above method are implemented when the computer program/instruction is executed by a processor.

根据本说明书实施例的第七方面,提供了一种计算机程序产品,包括计算机程序/指令,该计算机程序/指令被处理器执行时实现上述方法的步骤。According to a seventh aspect of the embodiments of this specification, a computer program product is provided, comprising a computer program/instruction, which implements the steps of the above method when executed by a processor.

本说明书一个实施例提供了一种图像处理方法,包括:接收图像处理任务,其中,所述图像处理任务携带三维脑图像;将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。An embodiment of the present specification provides an image processing method, comprising: receiving an image processing task, wherein the image processing task carries a three-dimensional brain image; inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image; generating a target cerebral cortex according to the gray matter information, and determining curvature feature information of at least one vertex on the target cerebral cortex; determining a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and the curvature feature information of at least one vertex on the target cerebral cortex; and determining a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information, and the subcortical structure information.

通过本申请实施例提供的图像处理方法,采用从粗粒度到细粒度的两步分割方法,将三维脑图像输入到解剖分割模型中,解剖分割模型仅根据对比度信息分割出灰质信息、白质信息和皮层下结构信息。再进一步的根据灰质信息生成大脑皮层,并计算大脑皮层上各顶点的曲率特征,曲率特征代表了大脑皮层的结构信息,将代表结构信息的曲率特征灰质信息确定大脑皮层对应的皮层表面分割结果,解决了大脑分割准确率低的问题。最后根据白质信息、皮层下结构信息和皮层表面分割结果,确定三维脑图像对应的全脑分割结果,既节省了模型处理时间,加快了数据处理速度,又能提升图像处理的准确度。The image processing method provided by the embodiment of the present application adopts a two-step segmentation method from coarse-grained to fine-grained, and the three-dimensional brain image is input into the anatomical segmentation model. The anatomical segmentation model only segments the gray matter information, white matter information and subcortical structure information based on the contrast information. Further, the cerebral cortex is generated based on the gray matter information, and the curvature characteristics of each vertex on the cerebral cortex are calculated. The curvature characteristics represent the structural information of the cerebral cortex. The gray matter information representing the curvature characteristics of the structural information is used to determine the cortical surface segmentation result corresponding to the cerebral cortex, which solves the problem of low brain segmentation accuracy. Finally, based on the white matter information, subcortical structure information and cortical surface segmentation results, the whole brain segmentation result corresponding to the three-dimensional brain image is determined, which not only saves model processing time, speeds up data processing speed, but also improves the accuracy of image processing.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1是本说明书一个实施例提供的一种图像处理方法的流程图;FIG1 is a flow chart of an image processing method provided by an embodiment of the present specification;

图2是本说明书一个实施例提供的一种图像处理方法的处理结构示意图;FIG2 is a schematic diagram of a processing structure of an image processing method provided by an embodiment of this specification;

图3是本说明书一个实施例提供的一种图像处理系统的架构图;FIG3 is an architecture diagram of an image processing system provided by an embodiment of the present specification;

图4是本说明书一个实施例提供的应用于云侧设备的图像处理方法的流程示意图;FIG4 is a schematic diagram of a flow chart of an image processing method applied to a cloud-side device provided by an embodiment of the present specification;

图5是本说明书一个实施例提供的一种脑部疾病的计算机辅助诊断方法的流程示意图;FIG5 is a flow chart of a computer-aided diagnosis method for brain diseases provided in one embodiment of the present specification;

图6是本说明书一个实施例提供的一种图像处理装置的结构示意图;FIG6 is a schematic diagram of the structure of an image processing device provided by an embodiment of this specification;

图7是本说明书一个实施例提供的一种计算设备的结构框图。FIG. 7 is a structural block diagram of a computing device provided by an embodiment of the present specification.

具体实施方式DETAILED DESCRIPTION

在下面的描述中阐述了很多具体细节以便于充分理解本说明书。但是本说明书能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本说明书内涵的情况下做类似推广,因此本说明书不受下面公开的具体实施的限制。Many specific details are described in the following description to facilitate a full understanding of this specification. However, this specification can be implemented in many other ways than those described herein, and those skilled in the art can make similar generalizations without violating the connotation of this specification, so this specification is not limited to the specific implementation disclosed below.

在本说明书一个或多个实施例中使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本说明书一个或多个实施例。在本说明书一个或多个实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本说明书一个或多个实施例中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。The terms used in one or more embodiments of this specification are only for the purpose of describing specific embodiments, and are not intended to limit one or more embodiments of this specification. The singular forms of "a", "said" and "the" used in one or more embodiments of this specification and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings. It should also be understood that the term "and/or" used in one or more embodiments of this specification refers to and includes any or all possible combinations of one or more associated listed items.

应当理解,尽管在本说明书一个或多个实施例中可能采用术语第一、第二等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本说明书一个或多个实施例范围的情况下,第一也可以被称为第二,类似地,第二也可以被称为第一。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。It should be understood that although the terms first, second, etc. may be used to describe various information in one or more embodiments of this specification, this information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of one or more embodiments of this specification, the first may also be referred to as the second, and similarly, the second may also be referred to as the first. Depending on the context, the word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

需要说明的是,本说明书所涉及的用户信息(包括但不限于用户设备信息、用户个人信息等)和数据(包括但不限于用于分析的数据、存储的数据、展示的数据等),均为经用户授权或者经过各方充分授权的信息和数据,并且相关数据的收集、使用和处理需要遵守相关地区的相关法律法规和标准,并提供有相应的操作入口,供用户选择授权或者拒绝。It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this manual are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the relevant laws, regulations and standards of the relevant regions, and provide corresponding operation entrances for users to choose to authorize or refuse.

首先,对本说明书一个或多个实施例涉及的名词术语进行解释。First, the terms involved in one or more embodiments of this specification are explained.

MRI:Magneticresonanceimaging,磁共振成像技术或常指磁共振成像技术采集的磁共振图像。MRI: Magnetic resonance imaging, magnetic resonance imaging technology or often refers to magnetic resonance images collected by magnetic resonance imaging technology.

WM:Whitematter,白质,指中枢神经系统中主要由被髓鞘包覆的神经轴突(可形成神经束)所组成的区域,在结构MRI中通常呈现为白色。WM: White matter, refers to the area in the central nervous system that is mainly composed of nerve axons covered by myelin (which can form nerve bundles), which usually appears white in structural MRI.

GM:Graymatter,又称皮层(Cortex),是中枢神经系统的重要组成部分,在结构MRI中通常呈现为灰色。GM: Graymatter, also known as the cortex, is an important part of the central nervous system and usually appears gray in structural MRI.

CorticalParcellation:皮层分区,即将皮层分成不同的区域,是脑影像分析的先决步骤。CorticalParcellation: Cortical parcellation, which is the division of the cortex into different regions, is a prerequisite step in brain imaging analysis.

全脑分割:全脑分割是一项重要的神经成像任务,即将大脑MRI图像细分成不同的区域,包括WM、GM、小脑、脑干,以及更精细的皮层分区等,是脑影像分析的先决步骤。Whole-brain segmentation: Whole-brain segmentation is an important neuroimaging task, which is to subdivide brain MRI images into different regions, including WM, GM, cerebellum, brainstem, and more detailed cortical partitions, etc. It is a prerequisite for brain image analysis.

在医学领域的脑部磁共振成像(Magnetic Resonance Imaging,MRI)分析中,全脑结构的分割是一个基本且关键的步骤。大脑结构磁共振图像的全脑分割将整个大脑划分为多个解剖相关的感兴趣区域,是后续定量分析不同脑结构的先决步骤。其为每个体素分配一个与特定神经解剖结构相对应的语义标签,从而可以获得诸如体积、厚度、各种感兴趣区域的面积等定量属性值,这些属性值可以用于脑形态分析、手术规划和术后评估等。In the analysis of brain magnetic resonance imaging (MRI) in the medical field, the segmentation of whole brain structure is a basic and critical step. Whole brain segmentation of brain structure MRI images divides the entire brain into multiple anatomically related regions of interest, which is a prerequisite for the subsequent quantitative analysis of different brain structures. It assigns a semantic label corresponding to a specific neuroanatomical structure to each voxel, so that quantitative attribute values such as volume, thickness, and area of various regions of interest can be obtained. These attribute values can be used for brain morphology analysis, surgical planning, and postoperative evaluation.

传统的全脑分割算法通常包括两类:Traditional whole-brain segmentation algorithms usually include two categories:

第一类是基于脑图谱的方法,即用可变形配准技术将手动分割的脑图谱标签配准迁移到每个图像中,例如FreeSurfer和BrainSuite,虽然该方法在分割高对比度的组织类型(如白质、灰质和脑脊液)以及皮层下结构(如脑室、小脑和杏仁核)时比较准确,但是将灰质(即皮层)细分为更小区域时,一般表现不佳。这是因为大脑皮层是一个高度褶皱且薄的结构,含有的组织密度又相似,传统的基于灰度值的配准算法很难将不同区域对齐好。而且,这类算法的复杂性和耗时特性更加阻碍了它们在临床应用中的常规使用。The first category is the brain atlas-based method, which uses deformable registration technology to migrate the manually segmented brain atlas labels to each image, such as FreeSurfer and BrainSuite. Although this method is relatively accurate in segmenting high-contrast tissue types (such as white matter, gray matter, and cerebrospinal fluid) and subcortical structures (such as ventricles, cerebellum, and amygdala), it generally performs poorly when subdividing gray matter (ie, cortex) into smaller areas. This is because the cerebral cortex is a highly wrinkled and thin structure, and the tissue density it contains is similar. Traditional grayscale-based registration algorithms have difficulty aligning different regions. Moreover, the complexity and time-consuming nature of such algorithms further hinder their routine use in clinical applications.

第二类全脑分割方法则利用深度学习技术开发端到端的分割网络。这类方法直接以整个大脑图像作为输入,并输出每个体素的标签分类结果。这些网络成功解决了第一类算法的速度问题,满足了临床实时应用的要求。然而,这些端到端的深度学习网络仅依赖输入图像的灰度特征,忽略了大脑解剖结构的复杂几何细节,因此经常会牺牲准确性。The second type of whole-brain segmentation methods use deep learning technology to develop end-to-end segmentation networks. This type of method directly takes the whole brain image as input and outputs the label classification results for each voxel. These networks successfully solve the speed problem of the first type of algorithms and meet the requirements of clinical real-time applications. However, these end-to-end deep learning networks only rely on the grayscale features of the input image, ignoring the complex geometric details of the brain's anatomical structure, and therefore often sacrifice accuracy.

基于此,在本说明书中,提供了一种图像处理方法,本说明书同时涉及一种脑部疾病的计算机辅助诊断方法、一种图像处理装置,一种计算设备,一种计算机可读存储介质以及一种计算机程序产品,在下面的实施例中逐一进行详细说明。Based on this, in this specification, an image processing method is provided. This specification also involves a computer-aided diagnosis method for brain diseases, an image processing device, a computing device, a computer-readable storage medium and a computer program product, which are described in detail one by one in the following embodiments.

参见图1,图1示出了根据本说明书一个实施例提供的一种图像处理方法的流程图,具体包括以下步骤。Referring to FIG. 1 , FIG. 1 shows a flow chart of an image processing method provided according to an embodiment of the present specification, which specifically includes the following steps.

步骤102:接收图像处理任务,其中,所述图像处理任务携带三维脑图像。Step 102: receiving an image processing task, wherein the image processing task carries a three-dimensional brain image.

在实际应用中,本说明书实施例提供的图像处理方法,可以应用于服务器,也可以应用于客户端。在本实施方式中,对图像处理方法的执行主体不做限定。In practical applications, the image processing method provided in the embodiments of this specification can be applied to a server or a client. In this implementation, the execution subject of the image processing method is not limited.

具体的,图像处理任务可以理解为针对三维脑图像的全脑分割任务,其目的是对三维脑图像进行分割,在图像处理任务中携带有三维脑图像,三维脑图像可以理解为通过脑部磁共振成像技术生成的三维图像。Specifically, the image processing task can be understood as a whole-brain segmentation task for three-dimensional brain images, and its purpose is to segment the three-dimensional brain images. The image processing task carries three-dimensional brain images, and the three-dimensional brain images can be understood as three-dimensional images generated by brain magnetic resonance imaging technology.

在实际应用中,某人通过核磁共振仪器拍摄了脑部的三维脑图像,需要针对该三维脑图像进行全脑分割任务。基于该三维脑图像触发图像处理任务,执行终端接收到该图像处理任务,在图像处理任务中携带有三维脑图像。In practical applications, a person takes a three-dimensional brain image of the brain through a nuclear magnetic resonance instrument, and needs to perform a whole brain segmentation task on the three-dimensional brain image. An image processing task is triggered based on the three-dimensional brain image, and the execution terminal receives the image processing task, which carries the three-dimensional brain image.

更进一步的,当本说明书提供的方法是应用于客户端的情况下,用户可以在客户端中基于三维脑图像触发图像处理任务,客户端即可接收到该图像处理任务。当本说明书提供的方法是应用于服务器的情况下,用户可以在客户端中基于三维脑图像生成图像处理任务,并将该图像处理任务发送指服务器,服务器即可接收到该图像处理任务。Furthermore, when the method provided in this specification is applied to a client, the user can trigger an image processing task based on a three-dimensional brain image in the client, and the client can receive the image processing task. When the method provided in this specification is applied to a server, the user can generate an image processing task based on a three-dimensional brain image in the client, and send the image processing task to the server, and the server can receive the image processing task.

步骤104:将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息。Step 104: input the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image.

其中,解剖分割模型可以理解为被预先训练好的机器学习模型,其被训练于基于MRI灰度值和对比度来区分三维脑图像中不同的解剖结构。Among them, the anatomical segmentation model can be understood as a pre-trained machine learning model, which is trained to distinguish different anatomical structures in three-dimensional brain images based on MRI grayscale values and contrast.

在实际应用中,磁共振图像是一种利用核磁共振原理对人体组织进行断层图像成像的医学诊断技术,在MRI图像中,灰度值是根据不同组织的信号强度来表示的。磁共振图像的检查结果通常有两个信号,分别是t1加权像信号和t2加权像信号。在典型的MRI图像中,由于白质、灰质和皮层下结构有相对较高的对比度,可以通过基于对比度信息实现对灰质信息、白质信息、皮层下结构信息的粗分割。In practical applications, magnetic resonance imaging is a medical diagnostic technology that uses the principle of nuclear magnetic resonance to perform tomographic imaging of human tissue. In MRI images, grayscale values are represented according to the signal intensity of different tissues. The examination results of magnetic resonance imaging usually have two signals, namely T1-weighted image signals and T2-weighted image signals. In typical MRI images, since white matter, gray matter and subcortical structures have relatively high contrast, rough segmentation of gray matter information, white matter information and subcortical structure information can be achieved based on contrast information.

在本实施方式中,皮层下结构信息包括皮层下各个脑组织结构的掩膜,皮层下结构包括有脑室、小脑、丘脑、尾状核、壳状核等。In this embodiment, the subcortical structure information includes masks of various subcortical brain tissue structures, and the subcortical structures include ventricles, cerebellum, thalamus, caudate nucleus, putamen, and the like.

在实际应用中,将三维脑图像输入到解剖分割模型中进行初步分割,将三维脑图像粗分割为白质、灰质和皮层下结构(如脑室、小脑、丘脑、尾状核、壳状核等)。在解剖分割模型中仅强调学习灰度值和对比度特征,而无需考虑脑结构中的几何结构特征,可以加快解剖分割模型处理三维脑图像的数据处理速度,进而满足了临床实施应用的要求。In practical applications, the 3D brain image is input into the anatomical segmentation model for preliminary segmentation, and the 3D brain image is roughly segmented into white matter, gray matter and subcortical structures (such as ventricles, cerebellum, thalamus, caudate nucleus, putamen, etc.). In the anatomical segmentation model, only the grayscale value and contrast features are emphasized, without considering the geometric structure features in the brain structure, which can speed up the data processing speed of the anatomical segmentation model to process 3D brain images, thereby meeting the requirements of clinical implementation and application.

在本说明书提供的一具体实施方式中,解剖分割模型通过下述步骤训练获得:In a specific embodiment provided in this specification, the anatomical segmentation model is obtained by training through the following steps:

获取样本三维脑图像和所述样本三维脑图像对应的样本灰质信息、样本白质信息和样本皮层下结构信息;Acquire a sample three-dimensional brain image and sample gray matter information, sample white matter information, and sample subcortical structure information corresponding to the sample three-dimensional brain image;

将所述样本三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于所述样本三维脑图像中的对比度信息输出的预测灰质信息、预测白质信息和预测皮层下结构信息;Inputting the sample three-dimensional brain image into an anatomical segmentation model to obtain predicted gray matter information, predicted white matter information, and predicted subcortical structure information output by the anatomical segmentation model based on contrast information in the sample three-dimensional brain image;

根据所述样本灰质信息和所述预测灰质信息计算灰质损失值,根据所述样本白质信息和所述预测白质信息计算白质损失值,根据所述样本皮层下结构信息和所述预测皮层下结构信息计算皮层下结构损失值;Calculating a gray matter loss value according to the sample gray matter information and the predicted gray matter information, calculating a white matter loss value according to the sample white matter information and the predicted white matter information, and calculating a subcortical structure loss value according to the sample subcortical structure information and the predicted subcortical structure information;

根据所述灰质损失值、所述白质损失值和所述皮层下结构损失值调整所述解剖分割模型的模型参数,并继续训练所述解剖分割模型,直至达到模型训练停止条件。The model parameters of the anatomical segmentation model are adjusted according to the gray matter loss value, the white matter loss value and the subcortical structure loss value, and the anatomical segmentation model is continuously trained until a model training stop condition is reached.

在本说明书实施例提供的方法中,解剖分割模型通过有监督的方式进行训练,其包括有多个训练样本对,在训练样本对中包括有样本三维脑图像和样本三维脑图像对应的样本标签。更进一步的,样本标签具体包括样本灰质信息、样本白质信息和样本皮层下结构。In the method provided in the embodiment of the present specification, the anatomical segmentation model is trained in a supervised manner, which includes a plurality of training sample pairs, wherein the training sample pairs include a sample three-dimensional brain image and a sample label corresponding to the sample three-dimensional brain image. Furthermore, the sample label specifically includes sample gray matter information, sample white matter information, and sample subcortical structure.

在实际应用中,样本三维脑图像和样本三维脑图像可以由人工进行标注,在经过授权后,可以获得样本三维脑图像,由技术人员在样本三维脑图像上进行标注,或将样本三维脑图像输入到样本标注模型中进行处理,可以获得样本标注模型输出的样本标注信息。将样本三维脑图像和样本标注信息组成训练样本对,用于后续对解剖分割模型的模型训练。In practical applications, the sample 3D brain image and the sample 3D brain image can be manually annotated. After authorization, the sample 3D brain image can be obtained, and the technicians can annotate the sample 3D brain image, or the sample 3D brain image can be input into the sample annotation model for processing, and the sample annotation information output by the sample annotation model can be obtained. The sample 3D brain image and the sample annotation information form a training sample pair for subsequent model training of the anatomical segmentation model.

具体的,可以将皮层感兴趣区域标签合并为样本灰质信息。更进一步的,样本灰质信息根据在三维脑图像中的位置关系,还可以进一步分为样本左侧灰质信息和样本右侧灰质信息,在实际应用中,可以根据实际情况来确认样本灰质信息是否需要拆分为样本左侧灰质信息和样本右侧灰质信息。在本实施方式中对此不做限定。Specifically, the cortical region of interest labels can be merged into the sample gray matter information. Furthermore, the sample gray matter information can be further divided into the sample left gray matter information and the sample right gray matter information according to the positional relationship in the three-dimensional brain image. In practical applications, it can be determined whether the sample gray matter information needs to be split into the sample left gray matter information and the sample right gray matter information according to the actual situation. This is not limited in the present embodiment.

将样本三维脑图像输入至解剖分割模型中进行处理,解剖分割模型基于对比度信息、灰度值对样本三维脑图像中进行图像识别,由于灰质、白质和皮层下结构在MRI图像中表现出的较高对比度,样本三维脑图像可以被分为预测灰质信息、预测白质信息和预测皮层下结构信息。The sample three-dimensional brain image is input into the anatomical segmentation model for processing. The anatomical segmentation model performs image recognition on the sample three-dimensional brain image based on contrast information and grayscale values. Due to the high contrast of gray matter, white matter and subcortical structures in MRI images, the sample three-dimensional brain image can be divided into predicted gray matter information, predicted white matter information and predicted subcortical structure information.

此时的解剖分割模型是还未训练好的机器学习模型,需要根据预测出的结果与样本对中的样本标签进行比对,判断解剖分割模型的预测结果与实际结果之间的差距。具体的,预测结果中包括预测灰质信息、预测白质信息和预测皮层下结构信息;样本标签中包括样本灰质信息、样本白质信息和样本皮层下结构信息。则可以根据样本灰质信息和所述预测灰质信息计算灰质损失值,根据样本白质信息和预测白质信息计算白质损失值,根据样本皮层下结构信息和预测皮层下结构信息计算皮层下结构损失值。The anatomical segmentation model at this time is a machine learning model that has not been trained yet. It is necessary to compare the predicted results with the sample labels in the sample pair to determine the gap between the predicted results of the anatomical segmentation model and the actual results. Specifically, the predicted results include predicted gray matter information, predicted white matter information, and predicted subcortical structure information; the sample labels include sample gray matter information, sample white matter information, and sample subcortical structure information. The gray matter loss value can be calculated based on the sample gray matter information and the predicted gray matter information, the white matter loss value can be calculated based on the sample white matter information and the predicted white matter information, and the subcortical structure loss value can be calculated based on the sample subcortical structure information and the predicted subcortical structure information.

计算模型损失值的方法有很多,例如交叉熵损失函数、最大损失函数、平均值损失函数等等,在本说明书提供的实施方式中,对损失函数的具体方式不做限定,以实际应用为准。There are many methods for calculating the model loss value, such as cross entropy loss function, maximum loss function, average loss function, etc. In the implementation method provided in this specification, the specific method of the loss function is not limited and is subject to actual application.

在计算获得了灰质损失值、白质损失值和皮层下结构损失值之后,即可基于这三个损失值对解剖分割模型的模型参数进行调整,具体的,可以是将三个损失值进行拼接融合之后,获得模型损失值,并将模型损失值反向传播,更新解剖分割模型中的模型参数。After calculating the gray matter loss value, white matter loss value and subcortical structure loss value, the model parameters of the anatomical segmentation model can be adjusted based on these three loss values. Specifically, the three loss values can be spliced and fused to obtain the model loss value, and the model loss value can be back-propagated to update the model parameters in the anatomical segmentation model.

在调整完解剖分割模型的模型参数之后,可以重复上述的步骤,继续对解剖分割模型进行训练,直至达到模型训练停止条件。在实际应用中,解剖分割模型的训练停止条件包括:After adjusting the model parameters of the anatomical segmentation model, the above steps can be repeated to continue training the anatomical segmentation model until the model training stop condition is reached. In practical applications, the training stop condition of the anatomical segmentation model includes:

模型损失值小于预设阈值;和/或The model loss value is less than a preset threshold; and/or

训练轮次达到预设的训练轮次。The training rounds reach the preset training rounds.

具体的,对解剖分割模型进行训练的过程中,可以将模型的训练停止条件设置为模型损失值小于预设阈值。具体的,可以是灰质损失值、白质损失值和皮层下结构损失值分别小于各自对应的预设阈值,还可以是三者相加之和小于预设阈值,在本说明书提供的方法中对此不做限定。Specifically, during the training of the anatomical segmentation model, the training stop condition of the model can be set to that the model loss value is less than a preset threshold. Specifically, the gray matter loss value, the white matter loss value, and the subcortical structure loss value can be less than their respective corresponding preset thresholds, or the sum of the three can be less than the preset threshold, which is not limited in the method provided in this specification.

模型的训练停止条件还可以设置为训练轮次为预设的训练轮次,例如,模型训练为20轮。当模型训练轮次达到20次之后,即可认为达到模型训练停止条件。The training stop condition of the model can also be set to a preset training round number, for example, the model training is 20 rounds. When the model training rounds reach 20 times, it can be considered that the model training stop condition is met.

在实际应用中,可以将损失值作为训练停止条件,也可以将训练轮次作为训练停止条件,也可以将损失值和训练轮次同时作为训练停止条件。In practical applications, the loss value can be used as the training stop condition, the training round can be used as the training stop condition, or the loss value and the training round can be used as the training stop condition at the same time.

本说明书实施例提供的方法,训练解剖分割模型,使得解剖分割模型具备仅根据三维脑图像的对比度信息即可完成灰质、白质和皮层下结构的分割,仅利用图像中的灰度信息即可完成分割,获得灰质信息、白质信息,以及脑室、小脑、丘脑、尾状核、壳状核等皮层下结构信息。利用了机器学习模型的高效性,提高了全脑分割的处理速度。The method provided in the embodiments of this specification trains the anatomical segmentation model so that the anatomical segmentation model can complete the segmentation of gray matter, white matter and subcortical structure based only on the contrast information of the three-dimensional brain image, and can complete the segmentation using only the grayscale information in the image, and obtain gray matter information, white matter information, and subcortical structure information such as ventricles, cerebellum, thalamus, caudate nucleus, putamen, etc. The high efficiency of the machine learning model is utilized to improve the processing speed of whole brain segmentation.

步骤106:根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息。Step 106: Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex.

在经过上述步骤,完成粗糙的解剖结构分割之后,获得了灰质信息、白质信息和各个皮层下结构的掩膜。在后续处理过程中,需要用到大脑皮层的几何特征进行更精细的皮层划分。更进一步,是计算大脑皮层上定点的曲率特征信息。After the above steps, the rough anatomical structure segmentation is completed, and the gray matter information, white matter information and masks of each subcortical structure are obtained. In the subsequent processing, the geometric features of the cerebral cortex are needed to perform more detailed cortical division. Furthermore, the curvature feature information of fixed points on the cerebral cortex is calculated.

在本说明书实施例提供的方法中,仅保留了计算曲率做必须的关键步骤,即重建大脑皮层。具体的,在上述处理之后,获得了灰质信息,根据灰质信息可以生成三维脑图像对应的目标大脑皮层。In the method provided in the embodiment of this specification, only the necessary key step of calculating curvature is retained, that is, reconstructing the cerebral cortex. Specifically, after the above processing, gray matter information is obtained, and the target cerebral cortex corresponding to the three-dimensional brain image can be generated according to the gray matter information.

目标大脑皮层是基于三维重建基础生成的。在目标大脑皮层上有很多的沟壑,三维重建后的目标大脑皮层由很多个三角面片拼接而成,在本说明书提供的方法中,要根据大脑皮层的结构信息来进一步进行大脑皮层的分割,其中,大脑皮层的结构信息可以理解为曲率特征信息。因此,在本说明书实施例提供的方法中,根据灰质信息生成目标大脑皮层,并确定目标大脑皮层上的各个顶点的曲率信息,用于后续进行大脑皮层的分割。The target cerebral cortex is generated based on three-dimensional reconstruction. There are many gullies on the target cerebral cortex, and the target cerebral cortex after three-dimensional reconstruction is composed of many triangular facets. In the method provided in this specification, the cerebral cortex is further segmented according to the structural information of the cerebral cortex, wherein the structural information of the cerebral cortex can be understood as curvature feature information. Therefore, in the method provided in the embodiment of this specification, the target cerebral cortex is generated according to the gray matter information, and the curvature information of each vertex on the target cerebral cortex is determined for subsequent segmentation of the cerebral cortex.

在本说明书提供的一具体实施方式中,根据所述灰质信息生成目标大脑皮层,包括:In a specific embodiment provided in this specification, generating a target cerebral cortex according to the gray matter information includes:

基于所述灰质信息进行三维重建,生成初始大脑皮层;Perform three-dimensional reconstruction based on the gray matter information to generate an initial cerebral cortex;

针对所述初始大脑皮层进行迭代平滑处理,生成目标大脑皮层。An iterative smoothing process is performed on the initial cerebral cortex to generate a target cerebral cortex.

在本说明书实施例提供的方法中,为了加快整体流程的处理速度,抛弃了传统处理方式中拓扑修正、球面映射、表面配准等步骤,仅保留了大脑皮层重建的步骤。具体的,先移除小脑、脑干等脑组织结构,并填充其他皮层下区域。基于灰质信息进行三维重建,生成初始大脑皮层。In the method provided in the embodiment of this specification, in order to speed up the processing speed of the overall process, the steps of topological correction, spherical mapping, surface registration, etc. in the traditional processing method are abandoned, and only the step of cerebral cortex reconstruction is retained. Specifically, brain tissue structures such as cerebellum and brainstem are first removed, and other subcortical areas are filled. Three-dimensional reconstruction is performed based on gray matter information to generate the initial cerebral cortex.

快速行进立方体方法(Marching Cubes Algorithm)是一种计算机图形学算法,用于从三维体积数据中抽取等值面,常用于医学成像、地质勘探数据、流体动力学模拟等领域中的三维形状重建。The Marching Cubes Algorithm is a computer graphics algorithm used to extract isosurfaces from three-dimensional volume data. It is often used for three-dimensional shape reconstruction in fields such as medical imaging, geological exploration data, and fluid dynamics simulation.

在灰质信息中包括有各个体素对应的灰度值,体素是体积元素的简称,是三维空间分割上的最小单位。包含体素的立体可以通过立体渲染或者提取给定阈值轮廓的多边形等值面表现出来。体素用于三维成像、科学数据与医学影响等领域。其概念上类似二维空间的最小单位像素。The gray matter information includes the grayscale value corresponding to each voxel. Voxel is the abbreviation of volume element and is the smallest unit in three-dimensional space segmentation. The volume containing voxels can be expressed by stereo rendering or extracting polygonal isosurfaces with given threshold contours. Voxels are used in the fields of three-dimensional imaging, scientific data and medical impact. Its concept is similar to the smallest unit pixel in two-dimensional space.

从各个体素的灰度值中提取出等于或接近预设阈值的表面,这个表面被称为等值面。算法将整个体积空间划分为多个立方体单元,每个立方体中心对应一个采样点,立方体的八个顶点根据邻近采样点的值通过插值确定其标量值。对于每个立方体,算法检查哪些边上的提素质跨过了等值阈值,通过插值确定这些边上的交点,从而形成等值面穿过立方体的轮廓线。在所有交点被确定的情况下,算法会使用预先配置好的连接模式将这些交点连接成一系列的三角面片,形成一个平滑的表面,这些平滑的表面构成了初始大脑皮层。A surface that is equal to or close to a preset threshold is extracted from the grayscale value of each voxel. This surface is called an isosurface. The algorithm divides the entire volume space into multiple cubic units. The center of each cube corresponds to a sampling point, and the eight vertices of the cube determine their scalar values through interpolation based on the values of the adjacent sampling points. For each cube, the algorithm checks which edges have the quality that crosses the isovalue threshold, and determines the intersection points of these edges through interpolation, thereby forming the contour line of the isosurface passing through the cube. When all intersection points are determined, the algorithm uses a pre-configured connection pattern to connect these intersection points into a series of triangular facets to form a smooth surface. These smooth surfaces constitute the initial cerebral cortex.

此时生成的初始大脑皮层,皮层表面还凹凸不平,还可以进一步对初始大脑皮层进行迭代平滑处理,具体的,迭代平滑是用于减少形成的初始大脑皮层的棱角感,使其看起来更加平滑自然。在经过迭代平滑处理之后,即生成了目标大脑皮层。The initial cerebral cortex generated at this time still has an uneven surface, and the initial cerebral cortex can be further iteratively smoothed. Specifically, iterative smoothing is used to reduce the angularity of the initial cerebral cortex, making it look smoother and more natural. After iterative smoothing, the target cerebral cortex is generated.

通过本说明书实施例提供的处理方法,使得整个重建的处理流程耗时不到1秒,简化了大脑皮层的重建流程,极大缩小了重建大脑皮层的处理时间,Through the processing method provided in the embodiment of this specification, the entire reconstruction process takes less than 1 second, which simplifies the reconstruction process of the cerebral cortex and greatly reduces the processing time of reconstructing the cerebral cortex.

目标大脑皮层是由多个三角面片构成的,每个三角面片均包括有3个顶点,每个顶点对应的曲率信息代表了目标大脑皮层的结构特征信息,在本说明书提供的一具体实施方式中,确定所述目标大脑皮层上至少一个顶点的曲率特征信息,包括:The target cerebral cortex is composed of a plurality of triangular facets, each of which includes three vertices, and the curvature information corresponding to each vertex represents the structural characteristic information of the target cerebral cortex. In a specific embodiment provided in this specification, determining the curvature characteristic information of at least one vertex on the target cerebral cortex includes:

计算各顶点与各顶点对应相邻顶点之间的相邻顶点曲率;Calculate the adjacent vertex curvature between each vertex and the corresponding adjacent vertex;

根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率;Calculate the vertex curvature corresponding to each vertex according to the curvature of the adjacent vertices corresponding to each vertex;

将各顶点对应的顶点曲率映射到三维体素空间,确定各顶点对应的曲率特征信息。The vertex curvature corresponding to each vertex is mapped to the three-dimensional voxel space to determine the curvature feature information corresponding to each vertex.

传统的曲率计算方式中,需要用如多项式拟合等方法对曲面建模,再根据一阶和二阶导数计算出该曲面的最大曲率和最小曲率,最后平均最大曲率和最小曲率得到平均曲率,这个过程一般需要2-3分钟,耗时耗力。In the traditional curvature calculation method, it is necessary to use methods such as polynomial fitting to model the surface, and then calculate the maximum curvature and minimum curvature of the surface based on the first-order and second-order derivatives. Finally, the maximum curvature and minimum curvature are averaged to obtain the average curvature. This process generally takes 2-3 minutes and is time-consuming and labor-intensive.

在本说明书实施例提供的方法中,简化了传统的计算方式,首先计算各顶点与相邻顶点之间的相邻顶点曲率,通过各相邻顶点曲率计算各顶点的顶点曲率。In the method provided in the embodiment of the present specification, the traditional calculation method is simplified. First, the adjacent vertex curvature between each vertex and the adjacent vertex is calculated, and the vertex curvature of each vertex is calculated through the adjacent vertex curvature.

此时计算获得各顶点的顶点曲率是三维可视化模型上的顶点,在三维体素空间中会有各个顶点对应的体素,为了后续对皮层进行更精细的分区,还需要将各顶点的曲率信息映射到三维体素空间中,具体的,根据三维可视化模型和三维体素空间中各体素之间的对应关系,将各顶点对应的曲率添加到各顶点对应体素的体素特征信息中,用以表示各体素的曲率特征信息。At this time, the vertex curvature calculated for each vertex is the vertex on the three-dimensional visualization model. There will be voxels corresponding to each vertex in the three-dimensional voxel space. In order to perform a more refined partition of the cortex later, it is also necessary to map the curvature information of each vertex to the three-dimensional voxel space. Specifically, according to the correspondence between the three-dimensional visualization model and the voxels in the three-dimensional voxel space, the curvature corresponding to each vertex is added to the voxel feature information of the voxel corresponding to each vertex to represent the curvature feature information of each voxel.

在本说明书实施例提供的方法中,采用了更精简的曲率计算方法,具体的,计算各顶点与各顶点对应相邻顶点之间的相邻顶点曲率,包括:In the method provided in the embodiment of this specification, a more simplified curvature calculation method is adopted. Specifically, the curvature of adjacent vertices between each vertex and the adjacent vertices corresponding to each vertex is calculated, including:

确定目标顶点和所述目标顶点对应的相邻顶点,其中,所述目标顶点为至少一个顶点中的任意一个;Determine a target vertex and an adjacent vertex corresponding to the target vertex, wherein the target vertex is any one of the at least one vertex;

计算所述目标顶点与各相邻顶点对应的相邻顶点曲率;Calculate the adjacent vertex curvatures corresponding to the target vertex and each adjacent vertex;

相应的,根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率,包括:Accordingly, the vertex curvature corresponding to each vertex is calculated according to the curvature of the adjacent vertices corresponding to each vertex, including:

根据各相邻顶点曲率的均值确定所述目标顶点对应的目标顶点曲率。The target vertex curvature corresponding to the target vertex is determined according to the average of the curvatures of the adjacent vertices.

在实际应用中,在多个顶点中确定目标顶点,以及目标顶点对应的相邻顶点,其中,相邻顶点可以理解为目标顶点所在三角面片上的其他顶点。In practical applications, a target vertex and adjacent vertices corresponding to the target vertex are determined among multiple vertices, wherein the adjacent vertices can be understood as other vertices on the triangle where the target vertex is located.

分别计算目标顶点与各相邻顶点之间的相邻顶点曲率,并根据各相邻顶点曲率的均值确定为目标顶点的顶点曲率。The adjacent vertex curvatures between the target vertex and each adjacent vertex are calculated respectively, and the vertex curvature of the target vertex is determined according to the average value of the curvatures of each adjacent vertex.

例如,某个目标顶点有6个相邻顶点,则目标顶点和各相邻顶点之间有6个相邻顶点曲率,对6个相邻顶点曲率计算均值,并将该均值作为目标顶点对应的目标顶点曲率。For example, if a target vertex has 6 adjacent vertices, there are 6 adjacent vertex curvatures between the target vertex and each adjacent vertex. The average of the 6 adjacent vertex curvatures is calculated and used as the target vertex curvature corresponding to the target vertex.

具体的,计算所述目标顶点与各相邻顶点对应的相邻顶点曲率,包括:Specifically, calculating the adjacent vertex curvatures corresponding to the target vertex and each adjacent vertex includes:

确定所述目标顶点的法向量,其中,所述法向量根据所述目标顶点对应的邻接三角面片的法向量确定;Determine a normal vector of the target vertex, wherein the normal vector is determined according to normal vectors of adjacent triangles corresponding to the target vertex;

根据所述目标顶点和目标相邻顶点的中垂线和所述法向量确定所述目标相邻顶点对应的目标圆心,其中,所述目标相邻顶点为所述目标顶点对应的相邻顶点中的任一个;Determine a target circle center corresponding to the target adjacent vertex according to a perpendicular bisector between the target vertex and the target adjacent vertex and the normal vector, wherein the target adjacent vertex is any one of the adjacent vertices corresponding to the target vertex;

根据所述目标圆心和所述目标顶点确定所述目标相邻顶点对应的相邻顶点曲率。The adjacent vertex curvature corresponding to the target adjacent vertex is determined according to the target circle center and the target vertex.

在实际应用中,以其中一个目标顶点为例进行解释说明,首先确定目标顶点的法向量。目标顶点的法向量有其邻接三角面片的法向量计算获得,在实际应用中,先确定目标顶点对应的邻接三角面片,根据各邻接三角面片的任意两条边的向量叉乘之积,再对结果进行单位正则化处理,获得目标顶点的法向量。In practical applications, we take one of the target vertices as an example to explain, first determine the normal vector of the target vertex. The normal vector of the target vertex is calculated from the normal vectors of its adjacent triangles. In practical applications, first determine the adjacent triangles corresponding to the target vertex, and then perform unit normalization on the result based on the cross product of any two edges of each adjacent triangle to obtain the normal vector of the target vertex.

在确定每个顶点的法向量之后,连接目标顶点与各相邻顶点,就可以得到此目标顶点处经过其每个相邻顶点的局部二维曲线,做目标顶点与目标相邻顶点连线的中垂线,中垂线与法向量的交点即为经过该目标相邻顶点与目标顶点的目标圆心。根据该目标圆心和目标顶点即可计算出该目标圆心与目标顶点对应的半径r,进而可以确定出该目标相邻顶点对应的相邻顶点曲率1/r。After determining the normal vector of each vertex, connect the target vertex with each adjacent vertex to obtain a local two-dimensional curve at the target vertex passing through each of its adjacent vertices. Draw the perpendicular bisector of the line connecting the target vertex and the target adjacent vertex. The intersection of the perpendicular bisector and the normal vector is the target circle center passing through the target adjacent vertex and the target vertex. Based on the target circle center and the target vertex, the radius r corresponding to the target circle center and the target vertex can be calculated, and then the adjacent vertex curvature 1/r corresponding to the target adjacent vertex can be determined.

便利目标顶点对应的所有相邻顶点,并通过上述方法计算出各相邻顶点对应的相邻顶点曲率。最后,再取各相邻顶点曲率的均值,得到目标顶点对应的目标顶点曲率。All adjacent vertices corresponding to the target vertex are conveniently obtained, and the adjacent vertex curvature corresponding to each adjacent vertex is calculated by the above method. Finally, the average of the curvatures of each adjacent vertex is taken to obtain the target vertex curvature corresponding to the target vertex.

在本说明书实施例提供的曲率计算方式中,简化了曲率计算方式。无需根据一阶、二阶导数计算曲面的最大曲率和最小曲率。仅需计算目标顶点与相邻顶点之间的相邻顶点曲率,再根据相邻顶点曲率求均值确定目标顶点的目标顶点曲率。经过实验,该处理方法耗时仅需1.5秒左右,极大减少了计算曲率的计算耗时,节省了计算资源,进一步满足了临床实时应用的要求。In the curvature calculation method provided in the embodiment of this specification, the curvature calculation method is simplified. There is no need to calculate the maximum curvature and minimum curvature of the surface based on the first-order and second-order derivatives. It is only necessary to calculate the adjacent vertex curvature between the target vertex and the adjacent vertex, and then determine the target vertex curvature of the target vertex based on the average of the adjacent vertex curvature. After experiments, this processing method only takes about 1.5 seconds, which greatly reduces the calculation time of curvature, saves computing resources, and further meets the requirements of clinical real-time applications.

步骤108:根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果。Step 108: Determine a cortical surface segmentation result corresponding to the target cerebral cortex based on the gray matter information and the curvature feature information of at least one vertex on the target cerebral cortex.

在获得了各体素的曲率特征信息后,可以将计算得到的曲率特征信息与原始输入的灰度图像在体素空间中进行拼接,通过曲率特征信息为大脑皮层的皮层表面分割提供几何特征的参考。从而根据灰质信息和目标大脑皮层上各顶点的曲率特征信息确定皮层表面分割结果。After obtaining the curvature feature information of each voxel, the calculated curvature feature information can be spliced with the original input grayscale image in the voxel space, and the curvature feature information can be used to provide a reference for geometric features for the cortical surface segmentation of the cerebral cortex. Thus, the cortical surface segmentation result is determined based on the gray matter information and the curvature feature information of each vertex on the target cerebral cortex.

在本说明书提供的一具体实施方式中,根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果,包括:In a specific embodiment provided in this specification, determining a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and the curvature characteristic information of at least one vertex on the target cerebral cortex includes:

根据所述灰质信息对所述三维脑图像进行裁剪,获得三维灰质脑图像;Cropping the three-dimensional brain image according to the gray matter information to obtain a three-dimensional gray matter brain image;

将所述三维灰质脑图像和所述目标大脑皮层上至少一个顶点的曲率特征信息输入至皮层分割模型,获得所述皮层分割模型输出的皮层表面分割结果。The three-dimensional gray matter brain image and curvature feature information of at least one vertex on the target cerebral cortex are input into a cortical segmentation model to obtain a cortical surface segmentation result output by the cortical segmentation model.

在实际应用中,为了更好的对大脑皮层进行分割,排除前一阶段中已经分割出来的白质和皮层下结构的影响。在本实施方式中,使用灰质信息对三维脑图像进行裁剪,保留三维灰质脑图像。In practical applications, in order to better segment the cerebral cortex and eliminate the influence of the white matter and subcortical structures that have been segmented in the previous stage, in this embodiment, the three-dimensional brain image is cropped using the gray matter information to retain the three-dimensional gray matter brain image.

将三维灰质脑图像和各顶点的曲率特征信息进行拼接之后,输入至皮层分割模型进行处理,获得皮层分割模型输出的皮层表面分割结果。After the three-dimensional gray matter brain image and the curvature feature information of each vertex are spliced, they are input into the cortical segmentation model for processing to obtain the cortical surface segmentation result output by the cortical segmentation model.

在三维灰质脑图像中包括各体素点对应的灰度值,同时各顶点的曲率特征信息代表各体素点对应的曲率,将三维灰质脑图像和各顶点的曲率特征信息进行拼接,获得各体素点对应的拼接特征信息,将各体素点对应的拼接特征信息输入至皮层分割模型,获得所述皮层分割模型输出的皮层表面分割结果。The three-dimensional gray matter brain image includes the grayscale value corresponding to each voxel point, and the curvature feature information of each vertex represents the curvature corresponding to each voxel point. The three-dimensional gray matter brain image and the curvature feature information of each vertex are spliced to obtain the splicing feature information corresponding to each voxel point, and the splicing feature information corresponding to each voxel point is input into the cortical segmentation model to obtain the cortical surface segmentation result output by the cortical segmentation model.

皮层分割模型被训练于根据体素点的拼接特征信息预测出各体素点对应的脑结构区域表示,根据各体素点对应的脑结构区域标识可以将大脑皮层分割成不同的皮层表面分割结果。在拼接特征信息中即保留了各体素点的灰度信息,又保留了各体素点对应的几何结构特征,使得皮层分割模型可以同时关注到结构曲率特征,解决了以往深度学习模型中识别准确率低的问题。The cortical segmentation model is trained to predict the brain structure region representation corresponding to each voxel point based on the splicing feature information of the voxel points. According to the brain structure region identification corresponding to each voxel point, the cerebral cortex can be segmented into different cortical surface segmentation results. The splicing feature information retains both the grayscale information of each voxel point and the geometric structure features corresponding to each voxel point, so that the cortical segmentation model can pay attention to the structural curvature features at the same time, solving the problem of low recognition accuracy in previous deep learning models.

在本说明书提供的一具体实施方式中,皮层分割模型通过下述步骤训练获得:In a specific embodiment provided in this specification, the cortical segmentation model is obtained by training through the following steps:

获取样本三维脑图像、所述样本三维脑图像对应的样本灰质信息和所述样本三维脑图像对应的样本皮层表面分割结果;Acquire a sample three-dimensional brain image, sample gray matter information corresponding to the sample three-dimensional brain image, and a sample cortical surface segmentation result corresponding to the sample three-dimensional brain image;

根据所述样本灰质信息生成样本大脑皮层,并确定所述样本大脑皮层上至少一个顶点的样本曲率特征信息;Generate a sample cerebral cortex according to the sample gray matter information, and determine sample curvature characteristic information of at least one vertex on the sample cerebral cortex;

根据所述灰质信息对所述样本三维脑图像进行裁剪,获得样本三维灰质脑图像;Cropping the sample three-dimensional brain image according to the gray matter information to obtain a sample three-dimensional gray matter brain image;

将所述样本三维灰质脑图像和所述样本大脑皮层上至少一个顶点的样本曲率特征信息输入至皮层分割模型,获得所述皮层分割模型输出的预测皮层表面分割结果;Inputting the sample three-dimensional gray matter brain image and sample curvature feature information of at least one vertex on the sample cerebral cortex into a cortical segmentation model to obtain a predicted cortical surface segmentation result output by the cortical segmentation model;

根据所述样本皮层表面分割结果和所述预测皮层表面分割结果计算模型损失值;Calculating a model loss value according to the sample cortical surface segmentation result and the predicted cortical surface segmentation result;

根据所述模型损失值调整所述皮层分割模型的模型参数,并继续训练所述皮层分割模型,直至达到模型训练停止条件。The model parameters of the cortical segmentation model are adjusted according to the model loss value, and the cortical segmentation model is continuously trained until a model training stop condition is reached.

在本说明书实施例提供的方法中,皮层分割模型也是通过有监督的方式进行训练,其包括有样本三维脑图像、样本三维脑图像对应的样本灰质信息和样本皮层表面分割结果。In the method provided in the embodiments of this specification, the cortical segmentation model is also trained in a supervised manner, which includes a sample three-dimensional brain image, sample gray matter information corresponding to the sample three-dimensional brain image, and a sample cortical surface segmentation result.

先根据样本灰质信息进行三维重建,生成样本大脑皮层,并计算大脑皮层上各顶点的样本曲率特征信息。根据样本灰质信息进行三维重建,生成样本大脑皮层,并计算各顶点的样本曲率特征信息的方式,参见上文相关描述,在此不再赘述。First, a three-dimensional reconstruction is performed based on the sample gray matter information to generate a sample cerebral cortex, and the sample curvature characteristic information of each vertex on the cerebral cortex is calculated. The method of performing three-dimensional reconstruction based on the sample gray matter information to generate a sample cerebral cortex and calculate the sample curvature characteristic information of each vertex is described above and will not be repeated here.

根据样本三维灰质脑图像和样本大脑皮层上各顶点对应的样本曲率特征信息输入到皮层分割模型中处理,获得预测皮层表面分割结果的实现方式,也烦请参见上文中的相关描述。Please refer to the above description for the implementation method of inputting the sample curvature feature information corresponding to each vertex on the sample cerebral cortex into the cortical segmentation model for processing to obtain the predicted cortical surface segmentation result.

此时的皮层分割模型是还未训练好的模型,可以根据模型输出的预测皮层表面分割结果和样本皮层表面分割结果计算模型损失值,具体的计算模型损失值的方法可以交叉熵损失函数、最大损失函数、平均值损失函数等等,在本说明书提供的实施方式中,对损失函数的具体方式不做限定,以实际应用为准。At this time, the cortical segmentation model is an untrained model. The model loss value can be calculated based on the predicted cortical surface segmentation results and the sample cortical surface segmentation results output by the model. The specific method for calculating the model loss value can be cross entropy loss function, maximum loss function, average loss function, etc. In the implementation mode provided in this specification, the specific method of the loss function is not limited and is subject to actual application.

在计算了模型损失值之后,可以根据模型损失值调整皮层分割模型的模型参数,并继续用样本数据训练皮层分割模型,直至达到模型训练停止条件。After calculating the model loss value, the model parameters of the cortical segmentation model can be adjusted according to the model loss value, and the cortical segmentation model can continue to be trained with the sample data until the model training stop condition is reached.

在本说明书实施例提供的方法中,在精细的皮层细分阶段,将代表结构特征的曲率特征信息输入到深度神经网络中,使得深度神经网络模型可以参考结构特征,解决了之前深度学习方案中识别准确率低的问题。提升了模型识别的准确率。In the method provided in the embodiment of this specification, in the fine cortical subdivision stage, the curvature feature information representing the structural features is input into the deep neural network, so that the deep neural network model can refer to the structural features, solving the problem of low recognition accuracy in previous deep learning solutions and improving the accuracy of model recognition.

步骤110:根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。Step 110: Determine a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

皮层表面分割结果可以理解为大脑皮层的感兴趣区域标签,使用皮层分割模型输出的灰质掩码对大脑皮层进行裁剪,为每个裁剪区域分配对应的标签。再与之前已经预测出来的白质和皮层下结构进行结合,获得最终的全脑分割结果。The cortical surface segmentation results can be understood as the labels of the regions of interest of the cerebral cortex. The cerebral cortex is cropped using the gray matter mask output by the cortical segmentation model, and a corresponding label is assigned to each cropped region. This is then combined with the previously predicted white matter and subcortical structures to obtain the final whole-brain segmentation results.

本说明书实施例提供的方法,提供了两个分割阶段,在第一个分割阶段进行粗糙的解剖结构的分割,仅根据三维脑图像的对比度信息进行粗滤分割。确定出灰质信息、白质信息和皮层下结构信息。在第二个阶段进一步根据灰质信息生成大脑皮层,并计算大脑皮层上各顶点的曲率信息,通过代表结构特征的曲率信息输入到皮层分割模型中,使得皮层分割模型可以参考结构信息对大脑皮层进行进一步的分割,将大脑皮层划分为多个不同的细分区域,从而获得全脑分割结果。全操作流程仅需10秒左右就可以得到一个脑结构的全脑分割结果,提高处理效率,减少处理时间的情况下,同时提升了处理精度和准确度。在公开的数据集上的实验表明,本说明书实施例提供的方法中,比其他的全脑分割方法,在速度相当的情况下,整体分割准确率可以提升7%左右。The method provided in the embodiment of this specification provides two segmentation stages. In the first segmentation stage, a rough anatomical structure is segmented, and only a coarse filtering segmentation is performed based on the contrast information of the three-dimensional brain image. Gray matter information, white matter information and subcortical structure information are determined. In the second stage, the cerebral cortex is further generated based on the gray matter information, and the curvature information of each vertex on the cerebral cortex is calculated. The curvature information representing the structural features is input into the cortical segmentation model, so that the cortical segmentation model can refer to the structural information to further segment the cerebral cortex, and divide the cerebral cortex into multiple different subdivision areas, thereby obtaining the whole brain segmentation result. The whole operation process only takes about 10 seconds to obtain the whole brain segmentation result of a brain structure, which improves processing efficiency, reduces processing time, and improves processing precision and accuracy at the same time. Experiments on public data sets show that in the method provided in the embodiment of this specification, compared with other whole brain segmentation methods, the overall segmentation accuracy can be improved by about 7% at comparable speeds.

在本说明书提供的一具体实施方式中,还包括:In a specific implementation provided in this specification, it also includes:

在所述全脑分割结果中确定目标脑结构;determining a target brain structure in the whole-brain segmentation result;

针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果。Detection is performed on the target brain structure to obtain segmentation and detection results corresponding to the target brain structure.

在本说明书实施例提供的一具体实施方式中,在获得了全脑分割结果之后,可以根据全脑分割结果进行后续针对脑部结构的诊断。具体的,可以在全脑分割结果中确定目标脑结构,针对目标脑结构进行检测,获得目标脑结构对应的分割和检测结果,并获得相应的检测结论。In a specific implementation provided in the embodiments of this specification, after obtaining the whole brain segmentation result, a subsequent diagnosis of the brain structure can be performed based on the whole brain segmentation result. Specifically, the target brain structure can be determined in the whole brain segmentation result, and the target brain structure can be tested to obtain the segmentation and test results corresponding to the target brain structure, and obtain the corresponding test conclusion.

例如,以目标脑结构为海马体为例进行解释说明,在根据某个用户的磁共振图像获得了全脑分割结果之后,可以从全脑分割结果中确定出该用户的海马体结构,通过对该海马体结构进行检测,判断其是否发生萎缩,从而判断该用户是否患有阿尔兹海默症。For example, taking the target brain structure as the hippocampus as an example, after obtaining the whole-brain segmentation result based on the magnetic resonance image of a certain user, the hippocampus structure of the user can be determined from the whole-brain segmentation result, and by testing the hippocampus structure, it can be determined whether it has atrophied, thereby determining whether the user suffers from Alzheimer's disease.

又例如,某个用户出现步态不稳的症状,可以通过拍摄磁共振图像,并对其磁共振图像进行全脑分割后,识别出该用户的小脑结构,判断该小脑结构是否出现萎缩,从而判断该用户出现的步态不稳的症状是否是小脑萎缩引起的。For another example, if a user has symptoms of unstable gait, we can take a magnetic resonance image and perform whole-brain segmentation on the magnetic resonance image to identify the user's cerebellum structure and determine whether the cerebellum structure has atrophied, thereby determining whether the user's unstable gait symptoms are caused by cerebellar atrophy.

通过本说明书实施例提供的方法,提供了从粗粒度到细粒度的两阶段模型处理方法,将不同类型的分割任务分在两个不同的阶段,使用不同的模型分别对其对应的任务进行处理,提高全脑分割的处理速度,同时还可以提升每个模型的准确率。Through the method provided in the embodiments of this specification, a two-stage model processing method from coarse-grained to fine-grained is provided, which divides different types of segmentation tasks into two different stages, and uses different models to process their corresponding tasks respectively, thereby improving the processing speed of whole-brain segmentation and at the same time improving the accuracy of each model.

在第二个阶段对大脑皮层进行精细分割的过程中,结合了体素点的曲率特征信息和灰度值,从本质上提高了模型对皮层分割任务的理解能力,获得更好的皮层分割结果。In the second stage of fine segmentation of the cerebral cortex, the curvature feature information and grayscale value of the voxel points are combined, which essentially improves the model's ability to understand the cortical segmentation task and obtains better cortical segmentation results.

图2示出了本说明书一实施例提供的一种图像处理方法的处理结构示意图,如图2所示,在本说明书实施例提供的方法中,提供了两个处理阶段,分别为粗分割阶段和精细分割阶段。FIG2 shows a schematic diagram of a processing structure of an image processing method provided in an embodiment of the present specification. As shown in FIG2 , in the method provided in the embodiment of the present specification, two processing stages are provided, namely a coarse segmentation stage and a fine segmentation stage.

在粗分割阶段,将三维脑图像输入至解剖分割模型中进行处理,解剖分割模型根据三维脑图像中各个点的灰度值和对比度信息进行分割,获得解剖分割结果。In the rough segmentation stage, the three-dimensional brain image is input into the anatomical segmentation model for processing. The anatomical segmentation model performs segmentation according to the gray value and contrast information of each point in the three-dimensional brain image to obtain the anatomical segmentation result.

解剖分割结果中包括灰质信息、白质信息和皮层下结构信息。The anatomical segmentation results include gray matter information, white matter information, and subcortical structure information.

在精细分割阶段,根据灰质信息先进行大脑皮层重建,在重建完成之后,计算大脑皮层上每个顶点的曲率信息,再将各顶点的曲率信息映射回体素空间中。将体素空间中的曲率特征信息和灰度信息进行拼接后,输入到皮层分割模型中进行处理,获得皮层表面分割结果。In the fine segmentation stage, the cerebral cortex is reconstructed based on the gray matter information. After the reconstruction is completed, the curvature information of each vertex on the cerebral cortex is calculated, and then the curvature information of each vertex is mapped back to the voxel space. After the curvature feature information and grayscale information in the voxel space are spliced, they are input into the cortical segmentation model for processing to obtain the cortical surface segmentation result.

最后根据皮层表面分割结果和粗分割阶段中生成的白质信息、皮层下结构信息生成全脑分割结果。Finally, the whole brain segmentation result is generated based on the cortical surface segmentation result and the white matter information and subcortical structure information generated in the rough segmentation stage.

本说明书实施例提供的方法,提供了两个分割阶段,在第一个分割阶段进行粗糙的解剖结构的分割,仅根据三维脑图像的对比度信息进行粗滤分割。确定出灰质信息、白质信息和皮层下结构信息。在第二个阶段进一步根据灰质信息生成大脑皮层,并计算大脑皮层上各顶点的曲率信息,通过代表结构特征的曲率信息输入到皮层分割模型中,使得皮层分割模型可以参考结构信息对大脑皮层进行进一步的分割,将大脑皮层划分为多个不同的细分区域,从而获得全脑分割结果。全操作流程仅需10秒左右就可以得到一个脑结构的全脑分割结果,提高处理效率,减少处理时间的情况下,同时提升了处理精度和准确度。The method provided in the embodiment of this specification provides two segmentation stages. In the first segmentation stage, a rough segmentation of the anatomical structure is performed, and a coarse filtering segmentation is performed only based on the contrast information of the three-dimensional brain image. Gray matter information, white matter information and subcortical structure information are determined. In the second stage, the cerebral cortex is further generated based on the gray matter information, and the curvature information of each vertex on the cerebral cortex is calculated. The curvature information representing the structural characteristics is input into the cortical segmentation model, so that the cortical segmentation model can refer to the structural information to further segment the cerebral cortex, and divide the cerebral cortex into multiple different subdivision areas, thereby obtaining the whole brain segmentation result. The whole operation process only takes about 10 seconds to obtain the whole brain segmentation result of a brain structure, which improves processing efficiency, reduces processing time, and improves processing precision and accuracy at the same time.

图3示出了本说明书一实施例提供的一种图像处理系统的架构图,图像处理系统可以包括客户端100和服务端200;FIG3 shows an architecture diagram of an image processing system provided by an embodiment of the present specification. The image processing system may include a client 100 and a server 200;

客户端100,用于向服务端200发送图像处理任务,其中,所述图像处理任务携带三维脑图像;The client 100 is used to send an image processing task to the server 200, wherein the image processing task carries a three-dimensional brain image;

服务端200,用于将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息;根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果;向客户端100发送全脑分割结果;The server 200 is used to input the three-dimensional brain image into the anatomical segmentation model, obtain the gray matter information, white matter information and subcortical structure information output by the anatomical segmentation model based on the contrast information of the three-dimensional brain image; generate a target cerebral cortex according to the gray matter information, and determine the curvature feature information of at least one vertex on the target cerebral cortex; determine the cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and the curvature feature information of at least one vertex on the target cerebral cortex; determine the whole brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information; and send the whole brain segmentation result to the client 100;

客户端100,还用于接收服务端200发送的全脑分割结果。The client 100 is also used to receive the whole brain segmentation result sent by the server 200.

图像处理系统可以包括多个客户端100以及服务端200,其中,客户端100可以称为端侧设备,服务端200可以称为云侧设备。多个客户端100之间通过服务端200可以建立通信连接,在图像处理场景中,服务端200即用来在多个客户端100之间提供图像处理服务,多个客户端100可以分别作为发送端或接收端,通过服务端200实现通信。The image processing system may include multiple clients 100 and a server 200, wherein the client 100 may be referred to as a client-side device and the server 200 may be referred to as a cloud-side device. A communication connection may be established between multiple clients 100 through the server 200. In the image processing scenario, the server 200 is used to provide image processing services between multiple clients 100. Multiple clients 100 may serve as a sender or a receiver respectively, and realize communication through the server 200.

用户通过客户端100可与服务端200进行交互以接收其它客户端100发送的数据,或将数据发送至其它客户端100等。在图像处理场景中,可以是用户通过客户端100向服务端200发布数据流,服务端200根据该数据流生成全脑分割结果,并将全脑分割结果推送至其他建立通信的客户端中。The user can interact with the server 200 through the client 100 to receive data sent by other clients 100, or send data to other clients 100, etc. In the image processing scenario, the user can publish a data stream to the server 200 through the client 100, and the server 200 generates a whole brain segmentation result based on the data stream, and pushes the whole brain segmentation result to other clients that have established communication.

其中,客户端100与服务端200之间通过网络建立连接。网络为客户端100与服务端200之间提供了通信链路的介质。网络可以包括各种连接类型,例如有线、无线通信链路或者光纤电缆等等。客户端100所传输的数据可能需要经过编码、转码、压缩等处理之后才发布至服务端200。The client 100 and the server 200 are connected via a network. The network provides a medium for a communication link between the client 100 and the server 200. The network may include various connection types, such as wired or wireless communication links or optical fiber cables, etc. The data transmitted by the client 100 may need to be encoded, transcoded, compressed, etc. before being released to the server 200.

客户端100可以为浏览器、APP(Application,应用程序)、或网页应用如H5(HyperText Markup Language5,超文本标记语言第5版)应用、或轻应用(也被称为小程序,一种轻量级应用程序)或云应用等,客户端100可以基于服务端200提供的相应服务的软件开发工具包(SDK,Software Development Kit),如基于实时通信(RTC,Real TimeCommunication)SDK开发获得等。客户端100可以部署在电子设备中,需要依赖设备运行或者设备中的某些APP而运行等。电子设备例如可以具有显示屏并支持信息浏览等,如可以是个人移动终端如手机、平板电脑、个人计算机等。在电子设备中通常还可以配置各种其它类应用,例如人机对话类应用、模型训练类应用、文本处理类应用、网页浏览器应用、购物类应用、搜索类应用、即时通信工具、邮箱客户端、社交平台软件等。The client 100 can be a browser, an APP (Application), or a web application such as an H5 (HyperText Markup Language 5, Hypertext Markup Language Version 5) application, or a light application (also known as a mini-program, a lightweight application) or a cloud application, etc. The client 100 can be based on the software development kit (SDK, Software Development Kit) of the corresponding service provided by the server 200, such as based on the real-time communication (RTC, Real Time Communication) SDK development and acquisition. The client 100 can be deployed in an electronic device and needs to rely on the device to run or some APPs in the device to run. For example, the electronic device can have a display screen and support information browsing, such as a personal mobile terminal such as a mobile phone, a tablet computer, a personal computer, etc. Various other types of applications can also be configured in the electronic device, such as human-computer dialogue applications, model training applications, text processing applications, web browser applications, shopping applications, search applications, instant messaging tools, email clients, social platform software, etc.

服务端200可以包括提供各种服务的服务器,例如为多个客户端提供通信服务的服务器,又如为客户端上使用的模型提供支持的用于后台训练的服务器,又如对客户端发送的数据进行处理的服务器等。需要说明的是,服务端200可以实现成多个服务器组成的分布式服务器集群,也可以实现成单个服务器。服务器也可以为分布式系统的服务器,或者是结合了区块链的服务器。服务器也可以是云服务、云数据库、云计算、云函数、云存储、网络服务、云通信、中间件服务、域名服务、安全服务、内容分发网络(CDN,Content DeliveryNetwork)以及大数据和人工智能平台等基础云计算服务的云服务器,或者是带人工智能技术的智能云计算服务器或智能云主机。The server 200 may include servers that provide various services, such as servers that provide communication services to multiple clients, servers for background training that support models used on clients, and servers that process data sent by clients. It should be noted that the server 200 can be implemented as a distributed server cluster consisting of multiple servers, or as a single server. The server can also be a server of a distributed system, or a server combined with a blockchain. The server can also be a cloud server for basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, content distribution networks (CDN, Content Delivery Network), and big data and artificial intelligence platforms, or an intelligent cloud computing server or intelligent cloud host with artificial intelligence technology.

值得说明的是,本说明书实施例中提供的图像处理方法一般由服务端执行,但是,在本说明书的其它实施例中,客户端也可以与服务端具有相似的功能,从而执行本说明书实施例所提供的图像处理方法。在其它实施例中,本说明书实施例所提供的图像处理方法还可以是由客户端与服务端共同执行。It is worth noting that the image processing method provided in the embodiments of this specification is generally executed by the server, but in other embodiments of this specification, the client may also have similar functions to the server, thereby executing the image processing method provided in the embodiments of this specification. In other embodiments, the image processing method provided in the embodiments of this specification may also be executed jointly by the client and the server.

图4是本说明书一实施例提供的应用于云侧设备的图像处理方法的流程示意图,具体包括:FIG4 is a flow chart of an image processing method applied to a cloud-side device provided in an embodiment of the present specification, which specifically includes:

步骤402:接收端侧发送的图像处理任务,其中,所述图像处理任务携带三维脑图像。Step 402: receiving an image processing task sent by the end side, wherein the image processing task carries a three-dimensional brain image.

步骤404:将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息。Step 404: input the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image.

步骤406:根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息。Step 406: Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex.

步骤408:根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果。Step 408: Determine a cortical surface segmentation result corresponding to the target cerebral cortex based on the gray matter information and the curvature feature information of at least one vertex on the target cerebral cortex.

步骤410:根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。Step 410: Determine a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

步骤412:向所述端侧设备发送所述全脑分割结果。Step 412: Send the whole-brain segmentation result to the client-side device.

需要说明的是,步骤402至步骤410的视线方式,与上述步骤102至步骤110的视线方式相同,在本说明书中不再赘述。It should be noted that the line of sight method of step 402 to step 410 is the same as the line of sight method of step 102 to step 110 described above, and will not be described in detail in this specification.

可选的,所述方法还包括:Optionally, the method further includes:

接收端侧发送的目标检测任务,其中,所述目标检测任务中携带有目标脑结构标识;A target detection task sent by a receiving end, wherein the target detection task carries a target brain structure identifier;

根据所述目标脑结构标识在所述全脑分割结果中确定目标脑结构;Determining a target brain structure in the whole-brain segmentation result according to the target brain structure identifier;

针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果;Performing detection on the target brain structure to obtain segmentation and detection results corresponding to the target brain structure;

向所述端侧设备发送目标脑结构对应的分割和检测结果。The segmentation and detection results corresponding to the target brain structure are sent to the end-side device.

在实际应用中,云服务器可以将神经网络模型部署在云端,并为端侧设备提供API调用接口,用户可以通过AIP调用接口向云侧发送图像处理任务。在云侧设备中执行上述图像处理方法,并将最终的全脑分割结果返回给端侧设备。本方法无需在端侧部署神经网络模型,实现了端侧设备的轻量化部署,减少了对端侧设备的硬件要求。提升了用户的使用体验。In actual applications, the cloud server can deploy the neural network model on the cloud and provide an API calling interface for the end-side device. Users can send image processing tasks to the cloud side through the AIP calling interface. The above image processing method is executed in the cloud-side device, and the final whole-brain segmentation result is returned to the end-side device. This method does not require the deployment of a neural network model on the end-side, realizes lightweight deployment of the end-side device, and reduces the hardware requirements for the end-side device. It improves the user experience.

图5是本说明书一实施例提供的一种脑部疾病的计算机辅助诊断方法的流程示意图,具体包括:FIG5 is a flowchart of a computer-aided diagnosis method for brain diseases provided in an embodiment of this specification, which specifically includes:

步骤502:接收脑部疾病检测任务,其中,所述脑部疾病检测任务携带三维脑图像和目标脑结构标识。Step 502: Receive a brain disease detection task, wherein the brain disease detection task carries a three-dimensional brain image and a target brain structure identifier.

步骤504:将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息。Step 504: input the three-dimensional brain image into the anatomical segmentation model to obtain gray matter information, white matter information and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image.

步骤506:根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息。Step 506: Generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex.

步骤508:根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果。Step 508: Determine a cortical surface segmentation result corresponding to the target cerebral cortex based on the gray matter information and the curvature feature information of at least one vertex on the target cerebral cortex.

步骤510:根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。Step 510: Determine a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

步骤512:根据所述目标脑结构标识在所述全脑分割结果中确定目标脑结构,针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果。Step 512: Determine the target brain structure in the whole-brain segmentation result according to the target brain structure identifier, perform detection on the target brain structure, and obtain the segmentation and detection results corresponding to the target brain structure.

本说明书实施例提供的脑部疾病的检测方法,适用于检测由于脑部结构发生病变,需要对脑部结构进行检查的脑部疾病,脑部疾病包括神经退行性疾病的检测,例如阿尔兹海默症(Alzheimer’s disease,AD)、帕金森症(Parkinson disease,PD)小脑萎缩等;脑部疾病还包括如精神类疾病,如精神分裂症、抑郁症等;脑部疾病还包括神经发育类疾病,如自闭症等。通过精确的对三维脑图像进行全脑分割,获得全脑分割结构后,从全脑分割结果中确定目标脑结构,再进一步对目标脑结构进行检测,从而确定用户是否患有相应的脑部疾病。The brain disease detection method provided in the embodiments of this specification is suitable for detecting brain diseases that require examination of brain structures due to pathological changes in brain structures. Brain diseases include the detection of neurodegenerative diseases, such as Alzheimer’s disease (AD), Parkinson’s disease (PD), cerebellar atrophy, etc.; brain diseases also include mental illnesses, such as schizophrenia, depression, etc.; brain diseases also include neurodevelopmental diseases, such as autism, etc. After obtaining the whole-brain segmentation structure by accurately performing whole-brain segmentation on the three-dimensional brain image, the target brain structure is determined from the whole-brain segmentation result, and the target brain structure is further detected to determine whether the user suffers from the corresponding brain disease.

本说明书实施例提供的方法,提供了两个分割阶段,在第一个分割阶段进行粗糙的解剖结构的分割,仅根据三维脑图像的对比度信息进行粗滤分割。确定出灰质信息、白质信息和皮层下结构信息。在第二个阶段进一步根据灰质信息生成大脑皮层,并计算大脑皮层上各顶点的曲率信息,通过代表结构特征的曲率信息输入到皮层分割模型中,使得皮层分割模型可以参考结构信息对大脑皮层进行进一步的分割,将大脑皮层划分为多个不同的细分区域,从而获得全脑分割结果。全操作流程仅需10秒左右就可以得到一个脑结构的全脑分割结果,提高处理效率,减少处理时间的情况下,同时提升了处理精度和准确度。The method provided in the embodiment of this specification provides two segmentation stages. In the first segmentation stage, a rough segmentation of the anatomical structure is performed, and a coarse filtering segmentation is performed only based on the contrast information of the three-dimensional brain image. Gray matter information, white matter information and subcortical structure information are determined. In the second stage, the cerebral cortex is further generated based on the gray matter information, and the curvature information of each vertex on the cerebral cortex is calculated. The curvature information representing the structural characteristics is input into the cortical segmentation model, so that the cortical segmentation model can refer to the structural information to further segment the cerebral cortex, and divide the cerebral cortex into multiple different subdivision areas, thereby obtaining the whole brain segmentation result. The whole operation process only takes about 10 seconds to obtain the whole brain segmentation result of a brain structure, which improves processing efficiency, reduces processing time, and improves processing precision and accuracy at the same time.

与上述方法实施例相对应,本说明书还提供了图像处理装置实施例,图6示出了本说明书一个实施例提供的一种图像处理装置的结构示意图。如图6所示,该装置包括:Corresponding to the above method embodiment, this specification also provides an image processing device embodiment. FIG6 shows a schematic diagram of the structure of an image processing device provided by an embodiment of this specification. As shown in FIG6, the device includes:

接收模块602,被配置为接收图像处理任务,其中,所述图像处理任务携带三维脑图像。The receiving module 602 is configured to receive an image processing task, wherein the image processing task carries a three-dimensional brain image.

第一分割模块604,被配置为将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息。The first segmentation module 604 is configured to input the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image.

生成模块606,被配置为根据所述灰质信息生成目标大脑皮层,并确定所述目标大脑皮层上至少一个顶点的曲率特征信息。The generation module 606 is configured to generate a target cerebral cortex according to the gray matter information, and determine curvature characteristic information of at least one vertex on the target cerebral cortex.

第二分割模块608,被配置为根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果。The second segmentation module 608 is configured to determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex.

确定模块610,被配置为根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。The determination module 610 is configured to determine the whole brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information.

本说明书实施例提供的装置,提供了两个分割阶段,在第一个分割阶段进行粗糙的解剖结构的分割,仅根据三维脑图像的对比度信息进行粗滤分割。确定出灰质信息、白质信息和皮层下结构信息。在第二个阶段进一步根据灰质信息生成大脑皮层,并计算大脑皮层上各顶点的曲率信息,通过代表结构特征的曲率信息输入到皮层分割模型中,使得皮层分割模型可以参考结构信息对大脑皮层进行进一步的分割,将大脑皮层划分为多个不同的细分区域,从而获得全脑分割结果。全操作流程仅需10秒左右就可以得到一个脑结构的全脑分割结果,提高处理效率,减少处理时间的情况下,同时提升了处理精度和准确度。The device provided in the embodiment of this specification provides two segmentation stages. In the first segmentation stage, a rough anatomical structure segmentation is performed, and a coarse filtering segmentation is performed only based on the contrast information of the three-dimensional brain image. Gray matter information, white matter information and subcortical structure information are determined. In the second stage, the cerebral cortex is further generated based on the gray matter information, and the curvature information of each vertex on the cerebral cortex is calculated. The curvature information representing the structural characteristics is input into the cortical segmentation model, so that the cortical segmentation model can refer to the structural information to further segment the cerebral cortex, and divide the cerebral cortex into multiple different subdivision areas, thereby obtaining the whole brain segmentation result. The whole operation process only takes about 10 seconds to obtain the whole brain segmentation result of a brain structure, which improves processing efficiency, reduces processing time, and improves processing precision and accuracy at the same time.

上述为本实施例的一种图像处理装置的示意性方案。需要说明的是,该图像处理装置的技术方案与上述的图像处理方法的技术方案属于同一构思,图像处理装置的技术方案未详细描述的细节内容,均可以参见上述图像处理方法的技术方案的描述。The above is a schematic scheme of an image processing device of this embodiment. It should be noted that the technical scheme of the image processing device and the technical scheme of the above-mentioned image processing method belong to the same concept, and the details not described in detail in the technical scheme of the image processing device can be referred to the description of the technical scheme of the above-mentioned image processing method.

图7示出了根据本说明书一个实施例提供的一种计算设备700的结构框图。该计算设备700的部件包括但不限于存储器710和处理器720。处理器720与存储器710通过总线730相连接,数据库750用于保存数据。Fig. 7 shows a block diagram of a computing device 700 according to an embodiment of the present specification. The components of the computing device 700 include but are not limited to a memory 710 and a processor 720. The processor 720 is connected to the memory 710 via a bus 730, and the database 750 is used to store data.

计算设备700还包括接入设备740,接入设备740使得计算设备700能够经由一个或多个网络760通信。这些网络的示例包括公用交换电话网(PSTN,Public SwitchedTelephone Network)、局域网(LAN,Local Area Network)、广域网(WAN,Wide AreaNetwork)、个域网(PAN,Personal Area Network)或诸如因特网的通信网络的组合。接入设备740可以包括有线或无线的任何类型的网络接口(例如,网络接口卡(NIC,networkinterface controller))中的一个或多个,诸如IEEE802.11无线局域网(WLAN,WirelessLocal Area Network)无线接口、全球微波互联接入(Wi-MAX,WorldwideInteroperability for Microwave Access)接口、以太网接口、通用串行总线(USB,Universal Serial Bus)接口、蜂窝网络接口、蓝牙接口、近场通信(NFC,Near FieldCommunication)。The computing device 700 also includes an access device 740 that enables the computing device 700 to communicate via one or more networks 760. Examples of these networks include a public switched telephone network (PSTN), a local area network (LAN), a wide area network (WAN), a personal area network (PAN), or a combination of communication networks such as the Internet. The access device 740 may include one or more of any type of network interface (e.g., a network interface card (NIC)) of wired or wireless, such as an IEEE 802.11 wireless local area network (WLAN) wireless interface, a world-wide interoperability for microwave access (Wi-MAX) interface, an Ethernet interface, a universal serial bus (USB) interface, a cellular network interface, a Bluetooth interface, and a near field communication (NFC).

在本说明书的一个实施例中,计算设备700的上述部件以及图7中未示出的其他部件也可以彼此相连接,例如通过总线。应当理解,图7所示的计算设备结构框图仅仅是出于示例的目的,而不是对本说明书范围的限制。本领域技术人员可以根据需要,增添或替换其他部件。In one embodiment of the present specification, the above components of the computing device 700 and other components not shown in FIG. 7 may also be connected to each other, for example, through a bus. It should be understood that the computing device structure block diagram shown in FIG. 7 is only for illustrative purposes and is not intended to limit the scope of the present specification. Those skilled in the art may add or replace other components as needed.

计算设备700可以是任何类型的静止或移动计算设备,包括移动计算机或移动计算设备(例如,平板计算机、个人数字助理、膝上型计算机、笔记本计算机、上网本等)、移动电话(例如,智能手机)、可佩戴的计算设备(例如,智能手表、智能眼镜等)或其他类型的移动设备,或者诸如台式计算机或个人计算机(PC,Personal Computer)的静止计算设备。计算设备700还可以是移动式或静止式的服务器。The computing device 700 may be any type of stationary or mobile computing device, including a mobile computer or mobile computing device (e.g., a tablet computer, a personal digital assistant, a laptop computer, a notebook computer, a netbook, etc.), a mobile phone (e.g., a smart phone), a wearable computing device (e.g., a smart watch, smart glasses, etc.), or other types of mobile devices, or a stationary computing device such as a desktop computer or a personal computer (PC). The computing device 700 may also be a mobile or stationary server.

其中,处理器720用于执行如下计算机程序/指令,该计算机程序/指令被处理器执行时实现上述图像处理方法的步骤。The processor 720 is used to execute the following computer program/instructions, which implement the steps of the above-mentioned image processing method when executed by the processor.

本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于计算设备实施例而言,由于其基本相似于图像处理方法实施例,所以描述的比较简单,相关之处参见图像处理方法实施例的部分说明即可。Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the computing device embodiment, since it is basically similar to the image processing method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the image processing method embodiment.

本说明书一实施例还提供一种计算机可读存储介质,其存储有计算机程序/指令,该计算机程序/指令被处理器执行时实现上述图像处理方法的步骤。An embodiment of the present specification further provides a computer-readable storage medium storing a computer program/instruction, which implements the steps of the above-mentioned image processing method when executed by a processor.

本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于计算机可读存储介质实施例而言,由于其基本相似于图像处理方法实施例,所以描述的比较简单,相关之处参见图像处理方法实施例的部分说明即可。Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the computer-readable storage medium embodiment, since it is basically similar to the image processing method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the image processing method embodiment.

本说明书一实施例还提供一种计算机程序产品,包括计算机程序/指令,该计算机程序/指令被处理器执行时实现上述图像处理方法的步骤。An embodiment of the present specification further provides a computer program product, including a computer program/instruction, which implements the steps of the above-mentioned image processing method when executed by a processor.

上述为本实施例的一种计算机程序产品的示意性方案。需要说明的是,该计算机程序产品的技术方案与上述的图像处理方法的技术方案属于同一构思,计算机程序产品的技术方案未详细描述的细节内容,均可以参见上述图像处理方法的技术方案的描述。The above is a schematic solution of a computer program product of this embodiment. It should be noted that the technical solution of the computer program product and the technical solution of the above-mentioned image processing method belong to the same concept, and the details not described in detail in the technical solution of the computer program product can be referred to the description of the technical solution of the above-mentioned image processing method.

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。The above is a description of a specific embodiment of the present specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

所述计算机指令包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据专利实践的要求进行适当的增减,例如在某些地区,根据专利实践,计算机可读介质不包括电载波信号和电信信号。The computer instructions include computer program codes, which may be in source code form, object code form, executable files or some intermediate forms, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium may be appropriately increased or decreased according to the requirements of patent practice. For example, in some regions, according to patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

需要说明的是,上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下, 在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定都是本说明书实施例所必须的。It should be noted that the above is a description of a specific embodiment of the present specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required by the embodiments of the present specification.

在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其它实施例的相关描述。In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

以上公开的本说明书优选实施例只是用于帮助阐述本说明书。可选实施例并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书实施例的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本说明书实施例的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本说明书。本说明书仅受权利要求书及其全部范围和等效物的限制。The preferred embodiments of this specification disclosed above are only used to help explain this specification. The optional embodiments do not describe all the details in detail, nor do they limit the invention to only the specific implementation methods described. Obviously, many modifications and changes can be made according to the content of the embodiments of this specification. This specification selects and specifically describes these embodiments in order to better explain the principles and practical applications of the embodiments of this specification, so that technicians in the relevant technical field can understand and use this specification well. This specification is only limited by the claims and their full scope and equivalents.

Claims (14)

1.一种图像处理方法,包括:1. An image processing method, comprising: 接收图像处理任务,其中,所述图像处理任务携带三维脑图像;receiving an image processing task, wherein the image processing task carries a three-dimensional brain image; 将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image; 根据所述灰质信息生成目标大脑皮层,并确定目标顶点和所述目标顶点对应的相邻顶点,其中,所述目标顶点为至少一个顶点中的任意一个;确定所述目标顶点的法向量,其中,所述法向量根据所述目标顶点对应的邻接三角面片的法向量确定;根据所述目标顶点和目标相邻顶点的中垂线和所述法向量确定所述目标相邻顶点对应的目标圆心,其中,所述目标相邻顶点为所述目标顶点对应的相邻顶点中的任一个;根据所述目标圆心和所述目标顶点确定所述目标相邻顶点对应的相邻顶点曲率;根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率;将各顶点对应的顶点曲率映射到三维体素空间,确定各顶点对应的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine a target vertex and an adjacent vertex corresponding to the target vertex, wherein the target vertex is any one of at least one vertex; determine a normal vector of the target vertex, wherein the normal vector is determined according to a normal vector of an adjacent triangle corresponding to the target vertex; determine a target center point corresponding to the target adjacent vertex according to a perpendicular bisector between the target vertex and the target adjacent vertex and the normal vector, wherein the target adjacent vertex is any one of the adjacent vertices corresponding to the target vertex; determine an adjacent vertex curvature corresponding to the target adjacent vertex according to the target center point and the target vertex; calculate a vertex curvature corresponding to each vertex according to the adjacent vertex curvature corresponding to each vertex; map the vertex curvature corresponding to each vertex to a three-dimensional voxel space, and determine curvature feature information corresponding to each vertex; 根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex; 根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。A whole-brain segmentation result corresponding to the three-dimensional brain image is determined according to the cortical surface segmentation result, the white matter information and the subcortical structure information. 2.如权利要求1所述的方法,根据所述灰质信息生成目标大脑皮层,包括:2. The method according to claim 1, generating a target cerebral cortex according to the gray matter information, comprising: 基于所述灰质信息进行三维重建,生成初始大脑皮层;Perform three-dimensional reconstruction based on the gray matter information to generate an initial cerebral cortex; 针对所述初始大脑皮层进行迭代平滑处理,生成目标大脑皮层。An iterative smoothing process is performed on the initial cerebral cortex to generate a target cerebral cortex. 3.如权利要求1所述的方法,根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率,包括:3. The method according to claim 1, wherein the vertex curvature corresponding to each vertex is calculated according to the curvature of the adjacent vertices corresponding to each vertex, comprising: 根据各相邻顶点曲率的均值确定所述目标顶点对应的目标顶点曲率。The target vertex curvature corresponding to the target vertex is determined according to the average of the curvatures of the adjacent vertices. 4.如权利要求1所述的方法,根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果,包括:4. The method according to claim 1, determining a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and the curvature characteristic information of at least one vertex on the target cerebral cortex, comprising: 根据所述灰质信息对所述三维脑图像进行裁剪,获得三维灰质脑图像;Cropping the three-dimensional brain image according to the gray matter information to obtain a three-dimensional gray matter brain image; 将所述三维灰质脑图像和所述目标大脑皮层上至少一个顶点的曲率特征信息输入至皮层分割模型,获得所述皮层分割模型输出的皮层表面分割结果。The three-dimensional gray matter brain image and curvature feature information of at least one vertex on the target cerebral cortex are input into a cortical segmentation model to obtain a cortical surface segmentation result output by the cortical segmentation model. 5.如权利要求1所述的方法,还包括:5. The method of claim 1, further comprising: 在所述全脑分割结果中确定目标脑结构;determining a target brain structure in the whole-brain segmentation result; 针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果。Detection is performed on the target brain structure to obtain segmentation and detection results corresponding to the target brain structure. 6.如权利要求1所述的方法,解剖分割模型通过下述步骤训练获得:6. The method according to claim 1, wherein the anatomical segmentation model is obtained by training through the following steps: 获取样本三维脑图像和所述样本三维脑图像对应的样本灰质信息、样本白质信息和样本皮层下结构信息;Acquire a sample three-dimensional brain image and sample gray matter information, sample white matter information, and sample subcortical structure information corresponding to the sample three-dimensional brain image; 将所述样本三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于所述样本三维脑图像中的对比度信息输出的预测灰质信息、预测白质信息和预测皮层下结构信息;Inputting the sample three-dimensional brain image into an anatomical segmentation model to obtain predicted gray matter information, predicted white matter information, and predicted subcortical structure information output by the anatomical segmentation model based on contrast information in the sample three-dimensional brain image; 根据所述样本灰质信息和所述预测灰质信息计算灰质损失值,根据所述样本白质信息和所述预测白质信息计算白质损失值,根据所述样本皮层下结构信息和所述预测皮层下结构信息计算皮层下结构损失值;Calculating a gray matter loss value according to the sample gray matter information and the predicted gray matter information, calculating a white matter loss value according to the sample white matter information and the predicted white matter information, and calculating a subcortical structure loss value according to the sample subcortical structure information and the predicted subcortical structure information; 根据所述灰质损失值、所述白质损失值和所述皮层下结构损失值调整所述解剖分割模型的模型参数,并继续训练所述解剖分割模型,直至达到模型训练停止条件。The model parameters of the anatomical segmentation model are adjusted according to the gray matter loss value, the white matter loss value and the subcortical structure loss value, and the anatomical segmentation model is continuously trained until a model training stop condition is reached. 7.如权利要求4所述的方法,皮层分割模型通过下述步骤训练获得:7. The method according to claim 4, wherein the cortical segmentation model is obtained by training through the following steps: 获取样本三维脑图像、所述样本三维脑图像对应的样本灰质信息和所述样本三维脑图像对应的样本皮层表面分割结果;Acquire a sample three-dimensional brain image, sample gray matter information corresponding to the sample three-dimensional brain image, and a sample cortical surface segmentation result corresponding to the sample three-dimensional brain image; 根据所述样本灰质信息生成样本大脑皮层,并确定所述样本大脑皮层上至少一个顶点的样本曲率特征信息;Generate a sample cerebral cortex according to the sample gray matter information, and determine sample curvature characteristic information of at least one vertex on the sample cerebral cortex; 根据所述灰质信息对所述样本三维脑图像进行裁剪,获得样本三维灰质脑图像;Cropping the sample three-dimensional brain image according to the gray matter information to obtain a sample three-dimensional gray matter brain image; 将所述样本三维灰质脑图像和所述样本大脑皮层上至少一个顶点的样本曲率特征信息输入至皮层分割模型,获得所述皮层分割模型输出的预测皮层表面分割结果;Inputting the sample three-dimensional gray matter brain image and sample curvature feature information of at least one vertex on the sample cerebral cortex into a cortical segmentation model to obtain a predicted cortical surface segmentation result output by the cortical segmentation model; 根据所述样本皮层表面分割结果和所述预测皮层表面分割结果计算模型损失值;Calculating a model loss value according to the sample cortical surface segmentation result and the predicted cortical surface segmentation result; 根据所述模型损失值调整所述皮层分割模型的模型参数,并继续训练所述皮层分割模型,直至达到模型训练停止条件。The model parameters of the cortical segmentation model are adjusted according to the model loss value, and the cortical segmentation model is continuously trained until a model training stop condition is reached. 8.一种图像处理方法,应用于云侧设备,包括:8. An image processing method, applied to a cloud-side device, comprising: 接收端侧发送的图像处理任务,其中,所述图像处理任务携带三维脑图像;An image processing task sent by a receiving end, wherein the image processing task carries a three-dimensional brain image; 将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image; 根据所述灰质信息生成目标大脑皮层,并确定目标顶点和所述目标顶点对应的相邻顶点,其中,所述目标顶点为至少一个顶点中的任意一个;确定所述目标顶点的法向量,其中,所述法向量根据所述目标顶点对应的邻接三角面片的法向量确定;根据所述目标顶点和目标相邻顶点的中垂线和所述法向量确定所述目标相邻顶点对应的目标圆心,其中,所述目标相邻顶点为所述目标顶点对应的相邻顶点中的任一个;根据所述目标圆心和所述目标顶点确定所述目标相邻顶点对应的相邻顶点曲率;根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率;将各顶点对应的顶点曲率映射到三维体素空间,确定各顶点对应的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine a target vertex and an adjacent vertex corresponding to the target vertex, wherein the target vertex is any one of at least one vertex; determine a normal vector of the target vertex, wherein the normal vector is determined according to a normal vector of an adjacent triangle corresponding to the target vertex; determine a target center point corresponding to the target adjacent vertex according to a perpendicular bisector between the target vertex and the target adjacent vertex and the normal vector, wherein the target adjacent vertex is any one of the adjacent vertices corresponding to the target vertex; determine an adjacent vertex curvature corresponding to the target adjacent vertex according to the target center point and the target vertex; calculate a vertex curvature corresponding to each vertex according to the adjacent vertex curvature corresponding to each vertex; map the vertex curvature corresponding to each vertex to a three-dimensional voxel space, and determine curvature feature information corresponding to each vertex; 根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex; 根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果;Determining a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information, and the subcortical structure information; 向所述端侧设备发送所述全脑分割结果。The whole-brain segmentation result is sent to the client-side device. 9.如权利要求8所述的方法,还包括:9. The method of claim 8, further comprising: 接收端侧发送的目标检测任务,其中,所述目标检测任务中携带有目标脑结构标识;A target detection task sent by a receiving end, wherein the target detection task carries a target brain structure identifier; 根据所述目标脑结构标识在所述全脑分割结果中确定目标脑结构;Determining a target brain structure in the whole-brain segmentation result according to the target brain structure identifier; 针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果;Performing detection on the target brain structure to obtain segmentation and detection results corresponding to the target brain structure; 向所述端侧设备发送目标脑结构对应的分割和检测结果。The segmentation and detection results corresponding to the target brain structure are sent to the end-side device. 10.一种脑部疾病的计算机辅助诊断方法,包括:10. A computer-aided diagnosis method for brain diseases, comprising: 接收脑部疾病检测任务,其中,所述脑部疾病检测任务携带三维脑图像和目标脑结构标识;receiving a brain disease detection task, wherein the brain disease detection task carries a three-dimensional brain image and a target brain structure identifier; 将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;Inputting the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image; 根据所述灰质信息生成目标大脑皮层,并确定目标顶点和所述目标顶点对应的相邻顶点,其中,所述目标顶点为至少一个顶点中的任意一个;确定所述目标顶点的法向量,其中,所述法向量根据所述目标顶点对应的邻接三角面片的法向量确定;根据所述目标顶点和目标相邻顶点的中垂线和所述法向量确定所述目标相邻顶点对应的目标圆心,其中,所述目标相邻顶点为所述目标顶点对应的相邻顶点中的任一个;根据所述目标圆心和所述目标顶点确定所述目标相邻顶点对应的相邻顶点曲率;根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率;将各顶点对应的顶点曲率映射到三维体素空间,确定各顶点对应的曲率特征信息;Generate a target cerebral cortex according to the gray matter information, and determine a target vertex and an adjacent vertex corresponding to the target vertex, wherein the target vertex is any one of at least one vertex; determine a normal vector of the target vertex, wherein the normal vector is determined according to a normal vector of an adjacent triangle corresponding to the target vertex; determine a target center point corresponding to the target adjacent vertex according to a perpendicular bisector between the target vertex and the target adjacent vertex and the normal vector, wherein the target adjacent vertex is any one of the adjacent vertices corresponding to the target vertex; determine an adjacent vertex curvature corresponding to the target adjacent vertex according to the target center point and the target vertex; calculate a vertex curvature corresponding to each vertex according to the adjacent vertex curvature corresponding to each vertex; map the vertex curvature corresponding to each vertex to a three-dimensional voxel space, and determine curvature feature information corresponding to each vertex; 根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;Determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex; 根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果;Determining a whole-brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information, and the subcortical structure information; 根据所述目标脑结构标识在所述全脑分割结果中确定目标脑结构,针对所述目标脑结构进行检测,获得所述目标脑结构对应的分割和检测结果。The target brain structure is determined in the whole-brain segmentation result according to the target brain structure identifier, and the target brain structure is detected to obtain the segmentation and detection results corresponding to the target brain structure. 11.一种图像处理装置,包括:11. An image processing device, comprising: 接收模块,被配置为接收图像处理任务,其中,所述图像处理任务携带三维脑图像;A receiving module, configured to receive an image processing task, wherein the image processing task carries a three-dimensional brain image; 第一分割模块,被配置为将所述三维脑图像输入至解剖分割模型,获得所述解剖分割模型基于三维脑图像的对比度信息输出的灰质信息、白质信息和皮层下结构信息;A first segmentation module is configured to input the three-dimensional brain image into an anatomical segmentation model to obtain gray matter information, white matter information, and subcortical structure information output by the anatomical segmentation model based on contrast information of the three-dimensional brain image; 生成模块,被配置为根据所述灰质信息生成目标大脑皮层,并确定目标顶点和所述目标顶点对应的相邻顶点,其中,所述目标顶点为至少一个顶点中的任意一个;确定所述目标顶点的法向量,其中,所述法向量根据所述目标顶点对应的邻接三角面片的法向量确定;根据所述目标顶点和目标相邻顶点的中垂线和所述法向量确定所述目标相邻顶点对应的目标圆心,其中,所述目标相邻顶点为所述目标顶点对应的相邻顶点中的任一个;根据所述目标圆心和所述目标顶点确定所述目标相邻顶点对应的相邻顶点曲率;根据各顶点对应的相邻顶点曲率计算各顶点对应的顶点曲率;将各顶点对应的顶点曲率映射到三维体素空间,确定各顶点对应的曲率特征信息;A generation module is configured to generate a target cerebral cortex according to the gray matter information, and determine a target vertex and an adjacent vertex corresponding to the target vertex, wherein the target vertex is any one of at least one vertex; determine a normal vector of the target vertex, wherein the normal vector is determined according to a normal vector of an adjacent triangle corresponding to the target vertex; determine a target center point corresponding to the target adjacent vertex according to a perpendicular bisector between the target vertex and the target adjacent vertex and the normal vector, wherein the target adjacent vertex is any one of the adjacent vertices corresponding to the target vertex; determine an adjacent vertex curvature corresponding to the target adjacent vertex according to the target center point and the target vertex; calculate a vertex curvature corresponding to each vertex according to the adjacent vertex curvature corresponding to each vertex; map the vertex curvature corresponding to each vertex to a three-dimensional voxel space, and determine curvature feature information corresponding to each vertex; 第二分割模块,被配置为根据所述灰质信息和所述目标大脑皮层上至少一个顶点的曲率特征信息,确定所述目标大脑皮层对应的皮层表面分割结果;A second segmentation module is configured to determine a cortical surface segmentation result corresponding to the target cerebral cortex according to the gray matter information and curvature feature information of at least one vertex on the target cerebral cortex; 确定模块,被配置为根据所述皮层表面分割结果、所述白质信息和所述皮层下结构信息,确定所述三维脑图像对应的全脑分割结果。The determination module is configured to determine the whole brain segmentation result corresponding to the three-dimensional brain image according to the cortical surface segmentation result, the white matter information and the subcortical structure information. 12.一种计算设备,包括:12. A computing device comprising: 存储器和处理器;Memory and processor; 所述存储器用于存储计算机程序/指令,所述处理器用于执行所述计算机程序/指令,该计算机程序/指令被处理器执行时实现权利要求1至10任意一项所述方法的步骤。The memory is used to store computer programs/instructions, and the processor is used to execute the computer programs/instructions. When the computer programs/instructions are executed by the processor, the steps of the method described in any one of claims 1 to 10 are implemented. 13.一种计算机可读存储介质,其存储有计算机程序/指令,该计算机程序/指令被处理器执行时实现权利要求1至10任意一项所述方法的步骤。13. A computer-readable storage medium storing a computer program/instruction, wherein the computer program/instruction, when executed by a processor, implements the steps of the method according to any one of claims 1 to 10. 14.一种计算机程序产品,包括计算机程序/指令,该计算机程序/指令被处理器执行时实现权利要求1至10任意一项所述方法的步骤。14. A computer program product, comprising a computer program/instruction, which, when executed by a processor, implements the steps of the method according to any one of claims 1 to 10.
CN202410766838.5A 2024-06-14 2024-06-14 Image processing method and device and computer-aided diagnosis method of brain diseases Active CN118365632B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410766838.5A CN118365632B (en) 2024-06-14 2024-06-14 Image processing method and device and computer-aided diagnosis method of brain diseases

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410766838.5A CN118365632B (en) 2024-06-14 2024-06-14 Image processing method and device and computer-aided diagnosis method of brain diseases

Publications (2)

Publication Number Publication Date
CN118365632A CN118365632A (en) 2024-07-19
CN118365632B true CN118365632B (en) 2024-09-13

Family

ID=91882089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410766838.5A Active CN118365632B (en) 2024-06-14 2024-06-14 Image processing method and device and computer-aided diagnosis method of brain diseases

Country Status (1)

Country Link
CN (1) CN118365632B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101515367A (en) * 2009-04-01 2009-08-26 西北工业大学 Method for segmenting sulus regions on surface of pallium of a three-dimensional cerebral magnetic resonance image
CN105590324A (en) * 2016-02-03 2016-05-18 上海联影医疗科技有限公司 Segmentation method and device of medical images
CN113139964A (en) * 2020-01-20 2021-07-20 上海微创医疗器械(集团)有限公司 Multi-modal image segmentation method and device, electronic equipment and storage medium

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8682045B2 (en) * 1997-02-25 2014-03-25 Wake Forest University Health Sciences Virtual endoscopy with improved image segmentation and lesion detection
JP4022587B2 (en) * 2004-02-23 2007-12-19 国立精神・神経センター総長 Diagnosis support method and apparatus for brain disease
JP6755406B2 (en) * 2017-08-29 2020-09-16 富士フイルム株式会社 Medical image display devices, methods and programs
CN115964594B (en) * 2023-01-10 2023-11-21 广东科学技术职业学院 Method and system for rapidly acquiring vertex curvature of grid model

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101515367A (en) * 2009-04-01 2009-08-26 西北工业大学 Method for segmenting sulus regions on surface of pallium of a three-dimensional cerebral magnetic resonance image
CN105590324A (en) * 2016-02-03 2016-05-18 上海联影医疗科技有限公司 Segmentation method and device of medical images
CN113139964A (en) * 2020-01-20 2021-07-20 上海微创医疗器械(集团)有限公司 Multi-modal image segmentation method and device, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN118365632A (en) 2024-07-19

Similar Documents

Publication Publication Date Title
Liu et al. Perception consistency ultrasound image super-resolution via self-supervised CycleGAN
Ueda et al. Technical and clinical overview of deep learning in radiology
CN111932529B (en) Image classification and segmentation method, device and system
Öztürk et al. Cell‐type based semantic segmentation of histopathological images using deep convolutional neural networks
CN114298997A (en) Method and device for detecting forged picture and storage medium
Garcia-Cantero et al. Neurotessmesh: a tool for the generation and visualization of neuron meshes and adaptive on-the-fly refinement
WO2021184195A1 (en) Medical image reconstruction method, and medical image reconstruction network training method and apparatus
US12303246B2 (en) Methods and systems for precise quantification of human sensory cortical areas
Abd El-Fattah et al. Deep-learning-based super-resolution and classification framework for skin disease detection applications
CN118691813A (en) A three-dimensional medical image segmentation method and medium
CN117408948A (en) Image processing method and training method of image classification segmentation model
Zhao et al. A deep network for joint registration and parcellation of cortical surfaces
Yang et al. BareSkinNet: De‐makeup and De‐lighting via 3D Face Reconstruction
Roullier et al. Automated visual quality assessment for virtual and augmented reality based digital twins
Wang et al. DFC-dehaze: an improved cycle-consistent generative adversarial network for unpaired image dehazing
Athey et al. Hidden Markov modeling for maximum probability neuron reconstruction
CN118365632B (en) Image processing method and device and computer-aided diagnosis method of brain diseases
CN117745725B (en) Image processing method, image processing model training method, three-dimensional medical image processing method, computing device and storage medium
CN114612618A (en) Image generation method, device, equipment, storage medium and computer program product
Tran et al. Fast and high-quality tetrahedral mesh generation from neuroanatomical scans
Tekawade et al. 3d autoencoders for feature extraction in X-ray tomography
Shibu et al. Multi-modal brain image fusion using multi feature guided fusion network
KR101948701B1 (en) Method for determining brain disorder of subject based on latent variables which describe brain structure thereof and apparatus using the same
CN115115900B (en) Image reconstruction model training method, device, equipment, medium and program product
Kumar et al. CycleGAN-based Data Augmentation to Improve Generalizability Alzheimer’s Diagnosis using Deep Learning

Legal Events

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