【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

相关文章

java中使用POI生成Excel并导出过程

《java中使用POI生成Excel并导出过程》:本文主要介绍java中使用POI生成Excel并导出过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录需求说明及实现方式需求完成通用代码版本1版本2结果展示type参数为atype参数为b总结注:本文章中代码均为

Spring Boot3虚拟线程的使用步骤详解

《SpringBoot3虚拟线程的使用步骤详解》虚拟线程是Java19中引入的一个新特性,旨在通过简化线程管理来提升应用程序的并发性能,:本文主要介绍SpringBoot3虚拟线程的使用步骤,... 目录问题根源分析解决方案验证验证实验实验1:未启用keep-alive实验2:启用keep-alive扩展建

使用Java实现通用树形结构构建工具类

《使用Java实现通用树形结构构建工具类》这篇文章主要为大家详细介绍了如何使用Java实现通用树形结构构建工具类,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录完整代码一、设计思想与核心功能二、核心实现原理1. 数据结构准备阶段2. 循环依赖检测算法3. 树形结构构建4. 搜索子

基于Flask框架添加多个AI模型的API并进行交互

《基于Flask框架添加多个AI模型的API并进行交互》:本文主要介绍如何基于Flask框架开发AI模型API管理系统,允许用户添加、删除不同AI模型的API密钥,感兴趣的可以了解下... 目录1. 概述2. 后端代码说明2.1 依赖库导入2.2 应用初始化2.3 API 存储字典2.4 路由函数2.5 应

GORM中Model和Table的区别及使用

《GORM中Model和Table的区别及使用》Model和Table是两种与数据库表交互的核心方法,但它们的用途和行为存在著差异,本文主要介绍了GORM中Model和Table的区别及使用,具有一... 目录1. Model 的作用与特点1.1 核心用途1.2 行为特点1.3 示例China编程代码2. Tab

SpringBoot使用OkHttp完成高效网络请求详解

《SpringBoot使用OkHttp完成高效网络请求详解》OkHttp是一个高效的HTTP客户端,支持同步和异步请求,且具备自动处理cookie、缓存和连接池等高级功能,下面我们来看看SpringB... 目录一、OkHttp 简介二、在 Spring Boot 中集成 OkHttp三、封装 OkHttp

使用Python实现获取网页指定内容

《使用Python实现获取网页指定内容》在当今互联网时代,网页数据抓取是一项非常重要的技能,本文将带你从零开始学习如何使用Python获取网页中的指定内容,希望对大家有所帮助... 目录引言1. 网页抓取的基本概念2. python中的网页抓取库3. 安装必要的库4. 发送HTTP请求并获取网页内容5. 解

使用Python实现网络设备配置备份与恢复

《使用Python实现网络设备配置备份与恢复》网络设备配置备份与恢复在网络安全管理中起着至关重要的作用,本文为大家介绍了如何通过Python实现网络设备配置备份与恢复,需要的可以参考下... 目录一、网络设备配置备份与恢复的概念与重要性二、网络设备配置备份与恢复的分类三、python网络设备配置备份与恢复实

C#中的 StreamReader/StreamWriter 使用示例详解

《C#中的StreamReader/StreamWriter使用示例详解》在C#开发中,StreamReader和StreamWriter是处理文本文件的核心类,属于System.IO命名空间,本... 目录前言一、什么是 StreamReader 和 StreamWriter?1. 定义2. 特点3. 用

Python使用date模块进行日期处理的终极指南

《Python使用date模块进行日期处理的终极指南》在处理与时间相关的数据时,Python的date模块是开发者最趁手的工具之一,本文将用通俗的语言,结合真实案例,带您掌握date模块的六大核心功能... 目录引言一、date模块的核心功能1.1 日期表示1.2 日期计算1.3 日期比较二、六大常用方法详