专利名称：一种由非平行断层图像序列重建三维目标对象表面的方法
技术领域：
本发明涉及一种图像重建和可视化的方法，尤其涉及一种由非平行断层图像序列重建三
维目标对象表面的方法。 背景技术：
在三维图像重建和可视化技术领域的某些应用场合，人们经由数据采集得到的是三维目 标对象的一系列非平行的断层图像序列，而观察三维目标对象最常用的方式之一是表面显 示，这就要求我们必须能够由非平行断层图像序列重建三维目标对象的表面。这里，基于轮 廓的表面重建是应用最为广泛的方法。以往的基于轮廓的表面重建方法通常是首先将非平行 的断层图像序列进行重组重排后得到一个规则的体数据空间，在此体数据空间中，对若干平 行的断层图像进行手动或半自动分割得到其上目标对象的轮廓，然后再插值获得这些平行断 层图像之间的其它断层图像上的目标对象的轮廓，最后进行基于轮廓的面绘制显示。以往方 法中的重组重排插值处理计算量较大、处理过程较为耗时、方法的工程实用性较差。而且， 以往方法至少需要两次插值计算，降低了方法的重建精度。
据査，目前为止也有少量的研究者对由非平行断层图像序列重建三维目标对象表面问题 进行了研究，如G. M. Treece等人在文章"Fast surface and volume estimation from non-parallel cross-sections, for freehand 3_D ultrasound, 1998，，禾口 "Surface interpolation from sparse cross-sections using region correspondence, 1999"中基 于形状插值方法实现了由非平行断层图像序列重建三维目标对象的表面，但是其重建过程计 算量大、耗时长。A. L. Bogush等人在文章"3D Object Reconstruction from Non-parallel Cross-sections, 2004"中采用三次样条插值方法解决此问题，他们的方法对几何外形较为 简单的三维目标对象能够取得较好的应用效果，但对于几何外形较为复杂的三维目标对象应 用效果却很不理想，因为他们的方法不能正确计算轮廓之间的空间几何关系。
发明内容
针对背景技术所述重建方法的耗时长、精度低、实用性差、难以适应几何外形复杂的三 维目标对象的缺点，本发明提出了一种由原始非平行断层图像序列直接重建三维目标对象表 面的方法，它直接由部分原始断层图像分割得到三维目标对象的轮廓，再插值得到基于轮廓 的面绘制显示所需要的全部轮廓。
一种由非平行断层图像序列重建三维目标对象表面的方法，步骤如下
51) 对部分断层图像进行分割，提取目标对象轮廓
先对部分断层图像进行分割，提取这些断层图像中目标对象的轮廓；在己提取出目标对 象轮廓的断层图像序列中，对每相邻的两帧断层图像都依次执行以下步骤S2)-S6);
52) 对断层图像进行距离变换，得到距离图像
将得到目标对象轮廓的相邻的两帧断层图像转化为二值图像；然后进行Chamfer距离变换， 选取3x3的模板，如图2示；距离变换完毕后，令轮廓外部点的距离值取原来的负数，轮廓 内部点的距离值不变，以区分轮廓内外点，得到的图像为断层图像的距离图像，其中目标对 象轮廓外部的每个点，其像素值表示该点到目标对象轮廓上最近一点距离的负值；内部的每 个点，其像素值表示该点到目标对象轮廓上最近一点的距离值；目标对象轮廓上点的像素值
为零53) 提取距离图像中目标对象轮廓内的极大圆盘集，并去除冗余
提取圆盘集通过计算得到一系列与目标对象轮廓相切且完全在目标对象轮廓内部的圆 盘集，其中的每个圆盘都不会被其余所有圆盘所完全覆盖；这些圆盘集代表了目标对象的骨 架及走向趋势；
首先找出距离图像中目标轮廓内所有像素值的极大值点，即圆盘的中心点对目标对象 轮廓内每个点P的像素值都要与其8个邻点的像素值进行比较判断，若满足
/7^/(户)>/7^/(《)-"则该点即为极大值点；其中点《是点； 上下左右及对角线方向上的8
个相邻点，/为点《到点/ 的相对距离值；对于点/7上下左右四个方向上的四个邻点，/ = 3;
对于点;7对角线方向上的四个邻点，？ = 4;距离图像中的每个点都要进行判断，直到找出所 有的极大值点；每个极大值点都对应一个圆盘，且圆盘中心位于极大值点处；
然后计算每个极大值点对应圆盘的半径，对每个极大值点分别进行以下相同处理；首先 对距离图像进行初始化，除当前处理的极大值点外，其余点像素值都初始化为0;然后对图
像中每个点进行与距离变换相似的两次扫描，取3x3的模板，如图2,每次扫描后得到的像
素值由下式确定
《r = max(《S—
其中v;7为图像上(/,力位置处在第柳m次扫描后的像素值，(t/)为模板相对其中心(o，o)的
位置坐标，而^,/)为模板在(A:,/)处与中心位置的相对距离值；若点(A:,/)位于中心位置的上 下左右四个方向上的四个邻点处时，《、/) = 3;若(、/)位于中心位置的对角线方向上的四
个邻点处时，"A:,/) = 4;经过两次扫描之后的图像，其中只有极大值点附近点的像素值是非
零值，且这些点的像素值小于极大值点的像素值，其它点的像素值是零值；找出其中的最小 非零值w，则该圆盘的半径为r-v-w，其中v为极大值点处的像素值，r即为该极大值点 对应圆盘的半径；重复上述操作直到每个极大值点对应的圆盘半径都被计算出；
对圆盘集去冗余:令M为原始的极大圆盘集，iV为去冗余的极大圆盘集；取出M中半
径最大的圆盘/m，将该圆盘放入iV中；然后对M中剩下的每一个圆盘/"都与圆盘/m进行
比较，取/为两圆盘圆心之间的距离，取/fl^ra^Ay(/w)-0.5xrW^(/")，若_/72</^2，
则把/w从M中去掉；然后重复上述操作，直到M为空；最终得到的W即为去冗余后的圆盘 集；
54) 在断层图像中选择用来插值的圆盘集
插值两相邻断层图像之间的像素值时，将其中一幅图像X中目标对象轮廓内每一个极大 圆盘sm依次向另一幅图像S进行投影(对于非平行断层图像，沿两相邻断层图像的平均法
线方向进行投影)，如图4,判断圆盘^2的投影与被投影图像中目标对象轮廓内包含的任意圆盘训是否有重叠部分，若有重叠部分，则将sm放入集合C,， s"放入集合q，若^"在q中
己经存在，就不用往里放；然后重复上述过程，直到处理完^中所有圆盘；最后C,与(^分
别为两相邻断层图像中插值所用的圆盘集；
S5)计算断层图像中每一点指向图像中目标对象的矢量值
根据步骤S4)得到的两相邻断层图像中用来插值的圆盘集，对两相邻断层图像中的每 一个点fl,进行如下计算，
<formula>formula see original document page 7</formula>4
其中，c表示断层图像中用来插值的圆盘集，y^是圆盘c/的中心，。是圆盘c/的半径，a,是
断层图像中的任意点，S为一矢量值，表示断层图像中点a,到目标对象的方向与距离；结果 即得到与两相邻断层图像相对应的两幅新图像，其中每一点的像素值为一矢量值，表示了该
点到图像中目标对象的方向与距离；
56) 确定插值断层图像之间空间点的距离值
把两相邻断层图像的距离图像中其中一幅j作为出发点，遍历其中每一个点fl，过该点 沿两断层图像的平均法线方向做直线交第二幅断层图像s于一点6 ，然后连接两断层图像中
。点与&点处对应的两个矢量3c—2，得到新矢量''，如图5所示；然后过a点沿!'的方向，做 直线交断层图像S于点c ，线段"c即为a点对应的插值线段；
根据"和C两点处的距离值对线段^上的点g进行余弦插值，公式如下<formula>formula see original document page 7</formula>
其中/，和/2分别为点g沿线段"C方向到点"与点C的距离，《、^分别为两距离图像中fl点 与C点处的距离值，《是插值得到的点g处的距离值，如图6;如果由fl点做直线与断层图
像5无交点时，则舍弃fl点，再对下一点重复如此计算，断层图像^中的每一点都遍历完毕 后，两相邻断层图像之间空间点的距离值即被正确的插值出来了；根据距离值的正负可确定 该点是在轮廓内部还是外部，而零值点即为轮廓上的点；
若已提取出目标对象轮廓的断层图像序列中每相邻两帧断层图像全部处理完毕，则执行 步骤S7)，否则继续执行步骤S2) -S6);
57) 基于轮廓的面绘制显示
对所有的空间点距离值提取零值点，并进行基于轮廓的面绘制显示，整个处理流程结束。其中步骤S2)中的Chamfer距离是距离变换中经典的变换方法，该名称取自木工艺中一 个术语，即倒角，用以表示计算距离时该类方法先后按相反两个方向的两遍扫描这一特点。
一种使用本发明重建方法的自由臂三维超声成像系统，包括微型计算机、超声扫描仪、 定位装置、图像采集设备，其特征在于超声扫描仪通过数字输出端口直接连接到微型计算 机；或者通过模拟接口连接到图像采集设备，图像采集设备通过PCI接口或者USB接口连接 到微型计算机；定位装置的可移动部分固定在超声扫描仪的探头上面，定位装置通过有线或 者无线的方式连接到微型计算机；微型计算机通过数字输入端口或者图像采集设备得到超声 扫描仪的图像信息，通过定位装置得到超声扫描仪探头的位置信息。 本发明成像系统的工作过程如下
101. 采集得到感兴趣人体器官的断层图像序列
操作者自由移动超声扫描仪的探头扫描得到感兴趣人体器官的断层图像序列；
102. 对部分断层图像进行分割提取目标对象轮廓
首先进行分割，提取出部分断层图像中目标对象的轮廓；这里，基于自适应活动轮廓模 型(即T-Snake模型)提取轮廓；
在已提取出目标对象轮廓的断层图像序列中，对每相邻的两帧断层图像都按照本发明步骤依 次执行以下步骤103-107;
103. 对断层图像进行距离变换，得到距离图像
104. 提取距离图像中目标对象轮廓内的极大圆盘集，并去除冗余
提取圆盘集可以计算得到一系列与目标对象轮廓相切且完全在目标对象轮廓内部的圆 盘集，其中的每个原盘都不会被其余所有原盘所覆盖；这些圆盘集代表了目标对象的骨架及 走向趋势；
105. 在断层图像中选择用来插值的圆盘集
106. 计算断层图像中每一点指向图像中目标对象的矢量值
107. 确定插值断层图像之间空间点的距离值
若已提取出目标对象轮廓的断层图像序列中每相邻两帧断层图像全部处理完毕，则执行步骤 108，否则继续执行103-107;
108. 基于轮廓的面绘制显示
对所有的空间点距离值提取零值点，并进行基于轮廓的面绘制显示，整个处理流程结束。 本发明实现了由非平行断层图像序列直接重建三维目标对象表面，重建精度高，速度快， 对几何形状复杂的三维目标对象进行表面重建也具有很好的效果。
图l是本发明方法的流程图。
图2是3X3的Chamfer距离变换掩码模板示意图，它分为前向扫描模板和后向扫描模板。
图3是步骤S3中最终求出的一种极大圆盘集的示意图。
图4是步骤S4中如何选择用来插值的极大圆盘的示意图。图中相邻断层图像中目标对 象轮廓内只有一个圆盘，其中Ksw)， K^)分别是圆盘^^， m的半径。s/是投影后两个圆
盘圆心之间的距离。d(湖)、d(一分别是在图像厶B上圆心连线与轮廓交点处的距离值。"为圆盘的投影方向，即两断层图像的平均法线方向。
图5是步骤S5中计算插值方向的示意图。"、6分别为直线"6与两断层平面的交点，
和c(2)分别是a点与6点处的矢量，/为c(1)、 c(2)两矢量终端的连线。
图6是步骤S6中由插值方向进行插值的示意图。点g为待插值点，^1)与^2)分别为
过g点沿插值方向的直线与两断层图像的交点，/(1)与/(2)分别为点g经插值方向所在直线
到两平面的距离。 具体实施例方式
下面结合附图和实施例对本发明做进一步说明，但不限于此。
实施例
51) 对部分断层图像进行分割，提取目标对象轮廓
先对部分断层图像进行分割，提取这些断层图像中目标对象的轮廓；在己提取出目标对 象轮廓的断层图像序列中，对每相邻的两帧断层图像都依次执行以下步骤S2)-S6);
52) 对断层图像进行距离变换，得到距离图像
将得到目标对象轮廓的相邻的两帧断层图像转化为二值图像；然后进行Chamfer距离变 换，选取3x3的模板，如图2示；距离变换完毕后，令轮廓外部点的距离值取原来的负数， 轮廓内部点的距离值不变，以区分轮廓内外点，得到的图像为断层图像的距离图像，其中目 标对象轮廓外部的每个点，其像素值表示该点到目标对象轮廓上最近一点距离的负值；内部 的每个点，其像素值表示该点到目标对象轮廓上最近一点的距离值；目标对象轮廓上点的像 素值为零值；
53) 提取距离图像中目标对象轮廓内的极大圆盘集，并去除冗余
提取圆盘集通过计算得到一系列与目标对象轮廓相切且完全在目标对象轮廓内部的圆 盘集，其中的每个圆盘都不会被其余所有圆盘所完全覆盖；这些圆盘集代表了目标对象的骨
架及走向趋势；
首先找出距离图像中目标轮廓内所有像素值的极大值点，即圆盘的中心点对目标对象 轮廓内每个点p的像素值都要与其8个邻点的像素值进行比较判断，若满足
/7^/07)>/ ^%)-/，则该点即为极大值点；其中点《是点/ 上下左右及对角线方向上的8 个相邻点，/为点g到点p的相对距离值；对于点p上下左右四个方向上的四个邻点，？ = 3;
对于点P对角线方向上的四个邻点，M4;距离图像中的每个点都要进行判断，直到找出所
有的极大值点；每个极大值点都对应一个圆盘，且圆盘中心位于极大值点处；
然后计算每个极大值点对应圆盘的半径，对每个极大值点分别进行以下相同处理；首先 对距离图像进行初始化，除当前处理的极大值点外，其余点像素值都初始化为0;然后对图
像中每个点进行与距离变换相似的两次扫描，取3x3的模板，如图2，每次扫描后得到的像
素值由下式确定《、max:(c:;:,-狀o)
其中v;j'为图像上(/,_/)位置处在第次扫描后的像素值，&，/)为模板相对其中心(o,o)的
位置坐标，而,(A:，/)为模板在(A:,/)处与中心位置的相对距离值；若点(A:,/)位于中心位置的上 下左右四个方向上的四个邻点处时，"A，/) = 3;若(A,/)位于中心位置的对角线方向上的四
个邻点处时，/(A，/) = 4;经过两次扫描之后的图像，其中只有极大值点附近点的像素值是非
零值，且这些点的像素值小于极大值点的像素值，其它点的像素值是零值；找出其中的最小 非零值w，则该圆盘的半径为^ = "w,其中v为极大值点处的像素值，r即为该极大值点 对应圆盘的半径；重复上述操作直到每个极大值点对应的圆盘半径都被计算出；
对圆盘集去冗余:令M为原始的极大圆盘集，W为去冗余的极大圆盘集；取出M中半
径最大的圆盘/m，将该圆盘放入iV中；然后对M中剩下的每一个圆盘/w都与圆盘/m进行 比较，取./7为两圆盘圆心之间的距离，取/t^rac/^(/m) —0.5xraA^(/"),若y72 </c/2 ，
则把/"从M中去掉；然后重复上述操作，直到M为空；最终得到的W即为去冗余后的圆盘
54) 在断层图像中选择用来插值的圆盘集
插值两相邻断层图像之间的像素值时，将其中一幅图像^中目标对象轮廓内每一个极大 圆盘sm依次向另一幅图像5进行投影(对于非平行断层图像，沿两相邻断层图像的平均法 线方向进行投影)，如图4，判断圆盘sm的投影与被投影图像中目标对象轮廓内包含的任意
圆盘训是否有重叠部分，若有重叠部分，则将^2放入集合C,， s"放入集合q，若训在G中
己经存在，就不用往里放；然后重复上述过程，直到处理完J中所有圆盘；最后C,与&分
别为两相邻断层图像中插值所用的圆盘集；
55) 计算断层图像中每一点指向图像中目标对象的矢量值
根据步骤S4)得到的两相邻断层图像中用来插值的圆盘集，对两相邻断层图像中的每 一个点a,进行如下计算，
二LV 4
《
4
其中，C表示断层图像中用来插值的圆盘集，、是圆盘J的中心，^是圆盘^/的半径，fl,是
断层图像中的任意点，5为一矢量值，表示断层图像中点",到目标对象的方向与距离；结果 即得到与两相邻断层图像相对应的两幅新图像，其中每一点的像素值为一矢量值，表示了该点到图像中目标对象的方向与距离；
S6)确定插值断层图像之间空间点的距离值
把两相邻断层图像的距离图像中其中一幅^作为出发点，遍历其中每一个点fl，过该点 沿两断层图像的平均法线方向做直线交第二幅断层图像B于一点6 ，然后连接两断层图像中
"点与6点处对应的两个矢量3c"2，得到一新矢量/如图5所示；然后过fl点沿Z'的方向，做
直线交断层图像5于点c，线段ac即为a点对应的插值线段；
根据"和C两点处的距离值对线段flC上的点g进行余弦插值，公式如下
<formula>formula see original document page 11</formula>
其中/,和/2分别为点g沿线段w方向到点a与点c的距离，《、《分别为两距离图像中a点 与c点处的距离值，《是插值得到的点g处的距离值，如图6;如果由"点做直线与断层图
像^无交点时，则舍弃"点，再对下一点重复如此计算，断层图像」中的每一点都遍历完毕 后，两相邻断层图像之间空间点的距离值即被正确的插值出来了；根据距离值的正负可确定 该点是在轮廓内部还是外部，而零值点即为轮廓上的点；
若己提取出目标对象轮廓的断层图像序列中每相邻两帧断层图像全部处理完毕，则执行 步骤S7),否则继续执行步骤S2) -S6);
S7)基于轮廓的面绘制显示
对所有的空间点距离值提取零值点，并进行基于轮廓的面绘制显示，整个处理流程结束。
权利要求
1、一种由非平行断层图像序列重建三维目标对象表面的方法，步骤如下S1)对部分断层图像进行分割，提取目标对象轮廓先对部分断层图像进行分割，提取这些断层图像中目标对象的轮廓；在已提取出目标对象轮廓的断层图像序列中，对每相邻的两帧断层图像都依次执行以下步骤S2)-S6)；S2)对断层图像进行距离变换，得到距离图像将得到目标对象轮廓的相邻的两帧断层图像转化为二值图像；然后进行Chamfer距离变换，选取3×3的模板，距离变换完毕后，令轮廓外部点的距离值取原来的负数，轮廓内部点的距离值不变，以区分轮廓内外点，得到的图像为断层图像的距离图像，其中目标对象轮廓外部的每个点，其像素值表示该点到目标对象轮廓上最近一点距离的负值；内部的每个点，其像素值表示该点到目标对象轮廓上最近一点的距离值；目标对象轮廓上点的像素值为零值；S3)提取距离图像中目标对象轮廓内的极大圆盘集，并去除冗余提取圆盘集通过计算得到一系列与目标对象轮廓相切且完全在目标对象轮廓内部的圆盘集，其中的每个圆盘都不会被其余所有圆盘所完全覆盖；这些圆盘集代表了目标对象的骨架及走向趋势；首先找出距离图像中目标轮廓内所有像素值的极大值点，即圆盘的中心点对目标对象轮廓内每个点p的像素值都要与其8个邻点的像素值进行比较判断，若满足pixel(p)>pixel(q)-t，则该点即为极大值点；其中点q是点p上下左右及对角线方向上的8个相邻点，t为点q到点p的相对距离值；对于点p上下左右四个方向上的四个邻点，t＝3；对于点p对角线方向上的四个邻点，t＝4；距离图像中的每个点都要进行判断，直到找出所有的极大值点；每个极大值点都对应一个圆盘，且圆盘中心位于极大值点处；然后计算每个极大值点对应圆盘的半径，对每个极大值点分别进行以下相同处理；首先对距离图像进行初始化，除当前处理的极大值点外，其余点像素值都初始化为0；然后对图像中每个点进行与距离变换相似的两次扫描，取3×3的模板，每次扫描后得到的像素值由下式确定其中为图像上(i，j)位置处在第num次扫描后的像素值，(k，l)为模板相对其中心(0，0)的位置坐标，而t(k，l)为模板在(k，l)处与中心位置的相对距离值；若点(k，l)位于中心位置的上下左右四个方向上的四个邻点处时，t(k，l)＝3；若(k，l)位于中心位置的对角线方向上的四个邻点处时，t(k，l)＝4；经过两次扫描之后的图像，其中只有极大值点附近点的像素值是非零值，且这些点的像素值小于极大值点的像素值，其它点的像素值是零值；找出其中的最小非零值w，则该圆盘的半径为r＝v-w，其中v为极大值点处的像素值，r即为该极大值点对应圆盘的半径；重复上述操作直到每个极大值点对应的圆盘半径都被计算出；对圆盘集去冗余令M为原始的极大圆盘集，N为去冗余的极大圆盘集；取出M中半径最大的圆盘fm，将该圆盘放入N中；然后对M中剩下的每一个圆盘fn都与圆盘fm进行比较，取fl为两圆盘圆心之间的距离，取fd＝radius(fm)-0.5×radius(fn)，若fl2<fd2，则把fn从M中去掉；然后重复上述操作，直到M为空；最终得到的N即为去冗余后的圆盘集；S4)在断层图像中选择用来插值的圆盘集插值两相邻断层图像之间的像素值时，将其中一幅图像A中目标对象轮廓内每一个极大圆盘sm依次向另一幅图像B进行投影，判断圆盘sm的投影与被投影图像中目标对象轮廓内包含的任意圆盘sn是否有重叠部分，若有重叠部分，则将sm放入集合C1，sn放入集合C2，若sn在C2中已经存在，就不用往里放；然后重复上述过程，直到处理完A中所有圆盘；最后C1与C2分别为两相邻断层图像中插值所用的圆盘集；S5)计算断层图像中每一点指向图像中目标对象的矢量值根据步骤S4)得到的两相邻断层图像中用来插值的圆盘集，对两相邻断层图像中的每一个点ai进行如下计算，其中，c表示断层图像中用来插值的圆盘集，kd是圆盘d的中心，rd是圆盘d的半径，ai是断层图像中的任意点，为一矢量值，表示断层图像中点ai到目标对象的方向与距离；结果即得到与两相邻断层图像相对应的两幅新图像，其中每一点的像素值为一矢量值，表示了该点到图像中目标对象的方向与距离；S6)确定插值断层图像之间空间点的距离值把两相邻断层图像的距离图像中其中一幅A作为出发点，遍历其中每一个点a，过该点沿两断层图像的平均法线方向做直线交第二幅断层图像B于一点b，然后连接两断层图像中a点与b点处对应的两个矢量得到新矢量i，然后过a点沿i的方向，做直线交断层图像B于点c，线段ac即为a点对应的插值线段；根据a和c两点处的距离值对线段ac上的点g进行余弦插值，公式如下其中l1和l2分别为点g沿线段ac方向到点a与点c的距离，d1、d2分别为两距离图像中a点与c点处的距离值，dg是插值得到的点g处的距离值；如果由a点做直线与断层图像B无交点时，则舍弃a点，再对下一点重复如此计算，断层图像A中的每一点都遍历完毕后，两相邻断层图像之间空间点的距离值即被正确的插值出来了；根据距离值的正负可确定该点是在轮廓内部还是外部，而零值点即为轮廓上的点；若已提取出目标对象轮廓的断层图像序列中每相邻两帧断层图像全部处理完毕，则执行步骤S7)，否则继续执行步骤S2)-S6)；S7)基于轮廓的面绘制显示对所有的空间点距离值提取零值点，并进行基于轮廓的面绘制显示，整个处理流程结束。
全文摘要
一种由非平行断层图像序列重建三维目标对象表面的方法，属三维图像重建和可视化技术领域，通过对断层图像分割提取目标对象轮廓，并对其进行距离变换、提取圆盘集并去除冗余找出目标对象轮廓内的极大圆盘集，在相邻断层图像中选择用来插值的圆盘集，计算断层图像中每一点指向图像中目标对象的矢量值，确定插值断层图像之间空间点的距离值，并进行基于轮廓的面绘制显示。本发明实现了由非平行断层图像序列直接重建三维目标对象表面，重建精度高，速度快，对几何形状复杂的三维目标对象进行表面重建也具有很好的效果。
文档编号G06T11/00GK101533518SQ200910020718
公开日2009年9月16日 申请日期2009年4月21日 优先权日2009年4月21日
发明者孙丰荣, 张明强, 王丽梅, 王庆浩, 王文明 申请人:山东大学