此为opencv中remap函数移植和分析,整理了双线性的插值部分的代码
尚未完全移植,但最困难的部分已经完成,而恰巧在这时,发现其实现并不是那么的令我满意,于是终止,改为自己实现。
考虑到以后可能会用到,暂且保存于此。

#include "./CCC.h"

#ifndef __DCC_H__
#define __DCC_H__

#define saturate_cast_uchar(v) (uchar)((unsigned)(v) <= UCHAR_MAX ? (v) : (v) > 0 ? UCHAR_MAX : 0)
#define saturate_cast_short(v) (short)((unsigned)((v) - SHRT_MIN) <= (unsigned)USHRT_MAX ? (v) : (v) > 0 ? SHRT_MAX : SHRT_MIN)

#define INTER_REMAP_COEF_BITS 15
#define INTER_REMAP_COEF_SCALE 32768 //(1 << INTER_REMAP_COEF_BITS)
#define SHIFT INTER_REMAP_COEF_BITS
#define DELTA (1 << (INTER_REMAP_COEF_BITS - 1))

 + ][][];//1024=tabsz2=INTER_TAB_SIZE2
][] = (][])alignPtr(BilinearTab_iC4_buf, );

struct RemapVec_8u
{
    int operator()(const Mat& _src, void* _dst, const short* XY, const ushort* FXY, const void* _wtab, int width) const
    {
        , sstep = (int)_src.step;

        );
         ? (][][];
        uchar* D = (uchar*)_dst;
        __m128i delta = _mm_set1_epi32(INTER_REMAP_COEF_SCALE / );
        __m128i xy2ofs = _mm_set1_epi32(cn + (sstep << ));
        __m128i z = _mm_setzero_si128();
        ) iofs0[], iofs1[];

        )
        {
            ; x += )
            {
                __m128i xy0 = _mm_loadu_si128(());
                __m128i xy1 = _mm_loadu_si128(( + ));
                __m128i v0, v1, v2, v3, a0, a1, b0, b1;
                unsigned i0, i1;

                xy0 = _mm_madd_epi16(xy0, xy2ofs);
                xy1 = _mm_madd_epi16(xy1, xy2ofs);
                _mm_store_si128((__m128i*)iofs0, xy0);
                _mm_store_si128((__m128i*)iofs1, xy1);

                i0 = *(]) + (*(]) << );
                i1 = *(]) + (*(]) << );
                v0 = _mm_unpacklo_epi32(_mm_cvtsi32_si128(i0), _mm_cvtsi32_si128(i1));
                i0 = *(]) + (*(]) << );
                i1 = *(]) + (*(]) << );
                v1 = _mm_unpacklo_epi32(_mm_cvtsi32_si128(i0), _mm_cvtsi32_si128(i1));
                v0 = _mm_unpacklo_epi8(v0, z);
                v1 = _mm_unpacklo_epi8(v1, z);

                a0 = _mm_unpacklo_epi32(_mm_loadl_epi64((__m128i*)(tabIJ + FXY[x] * )),
                    _mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )));
                a1 = _mm_unpacklo_epi32(_mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )),
                    _mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )));
                b0 = _mm_unpacklo_epi64(a0, a1);
                b1 = _mm_unpackhi_epi64(a0, a1);
                v0 = _mm_madd_epi16(v0, b0);
                v1 = _mm_madd_epi16(v1, b1);
                v0 = _mm_add_epi32(_mm_add_epi32(v0, v1), delta);

                i0 = *(]) + (*(]) << );
                i1 = *(]) + (*(]) << );
                v2 = _mm_unpacklo_epi32(_mm_cvtsi32_si128(i0), _mm_cvtsi32_si128(i1));
                i0 = *(]) + (*(]) << );
                i1 = *(]) + (*(]) << );
                v3 = _mm_unpacklo_epi32(_mm_cvtsi32_si128(i0), _mm_cvtsi32_si128(i1));
                v2 = _mm_unpacklo_epi8(v2, z);
                v3 = _mm_unpacklo_epi8(v3, z);

                a0 = _mm_unpacklo_epi32(_mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )),
                    _mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )));
                a1 = _mm_unpacklo_epi32(_mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )),
                    _mm_loadl_epi64((__m128i*)(tabIJ + FXY[x + ] * )));
                b0 = _mm_unpacklo_epi64(a0, a1);
                b1 = _mm_unpackhi_epi64(a0, a1);
                v2 = _mm_madd_epi16(v2, b0);
                v3 = _mm_madd_epi16(v3, b1);
                v2 = _mm_add_epi32(_mm_add_epi32(v2, v3), delta);

                v0 = _mm_srai_epi32(v0, INTER_REMAP_COEF_BITS);
                v2 = _mm_srai_epi32(v2, INTER_REMAP_COEF_BITS);
                v0 = _mm_packus_epi16(_mm_packs_epi32(v0, v2), z);
                _mm_storel_epi64((__m128i*)(D + x), v0);
            }
        }
        return x;
    }
};

void remapBilinear1(Mat_<uchar> src, Mat_<uchar> dst, Mat& mapxy, Mat& mapa, short* tabIJ)
{
    uchar* srcdata = (uchar*)src.data;
    uchar* dstdata = (uchar*)dst.data;
    vec2s* mapxydata = (vec2s*)mapxy.data;
    ushort* mapadata = (ushort*)mapa.data;
    , cols_ = cols - ;
    int i, j, io, jo, jprev, jprev_, joprev, joprev_, k;

    RemapVec_8u vecOp;

    , io = ; i < rows; ++i, io += cols)
    {
        uchar* D = dst.ptr<uchar>(i);

        bool prevInlier = false;

        jprev_ = ; joprev_ = io;
        , jo = io; j <= cols; ++j, ++jo)
        {
            ] < cols_ && (] < rows_ : !prevInlier;

            if (curInlier == prevInlier) continue;
            prevInlier = curInlier;

            jprev = jprev_; joprev = joprev_;
            jprev_ = j; joprev_ = jo;

            if (curInlier == false)
            {
                int len = vecOp(src, D, (short*)(mapxydata + joprev), mapadata + joprev, tabIJ, j - jprev);
                D += len;
                jprev += len;
                joprev += len;

                for (; jprev < j; jprev++, D++, ++joprev)
                {
                    ];
                    ];
                    ;
                    uchar* S = srcdata + sy*cols + sx;
                    ] * w[] + S[] * w[] + S[cols] * w[] + S[cols + ] * w[];
                    *D = saturate_cast_uchar((val + DELTA) >> SHIFT);
                }
            }
            else
            {
                for (; jprev < j; ++jprev, ++D, ++joprev)
                {
                    ];
                    ];
                     <  || sy >= rows || sy +  < ))
                    {
                        D[] = ;
                    }
                    else
                    {
                        int sx0, sx1, sy0, sy1;
                        uchar v0, v1, v2, v3;
                        ;

                        sx0 = borderInterpolate(sx, cols, BORDER_CONSTANT);
                        sx1 = borderInterpolate(sx + , cols, BORDER_CONSTANT);
                        sy0 = borderInterpolate(sy, rows, BORDER_CONSTANT);
                        sy1 = borderInterpolate(sy + , rows, BORDER_CONSTANT);
                        v0 = sx0 >=  && sy0 >=  ? srcdata[sy0*cols + sx0] : ;
                        v1 = sx1 >=  && sy0 >=  ? srcdata[sy0*cols + sx1] : ;
                        v2 = sx0 >=  && sy1 >=  ? srcdata[sy1*cols + sx0] : ;
                        v3 = sx1 >=  && sy1 >=  ? srcdata[sy1*cols + sx1] : ;

                        ] + v1*w[] + v2*w[] + v3*w[];
                        D[] = saturate_cast_uchar((val + DELTA) >> SHIFT);
                    }
                }
            }
        }
    }
}

void auremap(Mat_<uchar> src, Mat_<uchar> dst, Mat_<float> mapx, Mat_<float> mapy, int interpolation, int borderType = BORDER_CONSTANT, const Scalar& borderValue = Scalar())
{
    , cols_ = cols - ;
    , tabsz_= tabsz - , tabsz2 = ;//tabsz=INTER_TAB_SIZE, tabsz2=INTER_TAB_SIZE2
    .f / tabsz, vy, v, tmp2;
    ][][], *tabIJ = tab2D[][];
    ];//=8*tabsz

    //1.初化始化双线性插值表(在实际程序中应放到主循环外)
#if 1
    , io = -; i < tabsz; ++i)
    {
        tmp2 = i * scale;
        tab1D[++io] = .f - tmp2;    //tab0=1, tab2=31/32, tab4=30/32, ..., tab62=01/32
        tab1D[++io] = tmp2;            //tab1=0, tab3=01/32, tab5=02/32, ..., tab63=31/32
    }
    , io = ; i < tabsz; ++i, io += )
        , jo = ; j < tabsz; ++j, jo += , tabIJ += )
        {
            sum = ;
            , iio = ; ii < ; ++ii, iio += )    //tabi分别于tab0, tab1,...,tab63相乘后存放于tab2D, 存储顺序如(省略前缀tab1D和tab2D)
            {
                vy = tab1D[io + ii];                        //0: 0*0=0, 0*1=1; 0*2=5, 0*3=6; ...; 0*62=124, 0*63=125;
                , jjo = iio; jj < ; ++jj, ++jjo)//1: 1*0=2, 1*1=3; 1*2=7, 1*3=8; ...; 1*62=126, 1*63=127;
                {                                            //2: 2*0=128, 2*1=129; 2*2=132, 2*3=133; ...; 2*62=252, 2*63=253;
                    v = vy*tab1D[jo + jj];                    //3: 3*0=130, 3*1=131; 3*2=134, 3*3=135; ...; 3*62=254, 3*63=255;
                    v = v*INTER_REMAP_COEF_SCALE;
                    tmp1 = __round(v);
                    sum += tabIJ[jjo] = saturate_cast_short(tmp1);
                    ;
                }
            }

            if (sum != INTER_REMAP_COEF_SCALE)
            {
                diff = sum - INTER_REMAP_COEF_SCALE;
                mm = ; //int m1 = 1, m2 = 1;
                MM = ; //int M1 = 1, M2 = 1;
                , iio = ; ii < ; ++ii, iio +=)
                    , jjo = iio + jj; jj < ; ++jj)//寻找最大最小值
                        tabIJ[jjo] < tabIJ[mm] ? mm = jjo : (tabIJ[jjo] > tabIJ[MM] ? MM = jjo : ); //if (tabIJ[jjo] < tabIJ[m1*2 + m2]) {m1 = ii; m2 = jj;} else if (tabIJ[jjo] > tabIJ[M1*2 + M2]) {M1 = ii; M2 = jj; }

                diff <  ? tabIJ[MM] = (short)(tabIJ[MM] - diff): tabIJ[mm] = (short)(tabIJ[mm] - diff); //if (diff < 0) tabIJ[M1 * 2 + M2] = (short)(tabIJ[M1 * 2 + M2] - diff); else tabIJ[m1 * 2 + m2] = (short)(tabIJ[m1 * 2 + m2] - diff);
            }
        }
#endif

    //2.
    tabIJ = tab2D[][]; float *mapxdata = (float*)mapx.data, *mapydata = (float*)mapy.data;
    Mat_<Vec2s> mapxy(rows, cols); vec2s* mapxydata = (vec2s*)mapxy.data;
    Mat_<ushort> mapa(rows, cols); ushort* mapadata = (ushort*)mapa.data;
    , io = ; i < rows; ++i, io += cols)
    {
        , jo = io; j < cols; ++j, ++jo)
        {
            int sx = cvRound(mapxdata[jo] * tabsz);
            int sy = cvRound(mapydata[jo] * tabsz);
            int v = (sy & tabsz_)*tabsz + (sx & tabsz_);

            mapxydata[jo].] = saturate_cast_short(sx >> INTER_BITS);
            mapxydata[jo].] = saturate_cast_short(sy >> INTER_BITS);
            mapadata[jo] = (ushort)v;
        }
        remapBilinear1(src, dst, mapxy, mapa, tabIJ);
    }
}

#endif

图像插值:OpenCV_remap的更多相关文章

  1. OpenCV中的图像插值示例

    本文地址:http://www.cnblogs.com/QingHuan/p/7384433.html,转载请注明出处 ======================================== ...

  2. 图像上采样(图像插值)增取样(Upsampling)或内插(Interpolating)下采样(降采样),

    缩小图像(或称为下采样(subsampled)或降采样(downsampled))的主要目的有两个:1.使得图像符合显示区域的大小:2.生成对应图像的缩略图.放大图像(或称为上采样(upsamplin ...

  3. 数学之路-python计算实战(9)-机器视觉-图像插值仿射

    插值 Python: cv2.resize(src, dsize[, dst[, fx[, fy[, interpolation]]]]) → dst interpolation – interpol ...

  4. matlab彩色图像插值

    软件是MATLAB R2014b 使用的是matlab中已有的函数imresize(); 彩色图像分别对R.G.B三个通道进行插值,之后融合.[imresize函数不用,可以直接插值] clear;c ...

  5. python图像插值

    最近邻:选择离它所映射到的位置最近的输入像素的灰度值为插值结果. 最临近插值 图像的缩放很好理解,就是图像的放大和缩小.传统的绘画工具中,有一种叫做“放大尺”的绘画工具,画家常用它来放大图画.当然,在 ...

  6. Atitit 图像处理和计算机视觉的分类 三部分 图像处理 图像分析 计算机视觉

    Atitit 图像处理和计算机视觉的分类 三部分 图像处理 图像分析 计算机视觉 1.1. 按照当前流行的分类方法,可以分为以下三部分:三部分 图像处理 图像分析 计算机视觉1 1.2. 图像处理需要 ...

  7. dcraw源码解析

    dcraw源码解析 Author:Maddock Date:2015-04-22 转载请注明出处: 首先吐槽一点: 程序中使用了相当多的全局变量, 看的人头大.全局变量的坏处参看 http://wen ...

  8. ov5640摄像头设备驱动

    http://www.cnblogs.com/firege/p/5806121.html  (驱动大神) http://blog.csdn.net/yanbixing123/article/detai ...

  9. OpenCV 3.0函数库索引

    ================================== 前言:值的个数: cvCrossProduct:计算两个三维向量的向量积(叉积): cvCvtColor:将数组的通道从一个颜色空 ...

随机推荐

  1. 浏览器根对象window之screen

    1. screen 1.1 availHeight/Width screen.availWidth返回浏览器窗口可占用的水平宽度(单位:像素). screen.availHeight返回浏览器窗口在屏 ...

  2. 如何用Fritzing实现元器件自定义接线图

    在用Micropython开发板完成小实验时,很多朋友反应对照接线图实际接线有时会有一些困扰.今天给大家介绍一款画图软件Fritzing   看看是怎么自定义制作接线图的. 前提条件 1.准备好元器件 ...

  3. [笔记]《JavaScript高级程序设计》- 最佳实践

    一.可维护性 1 什么是可维护的代码 可理解性--其他人可以接受代码并理解它的意图和一般途径,而无需原开发人员的完整解释. 直观性--代码中的东西一看就能明白,不管其操作过程多么复杂. 可适应性--代 ...

  4. SAP HANA中的SLT简介

    在以SAP系统作为主要ERP的企业中,不同系统之间的数据库数据同步是个重要的工作.对于这种需求,除了开发ABAP接口之外,也有高效的工具可用.SLT就是其中之一. SLT是SAP的第一个ETL(Ext ...

  5. java web 之 listen 与 filter

    一.Listener监听器 Javaweb开发中的监听器,是用于监听web常见对象 HttpServletRequest HttpSession ServletContext 监听它们的创建与销毁.属 ...

  6. iOS 开发 右滑返回上一级控制器

    #import <objc/runtime.h> @interface UINavigationController (Transition)<UIGestureRecognizer ...

  7. debian9 配置数据库

    //安装数据库, 在debian上安装mysql,会默认安装mariadb sudo apt install mariadb-server // 这时普通用户无法登录,我们切换到root用户下 sud ...

  8. 使用js对form表单base64加密

    利用js可以对前台数据加密,以防止被恶意获取,以下代码,用base64对数据进行加密,可以在后台进行解密. 引入的js <script type="text/javascript&qu ...

  9. Mvc 模板化的Razor引擎委托

    最近在研究NopCommerce,它后台用的富文本编辑器可根据语言库加载不同语言的编辑器,其中用到了模板化Razor引擎委托,参考这儿 废话不多说,直接上代码. public static class ...

  10. CTSC2017 && APIO2017 && THUSC2017 游记

    一去北京就是近20天,还是回来写写游记吧. 5.6 坐飞机到天津转动车到北京. 5.7 在天坛公园逛了一圈就去报到了. 下午试机,好像没发生什么. 5.8 CTSC一试 T1签到,开个桶打个标记就好了 ...