C#开发PACS医学影像三维重建(十三):基于人体CT值从皮肤渐变到骨骼的梯度透明思路

已于 2023-06-07 14:41:26 修改
阅读量1.3k 收藏 9
点赞数 2
分类专栏： C# PACS影像处理系统 文章标签： 3d c# pacs dicom wpf
于 2023-05-05 01:47:43 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/after5120/article/details/130498377
版权

当我们将CT切片重建为三维体之后，通常会消除一些不必要的外部组织来观察内部病灶，

一般思路是根据人体常见CT值范围来使得部分组织透明来达到效果，

但这是非黑即白的，即，要么显示皮肤，要么显示神经，要么显示骨骼…

然而，当病灶介于两个CT值范围中间或重合时，要观察病灶这将是一个非常困难的工作，需要手动旋转、切割等一系列操作，

严重降低判断准确性和工作效率，故此采用梯度透明思想来展现渐变过程。

准备一张人体CT值范围表参考：

当我们对三维体使用鼠标操作时，将CT值从最小值到最大值映射到三维体上面，并将透明度从0到100渐变，

如果能自己用代码实现效果最好，否则也可以用VTK进行梯度透明，部分C++代码如下：

vtkSmartPointer<vtkPiecewiseFunction> GradualOp = vtkSmartPointer<vtkPiecewiseFunction>::New();
GradualOp->AddPoint(10, 0.0);
GradualOp->AddPoint(20, 0.5);
GradualOp->AddPoint(30, 1.0);
volumeProperty->SetGradientOpacity(GradualOp);

梯度小于10点的不透明度乘子设为0，完全透明。

梯度大小为10时，不透明度为0；

梯度大小为20时，不透明度乘子为0.5；

梯度大小在10-20时，不透明度乘子通过线性映射至0-0.5之间的数值；

当梯度大小在20-30之间时，不透明度乘子通过线性映射至0.5-1.0之间；

梯度值大于30时，不透明度乘子为1
 

同样最后还是用观察者模式展现体重建

//设置相机
    vtkCamera aCamera = new vtkCamera();
 
    aCamera.SetViewUp(0, 0, -1);
 
    aCamera.SetPosition(0, 1, 0);
 
    aCamera.SetFocalPoint(0, 0, 0);
 
    aCamera.ComputeViewPlaneNormal();
 
 
    //设置Actor
    vtkActor coneActor = new vtkActor();
 
    coneActor.SetMapper(skinMapper);
 
    coneActor.GetProperty().SetAmbient(0.5);
 
    coneActor.GetProperty().SetDiffuse(1);
 
    coneActor.GetProperty().SetSpecular(0.6);
 
    //显示类
    vtkRenderer renderer = renderWindowControl1.RenderWindow.GetRenderers().GetFirstRenderer();
 
    renderer.AddActor(coneActor);//添加coneActor对象
 
    //renderer.AddActor2D(new vtkProp());//添加textActor对象
 
    renderer.SetBackground(0, 0, 0);
 
    renderer.SetActiveCamera(aCamera);//添加相机
 
    renderer.ResetCamera();
 
    vtkRenderWindow renWin = renderWindowControl1.RenderWindow;//设置绘图窗口renWin->AddRenderer(renderer);//装载绘图类
 
    vtkWin32RenderWindowInteractor iren = new vtkWin32RenderWindowInteractor();
 
    iren.SetRenderWindow(renWin);//装载绘图窗口

依照上述思想，实现最终效果，

案例①：从面部重建渐变到口腔牙齿，显露出牙床钢钉：

更为直观、方便的观察牙床钢钉及口腔其他情况：

案例②：颅脑及上胸腔平扫，从皮肤渐变到肌肉、神经、骨骼

渐变过程中观察骨骼、神经、血管等病灶：

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

菜鸟入门篇

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影像放大镜

医学影像三维篇

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#开发PACS医学影像三维重建(十三):基于人体CT值从皮肤渐变到骨骼的梯度透明思路

C#开发PACS医学影像三维重建(十四):基于能量模型算法将曲面牙床展开至二维平面

熟手进阶篇

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

C#处理医学影像(二):基于Hessian矩阵的医学影像增强与窗宽窗位

C#处理医学影像(三):基于漫水边界自动选取病灶范围的实现思路

C#处理医学影像(四):基于Stitcher算法拼接人体全景脊柱骨骼影像

胶片打印：

C#开发医学影像胶片打印系统(一)：万能花式布局的实现思路

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

C#开发医学影像胶片打印系统(三)：Pacs二维功能在排版中的应用

登峰造极篇

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

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

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

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

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

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

免费下载

免费下载使用本教程PACS软件

博客
C#处理医学影像(四):基于Stitcher算法拼接人体全景脊柱骨骼影像
05-05 1186
在拍摄脊柱或胸片时，经常会遇到因设备高度不够需要分段拍摄的情况，对于影像科诊断查阅影像时希望将分段影像合并成一张影像，有助于更直观的观察病灶，以下图为例的两个分段影像：我们使用OpenCVSharp中的Stitcher类的Stitch方法，导入两张图像并拼接：但结果却失败了，返回错误结果：ERR_NEED_MORE_IMGS，是由于医学影像的特征点匹配不够，导致无法确定对接点。一幅图中总存在着一些独特的像素点，这些点我们可以认为就是这幅图的特征，即为特征点。
博客
C#开发PACS医学影像三维重建(十四):基于能量模型算法将曲面牙床展开至二维平面
05-05 1084
在医学影像领域中，将三维重建中的人体组织展开平铺至二维，用来研判病灶和制定治疗方案的重要手段之一，众所周知，人体牙床正常情况下是有弧度的，无论是从俯视位还是冠状位观察都是不能直观的了解牙齿状况，同理，还可以应用到血管和骨骼的平铺展开，更好的观察血管阻塞和破裂、骨折骨裂等情况。牙床全景图、平铺血管、骨骼二维化展开(肋骨平铺)。③迭代应变能量释放法。
博客
C#处理医学影像(三):基于漫水边界自动选取病灶范围的实现思路
05-05 506
通过给定指定的中心坐标，和指定染色的颜色值，向四周扩撒，遇到一样或近似值将其包含在范围内，网上源码很多，很容易实现。其中直观区别是canny算子计算的结果清晰，但不连续，容易受噪点影响，而sobel算子线条相对柔和，连续性强。医生在实际使用PACS软件观察病灶时，经常会测量不规则病灶的周长和面积，使用画笔工具勾勒比较耗时且准度欠佳，这些点是不连续的，不首尾相连的，我们需要将非边缘坐标删除，并将边缘坐标按顺时针或逆时针排序；将漫水填充范围的最大外接矩形裁剪出来，提高后续计算效率。
博客
C#开发医学影像胶片打印系统(一)：万能花式布局的实现思路
05-05 1063
本篇文章将介绍开发医学影像胶片打印系统(printscu模式)遇到不规则排版时的一种思路，一般来讲，医院打印胶片时都是整张胶片打印，但有时需要将多个病人或一个病人的多个检查打印在同一张胶片上，这时候就需要不规则排版来满足打印需求，使胶片利用率最大化。
博客
C#处理医学图像(二):基于Hessian矩阵的医学图像增强与窗宽窗位
05-05 941
在Slider控件鼠标抬起的事件中，先将目标单元格内WPF图像转为base64，发送给我们生成的C++接口，再将返回的base64转为WPF图像。乳腺影像由原始dcm显示的絮状腺体在图像增强下变成丝状，对于乳腺中存在肿瘤或其他病症的显示更为明显，亦可自由调整参数达到自己想要的效果。考虑到性能问题，因为是base64转码和解码，所以不推荐用ValueChanged，它的触发频率要高得多，由原理中得知，求得特征值和特征向量反应出的变化上的各向异性
博客
C#处理医学图像(一):基于Hessian矩阵的血管肺纹理骨骼增强对比
05-05 810
如果不想用C++，直接用C#里面的opencv库也可以，直接用nuget搜索EmguCV，需要自己将MatLab代码或C++代码翻译成C#导出32位dll，复制到C#debug目录下，C#调用：将目标图像转为base64，发送给C++，返回处理后的base64，在转为图像。将其原理翻译写成C++类库，在C++中使用Opencv对于矩阵操作比较方便，导出dll后再由C#调用，在医院实际环境中，经常遇到有问题的患者，对于一些特殊的场景，比如骨折，肺结节，心脑血管问题。
博客
C#开发PACS医学影像三维重建(一):使用VTK重建3D影像
05-05 5861
VTK是一个开源的免费软件系统，主要用于三维计算机图形学、图像处理和可视化。Vtk是在面向对象原理的基础上设计和实现的，它的内核是用C++构建的。//血管CT值为200-300左右。因为使用C#语言开发，而VTK是C++的，所以推荐使用VTK的.Net开发库：ActiViz。//建立算法对象，从CT切片数据中提取出皮肤。最大密度投影（MIP）表面阴影遮盖（SSD）多层面重建（MPR）容积漫游技术（VR）虚拟内镜技术（VE）
博客
C#开发PACS医学影像处理系统(十九)：Dicom影像放大镜
05-05 846
在XAML代码设计器中，添加canvas画布与圆形几何对象，利用VisualBrush笔刷来复制画面内容到指定容器
博客
C#开发PACS医学影像处理系统(十八)：Dicom使用LUT色彩增强和反色
05-05 1106
在医生阅片确诊的过程中，当发线疑似病灶时在灰度显示下有时并不清晰，这时候就需要色彩增强效果来使灰度图像变为彩色图像。在有些类型的影像中，例如牙片(PX),胸片等，需要黑白反色的情况下观察病灶。//实例化文件处理对象并打开文件。//获取dicom图像对象。或者使用LUT逆灰度编码也能达到同样的效果。
博客
C#开发PACS医学影像处理系统(十七)：2D处理之影像旋转和翻转
05-05 722
根据鼠标按下时的位置和图片中心位置，利用WPF中Vector类计算向量夹角，得到角度，再设置元素的RotateTransform旋转属性变换。在XAML设计器中，设置RotateTransform属性。）文中提到的ScaleTransform属性中，当 st.ScaleX的值为-1时则水平翻转，当 st.ScaleY的值为-1时则垂直翻转。同样还是先获取变换对象。
博客
C#开发PACS医学影像处理系统(十六)：2D处理之影像平移和缩放
05-05 792
【代码】C#开发PACS医学影像处理系统(十六)：2D处理之影像平移和缩放。
博客
C#开发PACS医学影像处理系统(十五)：Dicom影像交叉定位线算法
05-05 1162
而左边图像是从左手到右手的扫描切片，那么右边图像的位置就恰好在左边图像的眼睛的位置，用红线标出的位置则为定位线，一般用于医生参考病灶在矢状面、冠状面和横断面的具体方位。举个例子：右边的是MR(核磁共振)的某一帧切片，这是从头开始扫描，扫描到眼睛这个位置，得到坐标之后就可以利用绘图操作类(参考本系列教程之图形标记)来自己绘制定位线。在笛卡尔空间直角坐标系中，Y 右肩膀到左肩膀，X 后背到前胸，Z 足到头。
博客
C#开发PACS医学影像处理系统(十四)：处理Dicom影像窗宽窗位
05-05 3292
反之，如果使用宽的窗宽，则显示的CT值范围大，每一灰阶代表的CT 值幅度大，则图像对比度差，但密度均匀，适于观察密度差别大的结构。则以窗位为中心（0Hu），向上包括+50Hu，向下包括-50Hu，凡是在这个100Hu 范围内的组织均可显示出来并为人眼所识别。纵隔内淋巴结的大小,纵隔内肿块及这些结构的比邻关系,设定纵隔窗可用窗宽300 Hu~500 Hu,窗位30 Hu~50 Hu。喉颈部、鼻咽、咽喉部的窗宽和窗位常设在300 Hu~350 Hu和30 Hu~50 Hu,能满足该部位的解剖和病灶显示,
博客
C#开发PACS医学影像处理系统(十三)：绘图处理之病灶测量
05-05 797
接上一篇文章，当我们可以绘制图形标记后，就可以在此操作类上面进行扩展，一些角度类型的测量计算相对复杂一点，需要计算角度大小和绘制弧线。计算两点之间的距离：其中dpi是缩放比例。
博客
C#开发PACS医学影像处理系统(十二)：绘图处理之图形标记
05-05 966
那么以此类推，通过命名空间System.Windows.Shapes，我们可以继续绘制其他图形了。在医生实际使用过程中，对于有病灶的影像需要一些2D绘图操作，例如对于病灶的标记和测量，最后，关于撤销上个标记和清除所有标记，只需要将图形对象保存在LIst集合中，绘制控制点：一个白色的圆形内嵌一个黑色的圆形，作为绘图完成之后的控制点。需要删除或撤销的时候从集合中移除，并在画布中移除即可。在WPF中， 推荐使用InkCanvas控件。在此类中创建绘制方法。
博客
C#开发PACS医学影像处理系统(十一)：Dicom影像挂片协议
05-05 1254
通俗点说，挂片协议可以看作整个系统的一个相对复杂一点的配置文件，可以用JSON或XML格式来读取与保存，另外，可以制作一个独立的exe配置程序来管理这些挂片协议。
博客
C#开发PACS医学影像处理系统(十)：Dicom影像下载策略与算法
09-15 1632
1.传输协议当我们PACS客户端接收到下载命令时需要根据dcm文件的url来下载，这里推荐使用http协议，支持大文件的断点续传。2.下载策略假设有10个序列，每个序列有100张图像，①广度优先：将任务池中的下载任务平均分配到10个序列，每个序列都同时下载，优点：是每个序列都能分布加载，当用户切换序列到单元格时不需要等待。缺点：当图像特别多的时候，比如CT薄层，序列切换图像滚动时显示不完整。②深度优先：将任务池中的下载任务根据先进先出原则分配到排序在前的序列，下载完一个后再开始下
博客
C#开发PACS医学影像处理系统(九)：序列控件与拖拽
09-15 1075
1.先看结构：创建WPF用户控件：YourTab创建WPF用户控件：YourItem创建选项卡时循环添加item，并设置序列缩略图到控件和异步下载的进度条，1个病人1个或多个Study检查，1个Study检查1个选项卡，并为每个选项卡添加自定义关闭按钮Tab自定义关闭按钮的XAML样式如下：<Style TargetType="{x:Type TabItem}"> <Setter Property="Templat...
博客
C#开发PACS医学影像处理系统(八)：单元格变换
09-15 1252
再看结构，Viewer(影像控件)分为N个Cell(单元格)，每个单元格分为N个Box(影像容器)，当默认情况为2x2的行列时，结构如下图：其中，我们在WPF设计器中创建一个用户控件叫Viewer，使用Grid容器，默认2x2，再根据自己的业务需要创建用户控件Cell和用户控件Box，在Viewer中声明变量集合来记录Cell，同理，在Cell中声明变量集合来记录Box：
博客
C#开发PACS医学影像处理系统(七)：读取影像Dicom信息
09-14 4141
Dicom的文件信息一般是由设备信息、图像信息、病人信息等Tag标记组成的，用字节长度和偏移量来区分，具体数据结构可以自行去参考Dicom文件解析，内容比较多，这里不赘述。对于常用的Tag，可以简单的理解为两层结构，第一层是组，第二层是组包含的元素，每个组号和元素号用2个字节表示(Uint16),例如常用的组有 0002组描述设备通讯，0008组描述特征参数，0010组描述患者信息，0028组描述图像信息参数，例如 0010组里面又包含病人姓名、年龄、性别等。同样的，按照上一..
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值