【ITK库学习】使用itk库进行图像配准:内插器(插值)

2023-12-23 15:12

本文主要是介绍【ITK库学习】使用itk库进行图像配准:内插器(插值),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

目录

  • 1、itkNearestNeighborInterpolateImageFunction 最近点插值
  • 2、itkLinearInterpolateImageFunction 线性插值
  • 3、itkBSplineInterpolateImageFunction B样条插值
  • 4、itkWindowedSincInterpolateImageFunction 窗口化Sinc插值
  • 5、itkRayCastInterpolateImageFunction 投射插值
  • 6、itkGaussianInterpolateImageFunction 高斯插值

配准时,参考图像的网格位置没有映射到待配准图像的网格位置时,就需要用差补法来计算待配准图像中非网格位置上的亮度值。在ITK中主要用interpolator类实现,能够插入配准函数中。

配准计算中,差补法会影响最后花的平滑度和总计算时间,因此选择具体的差补方法时要考虑平滑度和计算简化两个方面。

  • ImageFunction::IsInsideBuffer(): 可用于在调用方法之前检查边界,判断物理点是否落在图像像素的空间范围之内。
  • GetRadius():获取插值所需的半径,定义了在给定点处插值所需的周围像素的数量
  • Evaluate(const PointType & point):在点位置对图像进行插值,返回插值图像强度,指定点位置,不进行边界检查,假设该点位于图像缓冲区内
  • EvaluateAtContinuousIndex(const ContinuousIndexType & index):在连续索引位置对图像进行插值,返回指定索引位置处的插值图像强度,不进行边界检查,假设该点位于图像缓冲区内,子类必须重写此方法
  • EvaluateAtIndex(const IndexType& index):在索引位置插入图像,返回指定索引位置处的图像值,不进行边界检查,假设该点位于图像缓冲区内

1、itkNearestNeighborInterpolateImageFunction 最近点插值

该类是标量图像的最近邻像素点插值。

通过复制最近领域的强度来在非整数像素位置插值图像强度,此类根据输入图像类型和坐标表示类型(例如 float 或 double)进行模板化。

该差补法较简单。

模板设置:

typedef itk::NearestNeighborInterpolateImageFunction<InputImageType, double> InterpolatorType ;

2、itkLinearInterpolateImageFunction 线性插值

该类在指定位置线性插值图像。

在非整数像素位置线性插值图像强度,此类根据输入图像类型和坐标表示类型(例如 float 或 double)进行模板化。它的插值亮度是连续的,但是亮度梯度在网格点上是不连续的。

此函数适用于标量和矢量像素类型的图像以及 VectorImage 类型的图像。

模板设置:

typedef  itk::LinearInterpolateImageFunction<InputImageType, double> InterpolatorType ;

3、itkBSplineInterpolateImageFunction B样条插值

该类评估图像的B样条插值。

样条阶数范围是0~5,0:表示几乎和最近插补相同,1:表示和线性插补相同,>1:插补值和它的微分是空间连续的。

B-Spline的系数用递归滤波器BSplineDecompositionImageFilter计算,非网格位置的亮度值 = 亮度*已经改变了B-Spline kernels的系数来计算。

限制:样条线阶数必须在设置图像之前设置样条线阶数,使用镜像边界条件, 每个尺寸都需要相同的样条线顺序,样条线是在所有尺寸中确定的,不能选择性地选择尺寸来计算样条线。

注意:插补值肯定会超出输入图像亮度范围。在处理unsigned的数据时,插补值可能是负的。

  • SetInputImage():设置输入图像
  • Set/GetSplineOrder():获取/设置样条线阶数,0~5之间, 默认为3阶样条
  • Set/GetUseImageDirection():UseImageDirection标志确定是相对于图像网格还是相对于物理空间来计算图像导数;当该标志为 ON 时,将相对于物理空间坐标系计算导数,区别在于是否考虑图像方向,On将考虑图像方向,与Off时执行的计算量相比,将导致额外的矩阵乘法;默认值为“开”
  • SetInterpolationWeights():确定值 x 插值的权重
  • SetDerivativeWeights():确定值 x 的导数部分的权重

模板设置:

typedef itk::BSplineInterpolateImageFunction<InputImageType, double>InterpolatorType ;

4、itkWindowedSincInterpolateImageFunction 窗口化Sinc插值

该类使用加窗 sinc 函数进行插值。

与线性插值相比,该函数旨在提供具有最小混叠伪影的插值函数,根据采样理论,无限支持 sinc 滤波器(其傅里叶变换是盒式滤波器)对于函数重采样是最佳的,实际上,无限支持 sinc 滤波器是使用有限支持“窗口”sinc 滤波器来近似的。

在二维中,位置 (x,y) 处的插值由以下表达式给出:
在这里插入图片描述
其中 m 是窗口的“半径”(3,4 是合理的数字),K(t) 是核函数,由 sinc 函数和几个可能的窗口函数之一组成:
在这里插入图片描述
itk::Function 命名空间中提供了几个窗口函数。 参考论文的结论建议在 m = 4,5 时使用 Welch、Cosine、Kaiser 和 Lanczos 窗,这些是基于旋转医学图像的误差线性插值法,在某些情况下,结果的准确性提高了 20 倍。
能用interpolaterwindows包括:

  • Cosinus windo
    在这里插入图片描述

  • Hamming window
    在这里插入图片描述

  • Welch window
    在这里插入图片描述

  • Lancos window
    在这里插入图片描述

  • Blackman window
    在这里插入图片描述

这种方法可以通过窗函数控制插值结果的平滑程度,可以获得比线性插值更平滑的插值结果,但是计算速度较慢,且窗函数的选择可能会影响插值结果。

在设置模板类型时:itk::WindowedSincInterpolateImageFunction< TInputImage, VRadius, TWindowFunction, TBoundaryCondition, TCoordRep >
其中:

  • TInputImage:是图像类型,这是标准的。
  • VRadius:为kernel半径,例如上面公式中的m;
  • TWindowFunction:是window的函数对象,根据上述五中不同函数选择,默认为Hamming window,根据引用的论文,它很常用,但不是最佳的;
  • TBoundaryCondition:确定落在图像边缘的像素值boundary condition class
typedef itk::ConstantBoundaryCondition<InputImageType> BoundaryConditionType ;
constexpr unsigned int WindowRadius = 5;
typedef itk::Function::HammingWindowFunction<WindowRadius> WindowFunctionType ;
typedef itk::WindowedSincInterpolateImageFunction<InputImageType,WindowRadius,WindowFunctionType,BoundaryConditionType,double> InterpolatorType;

5、itkRayCastInterpolateImageFunction 投射插值

在指定位置对图像进行投影插值,它在图像中的光线路径上进行插值。

RayCastInterpolateImageFunction 通过 3 维图像投射光线,并使用双线性插值来整合所遍历的每个体素平面。

这种方法可以模拟投影像素在图像中的路径,从而产生更精确的插值结果。但是计算速度较慢。

警告
该插值器仅适用于 3 维图像。

  • SetTransform():设置连接变换,此变换用于计算新的焦点位置
  • GetModifiableTransform():获取连接变换,此变换用于计算新的焦点位置
  • Set/GetThreshold():设置/获取阈值,高于该阈值,沿射线路径的体素将被积分
  • SetInterpolator():连接插补器
  • GetModifiableInterpolator():获取指向插值器的指针
  • Set/GetFocalPoint():设置/获取射线源的焦点或位置
  • IsInsideBuffer():检查图像缓冲区内是否有连续索引 或 检查索引是否位于图像缓冲区内,考虑到每个体素的中心位于整数坐标并延伸到下一个整数坐标的一半

模板设置:

typedef itk::RayCastInterpolateImageFunction<InputImageType, double> RayCastInterpolateFunctionType;

6、itkGaussianInterpolateImageFunction 高斯插值

该类为图像的高斯插值。

此类使用 vnl 误差函数定义 N 维高斯插值函数,与该函数相关的两个参数是:

  • Sigma:大小为 ImageDimension 的标量数组,确定插值函数的宽度。
  • Alpha:指定计算函数的截止距离的标量。

该类对噪声具有平滑作用,插值结果更接近真实图像,但计算复杂度较高。

  • Set/GetSigma():设置/获取Sigma值
  • Set/GetAlpha():设置/获取Alpha值
  • SetParameters(RealType* sigma, RealType alpha ):设置/获取Sigma和Alpha值
  • SetCutOffDistance():设置/获取截止距离
  • SetBoundingBoxStart():设置/获取边界框起点
  • SetBoundingBoxEnd():设置/获取边界框终点
typedef itkGaussianInterpolateImageFunction<InputImageType, double>

这篇关于【ITK库学习】使用itk库进行图像配准:内插器(插值)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/528478

相关文章

Spring Boot 集成 Quartz并使用Cron 表达式实现定时任务

《SpringBoot集成Quartz并使用Cron表达式实现定时任务》本篇文章介绍了如何在SpringBoot中集成Quartz进行定时任务调度,并通过Cron表达式控制任务... 目录前言1. 添加 Quartz 依赖2. 创建 Quartz 任务3. 配置 Quartz 任务调度4. 启动 Sprin

Linux下如何使用C++获取硬件信息

《Linux下如何使用C++获取硬件信息》这篇文章主要为大家详细介绍了如何使用C++实现获取CPU,主板,磁盘,BIOS信息等硬件信息,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录方法获取CPU信息:读取"/proc/cpuinfo"文件获取磁盘信息:读取"/proc/diskstats"文

Java使用SLF4J记录不同级别日志的示例详解

《Java使用SLF4J记录不同级别日志的示例详解》SLF4J是一个简单的日志门面,它允许在运行时选择不同的日志实现,这篇文章主要为大家详细介绍了如何使用SLF4J记录不同级别日志,感兴趣的可以了解下... 目录一、SLF4J简介二、添加依赖三、配置Logback四、记录不同级别的日志五、总结一、SLF4J

使用Python实现一个优雅的异步定时器

《使用Python实现一个优雅的异步定时器》在Python中实现定时器功能是一个常见需求,尤其是在需要周期性执行任务的场景下,本文给大家介绍了基于asyncio和threading模块,可扩展的异步定... 目录需求背景代码1. 单例事件循环的实现2. 事件循环的运行与关闭3. 定时器核心逻辑4. 启动与停

如何使用Nginx配置将80端口重定向到443端口

《如何使用Nginx配置将80端口重定向到443端口》这篇文章主要为大家详细介绍了如何将Nginx配置为将HTTP(80端口)请求重定向到HTTPS(443端口),文中的示例代码讲解详细,有需要的小伙... 目录1. 创建或编辑Nginx配置文件2. 配置HTTP重定向到HTTPS3. 配置HTTPS服务器

Java使用ANTLR4对Lua脚本语法校验详解

《Java使用ANTLR4对Lua脚本语法校验详解》ANTLR是一个强大的解析器生成器,用于读取、处理、执行或翻译结构化文本或二进制文件,下面就跟随小编一起看看Java如何使用ANTLR4对Lua脚本... 目录什么是ANTLR?第一个例子ANTLR4 的工作流程Lua脚本语法校验准备一个Lua Gramm

Java Optional的使用技巧与最佳实践

《JavaOptional的使用技巧与最佳实践》在Java中,Optional是用于优雅处理null的容器类,其核心目标是显式提醒开发者处理空值场景,避免NullPointerExce... 目录一、Optional 的核心用途二、使用技巧与最佳实践三、常见误区与反模式四、替代方案与扩展五、总结在 Java

使用Java将DOCX文档解析为Markdown文档的代码实现

《使用Java将DOCX文档解析为Markdown文档的代码实现》在现代文档处理中,Markdown(MD)因其简洁的语法和良好的可读性,逐渐成为开发者、技术写作者和内容创作者的首选格式,然而,许多文... 目录引言1. 工具和库介绍2. 安装依赖库3. 使用Apache POI解析DOCX文档4. 将解析

QT进行CSV文件初始化与读写操作

《QT进行CSV文件初始化与读写操作》这篇文章主要为大家详细介绍了在QT环境中如何进行CSV文件的初始化、写入和读取操作,本文为大家整理了相关的操作的多种方法,希望对大家有所帮助... 目录前言一、CSV文件初始化二、CSV写入三、CSV读取四、QT 逐行读取csv文件五、Qt如何将数据保存成CSV文件前言

Qt中QUndoView控件的具体使用

《Qt中QUndoView控件的具体使用》QUndoView是Qt框架中用于可视化显示QUndoStack内容的控件,本文主要介绍了Qt中QUndoView控件的具体使用,具有一定的参考价值,感兴趣的... 目录引言一、QUndoView 的用途二、工作原理三、 如何与 QUnDOStack 配合使用四、自