【正文】 平方和法 for (int i = 0。 i nr。 ++i) { uchar *subwin_ptr = uchar(i)。 uchar *subwind_ptr= uchar(i)。 for (int j = 0。 j nc。 j++) { sum += (float)(*(subwin_ptr+j)*(subwind_ptr+j))*(*(subwin_ptr+j)*(subwind_ptr+j))。 } } return sum。 case 5://差 絕對(duì)值和法 for (int i = 0。 i nr。 ++i) { uchar *subwin_ptr = uchar(i)。 uchar *subwind_ptr= uchar(i)。 for (int j = 0。 j nc。 j++) { sum += (float)abs(*(subwin_ptr+j) *(subwind_ptr+j))。 25 } } return sum。 default: break。 } return sum。 } 圖 2 相關(guān)函數(shù) 圖 3 協(xié)方差 圖 4 相關(guān)系數(shù) 圖 5 差平方和 圖 6 差絕對(duì)值 同名核線提取 include include include vector include string include fstream include iostream include QFileDialog using namespace std。 //構(gòu)造函數(shù) HXDialog::HXDialog(QWidget *parent): QDialog(parent), ui(new Ui::Dialog) { uisetupUi(this)。 } HXDialog::~HXDialog() { } void HXDialog::on_pushButton_clicked() { path1 = QFileDialog::getOpenFileName(this)。 26 uilineEditsetText(path1)。 } void HXDialog::on_pushButton_2_clicked() { path2= QFileDialog::getOpenFileName(this)。 uilineEdit_2setText(path2)。 } void HXDialog::on_pushButton_3_clicked() { tpath1= QFileDialog::getOpenFileName(this)。 uilineEdit_3setText(tpath1)。 } void HXDialog::on_pushButton_4_clicked() { tpath2= QFileDialog::getOpenFileName(this)。 uilineEdit_4setText(tpath2)。 } //讀取同名點(diǎn)數(shù)據(jù) int HXDialog::readxy(string p1,string p2) { int index = 1。 ifstream fil。 xy tmp。 int head。 //讀取圖像 1的像點(diǎn)坐標(biāo) (p1)。 if (!fil) { cout read file error endl。 return 1。 } fil head。 while (!()) { fil 。 cout read index : endl。 index++。 (tmp)。 } ()。 ()。 //讀取圖像 2的像點(diǎn)坐標(biāo) index = 1。 (p2)。 if (!fil) { cout read file error endl。 return 1。 } fil head。 while (!()) { fil 。 cout read index : endl。 index++。 (tmp)。 } cout read end endl。 ()。 cout close objfile... endl。 return 0。 } class findxy { public: findxy(int index){ ind = index。 } int ind。 27 bool operator () (HXDialog::xy amp。pt){ return ind == 。 } }。 //同名核線算法主函數(shù) void HXDialog::HX() { string p1 = ()。 string p2 = ()。 //讀取 同名點(diǎn) 坐標(biāo) readxy(p1,p2)。 //初始值 p1 = ()。 p2 = ()。 cv::Mat img1 = cv::imread(p1)。 cv::Mat img2 = cv::imread(p2)。 cv::Mat B。 cv::Mat L。 if (() || ()) { cout read image error endl。 return。 } for (int i = 0。 i ()。 ++i) { double x = pic1[i].x。 double y = pic1[i].y。 double x_ = pic2[i].x。 double y_ = pic2[i].y。 cv::Mat tmpB(cv::Size(8, 1), CV_64FC1)。 double(0) = 。 double(1) = x。 double(2) = y。 double(3) = x_。 double(4) = x_*x。 double(5) = x*y_。 double(6) = y*x_。 double(7) = y*y_。 cout double(0) double(1) double(2) double(3) double(4) double(5) double(6) double(7) endl。 (tmpB)。 (y y_)。 } cv::Mat X = (()*B).inv() * () *L。 double L1 = double(0)。 double L2 = double(1)。 double L3 = double(2)。 double L4 = double(3)。 double L5 = double(4)。 double L6 = double(5)。 double L7 = double(6)。 double L8 = double(7)。 vectorxy::iterator itr1 = find_if((), (), findxy(4507))。 vectorxy::iterator itr2 = find_if((), (), findxy(4507))。 double x = (*itr1).x。 double y = (*itr1).y。 double x0 = 0。 double x1 = 。 double x_0 = 0。 double x_1 = 。 double y_0 = ((1 L3)*y L1 L2*x L4*x_0 L5*x*x_0 L7*y*x_0) / (1 + L6 * x + L8*y)。 double y_1 = ((1 L3)*y L1 L2*x L4*x_1 L5*x*x_1 L7*y*x_1) / (1 + L6 * x + L8*y)。 double y0 = (y_0 + L1 + L2*x0 + L4*x_0 + L5*x1*x_0 + L6*x*y_0) / (1 L3 L7*x_0 L8*y_0)。 double y1 = (y_0 + L1 + L2*x1 + L4*x_0 + L5*x1*x_0 + L6*x*y_0) / (1 L3 L7*x_0 28 L8*y_0)。 cv::line(img1, cv::Point(x0, y0), cv::Point(x1, y1), CV_RGB(255,0,0), 5)。 cv::line(img2, cv::Point(x_0, y_0), cv::Point(x_1, y_1), CV_RGB(255, 0, 0), 5)。 cv::circle(img1, cv::Point(x, y), 10, CV_RGB(0, 255, 0), 5)。 cv::circle(img2, cv::Point((*itr2).x, (*itr2).y), 10, CV_RGB(0, 255, 0), 5)。 cv::resize(img1, dst1, cv::Size(1000, 1000 * / ))。 cv::resize(img2, dst2, cv::Size(1000, 1000 * / ))。 } void HXDialog::on_okButton_clicked() { HX()。 emit signelHX(dst1,dst2)。 } void HXDialog::on_cancelButton_clicked() { } 圖 7 左核線 圖 8 右核線 五、 心得體會(huì) 通過這次的實(shí)習(xí)，我學(xué)習(xí)了數(shù)字?jǐn)z影測(cè)量系統(tǒng) VirtuoZo和 GeoWay的操作方法，了解了通過這兩套系統(tǒng)生成 4D 產(chǎn)品的生產(chǎn)流程。在操作過程中，結(jié)合老師的講解，對(duì)于以前課堂上學(xué)過的有關(guān)知識(shí)也進(jìn)行了回顧。同時(shí)我體會(huì)到了 VirtuoZo 和GeoWay 這兩款數(shù)字?jǐn)z影測(cè)量系統(tǒng)的強(qiáng)大功能，操作的步驟也比較復(fù)雜，但是只要熟練掌握，就能夠快速生成產(chǎn)品。 這次的數(shù)字?jǐn)z影測(cè)量系統(tǒng)操作過程中也遇到了這樣那樣的問題，比如說內(nèi)定向失敗、相對(duì)定向失敗等，最后通過問老師同學(xué)都得到了解決。大部分的問題 歸結(jié)起來都是由于沒有嚴(yán)格按照流程來做，比如說漏掉哪個(gè)步驟或者哪個(gè)參數(shù)沒有設(shè)置等等。通過老師的講解，也學(xué)到了很多在實(shí)際生產(chǎn)中有用的方法。 通過這次的編程實(shí)習(xí)，我回顧了學(xué)過的相關(guān)攝影測(cè)量原理，并把他們通過程序?qū)崿F(xiàn)了。在這過程中提高了自己的編程能力以及加深對(duì)于相關(guān)原理的理解。同時(shí)也認(rèn)識(shí)到了一些實(shí)際的問題。比如說在做匹配的時(shí)候，五種方法得到的結(jié)果區(qū)別十分大，相關(guān)函數(shù)法基本不能夠成功匹配。 29 這次實(shí)習(xí)我感受到，即使是在計(jì)算機(jī)技術(shù)如此發(fā)達(dá)的時(shí)代，即使有如此完整的數(shù)字?jǐn)z影測(cè)量系統(tǒng)，但是整個(gè)生產(chǎn)過程也離不開人的干預(yù)。而且人的 干預(yù)直接影響了得到的結(jié)果。所以我們學(xué)習(xí)數(shù)字?jǐn)z影測(cè)量具有十分重要的意義，一方面我們只有理解了相關(guān)的原理，才能寫出相應(yīng)的算法，以推動(dòng)數(shù)字?jǐn)z影測(cè)量系統(tǒng)的發(fā)展，另一方面，我們只有了解了相關(guān)知識(shí)，才能在操作中準(zhǔn)確的進(jìn)行干預(yù)，得到最好的結(jié)果。