发明名称 一种在医学图像上显示矢量图形的方法
摘要 本发明提供一种在医学图像上显示矢量图形的方法,其根据医学图像和矢量图形边界矩形的相对变化,使矢量图形跟随医学图像变化。使用本发明的方法后,医生可以通过将医学DICOM图像和符合工业标准的DXF格式的人工关节工业模板文件显示在同一个坐标系下,并且通过对人工关节图形的移动和旋转等操作进行精确的型号和大小的匹配。这样既方便又快捷,大大方便了医生对人工关节的选型,对医生的手术起到了辅助和计划的作用。
申请公布号 CN102074026B 申请公布日期 2012.10.03
申请号 CN201010503356.9 申请日期 2010.09.29
申请人 深圳市蓝韵实业有限公司 发明人 程延俊;赵传氐;邓庆林
分类号 G06T11/00(2006.01)I;G06T3/00(2006.01)I;G06F19/00(2011.01)I 主分类号 G06T11/00(2006.01)I
代理机构 深圳市百瑞专利商标事务所(普通合伙) 44240 代理人 金辉
主权项 一种在医学图像上显示矢量图形的方法,其特征在于包括以下步骤:A.建立统一的坐标系;B.在坐标系中装载医学图像;C.在坐标系中装载矢量图形,矢量图形显示在医学图像之上;D.医学图像像素f(xi,yi)(0≤i<N)在坐标系中的坐标为(xi,yi),其中N是医学图像像素的个数;矢量图形在坐标系的图像像素f1(x1i,y1i)(0≤i<N1)的坐标为(x1i,y1i),其中N1是矢量图形在坐标系的图像像素的个数;根据医学图像像素的坐标,得到医学图像区域的边界矩形的大小,边界矩形有以坐标系为参考的图像区域的最左边left、最右边right、最顶端top与最底部bottom;该边界矩形的四个属性分别由以下公式得到:left=min{x0,x1,...,xN‑1};right=max{x0,x1,...,xN‑1};top=max{y0,y1,...,yN‑1};bottom=min{y0,y1,...yN‑1};医学图像的边界矩形为(left,right,top,bottom),矢量图形的边界矩形为(left1,right1,top1,bottom1);E.由医学图像的边界矩形得到该矩形的中心为: <mrow> <mo>(</mo> <mfrac> <mrow> <mo>(</mo> <mi>left</mi> <mo>+</mo> <mi>right</mi> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <mrow> <mi>top</mi> <mo>+</mo> <mi>bottom</mi> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow>由矢量图形的边界矩形得到该矩形的中心为: <mrow> <mo>(</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>left</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>right</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <mrow> <msub> <mi>top</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>bottom</mi> <mn>1</mn> </msub> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow>当医学图像发生移动后,矢量图形的移动方法步骤如下:E1.求出医学图像新的边界矩形的中心相对原矩形中心在坐标系中移动的距离为(x,y),E2.当移动的距离不为0时,矢量图形的边界矩形移动相同的距离,E3.在新的矢量图形的边界矩形内显示新的矢量图形;F.医学图像缩放时,医学图像边界矩形的中心不改变,医学图像缩放率为α,矢量图形的边界矩形与医学图像的边界矩形的中心点的距离为: <mrow> <mi>s</mi> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>left</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>right</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mi>left</mi> <mo>+</mo> <mi>right</mi> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>top</mi> <mo>+</mo> <mi>bottom</mi> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>top</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>bottom</mi> <mn>1</mn> </msub> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mrow> <mn>1</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mrow>缩放之后新的距离为:s1=α*s原矢量图形的中心与医学图像的中心的夹角β为: <mrow> <mi>&beta;</mi> <mo>=</mo> <mi>arctg</mi> <mrow> <mo>(</mo> <mo>|</mo> <mfrac> <mrow> <mfrac> <mrow> <mi>top</mi> <mo>+</mo> <mi>bottom</mi> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>top</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>bottom</mi> <mn>1</mn> </msub> </mrow> <mn>2</mn> </mfrac> </mrow> <mrow> <mfrac> <mrow> <msub> <mi>left</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>right</mi> <mn>1</mn> </msub> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>left</mi> <mo>+</mo> <mi>right</mi> </mrow> <mn>2</mn> </mfrac> </mrow> </mfrac> <mo>|</mo> <mo>)</mo> </mrow> <mo>,</mo> </mrow>则新的矢量图形的边界矩形的中心是: <mrow> <mo>(</mo> <mfrac> <mrow> <mo>(</mo> <mi>left</mi> <mo>+</mo> <mi>right</mi> <mo>)</mo> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>*</mo> <mi>cos</mi> <mi>&beta;</mi> <mo>,</mo> <mfrac> <mrow> <mi>top</mi> <mo>+</mo> <mi>bottom</mi> </mrow> <mn>2</mn> </mfrac> <mo>+</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>*</mo> <mi>sin</mi> <mi>&beta;</mi> <mo>)</mo> </mrow>原矢量图形边界矩形按照同样的方式计算后可以得到新的矢量图形的边界矩形,在新的矢量图形的边界矩形内绘制新的矢量图形。
地址 518000 广东省深圳市福田区景田北路81号碧景园E栋601