【正文】 0 = sqrt(double(0) / (2 * () 6))。 index++。 kappa+=double(5,0)。 phi+=double(3,0)。 Ys+=double(1,0)。 } //計(jì)算 dX = (() * A).inv() * () * L。 (Ai)。 cv::Mat Ai(2, 6, CV_64FC1, a)。 //合成矩陣 double a[12] = { a11, a12, a13, a14, a15, a16, a21, a22, a23, a24, a25, 20 a26 }。 a25 = fy*cos(kappa) y / fy*(x*sin(kappa) + y*cos(kappa))。 a16 = y。 a14 = y*sin(omega) (x / fx*(x*cos(kappa) y*sin(kappa)) + fx*cos(kappa))*cos(omega)。 a22 = 1 / Z_*(b2*fy + b3*y)。 a13 = 1 / Z_*(c1*fx + c3*x)。 //求矩陣 A各元素的值 a11 = 1 / Z_*(a1*fx + a3*x)。 //求矩陣 L各元素的值 lx = x x_。 //求 x， y的 近似值 double x_ = fx*X_ / Z_。 Y_ = a2*(X Xs) + b2*(Y Ys) + c2*(Z Zs)。 c3 = cos(phi)*cos(omega)。 c1 = sin(phi)*cos(kappa) + cos(phi)*sin(omega)*sin(kappa)。 b2 = cos(omega)*cos(kappa)。 a3 = sin(phi)*cos(omega)。 //求矩陣 R各元素的值 a1 = cos(phi)*cos(kappa) sin(phi)*sin(omega)*sin(kappa)。 //經(jīng)過畸變改正后的像點(diǎn)坐標(biāo)觀測值 x = (x x0) dx。 dx = (x x0)*(k1*r_r + k2*r_r*r_r) + p1*(r_r + 2 * (x x0)*(x x0)) 2 * p2*(x x0)*(y y0)。 continue。 char *arg2 = 沒有找到對(duì)應(yīng)的控制點(diǎn) 。 break。 Y = obj[j].Y。 j ()。 Z = 0。 X = 0。 ++i) { 19 x = pic[i].x。 do { for (size_t i = 0。 x0 = W / 2 + x0。 //求取主點(diǎn)坐標(biāo) double W = 。 //單位權(quán)中誤差 double m0。 //, A0, L0, A1。 double r_r, dx, dy。 double a21,a22,a23,a24,a25,a26。 double X_,Y_,Z_。 //定義 double x, y, X, Y, Z。 } //空間后方交會(huì)算法主體函數(shù) void HFJHDialog::HFJH() { //結(jié)果值 double Xs=Xs0, Ys=Ys0, Zs=Zs0。 } ()。 uitextBrowserappend(QString::fromLocal8Bit(log))。 (path)。 //path =F:\\我的課程 \\數(shù)字?jǐn)z影測量影像數(shù)據(jù) \\后方交會(huì) \\像點(diǎn)坐標(biāo) .txt。 return 0。 (tmp)。 uitextBrowserappend(QString(Point%1: X:%2 Y:%3 Z:%3).arg(index).arg().arg().arg())。 ++i) { filhead。 for(int i = 0。 } XYZ tmp。 uitextBrowserappend(QString::fromLocal8Bit(log))。 (path)。 } //讀取物方點(diǎn)坐標(biāo) int HFJHDialog::readXYZ(string path) { int index = 1。 } ()。 index++。 while(!()) { fil。 return 1。 if(!fil) { char *log = 打開像點(diǎn)文件失敗 。 ifstream fil。 } //讀取像點(diǎn)坐標(biāo) int HFJHDialog::readxy(string path) { int index = 1。 kappa0=。 //單位 mm phi0=。 Ys0=。 p2=。 k2=。 y0=。 fy=。 uitextBrowsersetText(HFJH)。 QString qtitle = QString::fromLocal8Bit(title)。 HFJHDialog::HFJHDialog(QWidget *parent) : QDialog(parent), ui(new Ui::HFJHDialog) { uisetupUi(this)。對(duì)于拓?fù)錁?gòu)建中被遺漏的要手工處理后再次構(gòu)建 圖 5 拓?fù)錁?gòu)建后的圖 GEOWAY 空三加密 1） 新建工程，數(shù)據(jù)準(zhǔn)備并檢查 依次點(diǎn)擊“新建工程 ?選擇數(shù)據(jù)所在的目錄初始構(gòu)網(wǎng)”然后點(diǎn)擊“生成成果 ?生成單 色點(diǎn)云”用 mashlab 打開 .ply 文件，查看 點(diǎn)云和 構(gòu)網(wǎng)結(jié)果 圖 6 點(diǎn)云圖 圖 7 網(wǎng)形圖 11 2） 控制網(wǎng)平差 “常規(guī)處理 ?光束法平差 ?刪除殘差大于域值的點(diǎn) ?自由網(wǎng)平差”，重復(fù)上述步驟直到最小殘差小于閾值 60（ PS：每次刪除的點(diǎn)不能太多，以防刪除過多正確的點(diǎn)） 3）精確匹配并查看匹配結(jié)果 點(diǎn)擊“常規(guī)處理 ?精確匹配”，進(jìn)一步匹配 圖 8 圖 9 4）刺點(diǎn) 常規(guī)處理光束法平差 ?常規(guī)處理 ?絕對(duì)定向 （ GPS） ?刺點(diǎn) 12 圖 10 刺點(diǎn)界面 5）空三結(jié)果 圖 11 立體測圖（無人機(jī)） 1）新建方案和數(shù)據(jù)準(zhǔn)備 影像管理 ?加載影 模型管理 ?新建模型 控制點(diǎn)編輯 ?導(dǎo)入控制點(diǎn) 像點(diǎn)導(dǎo)入 ?導(dǎo)入像點(diǎn) 外方為元素導(dǎo)入 ?導(dǎo)入外方位元素 2） 初步處理 內(nèi)定向 — 相對(duì)定向 — 絕對(duì)定向 — 核線重采樣 13 圖 12 相對(duì)定向 3）批處理 圖 13 4） 地物采集 新建加工工程，然后依次進(jìn)行高程點(diǎn)、道路、房屋、比高點(diǎn)、等高線的采集 圖 14 地物采集結(jié)果 5） DEM 編輯 點(diǎn)線處理 采集特征線處理 — 選擇特征線處理 面 處理 — 自適應(yīng)內(nèi) 插，量測點(diǎn)內(nèi)插局部置平 、 水域處理 14 圖 15 DEM 編輯前 圖 16 DEM 編輯后 6） TIN 網(wǎng)處理 圖 17 編輯前 圖 18 編輯后 15 圖 19 編輯后 7） 質(zhì)量檢查 設(shè)置質(zhì)檢方案 — 制作方案 — 加入工作方案 圖 20 檢查報(bào)告 圖 21 方案設(shè)置 新建檢查任務(wù)， 指定目錄到數(shù)據(jù)所在目錄 。 圖 4 節(jié)點(diǎn)編輯后的圖 面處理 。 懸掛點(diǎn)和偽節(jié)點(diǎn)處理 ， 工具 — 檢查 — 懸掛偽節(jié)點(diǎn) 進(jìn)行偽節(jié)點(diǎn)檢測。 點(diǎn)擊柵格圖層，加載 tfw 數(shù)據(jù)，加入柵格圖層，然后生成內(nèi)圖廓，格網(wǎng)線 圖 3 加載柵格圖層 勾選跨層輸入編輯 ， 設(shè)置激活圖層，然后點(diǎn)擊相應(yīng)地物 ?設(shè)為 輸入地物類 。至此，三角測量完成 。至此，交互式編輯完成 。 圖 18 選點(diǎn) 圖 19 使用 PATB 平差 然后加入控制點(diǎn)，只要在測區(qū)中加入三個(gè)控制點(diǎn)，然后 使用 PATB 進(jìn)行平差，就可以使用 “ 預(yù)測控制點(diǎn) ” 預(yù)測出其余的控制點(diǎn)的位置 ， 然后依次選出控制點(diǎn)坐標(biāo)。 然后進(jìn)行 “ 內(nèi)定向 ?航帶間偏移量 ” 圖 17 航帶間偏移量 在圖中選取航帶間的同名點(diǎn)，系統(tǒng)會(huì)自動(dòng)計(jì)算航帶間偏移量 設(shè)置 好航帶間 偏移量之后 ， 開始提取連接點(diǎn) ， 最后運(yùn)行 自動(dòng)挑點(diǎn) 程序。依次點(diǎn)擊“高級(jí) ?正射影像制作”彈出制作界面，選定 DEM 以及相應(yīng)的影像數(shù)據(jù)，然后執(zhí)行，結(jié)果如圖 圖 13 拼接后的正射影像 IGS 數(shù)字化測圖 點(diǎn)擊“測圖 ?IGS 數(shù)字化 測圖”進(jìn)入到數(shù)字化測圖界面，然后依次“新建VZV 矢量文件 ?裝載影像數(shù)據(jù)”，然后通過紅綠眼睛對(duì)測區(qū)進(jìn)行矢量化。點(diǎn)擊“高 級(jí) DEM 拼接”彈出 DEM 拼接界面，加載需要拼接的 DEM，用鼠標(biāo)劃定拼接范圍，接下來就可以執(zhí)行拼接。 圖 10 編輯后的 DEM 圖中的紅色部分就是我們編輯的范圍，經(jīng)過編輯后加載 DEM 顯示立體圖，可以看出其中的房屋高度已經(jīng)貼合地面了。 所以需要對(duì)上上述 DEM 進(jìn)行編輯，生成正真意義上的 DEM。 二、 VirtuoZo 實(shí)習(xí)過程 數(shù)據(jù)準(zhǔn)備以及初步處理 1) 新建測區(qū) ?設(shè)置測區(qū)參數(shù)，新建模型 ?設(shè)置模型參數(shù)，設(shè)置相機(jī)參數(shù)?設(shè)置控制點(diǎn)數(shù)據(jù) ?導(dǎo)入影像數(shù)據(jù)，至此，準(zhǔn)備工作完成； 圖 1 設(shè)置測區(qū)參數(shù) 圖 2 設(shè)置模型參數(shù) 圖 3 導(dǎo)入控制點(diǎn) 圖 4 導(dǎo)入相機(jī)參數(shù) 圖 5 導(dǎo)入影像 2) 內(nèi)定向 ?相對(duì)定向 ?絕對(duì)定向，至此影像定向完成； 2 圖 6 內(nèi)定向 圖 7 相對(duì)定向 相對(duì)定向做完之后系統(tǒng)會(huì)自動(dòng)進(jìn)行絕對(duì)定向，我們可以到主界面查看絕對(duì)定向的結(jié)果：下圖是模型 157156 的絕對(duì)定向結(jié)果圖： 圖 8 157156 絕對(duì)定向結(jié)果 3) 然