Inter IPP 绘图 ippi/ipps

IPP的资料网上比较少，主要还是参考Inter官网和文档

官方文档ipps.pdf主要是对数据做处理，包括加减乘除、FFT、DFT等

文档ippi.pdf只要是对图像做处理，包括通道转换、图片处理等

IPP最新的几个版本删掉了生成图片的函数，所以使用IPP加速计算，然后再用opencv生成本地图片

1、搭建VS+IPP+Opencv环境（前一篇文章）

2、IPP的数据Ipp32f 转换成 Opencv的Mat数据

#define READ_BYTES 32768
#define RC 8192
Ipp32f *specturmData = ippsMalloc_32f(READ_BYTES);//创建Ipp32f
//添加给specturmData 赋值之类的操作
//........
//........
int num=0;
Mat mat;
mat.create(READ_BYTES, RC, CV_32F);//创建mat数组,设定数组高度为32768，长度为8192，也就是画出一幅图32768*8192
while(num<RC)//把specturmData 赋给mat
{
    for (int i = 0; i < READ_BYTES ; i++)
    {
    mat.at<float>(i, num) = specturmData[i];//逐个赋值
    if (specturmData[i] > max)
        max = specturmData[i];
    if (specturmData[i] < min)
        min = specturmData[i];
    }
    num++;
}

3、将mat存成本地图片

//添加对mat的数据处理
//......
//......
imwrite("my.bmp", mat);//画图

4、用IPP 做FFT变换

新版的IPP删除了ippsFFTInitAlloc_R_32f等函数，如图文档附录给出了替代方案

int   CalculateSpectrumByReal(short* rData, int samplePoint, float* specturmData, int &outSpecturmDataLength,
    int winOpt, int fftOrder)/*计算频谱*/
{
    IppStatus status;
    int segement, fftPoint;
    IppsFFTSpec_R_32f* spec;

    fftPoint = (int)pow(2.0, fftOrder);        //FFT 点数
    outSpecturmDataLength = fftPoint / 2;       //输出数据长度 = FFT点数的一半

    Ipp32f *fData = ippsMalloc_32f(samplePoint);//设置浮点数组 fData
    Ipp32f *x = ippsMalloc_32f(fftPoint);//设置浮点数组 x
    Ipp32f *X = ippsMalloc_32f(fftPoint + 2);//设置浮点数组 X
    Ipp32f *mag = ippsMalloc_32f(outSpecturmDataLength);//求abs()后的幅度数据 一半的下标

    memset(fData, 0, sizeof(Ipp32f)*samplePoint); //初始化
    memset(x, 0, sizeof(Ipp32f)*fftPoint);//初始化
    memset(X, 0, sizeof(Ipp32f)*(fftPoint + 2));//初始化
    memset(mag, 0, sizeof(Ipp32f)*outSpecturmDataLength);//初始化
    memset(specturmData, 0, sizeof(Ipp32f)*outSpecturmDataLength);//初始化
    ippsConvert_16s32f(rData, fData, samplePoint);//转ipp结构

    int sizeSpec = 0;
    int sizeInit = 0;
    int sizeBuffer = 0;
　　 /*替代已经去掉的ippsFFTInitAlloc_R_32f函数*/
    status = ippsFFTGetSize_R_32f(fftOrder, IPP_FFT_DIV_INV_BY_N, ippAlgHintAccurate, &sizeSpec, &sizeInit, &sizeBuffer);
    Ipp8u *pMemSpec = ippsMalloc_8u(sizeSpec);
    Ipp8u *pMemInit = ippsMalloc_8u(sizeInit);
    Ipp8u *pMemBuffer = ippsMalloc_8u(sizeBuffer);
    status = ippsFFTInit_R_32f(&spec, fftOrder, IPP_FFT_DIV_FWD_BY_N, ippAlgHintAccurate, pMemSpec, pMemInit);
    segement = (int)(samplePoint / fftPoint);  //做FFT时的分段数 参数设定
    for (int i = 0; i<segement; i++)
    {
        memcpy(x, fData + i*fftPoint, fftPoint * sizeof(Ipp32f));
        status = ippsFFTFwd_RToCCS_32f(x, X, spec, pMemBuffer);  //do fft...
        ippsMagnitude_32fc((Ipp32fc*)X, mag, outSpecturmDataLength); //abs()...
        ippsAdd_32f_I(mag, specturmData, outSpecturmDataLength);
    }
    ippsDivC_32f_I(segement, specturmData, outSpecturmDataLength);//求平均
    ippsLn_32f_I(specturmData, outSpecturmDataLength);   //求Ln，因为IPP函数不支持浮点的10Log10计算这里通过Ln进行转换计算
    ippsDivC_32f_I(2.3026, specturmData, outSpecturmDataLength);// log(e,10)=2.3026
    ippsMulC_32f_I(20.0, specturmData, outSpecturmDataLength);  //20log10()

    ippsFree(spec);
    ippsFree(fData);
    ippsFree(x);
    ippsFree(X);
    ippsFree(mag);
    return status;
}

 https://www.cnblogs.com/dzzy/p/11242466.html