0

私の目標は、2D 画像の各行 (x 方向) の計算を計算することです。この例 (http://www.itk.org/Wiki/ITK/Examples/ImageProcessing/ConvolutionImageFilter) を見て、私の質問は、ガウス フィルターの値をカーネルに含める方法です。コードを試してみました以下のようにコンパイルするのですが、convolutionFilter を更新すると実行が止まらないようです…。私のコードは以下にあり、それほど多くはありません。何か助けがあれば、本当に感謝しています!!!!.

typedef itk::ConvolutionImageFilter <ImageType> ConvFilterType;

ImageType::Pointer kernel = ImageType::New();
CreateKernel(kernel, Gaussian , 256);

//Convolve image with kernel
ConvFilterType::Pointer convolutionFilter = ConvFilterType::New();
convolutionFilter->SetInput(Li_itk);
convolutionFilter->SetImageKernelInput(kernel);
convolutionFilter->Update();

void Reg_image_v2::CreateKernel(Reg_image_v2::ImageType::Pointer kernel, double *Gaussian, unsigned int width) {

ImageType::IndexType start;
start.Fill(0);

ImageType::SizeType size;
    size[0] = width;
size[1] = 1;

ImageType::RegionType region;
region.SetSize(size);
region.SetIndex(start);

kernel->SetRegions(region);
kernel->Allocate();

ImageType::IndexType pixelIndex;
for (int i = 0 ; i < width ; i ++){
    pixelIndex[0] = i ;
    pixelIndex[1] = 0 ;
    kernel->SetPixel(pixelIndex,Gaussian[i]);
}

itk::ImageRegionIterator<ImageType> imageIterator(kernel, region);

while(!imageIterator.IsAtEnd()) 

{
//imageIterator.Set??

++imageIterator; //operator++  increments the iterator one position in the positive direction
}


} // end createKernel

アントニオ

4

1 に答える 1

1

itkGaussianOperator を見たことがありますか? http://www.itk.org/Wiki/ITK/Examples/Operators/GaussianOperator

于 2011-12-12T00:38:18.800 に答える