【正文】 圖像上選點，保證每個藍框周圍至少有 3 個點就行。 7 圖 20 控制點 圖 21 被挑出的點 可能一次選擇的點精度不夠會有誤差超限的點，顯示在 如圖 21 中帶有星號，然后對帶有星號的點從新定位，這樣直到?jīng)]有星號的點。 3) 三角測量 打開 VirtuozoNT，打開在 AAT 中建立的測區(qū)，自動創(chuàng)建測區(qū)，之就可以使用批處理功能進行自動化處理。 圖 22 批處理界面 4) DEM 拼接和正射影像制作 按照前面的拼接步驟進行 DEM 拼接和正射影像的制作，結(jié)果如下： 圖 23 DEM 圖 24 正射影像 8 三、 GeoWay 實習流程 交互式編輯 1） 新建工 程，設(shè)置相關(guān)參數(shù) （ PS：對照表添加 [LAYER] [ENDSEC]，然后將 SYM 符號庫文件拷貝到 C:\WINDOWS\GeoSymbol） 圖 1 參數(shù)設(shè)置界面 2） 生成內(nèi)圖廓線 圖 2 生成圖廓線 9 3）交互式 地物編輯 加載柵格圖層。 點、單線、雙線編輯 點、線處理 。 批處理 — 消除懸掛點 和偽節(jié)點能夠消除大部分錯誤，不能消除的部分進行手工編輯。 新建房屋對照表，在拓撲構(gòu)建時啟用房屋對照表，然后 10 拓撲弧段預處理，最后進行拓撲構(gòu)建。 然后 自動檢查錯誤并查看錯誤，然后對矢量圖進行修改 16 圖 22 編輯前 圖 23 編輯后 圖 24 編輯前 圖 25 編輯后 17 四、編程內(nèi)容 單像空間后方交會 include include include vector include string include fstream include iostream include QFileDialog using namespace std。 char *title= 空間后方交會 。 setWindowTitle(qtitle)。 //焦距 單位 pixel fx=。 //內(nèi)方位元素 pixel x0=。 //畸變參數(shù) k1=。 p1=。 //初始值 Xs0=。 Zs0=。 omega0=。 //單位 rad } HFJHDialog::~HFJHDialog() { delete ui。 //path =F:\\我的課程 \\數(shù)字攝影測量影像數(shù)據(jù) \\后方交會 \\像點坐標 .txt。 (path)。 uitextBrowserappend(QString::fromLocal8Bit(log))。 } xy tmp。 uitextBrowserappend(QString(Point%1: x:%2 y:%3).arg(index).arg().arg())。 (tmp)。 return 0。 18 ifstream fil。 if(!fil) { char *log = 打開控制點文件失敗 。 return 1。 string head。 i 4。 } while(!()) { fil。 index++。 } ()。 } //讀取相機參數(shù)文件 int HFJHDialog::readcamera(string path) { int index = 1。 ifstream fil。 if(!fil) { char *log = 打開相機文件失敗 。 return 1。 return 0。 //定義 double phi=phi0, omega=omega0, kappa=kappa0。 //定義 X,Y,Z double a1,a2,a3,b1,b2,b3,c1,c2,c3。 double a11,a12,a13,a14,a15,a16。 double lx,ly。 cv::Mat A, L, dX。 //迭代終止閾值 double e = / 206265。 int index = 1。 double H = 。 y0 = H / 2 y0。 i != ()。 y = pic[i].y。 Y = 0。 for (size_t j = 0。 j++) { if (obj[j].index == pic[i].index) { X = obj[j].X。 Z = obj[j].Z。 } } if (abs(x) 1) { char *arg1 = 像點 。 uitextBrowserappend(QString(%1%2%3).arg(QString::fromLocal8Bit(arg1)).arg(pic[i].index).arg(QString::fromLocal8Bit(arg2)))。 } //求觀測值的改正 r_r = (x x0)*(x x0) + (y y0)*(y y0)。 dy = (y y0)*(k1*r_r + k2*r_r*r_r) + p2*(r_r + 2 * (y y0)*(y y0)) 2 * p1*(x x0)*(y y0)。 y = (y y0) dy。 a2 = cos(phi)*sin(kappa) sin(phi)*sin(omega)*cos(kappa)。 b1 = cos(omega)*sin(kappa)。 b3 = sin(omega)。 c2 = sin(phi)*sin(kappa) + cos(phi)*sin(omega)*cos(kappa)。 //求 X_， Y_， Z_ X_ = a1*(X Xs) + b1*(Y Ys) + c1*(Z Zs)。 Z_ = a3*(X Xs) + b3*(Y Ys) + c3*(Z Zs)。 double y_ = fy*Y_ / Z_。 ly = y y_。 a12 = 1 / Z_*(b1*fx + b3*x)。 a21 = 1 / Z_*(a2*fy + a3*y)。 a23 = 1 / Z_*(c2*fy + c3*y)。 a15 = fx*sin(kappa) x / fx*(x*sin(kappa) + y*cos(kappa))。 a24 = x*sin(omega) (y / fy*(x*cos(kappa) y*sin(kappa)) fy*sin(kappa))*cos(omega)。 a26 = x。 double l[2] = { lx, ly }。 cv::Mat Li(2, 1, CV_64FC1, l)。 (Li)。 Xs+=double(0,0)。 Zs+=double(2,0)。 omega+=double(4,0)。 uitextBrowserappend(QString(%1: %2 %3 %4 %5 %6 %7).arg(index).arg(tdouble(0)).arg(double(1)).arg(double(2)).arg(double(3)).arg(double(4)).arg(double(5)))。 cv::Mat M0 = (A*dX L).t()*(A*dX L)。 ()。 if (index 100) break。 uitextBrowserappend(QString(Result:))。 char *log = 單位權(quán)中誤差： 。 } void HFJHDialog::on_pushButton_clicked() { QString qpath = QFileDialog::getOpenFileName(this)。 if(readxy(path) == 1) { char *log =讀取像點坐標錯誤 XXXXXX。 return。 uitextBrowserappend(QString::fromLocal8Bit(log))。 } void HFJHDialog::on_pushButton_2_clicked() { QString qpath = QFileDialog::getOpenFileName(this)。 if(readXYZ(path) == 1) { char *log = 讀取控制點坐標錯誤 XXXXXX。 return。 uitextBrowserappend(QString::fromLocal8Bit(log))。 } void HFJHDialog::on_pushButton_4_clicked() { HFJH()。 21 } 圖 1 控件后方交會 特征匹配 include include include iostream include fstream using namespace std。 //構(gòu)造函數(shù) Mmatch::Mmatch(int length, int flag) { Mlength = length。 match()。src, Mat amp。 int xl = i length / 2。 int xr = i + length / 2。 int xdl = 0。 int xdr = length 1。 if (xl0) { xdl = 0 xl。 } if (yl 0) { ydl = 0 yl。 } if (xr = ) { 22 xdr = ( 1) (xr ())。 } if (yr = ) { ydr = () (yr ())。 } Mat sROI = src(Range(yl+1, yr+1), Range(xl+1,