不扭曲失真参数

我试图用OpenCV(通过OpenCVSharp的包装器)来消除扭曲的图像，我很难找出矩阵的正确格式，我正在通过CvFileStorage改进它们，如在未扭曲的样本中使用，但是我得到了真正意想不到的结果。

从我使用的另一个软件，我得到以下数据(这应该是正确的，因为该软件不扭曲图像正好与这些参数)

    <sensor id="0" label="RP_OV5647 (3.5976 mm)" type="frame">
      <resolution width="2592" height="1944"/>
      <property name="focal_length" value="3.5975999999999999e+000"/>
      <property name="fixed" value="false"/>
      <calibration type="frame" class="adjusted">
        <resolution width="2592" height="1944"/>
        <fx>2.5667779156627239e+003</fx>
        <fy>2.5654169232963172e+003</fy>
        <cx>1.2933905910821311e+003</cx>
        <cy>9.6694162774618212e+002</cy>
        <skew>2.9824815660771562e+000</skew>
        <k1>9.1041365324307497e-002</k1>
        <k2>-4.0485507081497402e-001</k2>
        <k3>3.3943759073230600e-001</k3>
        <p1>3.4409596859645629e-004</p1>
        <p2>-3.9472652058529605e-005</p2>
      </calibration>
    </sensor>

我使用以下代码来消除扭曲图像，固有矩阵很好(在应用不失真时给出预期结果)，但是失真矩阵给了我问题，从我对预期布局的理解(没有找到任何样本，但这是我从OpenCV源收集的)一个4X1矩阵k1k2p1p2应该工作，所以我在我的配置文件中有这个:

<?xml version="1.0"?>
<opencv_storage>
   <intrinsic type_id="opencv-matrix">
   <rows>3</rows>
   <cols>3</cols>
   <dt>f</dt>
   <data>
      1293.3905910821311 0. 966.94162774618212 0. 2566.7779156627239 2565.4169232963172 0. 0. 1.
   </data>
 </intrinsic>
 <distortion type_id="opencv-matrix">
   <rows>1</rows>
   <cols>5</cols>
   <dt>f</dt>
   <data>
  0.091041365324307497 -0.40485507081497402 0.00034409596859645629 -0.000039472652058529605 0.33943759073230600
   </data>
 </distortion>

代码:

            using (IplImage srcImg = new IplImage(img.Path, LoadMode.Color))
            using (IplImage dstImg = srcImg.Clone())
            {
                CvMat intrinsic, distortion;
                using (CvFileStorage fs = new CvFileStorage("camera.xml", null, FileStorageMode.Read))
                {
                    CvFileNode param = fs.GetFileNodeByName(null, "intrinsic");
                    intrinsic = fs.Read<CvMat>(param);
                    param = fs.GetFileNodeByName(null, "distortion");
                    distortion = fs.Read<CvMat>(param);
                }
                Cv.Undistort2(srcImg, dstImg, intrinsic, distortion);
                dstImg.SaveImage(new System.IO.FileInfo(img.Path).Name);
            }

然而，我得到完全扭曲的图像(漩涡)，我有矩阵顺序错误吗?是否有任何代码样本(我找不到任何)从命名参数构造的矩阵(所有我发现的是使用一个样本文件的样本，其中有一个矩阵，没有关于组件订单的文档)?

int main()
{
    // camera resolution
    cv::Mat input= cv::imread("../inputData/distorted1.jpg");

    cv::Mat distCoeff;
    distCoeff = cv::Mat::zeros(8,1,CV_64FC1);
    // indices: k1, k2, p1, p2, k3, k4, k5, k6 
    // your coefficients here!
    double k1 = 9.1041365324307497e-002;
    double k2 = -4.0485507081497402e-001;
    double p1 = 3.4409596859645629e-004;
    double p2 = -3.9472652058529605e-005;
    double k3 = 3.3943759073230600e-001;
    double k4 = 0;
    double k5 = 0;
    double k6 = 0;
    distCoeff.at<double>(0,0) = k1;
    distCoeff.at<double>(1,0) = k2;
    distCoeff.at<double>(2,0) = p1;
    distCoeff.at<double>(3,0) = p2;
    distCoeff.at<double>(4,0) = k3;
    distCoeff.at<double>(5,0) = k4;
    distCoeff.at<double>(6,0) = k5;
    distCoeff.at<double>(7,0) = k6;

    // camera intrinsics
    cv::Mat cam1;
    cam1 = cv::Mat::eye(3,3,CV_32FC1);
    cam1.at<float>(0,0) = 2.5667779156627239e+003;
    //cam1.at<float>(0,1) = 0;
    cam1.at<float>(0,1) = 2.9824815660771562e+000; // skew?
    cam1.at<float>(0,2) = 1.2933905910821311e+003;
    cam1.at<float>(1,0) = 0;
    cam1.at<float>(1,1) = 2.5654169232963172e+003;
    cam1.at<float>(1,2) = 9.6694162774618212e+002;
    // direct undistort:
    cv::Mat output;
    cv::undistort(input, output, cam1,distCoeff);

    cv::imshow("input", input);
    cv::imshow("distCorr-direct", output);
    cv::waitKey(-1);
    return 0;
}