cv::Matの計算について

Question

前提・実現したいこと

uv[a]がこのプログラム内で何を表しているのか，どなたか詳しい方お願いします。

該当のソースコード

//Point(x,y)という意味を持たせる
struct myclass {
    bool operator() (cv::Point pt1, cv::Point pt2) { return (pt1.x < pt2.x);}
} myobject;

  

// AX = B を解く. 両辺にAの逆行列を前掛けする
	matX = matA.inv()*matB;
	MatConstIterator_<float> iter = matX.begin<float>();
	while(iter != matX.end<float>()) printf("%.1f\n", *iter++);

	//切り上げで長方形座標をそれぞれ求める(赤点)

	std::vector<cv::Point> uv(5);
	while(a < 5){
		x = pts[a].x;
		y = pts[a].y;
		mother = matX.at<float>(6,0)*x + matX.at<float>(7,0)*y + 1;
		uv[a] =Point(ceil((matX.at<float>(0,0)*x + matX.at<float>(1,0)*y + matX.at<float>(2,0))/mother),ceil((matX.at<float>(3,0)*x + matX.at<float>(4,0)*y + matX.at<float>(5,0))/mother));
		a++;
	}

補足情報（FW/ツールのバージョンなど）

visual studio2012

Comments

[fana]

同次座標系を用いた何らかの座標変換処理っぽく見えますが，断片すぎて意味不明なので，これを書いた本人に訊いてください．

[tiitoi]

まずそのコードはどこから持ってきたものなのか (どこかのサイトならそのURL) を追記してください。

Answer 1

コードの雰囲気からの単なる想像です．

(x,y)を(u,v)に変換する何らかの式を，下記のような同次座標で表していて，

[u]   [ a b c ][x]
[v] = [ d e f ][y]
[1]   [ g h 1 ][1]

前半のAX = Bがどうのこうのいう箇所では，この変換のパラメタベクトルX=[ a,b,c,d,e,f,g,h ]^tを，
既知の何点かの(x,y)←→(u,v)対応データを用いて最小二乗法か何かで求めていて，

最後のuv[a]というのは，求めた変換パラメタで(x,y)座標系側の何らかのデータ点を(u,v)の側に変換した結果でしょう．

Comments

[fana]

これでBAか…
まぁ何かが解決するヒントになったのであれば良いのだけど．