附合導(dǎo)線平差程序設(shè)計報告(編輯修改稿)
2025-03-12 06:31 本頁面
【文章內(nèi)容簡介】 經(jīng)編好的 導(dǎo)線網(wǎng)間接平差程序，對邊角網(wǎng)、三角網(wǎng)的實例進行了平差，最終得到導(dǎo)線網(wǎng)的待定數(shù)據(jù)以及相關(guān)的精度評定，通過對這些平差結(jié)果的檢查，最后，認為此導(dǎo)線網(wǎng)的間接平差程序可以為測量的平差工作提供一些服務(wù)。雖然此程序可以完成導(dǎo)線網(wǎng)平差一些工作，但是它還有很多的不足地方需要改進以及完善。 平差程序的實現(xiàn)定義了五個結(jié)構(gòu)體（間接平差結(jié)構(gòu)、平面網(wǎng)、平面點、測站、觀測值），通過一系列的函數(shù)將這五個結(jié)構(gòu)體聯(lián)系起來，最終完成平差任務(wù)。 下面闡述一下導(dǎo)線網(wǎng)的間接平差程序?qū)崿F(xiàn)的基本流程： （ 1） 以 DAT 文件的形式輸入導(dǎo)線網(wǎng)中的已知控制點、測 站信息、觀測值、測角中誤差、測距加、乘常數(shù)等所有的已知數(shù)據(jù)，并且這些數(shù)據(jù)的輸入必須符合前面所提到的數(shù)據(jù)輸入格式； （ 2） 運用程序讀取輸入數(shù)據(jù)文件，通過相應(yīng)的函數(shù)求得待定點的近似坐標； （ 3） 通過相應(yīng)函數(shù)，求取間接平差的結(jié)構(gòu)：誤差方程系數(shù)陣 A、觀測權(quán)陣 P、誤差方程常數(shù)項 l； （ 4） 通過間接平差函數(shù)，實現(xiàn)對已經(jīng)得到的間接平差結(jié)構(gòu)的平差工作。 （ 5） 最后，以 DAT 文件的格式輸出并保存平差處理結(jié)果，為以后的測量工作提供參考。 六、部分源代碼： //度分秒 化成弧度 double DMStoRAD(double dms) { int d, m, f。 double s。 f = dms=0 ? 1 : 1。 // ＝ dms += f * 。 d = (int)dms。 dms = (dms d) * 。 m = (int)dms。 s = (dms m) * 。 return (d + m / + s / ) * _TORAD f * 。 地圖學(xué)課程設(shè)計報告 8 } //弧度轉(zhuǎn)“度分秒” double RADtoDMS(double rad) { int f = rad = 0 ? 1 : 1。 // 符號 ＋ － //加 （用弧度表示），化為度 rad = (rad + f * ) * _TODEG。 int d = (int)rad。 rad = (rad d) * 。 int m = (int)rad。 double s = (rad m) * 。 return d + m / + s / f * 。 } //讀入已知點的坐標 void ReadKnwData(FILE * in, DyArray * array) { char buffer[256]。 double dx,dy。 SurPnt pnt。 while ( !feof(in) ) { fgets(buffer, sizeof(buffer), in)。 if ( sscanf(buffer, %s %lf %lf, , amp。, amp。) == 3) { Add(array, pnt)。 } else break。 } dx = arrayelem[1].x arrayelem[0].x。 dy = arrayelem[1].y arrayelem[0].y。 =atan2(dy, dx) + (dy 0 ? 1 : 0) * _2PI。//alfa0 dx = arrayelem[3].x arrayelem[2].x。 dy = arrayelem[3].y arrayelem[2].y。 =atan2(dy, dx) + (dy 0 ? 1 : 0) * _2PI。//alfan =arrayelem[1].x。 =arrayelem[1].y。 地圖學(xué)課程設(shè)計報告 9 =arrayelem[2].x。 =arrayelem[2].y。 Destroy(array)。// } //觀測值信息寫到動態(tài)數(shù)組 void ReadObsValue(FILE * in, DyArray * array) { char buffer[256]。 double dir0,dir1,dist0,dist1。 SurPnt pnt。//用于臨時保存點的坐標 while ( !feof(in) ) {
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1