void ConvolveSeparable(CImageOf<T> src, CImageOf<T>& dst,
CFloatImage x_kernel, CFloatImage y_kernel,
int subsample)
{
// Allocate the result, if necessary
CShape dShape = src.Shape();
if (subsample > 1) {
dShape.width = (dShape.width + subsample - 1) / subsample;
dShape.height = (dShape.height + subsample - 1) / subsample;
}
dst.ReAllocate(dShape, false);
// Allocate the intermediate (or final) result
CImageOf<T> tmpImg;
if (subsample > 1)
tmpImg.ReAllocate(src.Shape());
CImageOf<T>& tmp = (subsample > 1) ? tmpImg : dst;
// Create the IPL images
IplImage* srcImage = IPLCreateImage(src);
IplImage* dstImage = IPLCreateImage(tmp);
srcImage->alphaChannel = 0; // convolve the A channel also
dstImage->alphaChannel = 0; // convolve the A channel also
// Call the convolution code
if (typeid(T) == typeid(float)) {
IplConvKernelFP* xKernel = IPLConvKernelFP(x_kernel, false);
IplConvKernelFP* yKernel = IPLConvKernelFP(y_kernel, true);
iplConvolveSep2DFP(srcImage, dstImage, xKernel, yKernel);
iplDeleteConvKernelFP(xKernel);
iplDeleteConvKernelFP(yKernel);
} else {
IplConvKernel* xKernel = IPLConvKernel(x_kernel, false);
IplConvKernel* yKernel = IPLConvKernel(y_kernel, true);
iplConvolveSep2D(srcImage, dstImage, xKernel, yKernel);
iplDeleteConvKernel(xKernel);
iplDeleteConvKernel(yKernel);
}
iplDeallocate(srcImage, IPL_IMAGE_HEADER);
iplDeallocate(dstImage, IPL_IMAGE_HEADER);
// Downsample if necessary
if (subsample > 1) {
for (int y = 0; y < dShape.height; y++) {
T* sPtr = &tmp.Pixel(0, y * subsample, 0);
T* dPtr = &dst.Pixel(0, y, 0);
int nB = dShape.nBands;
for (int x = 0; x < dShape.width; x++) {
for (int b = 0; b < nB; b++)
dPtr[b] = sPtr[b];
sPtr += subsample * nB;
dPtr += nB;
}
}
}
}
void Convolve(CImageOf<T> src, CImageOf<T>& dst,
CFloatImage kernel)
{
// Allocate the result, if necessary
dst.ReAllocate(src.Shape(), false);
// Create the IPL images
IplImage* srcImage = IPLCreateImage(src);
IplImage* dstImage = IPLCreateImage(dst);
// Call the convolution code
if (typeid(T) == typeid(float))
{
IplConvKernelFP* iplKernel = IPLConvKernelFP(kernel);
iplConvolve2DFP(srcImage, dstImage, &iplKernel, 1, IPL_SUM);
iplDeleteConvKernelFP(iplKernel);
}
else
{
IplConvKernel* iplKernel = IPLConvKernel(kernel);
iplConvolve2D(srcImage, dstImage, &iplKernel, 1, IPL_SUM);
iplDeleteConvKernel(iplKernel);
}
iplDeallocate(srcImage, IPL_IMAGE_HEADER);
iplDeallocate(dstImage, IPL_IMAGE_HEADER);
}
void YARPLpFirstOrderFlow::_free_kernels(void)
{
iplDeleteConvKernelFP (DerivX);
iplDeleteConvKernelFP (DerivY);
iplDeleteConvKernelFP (Gauss);
}
static int fcaScharr( void )
{
int KerX[15];
int KerY[15];
CvSize KerLens;
int dx,dy;
AtsRandState state;
AtsRandState sizegen;
double Error = 0;
IplImage* src8u;
IplImage* src8s;
IplImage* src32f;
IplImage* src16;
IplImage* dst16;
IplImage* test16;
IplImage* dst32f;
IplImage* test32f;
IplImage* buf;
IplConvKernel* KernelX;
IplConvKernel* KernelY;
IplConvKernelFP* KX;
IplConvKernelFP* KY;
/* Initialization global parameters */
if( !read_param )
{
read_param = 1;
/* Determining which test are needed to run */
trsCaseRead( &data_types,"/u/s/f/a", "a",
"u - unsigned char, s - signed char, f - float, a - all" );
/* Reading test-parameters */
trsiRead( &lMinImageSize, "16", "Image height" );
trsiRead( &lMaxImageSize, "256", "Image width" );
trsiRead( &lTestNum, "20", "Test count" );
}
if( data_types != 3 && data_types != 0 ) return TRS_UNDEF;
atsRandInit(&sizegen,lMinImageSize,lMaxImageSize,0);
/* Creating images for testing */
for(int i = 0; i < lTestNum; i++)
{
lImageWidth = atsRand32s(&sizegen);
lImageHeight = atsRand32s(&sizegen);
src8u= cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_8U, 1);
src8s= cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_8S, 1);
src32f = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_32F, 1);
src16 = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_16S, 1);
dst16 = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_16S, 1);
test16 = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_16S, 1);
buf = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_32F, 1);
dst32f = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_32F, 1);
test32f = cvCreateImage(cvSize(lImageWidth, lImageHeight), IPL_DEPTH_32F, 1);
atsRandInit(&state,0,64,32);
atsFillRandomImageEx(src8u, &state );
atsFillRandomImageEx(src8s, &state );
atsFillRandomImageEx(src32f, &state );
/* Calculating the kernels */
for(dx = 0; dx<=1; dx++)
{
dy=1-dx;
atsCalcKernel(cv8u,dx,dy,CV_SCHARR,(char*)KerX,(char*)KerY,&KerLens,CV_ORIGIN_TL);
KernelX = iplCreateConvKernel(KerLens.width,1,KerLens.width/2,0,KerX,0);
KernelY = iplCreateConvKernel(1,KerLens.height,0,KerLens.height/2,KerY,0);
/* Calculating the convolution */
atsConvert(src8u,src16);
iplSetBorderMode(src16,
IPL_BORDER_REPLICATE,
IPL_SIDE_TOP|IPL_SIDE_BOTTOM|IPL_SIDE_RIGHT|IPL_SIDE_LEFT,
0);
iplConvolveSep2D(src16,dst16,KernelX,KernelY);
cvSobel(src8u,test16,dx,dy,CV_SCHARR);
Error += iplNorm(dst16,test16,IPL_C);
atsConvert(src8s,src16);
iplSetBorderMode(src16,
IPL_BORDER_REPLICATE,
IPL_SIDE_TOP|IPL_SIDE_BOTTOM|IPL_SIDE_RIGHT|IPL_SIDE_LEFT,
0);
iplSetBorderMode(src16,
IPL_BORDER_REPLICATE,
IPL_SIDE_TOP|IPL_SIDE_BOTTOM|IPL_SIDE_RIGHT|IPL_SIDE_LEFT,
0);
iplConvolveSep2D(src16,dst16,KernelX,KernelY);
cvSobel(src8s,test16,dx,dy,CV_SCHARR);
Error += iplNorm(dst16,test16,IPL_C);
atsCalcKernel(IPL_DEPTH_32F,dx,dy,CV_SCHARR,(char*)KerX,(char*)KerY,&KerLens,CV_ORIGIN_TL);
KX = iplCreateConvKernelFP(KerLens.width,1,KerLens.width/2,0,(float*)KerX);
KY = iplCreateConvKernelFP(1,KerLens.height,0,KerLens.height/2,(float*)KerY);
iplSetBorderMode(src32f,
IPL_BORDER_REPLICATE,
IPL_SIDE_TOP|IPL_SIDE_BOTTOM|IPL_SIDE_RIGHT|IPL_SIDE_LEFT,
0);
/* Calculating the convolution */
iplConvolveSep2DFP(src32f,dst32f,KX,KY);
cvSobel(src32f,test32f,dx,dy,CV_SCHARR);
/*for(i = 0; i<lImageHeight; i++)
for(int j = 0; j<lImageWidth; j++)
{
float a = ((float*)(dst32f->imageData))[i*dst32f->widthStep/4+j];
float b = ((float*)(test32f->imageData))[i*test32f->widthStep/4+j];
}
*/
Error += iplNorm(test32f,dst32f,IPL_C);
iplDeleteConvKernel(KernelX);
iplDeleteConvKernel(KernelY);
iplDeleteConvKernelFP(KX);
iplDeleteConvKernelFP(KY);
}
/* Free Memory */
cvReleaseImage( &src8u );
cvReleaseImage( &src8s );
cvReleaseImage( &src32f );
cvReleaseImage( &src16 );
cvReleaseImage( &dst16 );
cvReleaseImage( &test16 );
cvReleaseImage( &dst32f );
cvReleaseImage( &test32f );
cvReleaseImage( &buf);
}
if(Error/lTestNum>=EPSILON)return TRS_FAIL;
return TRS_OK;
} /* fcaSobel8uC1R */
