ITK 生成sphere

itk::EllipseSpatialObject 画球

int main(int argc, char * argv[])
{
   

    using PixelType = double;
    using ImageType = itk::Image<PixelType, 3>;

    using EllipseType = itk::EllipseSpatialObject<3>;

    using SpatialObjectToImageFilterType = itk::SpatialObjectToImageFilter<EllipseType, RealImageType>;

    SpatialObjectToImageFilterType::Pointer imageFilter = SpatialObjectToImageFilterType::New();
    ImageType::SizeType size;
    size[0] = 50;
    size[1] = 50;
    size[2] = 50;
    imageFilter->SetSize(size);

    ImageType::SpacingType spacing;
    spacing[0] = 100.0 / size[0];
    spacing[1] = 100.0 / size[1];
    spacing[2] = 100.0 / size[2];
    imageFilter->SetSize(size);
    imageFilter->SetSpacing(spacing);

    EllipseType::Pointer  ellipse = EllipseType::New();
    EllipseType::ArrayType radiusArray;
    radiusArray[0] = 10;
    radiusArray[1] = 10;
    radiusArray[2] = 10;
    ellipse->SetRadiusInObjectSpace(radiusArray);

    using TransformType = EllipseType::TransformType;
    TransformType::Pointer transform = TransformType::New();
    transform->SetIdentity();
    TransformType::OutputVectorType translation;
    translation[0] = size[0] * spacing[0] / 2.0;
    translation[1] = size[1] * spacing[1] / 2.0;
    translation[2] = size[2] * spacing[2] / 2.0;
    transform->Translate(translation, false);
    ellipse->SetObjectToParentTransform(transform);
    ellipse->SetDefaultInsideValue(255);
    ellipse->SetDefaultOutsideValue(0);
    imageFilter->SetInput(ellipse);
    imageFilter->SetUseObjectValue(true);
    imageFilter->SetOutsideValue(0);
    using WriterType = itk::ImageFileWriter<ImageType>;
    WriterType::Pointer writer = WriterType::New();

    writer->SetFileName("C:/Users/Curaway/Desktop/hh.mha");
    writer->SetInput(imageFilter->GetOutput());

    try
    {
        imageFilter->Update();
        writer->Update();
    }
    catch (itk::ExceptionObject & excp)
    {
        std::cerr << excp << std::endl;
        return EXIT_FAILURE;
    }
    return EXIT_SUCCESS;
}