Chinaunix首页 | 论坛 | 博客

锦运_2012的ChinaUnix博客

首页　| 　博文目录　| 　关于我

锦运_2012

  • 博客访问： 263420
  • 博文数量： 28
  • 博客积分： 688
  • 博客等级： 上士
  • 技术积分： 365
  • 用 户 组： 普通用户
  • 注册时间： 2012-03-16 10:25
文章分类

全部博文（28）

文章存档

2012年（28）

我的朋友
最近访客
推荐博文
《学习 OpenCV》第六章 214页练习 3

分类：

2012-09-17 10:59:34

原文地址：《学习 OpenCV》第六章 214页练习 3 作者：duanlin

3. 你能从图6-5显示的滤波器里建立一个可分核吗？如果可以，显示一下？
   图6-5如下所示。
 
 
核的部分点有负数，想求出与该核效果等同的分核，则需要知道这些数与图像中的像素数据计算时相互之间是什么具体关系，计算是如何进行的。因此，解答此题目我认为要查看OpenCV源代码中cvSobel()函数的具体实现，至少要搞清楚当aperture_size等于CV_SCHARR时的计算方式。
于是找到cvSobel的具体实现位置：cvderiv.cpp
  1. CV_IMPL void
  2. cvSobel( const void* srcarr, void* dstarr, int dx, int dy, int aperture_size )
  3. {
  4.     CvSepFilter filter;
  5.     void* buffer = 0;
  6.     int local_alloc = 0;

  8.     CV_FUNCNAME( "cvSobel" );

  10.     __BEGIN__;

  12.     int origin = 0;
  13.     int src_type, dst_type;
  14.     CvMat srcstub, *src = (CvMat*)srcarr;
  15.     CvMat dststub, *dst = (CvMat*)dstarr;

  17.     if( !CV_IS_MAT(src) )
  18.         CV_CALL( src = cvGetMat( src, &srcstub ));
  19.     if( !CV_IS_MAT(dst) )
  20.         CV_CALL( dst = cvGetMat( dst, &dststub ));

  22.     if( CV_IS_IMAGE_HDR( srcarr ))
  23.         origin = ((IplImage*)srcarr)->origin;

  25.     src_type = CV_MAT_TYPE( src->type );
  26.     dst_type = CV_MAT_TYPE( dst->type );

  28.     if( !CV_ARE_SIZES_EQ( src, dst ))
  29.         CV_ERROR( CV_StsBadArg, "src and dst have different sizes" );

  31.     if( ((aperture_size == CV_SCHARR || aperture_size == 3 || aperture_size == 5) &&
  32.         dx <= 2 && dy <= 2 && dx + dy <= 2 && icvFilterSobelNegVertBorder_8u16s_C1R_p) &&
  33.         (src_type == CV_8UC1 && dst_type == CV_16SC1 ||
  34.         src_type == CV_32FC1 && dst_type == CV_32FC1) )
  35.     {
  36.         CvDerivGetBufSizeIPPFunc ipp_sobel_getbufsize_func = 0;
  37.         CvDerivIPPFunc_8u ipp_sobel_func_8u = 0;
  38.         CvDerivIPPFunc_32f ipp_sobel_func_32f = 0;

  40.         CvDeriv3x3GetBufSizeIPPFunc ipp_scharr_getbufsize_func = 0;
  41.         CvDeriv3x3IPPFunc_8u ipp_scharr_func_8u = 0;
  42.         CvDeriv3x3IPPFunc_32f ipp_scharr_func_32f = 0;

  44.         if( aperture_size == CV_SCHARR )
  45.         {
  46.             if( dx == 1 && dy == 0 )
  47.             {
  48.                 if( src_type == CV_8U )
  49.                     ipp_scharr_func_8u = icvFilterScharrVertBorder_8u16s_C1R_p,
  50.                     ipp_scharr_getbufsize_func = icvFilterScharrVertGetBufSize_8u16s_C1R_p;
  51.                 else
  52.                     ipp_scharr_func_32f = icvFilterScharrVertBorder_32f_C1R_p,
  53.                     ipp_scharr_getbufsize_func = icvFilterScharrVertGetBufSize_32f_C1R_p;
  54.             }
  55.             else if( dx == 0 && dy == 1 )
  56.             {
  57.                 if( src_type == CV_8U )
  58.                     ipp_scharr_func_8u = icvFilterScharrHorizBorder_8u16s_C1R_p,
  59.                     ipp_scharr_getbufsize_func = icvFilterScharrHorizGetBufSize_8u16s_C1R_p;
  60.                 else
  61.                     ipp_scharr_func_32f = icvFilterScharrHorizBorder_32f_C1R_p,
  62.                     ipp_scharr_getbufsize_func = icvFilterScharrHorizGetBufSize_32f_C1R_p;
  63.             }
  64.             else
  65.                 CV_ERROR( CV_StsBadArg, "Scharr filter can only be used to compute 1st image derivatives" );
  66.         }
  67.         else
  68.         {
  69.             if( dx == 1 && dy == 0 )
  70.             {
  71.                 if( src_type == CV_8U )
  72.                     ipp_sobel_func_8u = icvFilterSobelNegVertBorder_8u16s_C1R_p,
  73.                     ipp_sobel_getbufsize_func = icvFilterSobelNegVertGetBufSize_8u16s_C1R_p;
  74.                 else
  75.                     ipp_sobel_func_32f = icvFilterSobelNegVertBorder_32f_C1R_p,
  76.                     ipp_sobel_getbufsize_func = icvFilterSobelNegVertGetBufSize_32f_C1R_p;
  77.             }
  78.             else if( dx == 0 && dy == 1 )
  79.             {
  80.                 if( src_type == CV_8U )
  81.                     ipp_sobel_func_8u = icvFilterSobelHorizBorder_8u16s_C1R_p,
  82.                     ipp_sobel_getbufsize_func = icvFilterSobelHorizGetBufSize_8u16s_C1R_p;
  83.                 else
  84.                     ipp_sobel_func_32f = icvFilterSobelHorizBorder_32f_C1R_p,
  85.                     ipp_sobel_getbufsize_func = icvFilterSobelHorizGetBufSize_32f_C1R_p;
  86.             }
  87.             else if( dx == 2 && dy == 0 )
  88.             {
  89.                 if( src_type == CV_8U )
  90.                     ipp_sobel_func_8u = icvFilterSobelVertSecondBorder_8u16s_C1R_p,
  91.                     ipp_sobel_getbufsize_func = icvFilterSobelVertSecondGetBufSize_8u16s_C1R_p;
  92.                 else
  93.                     ipp_sobel_func_32f = icvFilterSobelVertSecondBorder_32f_C1R_p,
  94.                     ipp_sobel_getbufsize_func = icvFilterSobelVertSecondGetBufSize_32f_C1R_p;
  95.             }
  96.             else if( dx == 0 && dy == 2 )
  97.             {
  98.                 if( src_type == CV_8U )
  99.                     ipp_sobel_func_8u = icvFilterSobelHorizSecondBorder_8u16s_C1R_p,
  100.                     ipp_sobel_getbufsize_func = icvFilterSobelHorizSecondGetBufSize_8u16s_C1R_p;
  101.                 else
  102.                     ipp_sobel_func_32f = icvFilterSobelHorizSecondBorder_32f_C1R_p,
  103.                     ipp_sobel_getbufsize_func = icvFilterSobelHorizSecondGetBufSize_32f_C1R_p;
  104.             }
  105.             else if( dx == 1 && dy == 1 )
  106.             {
  107.                 if( src_type == CV_8U )
  108.                     ipp_sobel_func_8u = icvFilterSobelCrossBorder_8u16s_C1R_p,
  109.                     ipp_sobel_getbufsize_func = icvFilterSobelCrossGetBufSize_8u16s_C1R_p;
  110.                 else
  111.                     ipp_sobel_func_32f = icvFilterSobelCrossBorder_32f_C1R_p,
  112.                     ipp_sobel_getbufsize_func = icvFilterSobelCrossGetBufSize_32f_C1R_p;
  113.             }
  114.         }

  116.         if( (ipp_sobel_func_8u || ipp_sobel_func_32f) && ipp_sobel_getbufsize_func ||
  117.             (ipp_scharr_func_8u || ipp_scharr_func_32f) && ipp_scharr_getbufsize_func )
  118.         {
  119.             int bufsize = 0, masksize = aperture_size == 3 ? 33 : 55;
  120.             CvSize size = cvGetMatSize( src );
  121.             uchar* src_ptr = src->data.ptr;
  122.             uchar* dst_ptr = dst->data.ptr;
  123.             int src_step = src->step ? src->step : CV_STUB_STEP;
  124.             int dst_step = dst->step ? dst->step : CV_STUB_STEP;
  125.             const int bordertype = 1; // replication border
  126.             CvStatus status;

  128.             status = ipp_sobel_getbufsize_func ?
  129.                 ipp_sobel_getbufsize_func( size, masksize, &bufsize ) :
  130.                 ipp_scharr_getbufsize_func( size, &bufsize );

  132.             if( status >= 0 )
  133.             {
  134.                 if( bufsize <= CV_MAX_LOCAL_SIZE )
  135.                 {
  136.                     buffer = cvStackAlloc( bufsize );
  137.                     local_alloc = 1;
  138.                 }
  139.                 else
  140.                     CV_CALL( buffer = cvAlloc( bufsize ));
  141.             
  142.                 status =
  143.                     ipp_sobel_func_8u ? ipp_sobel_func_8u( src_ptr, src_step, dst_ptr, dst_step,
  144.                                                            size, masksize, bordertype, 0, buffer ) :
  145.                     ipp_sobel_func_32f ? ipp_sobel_func_32f( src_ptr, src_step, dst_ptr, dst_step,
  146.                                                              size, masksize, bordertype, 0, buffer ) :
  147.                     ipp_scharr_func_8u ? ipp_scharr_func_8u( src_ptr, src_step, dst_ptr, dst_step,
  148.                                                              size, bordertype, 0, buffer ) :
  149.                     ipp_scharr_func_32f ? ipp_scharr_func_32f( src_ptr, src_step, dst_ptr, dst_step,
  150.                                                                size, bordertype, 0, buffer ) :
  151.                         CV_NOTDEFINED_ERR;
  152.             }

  154.             if( status >= 0 &&
  155.                 (dx == 0 && dy == 1 && origin || dx == 1 && dy == 1 && !origin)) // negate the output
  156.                 cvSubRS( dst, cvScalarAll(0), dst );

  158.             if( status >= 0 )
  159.                 EXIT;
  160.         }
  161.     }

  163.     CV_CALL( filter.init_deriv( src->cols, src_type, dst_type, dx, dy,
  164.                 aperture_size, origin ? CvSepFilter::FLIP_KERNEL : 0));
  165.     CV_CALL( filter.process( src, dst ));

  167.     __END__;

  169.     if( buffer && !local_alloc )
  170.         cvFree( &buffer );
  171. }
其中分析关键部分
  1. if( aperture_size == CV_SCHARR )
  2.         {
  3.             if( dx == 1 && dy == 0 )
  4.             {
  5.                 if( src_type == CV_8U )
  6.                     ipp_scharr_func_8u = icvFilterScharrVertBorder_8u16s_C1R_p,
  7.                     ipp_scharr_getbufsize_func = icvFilterScharrVertGetBufSize_8u16s_C1R_p;
  8.                 else
  9.                     ipp_scharr_func_32f = icvFilterScharrVertBorder_32f_C1R_p,
  10.                     ipp_scharr_getbufsize_func = icvFilterScharrVertGetBufSize_32f_C1R_p;
  11.             }
  12.             else if( dx == 0 && dy == 1 )
  13.             {
  14.                 if( src_type == CV_8U )
  15.                     ipp_scharr_func_8u = icvFilterScharrHorizBorder_8u16s_C1R_p,
  16.                     ipp_scharr_getbufsize_func = icvFilterScharrHorizGetBufSize_8u16s_C1R_p;
  17.                 else
  18.                     ipp_scharr_func_32f = icvFilterScharrHorizBorder_32f_C1R_p,
  19.                     ipp_scharr_getbufsize_func = icvFilterScharrHorizGetBufSize_32f_C1R_p;
  20.             }
  21.             else
  22.                 CV_ERROR( CV_StsBadArg, "Scharr filter can only be used to compute 1st image derivatives" );
  23.         }
其中dx和dy是在x和y这两个方向上求导的阶数。通常为0和1，最多2。0表示这个方向上没有求导。如果调用cvSobel()并且传递CV_SCHARR后，就会进入该if条件。从此处可以看出，该版本opencv只能支持x和y方向上的0或1阶导数。进入该if条件后只对ipp_scharr_func_8u或ipp_scharr_func_32f和ipp_scharr_getbufsize_func这两个变量赋了一个值。真正的计算在后边。
 
  1. if( (ipp_sobel_func_8u || ipp_sobel_func_32f) && ipp_sobel_getbufsize_func ||
  2.             (ipp_scharr_func_8u || ipp_scharr_func_32f) && ipp_scharr_getbufsize_func )
  3.         {
  4.             int bufsize = 0, masksize = aperture_size == 3 ? 33 : 55;
  5.             CvSize size = cvGetMatSize( src );
  6.             uchar* src_ptr = src->data.ptr;
  7.             uchar* dst_ptr = dst->data.ptr;
  8.             int src_step = src->step ? src->step : CV_STUB_STEP;
  9.             int dst_step = dst->step ? dst->step : CV_STUB_STEP;
  10.             const int bordertype = 1; // replication border
  11.             CvStatus status;

  13.             status = ipp_sobel_getbufsize_func ?
  14.                 ipp_sobel_getbufsize_func( size, masksize, &bufsize ) :
  15.                 ipp_scharr_getbufsize_func( size, &bufsize );

  17.             if( status >= 0 )
  18.             {
  19.                 if( bufsize <= CV_MAX_LOCAL_SIZE )
  20.                 {
  21.                     buffer = cvStackAlloc( bufsize );
  22.                     local_alloc = 1;
  23.                 }
  24.                 else
  25.                     CV_CALL( buffer = cvAlloc( bufsize ));
  26.             
  27.                 status =
  28.                     ipp_sobel_func_8u ? ipp_sobel_func_8u( src_ptr, src_step, dst_ptr, dst_step,
  29.                                                            size, masksize, bordertype, 0, buffer ) :
  30.                     ipp_sobel_func_32f ? ipp_sobel_func_32f( src_ptr, src_step, dst_ptr, dst_step,
  31.                                                              size, masksize, bordertype, 0, buffer ) :
  32.                     ipp_scharr_func_8u ? ipp_scharr_func_8u( src_ptr, src_step, dst_ptr, dst_step,
  33.                                                              size, bordertype, 0, buffer ) :
  34.                     ipp_scharr_func_32f ? ipp_scharr_func_32f( src_ptr, src_step, dst_ptr, dst_step,
  35.                                                                size, bordertype, 0, buffer ) :
  36.                         CV_NOTDEFINED_ERR;
  37.             }

  39.             if( status >= 0 &&
  40.                 (dx == 0 && dy == 1 && origin || dx == 1 && dy == 1 && !origin)) // negate the output
  41.                 cvSubRS( dst, cvScalarAll(0), dst );

  43.             if( status >= 0 )
  44.                 EXIT;
  45.         }
由返回status值的ipp_scharr_func_8u或ipp_scharr_func_32f计算。这两个是函数指针变量，拿ipp_scharr_func_8u来讲，他被赋值为icvFilterScharrHorizBorder_8u16s_C1R_p，而icvFilterScharrHorizBorder_8u16s_C1R_p由icvFilterScharrHorizBorder_32f_C1R_t定义。最终通过[cxmisc.h]IPCVAPI_EX()转到[cxmisc.h]IPCVAPI(type,declspec,name,args)调用
extern name##_t name##_p。
------------
未完待续
------------
 
阅读(1856) | 评论(0) | 转发(0) |
0

上一篇：fatal error C1083: 无法打开包括文件:“iostream.h”:

下一篇：《学习 OpenCV》第六章 214页练习 4

给主人留下些什么吧！~~
关于我们 | 关于IT168 | 联系方式 | 广告合作 | 法律声明 | 免费注册

Copyright 2001-2010 ChinaUnix.net All Rights Reserved 北京皓辰网域网络信息技术有限公司. 版权所有

感谢所有关心和支持过ChinaUnix的朋友们

16024965号-6