Beispiel #1
0
    void RunRigid() {
        _initialTransform = InitialTransform::New();

        _initialTransform->SetFixedParameters(_centerOfRotation);

        OptiReporter::Pointer optiReporter = OptiReporter::New();
        Metric::Pointer metric = Metric::New();
        metric->SetFixedImage(_dst);
        metric->SetMovingImage(_src);
        metric->SetInterpolator(Interpolator::New());
        metric->SetTransform(_initialTransform);
        //metric->SetFixedImageRegion(_dst->GetBufferedRegion());
        metric->SetFixedImageIndexes(_labelIndexes);
//        metric->SetUseAllPixels(true);
        metric->Initialize();

        Optimizer::Pointer opti = Optimizer::New();
        opti->SetCostFunction(metric);
        Optimizer::ScalesType scales;
        scales.SetSize(InitialTransform::ParametersDimension);
        scales.Fill(1);
        //scales[0] = scales[1] = scales[2] = .1;
        
        opti->SetUseUnitLengthGradient(true);
        opti->SetScales(scales);
        opti->SetInitialPosition(_initialTransform->GetParameters());
        opti->AddObserver(itk::StartEvent(), optiReporter);
        opti->AddObserver(itk::IterationEvent(), optiReporter);
        opti->StartOptimization();
        cout << "Current Cost: " << opti->GetCurrentCost() << endl;
        _initialResult = opti->GetCurrentPosition();
        _initialTransform->SetParameters(opti->GetCurrentPosition());
    }
Beispiel #2
0
    void RunNine() {
        _finalTransform = NineTransform::New();
//        NineTransform::OutputVectorType initialTranslation;
//        initialTranslation[0] = _initialResult[0];
//        initialTranslation[1] = _initialResult[1];
//        initialTranslation[2] = _initialResult[2];

        // _finalTransform->SetTranslation(initialTranslation);
        _finalTransform->SetFixedParameters(_centerOfRotation);
        
        OptiReporter::Pointer optiReporter = OptiReporter::New();
        Metric::Pointer metric = Metric::New();
        metric->SetFixedImage(_dst);
        metric->SetMovingImage(_src);
        metric->SetInterpolator(InterpolatorNN::New());
        metric->SetTransform(_finalTransform);
        metric->SetFixedImageIndexes(_labelIndexes);
//        metric->SetUseAllPixels(true);
        metric->Initialize();

        Optimizer::Pointer opti = Optimizer::New();
        opti->SetCostFunction(metric);
        Optimizer::ScalesType scales;
        scales.SetSize(NineTransform::ParametersDimension);
        scales.Fill(1);
        scales[0] = scales[1] = scales[2] = 10;
        scales[6] = scales[7] = scales[8] = 100;

        opti->SetMaximumIteration(500);
        opti->SetUseUnitLengthGradient(true);
        opti->SetStepLength(0.5);
        opti->SetScales(scales);
        opti->SetInitialPosition(_finalTransform->GetParameters());
        opti->AddObserver(itk::StartEvent(), optiReporter);
        opti->AddObserver(itk::IterationEvent(), optiReporter);
        opti->StartOptimization();
        cout << "Current Cost: " << opti->GetCurrentCost() << endl;
        _finalResult = opti->GetCurrentPosition();
        _finalTransform->SetParameters(opti->GetCurrentPosition());
    }