volume_enhance_contrast.cpp File Reference

#include "volume_enhance_contrast.h"

Include dependency graph for volume_enhance_contrast.cpp:

Go to the source code of this file.

Macros
#define	COMPUTE_STATISTICS(V, N, k, i, j, avg, stddev, cubeSize)
#define	COMPUTE_MEGV(VxE, V_E, N, k, i, j, pixel)

Macro Definition Documentation

COMPUTE_MEGV

#define COMPUTE_MEGV	(		VxE,
			V_E,
			N,
			k,
			i,
			j,
			pixel
	)

Value:

{\
   int k0,kF,i0,iF,j0,jF; \
     k0=XMIPP_MAX(k-N,STARTINGZ(V_E)); \
     kF=XMIPP_MIN(k+N,FINISHINGZ(V_E)); \
     i0=XMIPP_MAX(i-N,STARTINGY(V_E)); \
     iF=XMIPP_MIN(i+N,FINISHINGY(V_E)); \
     j0=XMIPP_MAX(j-N,STARTINGX(V_E)); \
         jF=XMIPP_MIN(j+N,FINISHINGX(V_E)); \
         double sum1=0, sum2=0; \
         for (int kk=k0; kk<=kF; kk++) \
    for (int ii=i0; ii<=iF; ii++) \
     for (int jj=j0; jj<=jF; jj++) \
     { \
      double v1=A3D_ELEM(VxE,kk,ii,jj); \
      sum1=sum1+v1; \
      double v2=A3D_ELEM(V_E,kk,ii,jj); \
      sum2=sum2+v2; \
     } \
   pixel=sum1/sum2; \
  }

COMPUTE_STATISTICS

#define COMPUTE_STATISTICS	(		V,
			N,
			k,
			i,
			j,
			avg,
			stddev,
			cubeSize
	)

Value:

{\
   int k0,kF,i0,iF,j0,jF; \
     k0=XMIPP_MAX(k-N,STARTINGZ(V)); \
     kF=XMIPP_MIN(k+N,FINISHINGZ(V)); \
     i0=XMIPP_MAX(i-N,STARTINGY(V)); \
     iF=XMIPP_MIN(i+N,FINISHINGY(V)); \
     j0=XMIPP_MAX(j-N,STARTINGX(V)); \
         jF=XMIPP_MIN(j+N,FINISHINGX(V)); \
         double sum=0, sum2=0; \
         for (int kk=k0; kk<=kF; ++kk) \
    for (int ii=i0; ii<=iF; ++ii) \
     for (int jj=j0; jj<=jF; ++jj) \
     { \
      double v=A3D_ELEM(V,kk,ii,jj); \
      sum+=v; \
      sum2+=v*v; \
     } \
   cubeSize=(kF-k0+1)*(iF-i0+1)*(jF-j0+1); \
   avg=sum/cubeSize; \
   stddev=sum2/cubeSize-avg*avg; \
   stddev=sqrt(XMIPP_MAX(stddev,0)); \
  }