1.基于逆球极投影的网格球面保角参数化方法,其特征在于,包括如下步骤:步骤一:将大脑医学影像数据图形化,获取数据的网格坐标位置和厚度、所属解剖区域等属性信息,进而获得封闭的大脑网格曲面;
步骤二:对步骤一获取的网格曲面进行切割,切割成两片亏格为一的网格曲面C1、C2步骤三:对步骤二形成的两片网格分别进行平面参数化:以网格的固定边界为单位圆的调和映射,把切缝均匀映射到单位圆上,从而把两个亏格为一的网格分别保角参数化到两个单位圆盘c1,c2上;
步骤四,对步骤三获得的c2进行复反演,使c1,c2覆盖整个扩充复平面C;
步骤五,针对步骤四获得的扩充复平面C进行球极投影的逆变换,使扩充复平面上的点映射到黎曼球面;
步骤六,针对步骤五获得的黎曼球面进行莫比乌斯变换,完成网格球面参数化;从莫比乌斯群里找到使调和能量最小且满足0质心约束的莫比乌斯变换,以使参数化后的网格尽可能保持的原始网格特征。
2.根据权利要求1所述基于逆球极投影的网格球面保角参数化方法,其特征在于,所述的步骤二实现过程如下:
(2.1)最远点采样:对于网格中的点,采用最远点采样算法抽取3个相互距离最远的点N1,N2,N3;
(2.2)初始切缝形成:由最短路径算法获取N1到N2、N2到N3、N3到N1之间的最短路径,并将它们联接在一起形成一个闭合环路作为初始切缝;
(2.3)对切缝的处理:确保在切缝上的任意连续的三个点所构成的三角形不属于网格;
具体为:当连续的三个点所构成的三角形属于网格,且出现在切缝的联接点N1,N2,N3附近,则在切缝中删除三点中处于中间位置的点,然后重新遍历直至不存在连续的三个点所构成的三角形属于网格这种情况;
(2.4)网格切割:沿切缝对网格进行切割,将网格切割成切割成两片亏格为1的曲面C1、C2。
3.根据权利要求1所述基于逆球极投影的网格球面保角参数化方法,其特征在于,所述的步骤三实现过程如下:
(3.1)固定平面参数化的边界:统计切缝点的数量n,将单位圆盘c1,c2的边界即单位圆均分成n段,将切缝点一一映射到边界上;
(3.2)进行Tutte映射:将曲面C1、C2;分别进行固定边界的Tutte映射,边界为(3.1)所述的单位圆;
(3.3)计算调和能量 使调和能量最小,实现参数化的过程中保角,其中α、β、f(u)、f(v)参数关系满足图2所示位置关系;
其中,u、v是原网格上一条边的两个顶点,f(u)、f(v)是u、v经过映射后的点,α、β是原网格上由u、v所构成的边的两个对角;
(3.4)保存平面参数化结果{(Xn,Yn)},其中Xn为由平面参数化结果的横坐标组成的向量,Yn为由平面参数化结果的纵坐标组成的向量,n={1,2}表示两组平面参数化结果。
4.根据权利要求1所述基于逆球极投影的网格球面保角参数化方法,其特征在于,所述的步骤四实现过程如下:
(4.1)坐标转换:将参数化结果{(Xn,Yn)}转换到扩充复平面上,记为复数{Xn+iYn};其中表示虚数;
(4.2)复反演:将{X2+iY2}做关于单位圆的复反演 表示z的共轭复数;经过复反衍{(Xn,Yn)}将无重叠地连续地映射到整个扩充复平面;
其中 表示虚数,x+iy表示一个复数,x为实部,y为虚部;X2,Y2分别由圆盘c2经参数化后实部和虚部;{X2+iY2}为参数化后有序集合。
5.根据权利要求1所述基于逆球极投影的网格球面保角参数化方法,其特征在于,所述的步骤五实现过程如下:
(5.1)球极投影:在球极投影下,黎曼球面等价于扩充复平面;将两份平面参数化结果无重叠地连续地一一映射到球面上;
(5.2)求取球面坐标:由扩充重复平面与球面的对应关系 求得球面参数化的坐标。
6.根据权利要求1所述基于逆球极投影的网格球面保角参数化方法,其特征在于,所述的步骤六实现过程如下:
(6.1)莫比乌斯变换: 其中z,a,b,c,d为满足ad‑bc≠0的扩充复数;
(6.2)0质心约束: 其中σM网格模型M的面元, 是质心;
由于该机算是非线性的,在实际实施中用以下步骤代替:(6.2.1)计算M的质心
(6.2.2)对于所有点v∈M,(6.2.3)对于所有点v∈M,(6.3)计算调和能量,以(6.2)的网格为输入重复步骤二到步骤六,使调和能量最小。