📄 cxmeansdv.cpp
字号:
#define ICV_MEAN_SDV_ENTRY_BLOCK_COMMON( block_size ) \
int remaining = block_size; \
ICV_MEAN_SDV_ENTRY_COMMON()
#define ICV_MEAN_SDV_ENTRY_BLOCK_C1( sumtype, sqsumtype, \
worktype, sqworktype, block_size ) \
sumtype sum0 = 0; \
sqsumtype sqsum0 = 0; \
worktype s0 = 0; \
sqworktype sq0 = 0; \
ICV_MEAN_SDV_ENTRY_BLOCK_COMMON( block_size )
#define ICV_MEAN_SDV_ENTRY_BLOCK_C2( sumtype, sqsumtype, \
worktype, sqworktype, block_size ) \
sumtype sum0 = 0, sum1 = 0; \
sqsumtype sqsum0 = 0, sqsum1 = 0; \
worktype s0 = 0, s1 = 0; \
sqworktype sq0 = 0, sq1 = 0; \
ICV_MEAN_SDV_ENTRY_BLOCK_COMMON( block_size )
#define ICV_MEAN_SDV_ENTRY_BLOCK_C3( sumtype, sqsumtype, \
worktype, sqworktype, block_size ) \
sumtype sum0 = 0, sum1 = 0, sum2 = 0; \
sqsumtype sqsum0 = 0, sqsum1 = 0, sqsum2 = 0; \
worktype s0 = 0, s1 = 0, s2 = 0; \
sqworktype sq0 = 0, sq1 = 0, sq2 = 0; \
ICV_MEAN_SDV_ENTRY_BLOCK_COMMON( block_size )
#define ICV_MEAN_SDV_ENTRY_BLOCK_C4( sumtype, sqsumtype, \
worktype, sqworktype, block_size ) \
sumtype sum0 = 0, sum1 = 0, sum2 = 0, sum3 = 0; \
sqsumtype sqsum0 = 0, sqsum1 = 0, sqsum2 = 0, sqsum3 = 0; \
worktype s0 = 0, s1 = 0, s2 = 0, s3 = 0; \
sqworktype sq0 = 0, sq1 = 0, sq2 = 0, sq3 = 0; \
ICV_MEAN_SDV_ENTRY_BLOCK_COMMON( block_size )
/////////////////////////////////////// exit macros //////////////////////////////////////
#define ICV_MEAN_SDV_EXIT_COMMON() \
scale = pix ? 1./pix : 0
#define ICV_MEAN_SDV_EXIT_CN( total, sqtotal, idx ) \
ICV_MEAN_SDV_EXIT_COMMON(); \
mean[idx] = tmp = scale*(double)total##idx; \
tmp = scale*(double)sqtotal##idx - tmp*tmp; \
sdv[idx] = sqrt(MAX(tmp,0.))
#define ICV_MEAN_SDV_EXIT_C1( total, sqtotal ) \
ICV_MEAN_SDV_EXIT_COMMON(); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 0 )
#define ICV_MEAN_SDV_EXIT_C2( total, sqtotal ) \
ICV_MEAN_SDV_EXIT_COMMON(); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 0 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 1 )
#define ICV_MEAN_SDV_EXIT_C3( total, sqtotal ) \
ICV_MEAN_SDV_EXIT_COMMON(); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 0 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 1 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 2 )
#define ICV_MEAN_SDV_EXIT_C4( total, sqtotal ) \
ICV_MEAN_SDV_EXIT_COMMON(); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 0 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 1 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 2 ); \
ICV_MEAN_SDV_EXIT_CN( total, sqtotal, 3 )
////////////////////////////////////// update macros /////////////////////////////////////
#define ICV_MEAN_SDV_UPDATE_COMMON( block_size )\
remaining = block_size
#define ICV_MEAN_SDV_UPDATE_C1( block_size ) \
ICV_MEAN_SDV_UPDATE_COMMON( block_size ); \
sum0 += s0; sqsum0 += sq0; \
s0 = 0; sq0 = 0
#define ICV_MEAN_SDV_UPDATE_C2( block_size ) \
ICV_MEAN_SDV_UPDATE_COMMON( block_size ); \
sum0 += s0; sqsum0 += sq0; \
sum1 += s1; sqsum1 += sq1; \
s0 = s1 = 0; sq0 = sq1 = 0
#define ICV_MEAN_SDV_UPDATE_C3( block_size ) \
ICV_MEAN_SDV_UPDATE_COMMON( block_size ); \
sum0 += s0; sqsum0 += sq0; \
sum1 += s1; sqsum1 += sq1; \
sum2 += s2; sqsum2 += sq2; \
s0 = s1 = s2 = 0; sq0 = sq1 = sq2 = 0
#define ICV_MEAN_SDV_UPDATE_C4( block_size ) \
ICV_MEAN_SDV_UPDATE_COMMON( block_size ); \
sum0 += s0; sqsum0 += sq0; \
sum1 += s1; sqsum1 += sq1; \
sum2 += s2; sqsum2 += sq2; \
sum3 += s3; sqsum3 += sq3; \
s0 = s1 = s2 = s3 = 0; sq0 = sq1 = sq2 = sq3 = 0
#define ICV_DEF_MEAN_SDV_BLOCK_FUNC_2D( flavor, cn, arrtype, \
sumtype, sqsumtype, worktype, \
sqworktype, block_size, sqr_macro ) \
IPCVAPI_IMPL( CvStatus, icvMean_StdDev_##flavor##_C##cn##R, \
( const arrtype* src, int step, \
CvSize size, double* mean, double* sdv ), \
(src, step, size, mean, sdv) ) \
{ \
ICV_MEAN_SDV_ENTRY_BLOCK_C##cn( sumtype, sqsumtype, \
worktype, sqworktype, (block_size)*(cn) ); \
pix = size.width * size.height; \
size.width *= (cn); \
\
for( ; size.height--; src += step ) \
{ \
int x = 0; \
while( x < size.width ) \
{ \
int limit = MIN( remaining, size.width - x ); \
remaining -= limit; \
limit += x; \
ICV_MEAN_SDV_CASE_C##cn( worktype, sqworktype, \
sqr_macro, limit ); \
if( remaining == 0 ) \
{ \
ICV_MEAN_SDV_UPDATE_C##cn( (block_size)*(cn) ); \
} \
} \
} \
\
ICV_MEAN_SDV_UPDATE_C##cn(0); \
ICV_MEAN_SDV_EXIT_C##cn( sum, sqsum ); \
return CV_OK; \
}
#define ICV_DEF_MEAN_SDV_FUNC_2D( flavor, cn, arrtype, \
sumtype, sqsumtype, worktype ) \
IPCVAPI_IMPL( CvStatus, icvMean_StdDev_##flavor##_C##cn##R, \
( const arrtype* src, int step, \
CvSize size, double* mean, double* sdv ), \
(src, step, size, mean, sdv) ) \
{ \
ICV_MEAN_SDV_ENTRY_C##cn( sumtype, sqsumtype ); \
pix = size.width * size.height; \
size.width *= (cn); \
\
for( ; size.height--; src += step ) \
{ \
int x = 0; \
ICV_MEAN_SDV_CASE_C##cn( worktype, sqsumtype, \
CV_SQR, size.width ); \
} \
\
ICV_MEAN_SDV_EXIT_C##cn( s, sq ); \
return CV_OK; \
}
#define ICV_DEF_MEAN_SDV_BLOCK_FUNC_2D_COI( flavor, arrtype, \
sumtype, sqsumtype, worktype, \
sqworktype, block_size, sqr_macro ) \
static CvStatus CV_STDCALL icvMean_StdDev_##flavor##_CnCR \
( const arrtype* src, int step, \
CvSize size, int cn, int coi, \
double* mean, double* sdv ) \
{ \
ICV_MEAN_SDV_ENTRY_BLOCK_C1( sumtype, sqsumtype, \
worktype, sqworktype, (block_size)*(cn) ); \
pix = size.width * size.height; \
size.width *= (cn); \
src += coi - 1; \
\
for( ; size.height--; src += step ) \
{ \
int x = 0; \
while( x < size.width ) \
{ \
int limit = MIN( remaining, size.width - x ); \
remaining -= limit; \
limit += x; \
ICV_MEAN_SDV_COI_CASE( worktype, sqworktype, \
sqr_macro, limit, cn); \
if( remaining == 0 ) \
{ \
ICV_MEAN_SDV_UPDATE_C1( (block_size)*(cn) ); \
} \
} \
} \
\
ICV_MEAN_SDV_UPDATE_C1(0); \
ICV_MEAN_SDV_EXIT_C1( sum, sqsum ); \
return CV_OK; \
}
#define ICV_DEF_MEAN_SDV_FUNC_2D_COI( flavor, arrtype, \
sumtype, sqsumtype, worktype )\
static CvStatus CV_STDCALL icvMean_StdDev_##flavor##_CnCR \
( const arrtype* src, int step, CvSize size,\
int cn, int coi, double* mean, double* sdv )\
{ \
ICV_MEAN_SDV_ENTRY_C1( sumtype, sqsumtype ); \
pix = size.width * size.height; \
size.width *= (cn); \
src += coi - 1; \
\
for( ; size.height--; src += step ) \
{ \
int x = 0; \
ICV_MEAN_SDV_COI_CASE( worktype, sqsumtype, \
CV_SQR, size.width, cn ); \
} \
\
ICV_MEAN_SDV_EXIT_C1( s, sq ); \
return CV_OK; \
}
#define ICV_DEF_MEAN_SDV_MASK_BLOCK_FUNC_2D( flavor, cn, \
arrtype, sumtype, sqsumtype, worktype, \
sqworktype, block_size, sqr_macro ) \
IPCVAPI_IMPL( CvStatus, icvMean_StdDev_##flavor##_C##cn##MR, \
( const arrtype* src, int step, \
const uchar* mask, int maskstep, \
CvSize size, double* mean, double* sdv ), \
(src, step, mask, maskstep, size, mean, sdv))\
{ \
ICV_MEAN_SDV_ENTRY_BLOCK_C##cn( sumtype, sqsumtype, \
worktype, sqworktype, block_size ); \
pix = 0; \
\
for( ; size.height--; src += step, mask += maskstep ) \
{ \
int x = 0; \
while( x < size.width ) \
{ \
int limit = MIN( remaining, size.width - x ); \
remaining -= limit; \
limit += x; \
ICV_MEAN_SDV_MASK_CASE_C##cn( worktype, sqworktype, \
sqr_macro, limit ); \
if( remaining == 0 ) \
{ \
ICV_MEAN_SDV_UPDATE_C##cn( block_size ); \
} \
} \
} \
\
ICV_MEAN_SDV_UPDATE_C##cn(0); \
ICV_MEAN_SDV_EXIT_C##cn( sum, sqsum ); \
return CV_OK; \
}
#define ICV_DEF_MEAN_SDV_MASK_FUNC_2D( flavor, cn, arrtype, \
sumtype, sqsumtype, worktype)\
IPCVAPI_IMPL( CvStatus, icvMean_StdDev_##flavor##_C##cn##MR, \
( const arrtype* src, int step, \
const uchar* mask, int maskstep, \
CvSize size, double* mean, double* sdv ), \
(src, step, mask, maskstep, size, mean, sdv))\
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -