C#开发PACS医学影像处理系统(十三):绘图处理之病灶测量

接上一篇文章,当咱们能够绘制图形标记后,就能够在此操做类上面进行扩展,html

好比测量类工具,目前整理出的经常使用绘图和测量功能以下:git

测量工具类:(图形标记类请参考本系列文章:绘图处理之图形标记算法

功能 说明
标尺 线段长度测量
折尺 曲线长度测量 
心胸比 两根线段按比例测量 
交叉尺 两根线段互相垂直测量 
Cobb角 两根线段的垂线交点角度测量,通常用于脊柱曲率 
开角 两根线段的延长交点的角度测量
T型尺 两根线段相交的任意角度摇摆测量和长度测量
角度测量 角度测量
 圆形测量 圆形或椭圆面积测量,包括CT值,平均值,方差,最大最小值
矩形测量  矩形面积测量,包括CT值,平均值,方差,最大最小值
多边形测量  多边形面积

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

以线段测量为例,看效果:canvas

 

 

 在绘制图形的基础上,咱们建立一个Text Block来显示测量结果:微信

     //文本对象
        TextBlock txtMeasure;

        /// <summary>
        /// 建立测量结果文本
        /// </summary>
        /// <param name="point"></param>
        public void CreateMeasure(Point point)
        {
            if (!isMeasure)
            {
                return;
            }

            txtMeasure = new TextBlock();
            txtMeasure.Text = "0.0mm";
            txtMeasure.FontSize = ShapeManager.shapeMeasureFontSize;
            txtMeasure.Foreground = ShapeManager.shapeMeasureColor;
            txtMeasure.Height = 30;
            txtMeasure.Width = 100;
            txtMeasure.SetValue(Canvas.LeftProperty, point.X);
            txtMeasure.SetValue(Canvas.TopProperty, point.Y);
            canvas.Children.Add(txtMeasure);
            measureList.Add(txtMeasure);

        }

计算两点之间的距离:其中dpi是缩放比例架构

      /// <summary>
        /// 求平面中两点之间距离
        /// </summary>
        /// <param name="p1">点1</param>
        /// <param name="p2">点2</param>
        /// <returns></returns>
        public static double GetDistance(Point p1, Point p2)
        {
            double result = 0;
            result = Math.Sqrt((p1.X * dpiX - p2.X * dpiX) * (p1.X * dpiX - p2.X * dpiX) + (p1.Y * dpiY - p2.Y * dpiY) * (p1.Y * dpiY - p2.Y * dpiY));
            return result;
        }

 

一些角度类型的测量计算相对复杂一点,须要计算角度大小和绘制弧线并发

 

     /// <summary>
        /// 弧线画笔对象
        /// </summary>
        public Stroke circleStroke;


        /// <summary>
        /// 绘制角度弧线
        /// </summary>
        private void DrawAxiesCircle(Point point0, Point point1, Point point2)
        {

            if (inkCanvas.Strokes.Count > 0 && circleStroke != null)
            {
                if (inkCanvas.Strokes.Contains(circleStroke))
                {
                    inkCanvas.Strokes.Remove(circleStroke);
                }
            }

            //计算角度
            double a = Math.Sqrt((point1.X - point2.X) * (point1.X - point2.X) + (point1.Y - point2.Y) * (point1.Y - point2.Y));
            double b = Math.Sqrt((point1.X - point0.X) * (point1.X - point0.X) + (point1.Y - point0.Y) * (point1.Y - point0.Y));
            double c = Math.Sqrt((point2.X - point0.X) * (point2.X - point0.X) + (point2.Y - point0.Y) * (point2.Y - point0.Y));
            double cTheta = (a * a + b * b - c * c) / (2 * a * b);
            double theta = Math.Acos(cTheta) * 180 / Math.PI;

            //绘制弧线
            double r = 30;
            r = a > b ? b : a;
            double rMax = a;
            if (rMax > b)
            {
                rMax = b;
            }
            if (r > 0.5 * rMax)
            {
                r = 0.5 * rMax;
            }

            double theta0 = Math.Atan((point1.Y - point2.Y) / (point2.X - point1.X + 1e-10)) * 180 / Math.PI;
            if (point1.X > point2.X)
            {
                theta0 = 180 + theta0;
            }
            List<Point> pointList = new List<Point>();
            double sin_ab = ((point2.X - point1.X) * (point0.Y - point1.Y) - (point2.Y - point1.Y) * (point0.X - point1.X)) / (a * b); ;
            if (sin_ab <= 0)
            {

                if (theta < 1)
                {
                    for (double delta = 0.0001; delta <= theta;)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                        delta = delta + 0.0001;
                    }
                }
                else if (theta > 1 && theta < 20)
                {
                    for (double delta = 0.01; delta <= theta;)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                        delta = delta + 0.01;
                    }
                }
                else
                {
                    for (double delta = 0; delta <= theta; delta++)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                    }
                }
            }
            else
            {
                if (theta < 1)
                {
                    for (double delta = -theta; delta <= 0;)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                        delta = delta + 0.0001;
                    }
                }
                else if (theta > 1 && theta < 20)
                {
                    for (double delta = -theta; delta <= 0;)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                        delta = delta + 0.01;
                    }
                }
                else
                {
                    for (double delta = -theta; delta <= 0; delta++)
                    {
                        double th = delta + theta0;
                        pointList.Add(new Point(point1.X + r * Math.Cos(th * Math.PI / 180), point1.Y - r * Math.Sin(th * Math.PI / 180)));
                    }
                }
            }
            if (pointList.Count > 0)
            {
                StylusPointCollection point = new StylusPointCollection(pointList);
                circleStroke = new Stroke(point)
                {
                    DrawingAttributes = inkCanvas.DefaultDrawingAttributes.Clone(),
                };
                inkCanvas.Strokes.Add(circleStroke);
                txt.Text = theta.ToString(ShapeManager.measureDigit) + "°";
            }
            ReSetAnglePoint();
        }

在圆形和矩形的面积测量中,咱们可使用GetArea()方法来获取图形面积:机器学习

      /// <summary>
        /// 计算面积测量结果
        /// </summary>
        public void CalculateMeasure()
        {
            if (!isMeasure)
            {
                return;
            }
            Ellipse ellipse = (Ellipse)shape;
            ellipse.UpdateLayout();
            txtMeasure.Text = Math.Sqrt(ellipse.RenderedGeometry.GetArea()).ToString(ShapeManager.measureDigit) + "mm²";
            ReSetMeasurePoint();
        }

其余一些须要注意的细节:高并发

1.当单元格放大缩小时,图形也要从新计算各个控制点的位置来同步放大或缩小工具

 

监听画布大小变化事件:

private void ToolInkCanvas_SizeChanged(object sender, SizeChangedEventArgs e)

从新绘制元素:

      /// <summary>
        /// 从新绘制画布元素大小
        /// </summary>
        public void ReSetShapeSize()
        {
            for (int i = 0; i < shapeManager.shapeList.Count; i++)
            {
                sizeScaleX = ToolInkCanvas.ActualWidth / shapeManager.shapeList[i].cvsWidth;
                sizeScaleY = ToolInkCanvas.ActualHeight / shapeManager.shapeList[i].cvsHeight;

                shapeManager.shapeList[i].cvsWidth = ToolInkCanvas.ActualWidth;
                shapeManager.shapeList[i].cvsHeight = ToolInkCanvas.ActualHeight;

                if (shapeManager.shapeList[i] is TextInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "text", shapeManager.shapeList[i]);
                }
                else if (shapeManager.shapeList[i] is CrossRulerInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "crossRuler", shapeManager.shapeList[i]);
                }
                else if (shapeManager.shapeList[i] is AngleRulerInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "angle", shapeManager.shapeList[i]);
                }
                else if (shapeManager.shapeList[i] is TRulerInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "tRulerInfo", shapeManager.shapeList[i]);
                }
                else if (shapeManager.shapeList[i] is CobbAngleInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "cobbAngleInfo", shapeManager.shapeList[i]);
                }
                else if (shapeManager.shapeList[i] is OpenAngleInfo)
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY, "openAngleInfo", shapeManager.shapeList[i]);
                }
                else
                {
                    shapeManager.shapeList[i].ReSetScalePoints(sizeScaleX, sizeScaleY);
                }
            }
        }

2.文字可能会挡住图像,要能分离并拖动到其余位置,监听鼠标事件和从新设置位置便可。

 3.标记的复制与粘贴,思路是将List集合中用户所选定的标记插入到另外一个单元格的List集合中,并支持Ctrl C和Ctrl V 快捷键。

 

 

测量标记汇总(Demo):

 

C#开发PACS、RIS、3D医学影像处理系统系列教程目录整理:

菜鸟入门篇:

PACS客户端:

C#开发PACS医学影像处理系统(一):开发背景和功能预览

C#开发PACS医学影像处理系统(二):界面布局之菜单栏

C#开发PACS医学影像处理系统(三):界面布局之工具栏

C#开发PACS医学影像处理系统(四):界面布局之状态栏

C#开发PACS医学影像处理系统(五):查询病人信息列表

C#开发PACS医学影像处理系统(六):加载Dicom影像

C#开发PACS医学影像处理系统(七):读取影像Dicom信息

C#开发PACS医学影像处理系统(八):单元格变换

C#开发PACS医学影像处理系统(九):序列控件与拖拽

C#开发PACS医学影像处理系统(十):Dicom影像下载策略与算法

C#开发PACS医学影像处理系统(十一):Dicom影像挂片协议

C#开发PACS医学影像处理系统(十二):绘图处理之图形标记

C#开发PACS医学影像处理系统(十三):绘图处理之病灶测量

C#开发PACS医学影像处理系统(十四):处理Dicom影像窗宽窗位

C#开发PACS医学影像处理系统(十五):Dicom影像交叉定位线算法

C#开发PACS医学影像处理系统(十六):2D处理之影像平移和缩放

C#开发PACS医学影像处理系统(十七):2D处理之影像旋转和翻转

C#开发PACS医学影像处理系统(十八):Dicom使用LUT色彩加强和反色

C#开发PACS医学影像处理系统(十九):Dicom影像放大镜

 

PACS三维处理医学图像:

C#开发PACS医学影像三维重建(一):使用VTK重建3D影像

C#开发PACS医学影像三维重建(二):使用VTK进行体绘制

C#开发PACS医学影像三维重建(三):纹理映射与颜色传输

C#开发PACS医学影像三维重建(四):3D网格平滑效果

C#开发PACS医学影像三维重建(五):基于梯度透明的组织漫游

C#开发PACS医学影像三维重建(六):三维光源与阴影效果

C#开发PACS医学影像三维重建(七):空间测量与标注

C#开发PACS医学影像三维重建(八):VR体绘制

C#开发PACS医学影像三维重建(九):MPR三视图切面重建

C#开发PACS医学影像三维重建(十):MIP最小密度投影

C#开发PACS医学影像三维重建(十一):CPR曲面重建

C#开发PACS医学影像三维重建(十二):VE虚拟内镜技术

 

熟手进阶篇:

医学图像算法:

C#处理医学图像(一):基于Hessian矩阵的血管肺纹理骨骼加强对比

C#处理医学图像(二):图像锐化加强对比

 

PACS网页端 开发Web版本的PACS:

C#开发Web端PACS(一):基于PACS客户端思想重写Web端

C#开发Web端PACS(二):使用 .Net MVC 开发手机端PACS服务端

C#开发Web端PACS(三):使用HTML5和CSS3开发PACS手机端页面

C#开发Web端PACS(四):Web端与服务端的DICOM传输

C#开发Web端PACS(五):Web端的平移缩放旋转2D操做

C#开发Web端PACS(六):Web端的窗宽窗位调整

C#开发Web端PACS(七):将移动端接入微信公众号实现医院云胶片

 

登峰造极篇:

C#开发基于Python人工智能的肺结节自动检测

C#开发基于Python人工智能的脊柱侧弯曲率算法

C#开发基于Python机器学习的医学影像骨骼仿真动画

C#开发基于Python机器学习的术后恢复模拟

C#开发基于U3D的VR眼镜设备虚拟人体三维重建

C#开发基于全息投影的裸眼3D医学影像显示技术

 

周边附加篇:

胶片打印:

C#开发医学影像胶片打印系统(一):功能与胶片排版

C#开发医学影像胶片打印系统(二):胶片打印机通信

C#开发医学影像胶片打印系统(三):不规则排版打印

 

PACS服务端:

C#开发PACS医学影像处理系统服务端(一):医疗设备的链接与收图

C#开发PACS医学影像处理系统服务端(二):高并发架构

 

PACS与RIS系统的通讯与集成

在RIS系统中调起PACS并打开Dicom影像

 

云PACS与远程会诊

C#开发PACS医学影像处理系统之云PACS(区域PACS)(一):架构概述

C#开发PACS医学影像处理系统之云PACS(区域PACS)(二):远程会诊与双向转诊

 

科幻级视频特效:

使用Adobe After Effects 制做PACS影像处理系统宣传视频

 

QQ:1850969244

近10年开发经验,主攻C#、ASP MVC,HTML5,

B/S C/S 皆可,目前研究医疗领域医学影像相关技术,

任何技术问题欢迎加QQ交流。