Opencv 中的 PCA 项目和后台项目

PCA project and backproject in Opencv

本文关键字：项目 后台 PCA 中的 Opencv      更新时间：2023-10-16

我在Ubuntu Opencv中工作，我正在尝试对单个图像进行PCA分析。我采用 3 通道图像并将其更改为具有 3 列和 r*c 行数的单通道图像.r 和 c 是原始图像的行和列。当我在PCA上进行反向投影后尝试显示重建的图像时，它会给我一个绿色图像。这是我的代码

Mat pcaset=cvCreateMat(image->height*image->width,image->nChannels,CV_8UC1);
for(int i=0;i<image->height;i++)
    {
        for(int j=0;j<image->width;j++)
        {
            for(int k=0;k<image->nChannels;k++)
            (ptrpcaset+i*pcaset.step)[k]=((ptrimage+i*image->widthStep)[3*j+k]);
       }
    }
int nEigens=3;
    Mat databackprojected;
    PCA pca(pcaset,Mat(),CV_PCA_DATA_AS_ROW,nEigens);
    Mat dataprojected(pcaset.rows,nEigens,CV_8UC1);
    pca.project(pcaset,dataprojected);
    pca.backProject(dataprojected,databackprojected);
    Mat backprojectnorm;//(databackprojected.rows,nEigens,CV_8UC1);
    normalize(databackprojected,backprojectnorm,0,255,NORM_MINMAX,-1);
    Mat finaldataafterreshaping(image->height,image->width,CV_8UC3);
    uchar* finalptr=(uchar*)finaldataafterreshaping.data;
    uchar* ptrnorm=(uchar*)backprojectnorm.data;
    int x=0,y=0,i=0;
    while(i<backprojectnorm.rows)
    {
        while(x<image->height)
        {
            while(y<image->width)
            {
                for(int k=0;k<image->nChannels;k++)
                {
                    (finalptr+x*finaldataafterreshaping.step)[3*y+k]=(ptrnorm+i*backprojectnorm.step)[k];
                }
                y=y+1;i=i+1;
            }
            x=x+1;y=0;
        }
    }
imshow("Reconstructed data",finaldataafterreshaping);

您需要进行以下更改：

(ptrpcaset+(j + i*image->width)*pcaset.step)[k]=((ptrimage+i*image->widthStep)[3*j+k]);

因为在转换数据时没有考虑 j 坐标，因此最后只将图像的最后一行保存在新矩阵中。

重塑数据时，需要执行以下操作： 

float* val = (float*)&(ptrnorm+i*backprojectnorm.step)[(k*4)];
(finalptr+x*finaldataafterreshaping.step)[3*y+k]=*val;

因为您作为结果获得的矩阵是 float 类型而不是 uchar .所以你需要某种转换。我不确定，这样做是否是个好主意，但它有效。我建议你看看OpenCV 2的C++ API，它可以以更好的方式处理这件事。

此外，不需要整个while(i<backprojectnrom.rows)循环。

相关文章：