文件名称:算法和函数-r语言常用函数速查手册
文件大小:4.63MB
文件格式:PDF
更新时间:2024-07-05 15:11:33
医学图像
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 );