interpolate_macros.h File Reference

Go to the source code of this file.

Macros
#define	LI_CLAMP_3D( mijk, mijk_f, mijk_r, li_frac_1, li_frac_2, moving)
#define	LI_VALUE(m_val, fx1, fx2, fy1, fy2, fz1, fz2, mvf, m_img, moving)
#define	LI_VALUE_x(m_val, rx, fy1, fy2, fz1, fz2, mvf, m_img, moving)
#define	LI_VALUE_y(m_val, fx1, fx2, ry, fz1, fz2, mvf, m_img, moving)
#define	LI_VALUE_z(m_val, fx1, fx2, fy1, fy2, rz, mvf, m_img, moving)
#define	LI_VALUE_ALT(m_val, fx1, fx2, fy1, fy2, fz1, fz2, mvf, m_img, moving, plane)
#define	LI_VALUE_INTERLEAVED(m_val, fx1, fx2, fy1, fy2, fz1, fz2, mvf, m_img, moving, plane)

Macro Definition Documentation

LI_CLAMP_3D

#define LI_CLAMP_3D	(		mijk,
			mijk_f,
			mijk_r,
			li_frac_1,
			li_frac_2,
			moving
	)

Value:

    do {                                                                \
	li_clamp (mijk[0], moving->dim[0]-1,                           \
            &mijk_f[0], &mijk_r[0],                                     \
            &li_frac_1[0], &li_frac_2[0]);                              \
	li_clamp (mijk[1], moving->dim[1]-1,                           \
            &mijk_f[1], &mijk_r[1],                                     \
            &li_frac_1[1], &li_frac_2[1]);                              \
	li_clamp (mijk[2], moving->dim[2]-1,                           \
            &mijk_f[2], &mijk_r[2],                                     \
            &li_frac_1[2], &li_frac_2[2]);                              \
    } while (0)

LI_VALUE

#define LI_VALUE	(		m_val,
			fx1,
			fx2,
			fy1,
			fy2,
			fz1,
			fz2,
			mvf,
			m_img,
			moving
	)

Value:

    do {                                                                \
        float m_x1y1z1, m_x2y1z1, m_x1y2z1, m_x2y2z1;                   \
        float m_x1y1z2, m_x2y1z2, m_x1y2z2, m_x2y2z2;                   \
                                                                        \
        m_x1y1z1 = fx1 * fy1 * fz1 * m_img[mvf];                        \
        m_x2y1z1 = fx2 * fy1 * fz1 * m_img[mvf+1];                      \
        m_x1y2z1 = fx1 * fy2 * fz1 * m_img[mvf+moving->dim[0]];         \
        m_x2y2z1 = fx2 * fy2 * fz1 * m_img[mvf+moving->dim[0]+1];       \
        m_x1y1z2 = fx1 * fy1 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]]; \
        m_x2y1z2 = fx2 * fy1 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+1]; \
        m_x1y2z2 = fx1 * fy2 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]]; \
        m_x2y2z2 = fx2 * fy2 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]+1]; \
        m_val = m_x1y1z1 + m_x2y1z1 + m_x1y2z1 + m_x2y2z1               \
                + m_x1y1z2 + m_x2y1z2 + m_x1y2z2 + m_x2y2z2;            \
    } while (0)

LI_VALUE_ALT

#define LI_VALUE_ALT	(		m_val,
			fx1,
			fx2,
			fy1,
			fy2,
			fz1,
			fz2,
			mvf,
			m_img,
			moving,
			plane
	)

Value:

    do {                                                                \
        int idx_x1y1z1, idx_x2y1z1, idx_x1y2z1, idx_x2y2z1;             \
        int idx_x1y1z2, idx_x2y1z2, idx_x1y2z2, idx_x2y2z2;             \
        float m_x1y1z1, m_x2y1z1, m_x1y2z1, m_x2y2z1;                   \
        float m_x1y1z2, m_x2y1z2, m_x1y2z2, m_x2y2z2;                   \
                                                                        \
        idx_x1y1z1 = mvf;                                               \
        idx_x2y1z1 = mvf+1;                                             \
        idx_x1y2z1 = mvf+moving->dim[0];                                \
        idx_x2y2z1 = mvf+moving->dim[0+1];                              \
        idx_x1y1z2 = mvf+moving->dim[1]*moving->dim[0];                 \
        idx_x2y1z2 = mvf+moving->dim[1]*moving->dim[0]+1;               \
        idx_x1y2z2 = mvf+moving->dim[1]*moving->dim[0]+moving->dim[0];  \
        idx_x2y2z2 = mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]+1; \
                                                                        \
        m_x1y1z1 = fx1 * fy1 * fz1 * m_img[idx_x1y1z1];                 \
        m_x2y1z1 = fx2 * fy1 * fz1 * m_img[idx_x2y1z1];                 \
        m_x1y2z1 = fx1 * fy2 * fz1 * m_img[idx_x1y2z1];                 \
        m_x2y2z1 = fx2 * fy2 * fz1 * m_img[idx_x2y2z1];                 \
        m_x1y1z2 = fx1 * fy1 * fz2 * m_img[idx_x1y1z2];                 \
        m_x2y1z2 = fx2 * fy1 * fz2 * m_img[idx_x2y1z2];                 \
        m_x1y2z2 = fx1 * fy2 * fz2 * m_img[idx_x1y2z2];                 \
        m_x2y2z2 = fx2 * fy2 * fz2 * m_img[idx_x2y2z2];                 \
        m_val = m_x1y1z1 + m_x2y1z1 + m_x1y2z1 + m_x2y2z1               \
                + m_x1y1z2 + m_x2y1z2 + m_x1y2z2 + m_x2y2z2;            \
    } while (0)

LI_VALUE_INTERLEAVED

#define LI_VALUE_INTERLEAVED	(		m_val,
			fx1,
			fx2,
			fy1,
			fy2,
			fz1,
			fz2,
			mvf,
			m_img,
			moving,
			plane
	)

LI_VALUE_x

#define LI_VALUE_x	(		m_val,
			rx,
			fy1,
			fy2,
			fz1,
			fz2,
			mvf,
			m_img,
			moving
	)

Value:

    do {                                                                \
        float m_x1y1z1, m_x2y1z1, m_x1y2z1, m_x2y2z1;                   \
        float m_x1y1z2, m_x2y1z2, m_x1y2z2, m_x2y2z2;                   \
                                                                        \
        m_x1y1z1 = (1/rx) * fy1 * fz1 * m_img[mvf];                     \
        m_x2y1z1 = -(1/rx) * fy1 * fz1 * m_img[mvf+1];                  \
        m_x1y2z1 = (1/rx) * fy2 * fz1 * m_img[mvf+moving->dim[0]];              \
        m_x2y2z1 = -(1/rx) * fy2 * fz1 * m_img[mvf+moving->dim[0]+1];   \
        m_x1y1z2 = (1/rx) * fy1 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]]; \
        m_x2y1z2 = -(1/rx) * fy1 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+1]; \
        m_x1y2z2 = (1/rx) * fy2 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]]; \
        m_x2y2z2 = -(1/rx) * fy2 * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]+1]; \
        m_val = m_x1y1z1 + m_x2y1z1 + m_x1y2z1 + m_x2y2z1               \
                + m_x1y1z2 + m_x2y1z2 + m_x1y2z2 + m_x2y2z2;            \
    } while (0)

LI_VALUE_y

#define LI_VALUE_y	(		m_val,
			fx1,
			fx2,
			ry,
			fz1,
			fz2,
			mvf,
			m_img,
			moving
	)

Value:

    do {                                                                \
        float m_x1y1z1, m_x2y1z1, m_x1y2z1, m_x2y2z1;                   \
        float m_x1y1z2, m_x2y1z2, m_x1y2z2, m_x2y2z2;                   \
                                                                        \
        m_x1y1z1 = fx1 * (1/ry) * fz1 * m_img[mvf];                     \
        m_x2y1z1 = fx2 * (1/ry) * fz1 * m_img[mvf+1];                   \
        m_x1y2z1 = -fx1 * (1/ry) * fz1 * m_img[mvf+moving->dim[0]];             \
        m_x2y2z1 = -fx2 * (1/ry) * fz1 * m_img[mvf+moving->dim[0]+1];   \
        m_x1y1z2 = fx1 * (1/ry) * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]]; \
        m_x2y1z2 = fx2 * (1/ry) * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+1]; \
        m_x1y2z2 = -fx1 * (1/ry) * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]]; \
        m_x2y2z2 = -fx2 * (1/ry) * fz2 * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]+1]; \
        m_val = m_x1y1z1 + m_x2y1z1 + m_x1y2z1 + m_x2y2z1               \
                + m_x1y1z2 + m_x2y1z2 + m_x1y2z2 + m_x2y2z2;            \
    } while (0)

LI_VALUE_z

#define LI_VALUE_z	(		m_val,
			fx1,
			fx2,
			fy1,
			fy2,
			rz,
			mvf,
			m_img,
			moving
	)

Value:

    do {                                                                \
        float m_x1y1z1, m_x2y1z1, m_x1y2z1, m_x2y2z1;                   \
        float m_x1y1z2, m_x2y1z2, m_x1y2z2, m_x2y2z2;                   \
                                                                        \
        m_x1y1z1 = fx1 * fy1 * (1/rz) * m_img[mvf];                     \
        m_x2y1z1 = fx2 * fy1 * (1/rz) * m_img[mvf+1];                   \
        m_x1y2z1 = fx1 * fy2 * (1/rz) * m_img[mvf+moving->dim[0]];              \
        m_x2y2z1 = fx2 * fy2 * (1/rz) * m_img[mvf+moving->dim[0]+1];    \
        m_x1y1z2 = -fx1 * fy1 * (1/rz) * m_img[mvf+moving->dim[1]*moving->dim[0]]; \
        m_x2y1z2 = -fx2 * fy1 * (1/rz) * m_img[mvf+moving->dim[1]*moving->dim[0]+1]; \
        m_x1y2z2 = -fx1 * fy2 * (1/rz) * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]]; \
        m_x2y2z2 = -fx2 * fy2 * (1/rz) * m_img[mvf+moving->dim[1]*moving->dim[0]+moving->dim[0]+1]; \
        m_val = m_x1y1z1 + m_x2y1z1 + m_x1y2z1 + m_x2y2z1               \
              + m_x1y1z2 + m_x2y1z2 + m_x1y2z2 + m_x2y2z2;              \
    } while (0)