图像的降采样函数高效实现

最新推荐文章于 2024-06-12 15:59:52 发布

Kena_M

最新推荐文章于 2024-06-12 15:59:52 发布

阅读量7.6k

点赞数 2

分类专栏：图像处理文章标签：降采样图像处理图像金字塔

本文链接：https://blog.csdn.net/Kena_M/article/details/47134185

版权

图像处理专栏收录该内容

6 篇文章 0 订阅

订阅专栏

降采样函数的高效实现

所谓降采样，就是将图像的尺寸缩减为原来的一半，隔行隔列取值，降采样在构建图像金字塔时是需要反复使用的。
假如原图像尺寸为 $W\times H$ ，那么降采样之后尺寸变为 $\frac{W}{2}\times\frac{H}{2}$ ，长度和宽度是奇数还是偶数，都会自动向下取整。
另外之所以发这么一个小儿科功能的文章是因为本人患有重度强迫症，这么简单的功能OpenCV肯定有，但我就是不想用，我不知道它是如何实现的，对速度不放心。由于强迫症，我对于图像数据的访问，坚持以下几个原则：

能用一维向量就坚决不用二维的，因为连续存储的数据本来就是按照一维线性存储的，所谓二维是人类具象化的产物
能用指针就坚决不要难过数组下标，因为根据标准C，访问指针要比访问数组更高效，访问数组时涉及到下标的索引计算，这对我来说是绝不能容忍的
在使用指针一直出错烧脑时才用回数组

由此，编写了如下的函数，能看懂就尽量看吧，反正不好看，看不懂也直接可以拿来用，目前只处理了8位灰度图，太小儿科，在此露相纯粹是展示一下强迫症的魅力，就是为了追求处理的快速性，不惜脑细胞的码这种代码，反正我就认为这样的代码运行起来快速高效

unsigned char* downSampling(unsigned char *src, int width, int height, int *dwidth, int *dheight)
{
    int oddw = width % 2;
    int oddh = height % 2;
    int flag;
    *dwidth = width / 2;
    *dheight = height / 2;
    unsigned char *dst = Malloc(unsigned char, (*dwidth) * (*dheight));
    if(oddw + oddh == 0)flag = 0;
    else if(oddw + oddh == 2) flag = 3;
    else if(oddw == 0 && oddh == 1) flag = 1;
    else flag = 2;
    int i;
    unsigned char *pCur, *pEnd, *pdCur, *pdEnd;
    switch(flag)
    {
        case 0:
            for(i = 0, pCur = src, pEnd = src + width * height, pdCur = dst, pdEnd = dst + (*dwidth) * (*dheight); pdCur < pdEnd; pCur += 2, pdCur++)
            {
                i++;
                *pdCur=*pCur;
                if(i==*dwidth)
                {
                    pCur+=width;
                    i=0;
                }

            }
            break;
        case 1:
            for(i = 0, pCur = src, pEnd = src + width * height, pdCur = dst, pdEnd = dst + (*dwidth) * (*dheight); pdCur < pdEnd; pCur += 2, pdCur++)
            {
                i++;
                *pdCur=*pCur;
                if(i==*dwidth)
                {
                    pCur+=width;
                    i=0;
                }
            }
            break;
        case 2:
            for(i = 0, pCur = src, pEnd = src + width * height, pdCur = dst, pdEnd = dst + (*dwidth) * (*dheight); pdCur < pdEnd; pCur += 2, pdCur++)
            {
                i++;
                *pdCur=*pCur;
                if(i==*dwidth)
                {
                    pCur+=width+1;
                    i=0;
                }
            }
            break;
        case 3:
            for(i = 0, pCur = src, pEnd = src + width * height, pdCur = dst, pdEnd = dst + (*dwidth) * (*dheight); pdCur < pdEnd; pCur += 2, pdCur++)
            {
                i++;
                *pdCur=*pCur;
                if(i==*dwidth)
                {
                    pCur+=width+1;
                    i=0;
                }
            }
            break;
    }
    return dst;
}