/*========================================================================= * * Copyright NumFOCUS * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0.txt * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. * *=========================================================================*/ #include "itkEllipseSpatialObject.h" #include "itkPointSetToSpatialObjectDemonsRegistration.h" #include "itkRegularSphereMeshSource.h" int itkPointSetToSpatialObjectDemonsRegistrationTest(int, char *[]) { constexpr unsigned int Dimension = 3; using EllipseType = itk::EllipseSpatialObject; // Create a ellipse. EllipseType::Pointer ellipse = EllipseType::New(); // Set the radius ellipse->SetRadiusInObjectSpace(50); // Set its position EllipseType::TransformType::OffsetType offset; offset[0] = 50; offset[1] = 50; offset[2] = 50; ellipse->Update(); using PointSetType = itk::Mesh; using SphereType = itk::RegularSphereMeshSource; SphereType::Pointer sphereSource = SphereType::New(); sphereSource->Update(); using DemonsRegistrationType = itk::PointSetToSpatialObjectDemonsRegistration; DemonsRegistrationType::Pointer demonsRegistration = DemonsRegistrationType::New(); demonsRegistration->SetFixedPointSet(sphereSource->GetOutput()); demonsRegistration->SetMovingSpatialObject(ellipse); try { demonsRegistration->Update(); } catch (const itk::ExceptionObject & excp) { std::cerr << "Exception thrown during the registration process" << std::endl; std::cerr << excp << std::endl; } std::cout << "Test Succeed!" << std::endl; return EXIT_SUCCESS; }