这两天编程实现同态滤波，可实现的结果却令我大跌眼镜，滤波后的图像严重发生了错位，简直分辨不出图像的内容，检查程序没有发现错误，这让我一直很郁闷。今天早上在百度上搜到一个帖子，帖子是别人转的，但原帖子的作者也遇到过这样的错误。作者说宽度为偶数的图片不会发生这样的情况，但宽度为奇数的图片就会发生上面的错误。我也拿了几张图片试验了一下，正如作者所说。究其原因，原来是IplImage的widthstep在作怪，当width为偶数时，widthstep与width*nchannels相等，但当width为奇数时，就会填充一个单位像素的长度使widthstep成为偶数。看了帖子，我也找到了出错的原因，我就马上修改程序。呵呵，结果正常显示了。只要把读取数据的i*width+j改为i*widthstep+j就可以了。
width是图像宽度，可为任意值；widthstep是行字节数，应该是4的倍数，不一定等于width，nchannels为图像通道数。 #define WIDTHBYTES(bits) (((bits)+31)/32*4) 看看这个公式你就明白了。
对ROI和widthStep的补充
ROI和widthStep在实际工作中有很重要的作用，在很多情况下，使用它们会提高计算机视觉代码的执行速度。这是因为它们允许对图像的某一小部分进行操作，而不是对整个图像进行运算。在OpenCV中，普遍支持ROI和widthStep，函数的操作被限于感兴趣区域。要设置或取消ROI，就要使用cvSetImageROI()和cvResetImageROI()函数。如果想设置ROI，可以使用函数cvSetImageROI()，并为其传递一个图像指针和矩形。而取消ROI，只需要为函数cvResetImageROI()传递一个图像指针。
- void cvSetImageROI( IplImage* image, CvRect rect );
- void cvResetImageROI( IplImage* image );
为了解释ROI的用法，我们假设要加载一幅图像并修改一些区域，如例3-12的代码，读取了一幅图像，并设置了想要的ROI的x，y，width和height的值，最后将ROI区域中像素都加上一个整数。本例程中通过内联的cvRect()构造函数设置ROI。通过cvResetImageROI()函数释放ROI是非常重要的，否则，将忠实地只显示ROI区域。
例3-12：用imageROI来增加某范围的像素
- // roi_add <image> <x> <y> <width> <height> <add>
- #include <cv.h>
- #include <highgui.h>
- int main(int argc, char** argv)
- {
- IplImage* src;
- if( argc == 7 && ((src=cvLoadImage(argv[1],1)) != 0 ))
- {
- int x = atoi(argv[2]);
- int y = atoi(argv[3]);
- int width = atoi(argv[4]);
- int height = atoi(argv[5]);
- int add = atoi(argv[6]);
- cvSetImage ROI(src, cvRect(x,y,width,height));
- cvAddS(src, cvScalar(add),src);
- cvResetImageROI(src);
- cvNamedWindow( "Roi_Add", 1 );
- cvShowImage( "Roi_Add", src );
- cvWaitKey();
- }
- return 0;
- }
使用例3-12中的代码把ROI集中于一张猫的脸部，并将其蓝色通道增加150后的效果如图3-3所示。【45～46】
 |
| （点击查看大图）图3-3：在猫脸上用ROI增加150像素的效果 |
通过巧妙地使用widthStep，我们可以达到同样的效果。要做到这一点，我们创建另一个图像头，让它的width和height的值等于interest_rect的width和height的值。我们还需要按interest_rect起点设置图像起点(左上角或者左下角)。下一步，我们设置子图像的widthStep与较大的interest_img相同。这样，即可在子图像中逐行地步进到大图像里子区域中下一行开始处的合适位置。最后设置子图像的imageDate指针指向兴趣子区域的开始，如例3-13所示。
例3-13：利用其他widthStep方法把interest_img的所有像素值增加1
- // Assuming IplImage *interest_img; and
- // CvRect interest_rect;
- // Use widthStep to get a region of interest
- //
- // (Alternate method)
- //
- IplImage *sub_img = cvCreateImageHeader(
- cvSize(
- interest_rect.width,
- interest_rect.height
- ),
- interest_img->depth,
- interest_img->nChannels
- );
- sub_img->origin = interest_img->origin;
- sub_img->widthStep = interest_img->widthStep;
- sub_img->imageData = interest_img->imageData +
- interest_rect.y * interest_img->widthStep +
- interest_rect.x * interest_img->nChannels;
- cvAddS( sub_img, cvScalar(1), sub_img );
- cvReleaseImageHeader(&sub_img);
看起来设置和重置ROI更方便一些，为什么还要使用widthStep？原因在于有些时候在处理的过程中，想在操作过程中设置和保持一幅图像的多个子区域处于活动状态，但是ROI只能串行处理并且必须不断地设置和重置。
最后，我们要在此提到一个词-- 掩码或模板，在代码示例中cvAddS()函数允许第四个参数默认值为空：const CvArr* mask=NULL。这是一个8位单通道数组，它允许把操作限制到任意形状的非0像素的掩码区，如果ROI随着掩码或模板变化，进程将会被限制在ROI和掩码的交集区域。掩码或模板只能在指定了其图像的函数中使用。