전체 프레임 물고기 눈을 수학적으로 등각 투구로 변환

Question

전체 프레임 어안 이미지를 정사각형으로 변환하고 파노라마 스티칭에 사용할 수 있습니다. 나는이 목적을 위해 Hugin과 libpano를 사용했으며 잘 동작한다. 다음 이미지는 hugin에 의해 변환되었습니다. 나는 그러한 또는 유사한 변형이 성취 될 수있는 수학 방정식을 알고 싶다.

는 오리지널 이미지 사용

변환 화상

Answer 1

I 성공적으로 생성 한 hugin에 가까운 화상으로 변환 한 기준이 link. 더 나은 이미지 품질을 위해 여전히 서브 픽셀 샘플링이 필요합니다.

#include <iostream> 
#include <sstream> 
#include <time.h> 
#include <stdio.h> 
#include <opencv2/core/core.hpp> 
#include <opencv2/imgproc/imgproc.hpp> 
#include <opencv2/calib3d/calib3d.hpp> 
#include <opencv2/highgui/highgui.hpp> 
using namespace std; 
using namespace cv; 
#define PI 3.1415926536 
Point2f getInputPoint(int x, int y,int srcwidth, int srcheight) 
{ 
    Point2f pfish; 
    float theta,phi,r; 
    Point3f psph; 
    //float FOV = PI; 
    float FOV =(float)PI/180 * 150; 
    float FOV2 = (float)PI/180 * 120;; 
    float width = srcwidth; 
    float height = srcheight; 
    theta = PI * (x/width - 0.5); // -pi to pi 
    phi = PI * (y/height - 0.5); // -pi/2 to pi/2 
    psph.x = cos(phi) * sin(theta); 
    psph.y = cos(phi) * cos(theta); 
    psph.z = sin(phi); 
    theta = atan2(psph.z,psph.x); 
    phi = atan2(sqrt(psph.x*psph.x+psph.z*psph.z),psph.y); 
    r = width * phi/FOV; 
    float r2 = height * phi/FOV2; 
    pfish.x = 0.5 * width + r * cos(theta); 
    pfish.y = 0.5 * height + r2 * sin(theta); 
    return pfish; 
} 
int main(int argc, char **argv) 
{ 
    if(argc< 3) 
     return 0; 
    Mat orignalImage = imread(argv[1]); 
    if(orignalImage.empty()) 
    { 
     cout<<"Empty image\n"; 
     return 0; 
    } 
    Mat outImage(orignalImage.rows,orignalImage.cols,CV_8UC3); 
     //getInputPoint(0,5,10,10); 
    namedWindow("result",CV_WINDOW_NORMAL); 
    for(int i=0; i<outImage.cols; i++) 
    { 
     for(int j=0; j<outImage.rows; j++) 
     { 

      Point2f inP = getInputPoint(i,j,orignalImage.cols,orignalImage.rows); 
      //cout<<"in "<<i<<","<<j<<endl; 
      //cout<<"out "<<inP.x<<","<<inP.y<<endl; 
      Point inP2((int)inP.x,(int)inP.y); 
      if(inP2.x >= orignalImage.cols || inP2.y >= orignalImage.rows) 
       continue; 
      if(inP2.x < 0 || inP2.y < 0) 
       continue; 
      Vec3b color = orignalImage.at<Vec3b>(inP2); 
      outImage.at<Vec3b>(Point(i,j)) = color; 

     } 
    } 

    imshow("result",outImage); 
    imwrite(argv[2],outImage); 


}