cv::remap (in opencv) and interp2 (matlab)
그러나, 이들 두 가지 기능의 출력에 상당한 차이가 있음을 알았다. 여기에 내가
NaN NaN NaN
NaN NaN NaN
NaN NaN 0.1000
로 출력 V 행렬을 얻을 내 MATLAB 코드
U = [0.1 0.1 0.1; 0.2 0.2 0.2; 0.3 0.3 0.3];
X = [0 0 0; 0.5 0.5 0.5; 1 1 1];
Y = [0 0.5 1;0 0.5 1;0 0.5 1];
V = interp2(U,X,Y,'linear',NaN)
이 내 OpenCV의 코드입니다
#include "highgui.h"
#include "cv.h"
using namespace cv;
using namespace std;
int main(int argc, char** argv)
{
//generate flowmap model
CvFileStorage* fx = cvOpenFileStorage("result.txt", 0, CV_STORAGE_WRITE);//ask storage for save file
Mat xmesh = cvCreateMat(3, 3, 5);
Mat ymesh = cvCreateMat(3, 3, 5);
for (int i = 0; i < xmesh.rows; i++)
for (int j = 0; j < xmesh.cols; j++)
{
xmesh.at<float>(i, j) = i*0.5;
ymesh.at<float>(i, j) = j*0.5;
}
//generate optical flow folder
Mat u = cvCreateMat(3, 3, 5);
for (int i = 0; i <u.rows; i++)
for (int j = 0; j < u.cols; j++)
{
u.at<float>(i, j) = (i + 1)*0.1;
}
Mat v = Mat::zeros(u.size(), u.type());
remap(u, v, xmesh, ymesh, INTER_LINEAR, 0, cvScalarAll(0));
//convert mat to Iplimage
IplImage* xmesh_a = cvCloneImage(&(IplImage)xmesh);
IplImage* ymesh_a = cvCloneImage(&(IplImage)ymesh);
IplImage* u_a = cvCloneImage(&(IplImage)u);
IplImage* v_a = cvCloneImage(&(IplImage)v);
//save end to txt
cvWrite(fx, "xmesh", xmesh_a, cvAttrList());
cvWrite(fx, "ymesh", ymesh_a, cvAttrList());
cvWrite(fx, "u", u_a, cvAttrList());
cvWrite(fx, "v", v_a, cvAttrList());
cvReleaseFileStorage(&fx);
waitKey();
}
내가 얻을 출력 V 행렬은
1.00000001e-001, 1.50000006e-001, 2.00000003e-001,
1.00000001e-001, 1.50000006e-001, 2.00000003e-001,
1.00000001e-001, 1.50000006e-001, 2.00000003e-001
입니다
도움이 될 것입니다. 감사! 결과의
하지 이런 식으로 함수 호출 :'interp2 (X, Y, U, ...)'? – Irreducible
@Irereducible 답장을 보내 주셔서 감사합니다.하지만 그리드 보간에 아무런 차이가 없다고 생각합니다. 이 포스터 읽기 https://stackoverflow.com/questions/32235379/how-to-do-grid-interpolation-interp2-in-opencv – SimaGuanxing