【文件属性】:
文件名称:算法和函数-r语言常用函数速查手册
文件大小:4.63MB
文件格式:PDF
更新时间:2021-06-15 21:24:53
医学图像
10.2 算法和函数
在以前的部分,我们介绍了在ITK统计子系统的数据箱。我们也需要数据处理算法和用
这些数据箱来指导统计分析或统计分类的统计函数。这里,我们定义一个算法,这个算法是
例子中对度量向量集合的操作。一个函数是对单独一个度量向量的操作。例如,如果我们执
行一个类( itk::Statistics::EuclideanDistance)去计算在度量向量之间的欧几里德距离,我们调
用它作为一个函数;如果我们执行一个类(itk::Statistics::MeanCalculator)去计算一个样本的平
均值,我们调用它作为一个算法。
10.2.1 样本统计
我们将介绍如何用(itk::Statistics::Sample)类得到样本统计值,如平均值和协方差。统计
结果能告诉我们一个样本的特性。这样的样本统计对于统计分类是非常重要的。当我们知道
样本分配的形式和它们的参数的时候,我们能够使用贝叶斯分类方法。在ITK中,样本平均
值和协方差的计算算法被执行。每一个算法也有它的加权形式。在 大期望参数估计过程使
用的是加权形式。
平均值和协方差
这部分的源码在文件Examples/Statistics/SampleStatistics.cxx中。
我们包含itk::Vector类的头文件,itk::Vector类在我们的例子中将是我们的度量向量模
板:
#include "itkVector.h"
我们将用itk::Statistics::ListSample作为我们的样本模板。我们也包含这个类的头文件:
#include "itkListSample.h"
接下来的头文件是样本统计算法。
#include "itkMeanCalculator.h"
#include "itkCovarianceCalculator.h"
接下来的代码片断将会创建一个3结构浮点度量向量的ListSample对象,并将5个度量向
量参数放进ListSample对象中:
const unsigned int MeasurementVectorLength = 3;
typedef itk::Vector< float, MeasurementVectorLength > MeasurementVectorType;
typedef itk::Statistics::ListSample< MeasurementVectorType > SampleType;
SampleType::Pointer sample = SampleType::New( );
sample->SetMeasurementVectorSize( MeasurementVectorLength );