forward_zernike_subtomos.cpp File Reference

#include "forward_zernike_subtomos.h"
#include "core/transformations.h"
#include "core/xmipp_image_extension.h"
#include "core/xmipp_image_generic.h"
#include "data/projection.h"
#include "data/mask.h"

Include dependency graph for forward_zernike_subtomos.cpp:

Go to the source code of this file.

Functions
double	continuousZernikeSubtomoCost (double *x, void *_prm)

Function Documentation

continuousZernikeSubtomoCost()

double continuousZernikeSubtomoCost	(	double *	x,
		void *	_prm
	)

Definition at line 372 of file forward_zernike_subtomos.cpp.

{
    ProgForwardZernikeSubtomos *prm=(ProgForwardZernikeSubtomos *)_prm;
    int idx = 3*(prm->vecSize);
    // TODO: Optimize parameters for each image (not sharing)
    // prm->deltaDefocusU[0]=x[idx+6]; prm->deltaDefocusU[1]=x[idx+6]; 
    // prm->deltaDefocusV[0]=x[idx+7]; prm->deltaDefocusV[1]=x[idx+7]; 
    // prm->deltaDefocusAngle[0]=x[idx+8]; prm->deltaDefocusAngle[1]=x[idx+8];

    prm->deltaX = x[idx + 1];
    prm->deltaY = x[idx + 2];
    prm->deltaZ = x[idx + 3];
    prm->deltaRot = x[idx + 4];
    prm->deltaTilt = x[idx + 5];
    prm->deltaPsi = x[idx + 6];
    prm->deltaDefocusU=x[idx + 7];
    prm->deltaDefocusV=x[idx + 8];
    prm->deltaDefocusAngle=x[idx + 9];

    MAT_ELEM(prm->A, 0, 3) = prm->old_shiftX + prm->deltaX;
    MAT_ELEM(prm->A, 1, 3) = prm->old_shiftY + prm->deltaY;
    MAT_ELEM(prm->A, 2, 3) = prm->old_shiftZ + prm->deltaZ;
    MAT_ELEM(prm->A, 0, 0) = 1;
    MAT_ELEM(prm->A, 0, 1) = 0;
    MAT_ELEM(prm->A, 0, 2) = 0;
    MAT_ELEM(prm->A, 1, 0) = 0;
    MAT_ELEM(prm->A, 1, 1) = 1;
    MAT_ELEM(prm->A, 1, 2) = 0;
    MAT_ELEM(prm->A, 2, 0) = 0;
    MAT_ELEM(prm->A, 2, 1) = 0;
    MAT_ELEM(prm->A, 2, 2) = 1;
    // MAT_ELEM(prm->A, 3, 0) = 0;
    // MAT_ELEM(prm->A, 3, 1) = 0;
    // MAT_ELEM(prm->A, 3, 2) = 0;
    // MAT_ELEM(prm->A, 3, 3) = 1;

    return prm->transformImageSph(x);

    // return prm->transformImageSph(x,prm->old_rot+deltaRot, prm->old_tilt+deltaTilt, prm->old_psi+deltaPsi,
    //      prm->A, deltaDefocusU, deltaDefocusV, deltaDefocusAngle);
}