測繪類c#程序代碼

-.z.常用測量程序設(shè)計代碼〔1〕用全站儀在A點觀測了B點斜邊和垂直角，求A到B的高差?！蔡崾荆海珼--斜邊，--垂直角，--儀器高，--反光鏡高，--大氣折光系數(shù)〕usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入斜邊=");doubleD=double.Parse(Console.ReadLine());Console.Write("請輸入垂直角[ddd.mmss]=");doublea=DEG(double.Parse(Console.ReadLine()));Console.Write("請輸入儀器高=");doublei=double.Parse(Console.ReadLine());Console.Write("請輸入反光鏡高=");doublev=double.Parse(Console.ReadLine());doubleh=D*Math.Sin(a)+(1-0.13)*D/6371000.0*D/6371000.0*Math.Cos(a)*Math.Cos(a)/2.0+i-v;Console.WriteLine("高差為{0}",h);}//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}}}〔2〕如下圖，A點的坐標(biāo)及A點到B點的邊長及方位角，計算B點的坐標(biāo)。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入A點的*坐標(biāo)=");double*A=double.Parse(Console.ReadLine());Console.Write("請輸入A點的Y坐標(biāo)=");doubleYA=double.Parse(Console.ReadLine());Console.Write("請輸入A到B的方位角[ddd.mmss]=");doublea=DEG(double.Parse(Console.ReadLine()));Console.Write("請輸入A到B的水平距離=");doubleS=double.Parse(Console.ReadLine());double*B=*A+S*Math.Cos(a);doubleYB=YA+S*Math.Sin(a);Console.WriteLine("B點的坐標(biāo)〔{0},{1}〕",*B,YB);}//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}}}〔3〕如下圖，A點和B點的坐標(biāo)，計算A點到的邊長及方位角。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入A點的*坐標(biāo)=");double*A=double.Parse(Console.ReadLine());Console.Write("請輸入A點的Y坐標(biāo)=");doubleYA=double.Parse(Console.ReadLine());Console.Write("請輸入B點的*坐標(biāo)=");double*B=double.Parse(Console.ReadLine());Console.Write("請輸入B點的Y坐標(biāo)=");doubleYB=double.Parse(Console.ReadLine());doubleS=距離(*A,YA,*B,YB);doublea=方位角(*A,YA,*B,YB);Console.WriteLine("AB間的距離={0}，從A到B的方位角={1}",S,DMS(a));}//將弧度轉(zhuǎn)為ddd.mmssstaticpublicdoubleDMS(doubleang){ang+=1.0E-15;//加上一個小量，以保證進位intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang)*180.0/Math.PI;intd=(int)ang;ang=(ang-d)*60.0;intm=(int)ang;doubles=(ang-m)*60.0;return(d+m/100.0+s/10000.0)*fuhao;}//計算方位角,返回弧度值publicstaticdouble方位角(double*1,doubley1,double*2,doubley2){doubledelta*=*2-*1;doubledeltaY=y2-y1;doubleangle=Math.PI*0.5;if(Math.Abs(delta*)>0.000000001){angle=Math.Atan2(deltaY,delta*);}if(angle<0){angle+=Math.PI;}if(deltaY<0.0){angle+=Math.PI;}returnangle;}//計算距離publicstaticdouble距離(double*1,doubley1,double*2,doubley2){returnMath.Sqrt((*2-*1)*(*2-*1)+(y2-y1)*(y2-y1));}}}〔4〕在如下圖的支中導(dǎo)線，A點到M點的坐標(biāo)方位角及每個左角，求每條邊的坐標(biāo)方位角。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入A點到M點的坐標(biāo)方位角=");//將A點到M點的坐標(biāo)方位角換算為M點到A點的坐標(biāo)方位角，以便利用公式doublea0=DEG(double.Parse(Console.ReadLine())+180.0);if(a0>2*Math.PI){a0-=2*Math.PI;}List<double>導(dǎo)線轉(zhuǎn)角集合=newList<double>();inti=1;do{Console.Write("請輸入第{0}個轉(zhuǎn)角的水平角[左角為正，右角為負]<直接回車完畢輸入>=",i++);stringstr=Console.ReadLine();if(str!=""){導(dǎo)線轉(zhuǎn)角集合.Add(DEG(double.Parse(str)));}else{break;}}while(true);i=1;foreach(doubleain導(dǎo)線轉(zhuǎn)角集合){a0+=a+Math.PI;if(a0>2*Math.PI){a0-=2*Math.PI;}elseif(a0<0.0){a0+=2*Math.PI;}Console.WriteLine("第{0}條邊的方位角為{1}",i++,DMS(a0));}}//將弧度轉(zhuǎn)為ddd.mmssstaticpublicdoubleDMS(doubleang){ang+=1.0E-15;//加上一個小量，以保證進位intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang)*180.0/Math.PI;intd=(int)ang;ang=(ang-d)*60.0;intm=(int)ang;doubles=(ang-m)*60.0;return(d+m/100.0+s/10000.0)*fuhao;}//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}}}〔5〕在如下圖的支中導(dǎo)線，A點和M點的坐標(biāo)及每個左角和每條邊長，求每個點的坐標(biāo)。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入M點的*坐標(biāo)=");doubleM*=double.Parse(Console.ReadLine());Console.Write("請輸入M點的y坐標(biāo)=");doubleMy=double.Parse(Console.ReadLine());Console.Write("請輸入A點的*坐標(biāo)=");doubleA*=double.Parse(Console.ReadLine());Console.Write("請輸入A點的y坐標(biāo)=");doubleAy=double.Parse(Console.ReadLine());List<double>導(dǎo)線轉(zhuǎn)角集合=newList<double>();List<double>導(dǎo)線邊長集合=newList<double>();inti=1;do{Console.Write("請輸入第{0}個轉(zhuǎn)角的水平角[左角為正，右角為負]<直接回車完畢輸入>=",i);stringstr=Console.ReadLine();if(str!=""){導(dǎo)線轉(zhuǎn)角集合.Add(DEG(double.Parse(str)));}else{break;}Console.Write("請輸入第{0}條邊長值=",i++);導(dǎo)線邊長集合.Add(double.Parse(Console.ReadLine()));}while(true);//計算M到A的坐標(biāo)方位角doublea0=方位角(M*,My,A*,Ay);//計算每個點的坐標(biāo)double*0=A*;doubley0=Ay;for(intj=0;j<導(dǎo)線轉(zhuǎn)角集合.Count;j++){a0+=導(dǎo)線轉(zhuǎn)角集合[j]+Math.PI;if(a0>2*Math.PI)a0-=2*Math.PI;*0=*0+導(dǎo)線邊長集合[j]*Math.Cos(a0);y0=y0+導(dǎo)線邊長集合[j]*Math.Sin(a0);Console.WriteLine("P{0}點的坐標(biāo)是:{1},{2}",j+2,*0,y0);}}//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}//計算方位角,返回弧度值publicstaticdouble方位角(double*1,doubley1,double*2,doubley2){doubledelta*=*2-*1;doubledeltaY=y2-y1;doubleangle=Math.PI*0.5;if(Math.Abs(delta*)>0.000000001){angle=Math.Atan2(deltaY,delta*);}if(angle<0){angle+=Math.PI;}if(deltaY<0.0){angle+=Math.PI;}returnangle;}}}〔6〕在如下圖的單一附合水準路線中，A點和B點的高程及每段的長度和高差，求每個點的高程。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入A點的高程=");doubleHa=double.Parse(Console.ReadLine());Console.Write("請輸入B點的高程=");doubleHb=double.Parse(Console.ReadLine());List<double>路線段長度集合=newList<double>();List<double>路線段高差集合=newList<double>();inti=1;do{Console.Write("請輸入第{0}段長度[公里為單位]<直接回車完畢輸入>=",i);stringstr=Console.ReadLine();if(str!=""){路線段長度集合.Add(double.Parse(str));}else{break;}Console.Write("請輸入第{0}段高差=",i++);路線段高差集合.Add(double.Parse(Console.ReadLine()));}while(true);//計算水準路線閉合差doublew=0;double總長度=0;for(intj=0;j<路線段長度集合.Count;j++){w+=路線段高差集合[j];總長度+=路線段長度集合[j];}w=w-(Hb-Ha);//閉合差doublev=-w/總長度;//每公里改正數(shù)//計算每個點的高程doubleH=Ha;for(intj=0;j<路線段長度集合.Count;j++){H+=路線段高差集合[j]+v*路線段長度集合[j];Console.WriteLine("P{0}點的高程是:{1}",j+1,H);}Console.WriteLine("最后一點就是B點！");}}}〔7〕在如下圖的前方交會中，ABJ三點按逆時針方向排列，AB兩點的坐標(biāo)和兩個交會角度，求J點的坐標(biāo)。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){Console.Write("請輸入A點的*坐標(biāo)=");doubleA*=double.Parse(Console.ReadLine());Console.Write("請輸入A點的y坐標(biāo)=");doubleAy=double.Parse(Console.ReadLine());Console.Write("請輸入B點的*坐標(biāo)=");doubleB*=double.Parse(Console.ReadLine());Console.Write("請輸入B點的y坐標(biāo)=");doubleBy=double.Parse(Console.ReadLine());Console.Write("請輸入α的角度值=");doubleα=DEG(double.Parse(Console.ReadLine()));Console.Write("請輸入β的角度值=");doubleβ=DEG(double.Parse(Console.ReadLine()));//計算A到B的坐標(biāo)方位角doublea0=方位角(A*,Ay,B*,By);//計算A到J的坐標(biāo)方位角a0-=α;//計算A到B的邊長doubles=距離(A*,Ay,B*,By);//計算A到J的邊長s=s/Math.Sin(Math.PI-α-β)*Math.Sin(β);//計算J的坐標(biāo)doubleJ*=A*+s*Math.Cos(a0);doubleJy=Ay+s*Math.Sin(a0);Console.WriteLine("J點的坐標(biāo)={0},{1}",J*,Jy);}//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}//計算方位角,返回弧度值publicstaticdouble方位角(double*1,doubley1,double*2,doubley2){doubledelta*=*2-*1;doubledeltaY=y2-y1;doubleangle=Math.PI*0.5;if(Math.Abs(delta*)>0.000000001){angle=Math.Atan2(deltaY,delta*);}if(angle<0){angle+=Math.PI;}if(deltaY<0.0){angle+=Math.PI;}returnangle;}//計算距離publicstaticdouble距離(double*1,doubley1,double*2,doubley2){returnMath.Sqrt((*2-*1)*(*2-*1)+(y2-y1)*(y2-y1));}}}〔8〕在如下圖的單個多邊形中，每個點的坐標(biāo)，求多邊形的面積。提示面積的計算公式為：usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespaceConsoleApplication1{classApplication{staticvoidMain(string[]args){List<double>多邊形*坐標(biāo)集=newList<double>();List<double>多邊形Y坐標(biāo)集=newList<double>();inti=1;do{Console.Write("請輸入第{0}點的*坐標(biāo)<直接回車就完畢輸入>=",i);stringstr=Console.ReadLine();if(str!=""){多邊形*坐標(biāo)集.Add(double.Parse(str));}else{break;}Console.Write("請輸入第{0}點的y坐標(biāo)=",i++);多邊形Y坐標(biāo)集.Add(double.Parse(Console.ReadLine()));}while(true);//計算面積doublearea=0;//點的個數(shù)intCount=多邊形*坐標(biāo)集.Count;for(intj=0;j<Count-1;j++){area+=多邊形*坐標(biāo)集[j]*多邊形Y坐標(biāo)集[j+1]-多邊形*坐標(biāo)集[j+1]*多邊形Y坐標(biāo)集[j];}//計算最后一點到起點的項area+=多邊形*坐標(biāo)集[Count-1]*多邊形Y坐標(biāo)集[0]-多邊形*坐標(biāo)集[0]*多邊形Y坐標(biāo)集[Count-1];Console.WriteLine("多邊形的面積={0}",Math.Abs(area*0.5));}}}〔9〕在如下圖的僅有一個連接角的附合導(dǎo)線中，A點和B點的高程及每條邊的斜邊、垂直角、儀器高和反光鏡高，求每個點的高程?！蔡崾荆?，D--斜邊，--垂直角，--儀器高，--反光鏡高，--大氣折光系數(shù)〕usingSystem;usingSystem.Collections.Generic;usingSystem.Linq;usingSystem.Te*t;namespace導(dǎo)線{classProgram{staticvoidMain(string[]args){Console.Write("請輸入A點的高程=");doubleHA=double.Parse(Console.ReadLine());Console.Write("請輸入B點的高程=");doubleHB=double.Parse(Console.ReadLine());List<double>高差集合=newList<double>();List<double>水平距離平方集合=newList<double>();inti=1;double高差和=0;double水平距離平方和=0;do{Console.Write("請輸入{0}點到{1}點的斜距D<直接回車完畢輸入>=",i,i+1);stringstr=Console.ReadLine();if(str=="")break;doubleD=double.Parse(str);Console.Write("請輸入{0}點到{1}點的垂直角=",i,i+1);doubleα=導(dǎo)線.DEG(double.Parse(Console.ReadLine()));//并化為弧度Console.Write("請輸入{0}的儀器高=",i);double儀器高=double.Parse(Console.ReadLine());Console.Write("請輸入{0}的反光鏡高=",i+1);double反光鏡高=double.Parse(Console.ReadLine());i++;double水平距離;doubleh=導(dǎo)線.高差(D,α,儀器高,反光鏡高,out水平距離);水平距離平方集合.Add(水平距離*水平距離);高差集合.Add(h);水平距離平方和+=水平距離*水平距離;高差和+=h;}while(true);//計算高差閉合差double高差閉合差=HA+高差和-HB;//計算單位長度改正數(shù),三角高程的權(quán)與距離的平方成反比double單位長度改正數(shù)=-高差閉合差/水平距離平方和;//計算每點高程doubleH=HA;for(intj=0;j<高差集合.Count;j++){H+=高差集合[j]+單位長度改正數(shù)*水平距離平方集合[j];Console.WriteLine("第{0}點的高程H={1}",j+2,H);}Console.WriteLine("最后一個點的高程應(yīng)該與B點的高程一樣!");}}class導(dǎo)線{//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}//高差計算,返回高差publicstaticdouble高差(double斜距,double垂直角,double儀器高,double反光鏡高,outdouble水平距離){doubleh=斜距*Math.Sin(垂直角)+斜距*Math.Cos(垂直角)/6371000.0*斜距*Math.Cos(垂直角)*(1-0.13)+儀器高-反光鏡高;水平距離=斜距*Math.Cos(垂直角);returnh;}}}〔10〕在如下圖的僅有一個連接角的附合導(dǎo)線中，M點、A點和B點的坐標(biāo)，每個點的左角及每條邊的水平距離，求每個點的坐標(biāo)。usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespace導(dǎo)線{classProgram{staticvoidMain(string[]args){Console.Write("請輸入M點的*坐標(biāo)=");doubleM*=double.Parse(Console.ReadLine());Console.Write("請輸入M點的Y坐標(biāo)=");doubleMY=double.Parse(Console.ReadLine());Console.Write("請輸入A點的*坐標(biāo)=");doubleA*=double.Parse(Console.ReadLine());Console.Write("請輸入A點的Y坐標(biāo)=");doubleAY=double.Parse(Console.ReadLine());Console.Write("請輸入B點的*坐標(biāo)=");doubleB*=double.Parse(Console.ReadLine());Console.Write("請輸入B點的Y坐標(biāo)=");doubleBY=double.Parse(Console.ReadLine());//計算M到A的方位角doubleα0=導(dǎo)線.方位角(M*,MY,A*,AY);List<double>坐標(biāo)增量*集合=newList<double>();List<double>坐標(biāo)增量Y集合=newList<double>();inti=1;double坐標(biāo)增量*和=0;double坐標(biāo)增量Y和=0;double坐標(biāo)增量*絕對值和=0;double坐標(biāo)增量Y絕對值和=0;double水平距離和=0;do{Console.Write("請輸入{0}點到{1}點的水平距離S<直接回車完畢輸入>=",i,i+1);stringstr=Console.ReadLine();if(str=="")break;doubleS=double.Parse(str);Console.Write("請輸入{0}點的水平角<左角為正,右角為負>=",i);doubleβ=導(dǎo)線.DEG(double.Parse(Console.ReadLine()));i++;//計算坐標(biāo)增量doubled*,dy;α0=導(dǎo)線.坐標(biāo)(0,0,β,S,α0,outd*,outdy);坐標(biāo)增量*集合.Add(d*);坐標(biāo)增量Y集合.Add(dy);坐標(biāo)增量*和+=d*;;坐標(biāo)增量Y和+=dy;;坐標(biāo)增量*絕對值和+=Math.Abs(d*);坐標(biāo)增量Y絕對值和+=Math.Abs(dy);水平距離和+=S;}while(true);//計算閉合差double閉合差*=A*+坐標(biāo)增量*和-B*;double閉合差Y=AY+坐標(biāo)增量Y和-BY;//計算單位長度改正數(shù)double*坐標(biāo)單位長度改正數(shù)=-閉合差*/坐標(biāo)增量*絕對值和;doubleY坐標(biāo)單位長度改正數(shù)=-閉合差Y/坐標(biāo)增量Y絕對值和;//計算導(dǎo)線相對精度intf=(int)(水平距離和/Math.Sqrt(閉合差**閉合差*+閉合差Y*閉合差Y));Console.WriteLine("導(dǎo)線的相對精度為1/{0}",f);//計算每點坐標(biāo)double*=A*;doubleY=AY;for(intj=0;j<坐標(biāo)增量*集合.Count;j++){*+=坐標(biāo)增量*集合[j]+*坐標(biāo)單位長度改正數(shù)*Math.Abs(坐標(biāo)增量*集合[j]);Y+=坐標(biāo)增量Y集合[j]+Y坐標(biāo)單位長度改正數(shù)*Math.Abs(坐標(biāo)增量Y集合[j]);Console.WriteLine("第{0}點的*坐標(biāo)={1},Y坐標(biāo)={2}",j+2,*,Y);}Console.WriteLine("最后一個點的坐標(biāo)應(yīng)該與B點的坐標(biāo)一樣!");}}class導(dǎo)線{//將ddd.mmss轉(zhuǎn)為弧度staticpublicdoubleDEG(doubleang){intfuhao=(int)(ang/Math.Abs(ang));ang=Math.Abs(ang);intd=(int)ang;intm=((int)(ang*100))-d*100;doubles=ang*10000-m*100-d*10000;return((d+m/60.0+s/3600.0)*fuhao)/180.0*Math.PI;}//計算方位角,返回弧度值publicstaticdouble方位角(double*1,doubley1,double*2,doubley2){doubledelta*=*2-*1;doubledeltaY=y2-y1;doubleangle=Math.PI*0.5;if(Math.Abs(delta*)>0.000000001){angle=Math.Atan2(deltaY,delta*);}if(angle<0){angle+=Math.PI;}if(deltaY<0.0){angle+=Math.PI;}returnangle;}//計算坐標(biāo),返回點到計算點的方位角publicstaticdouble坐標(biāo)(double*0,doubley0,double左角,double水平距離,double方位角,outdouble*,outdoubley){double方位角=方位角+左角+Math.PI;//將方位角調(diào)整到0到2π之間if(方位角>=Math.PI*2)方位角-=Math.PI*2;if(方位角<0.0)方位角+=Math.PI*2;*=*0+水平距離*Math.Cos(方位角);y=y0+水平距離*Math.Sin(方位角);return方位角;}}}〔11〕在如下圖的僅有一個連接角的附合導(dǎo)線中，M點、A點和B點的坐標(biāo)，A點和B點的高程，每個點的左角及每條邊的斜邊、垂直角、儀器高和反光鏡高，求每個點的坐標(biāo)和高程?！蔡崾荆海珼--斜邊，--垂直角，--儀器高，--反光鏡高，--大氣折光系數(shù)〕usingSystem;usingSystem.Collections.Generic;usingSystem.Te*t;namespace導(dǎo)線{classProgram{staticvoidMain(string[]args){Console.Write("請輸入M點的*坐標(biāo)=");doubleM*=double.Parse(Console.ReadLine());Console.Write("請輸入M點的Y坐標(biāo)=");doubleMY=double.Parse(Console.ReadLine());Console.Write("請輸入A點的*坐標(biāo)=");doubleA*=double.Parse(Console.ReadLine());Console.Write("請輸入A點的Y坐標(biāo)=");doubleAY=double.Parse(Console.ReadLine());Console.Write("請輸入A點的高程=");doubleHA=double.Parse(Console.ReadLine());Console.Write("請輸入B點的*坐標(biāo)=");doubleB*=double.Parse(Console.ReadLine());Console.Write("請輸入B點的Y坐標(biāo)=");doubleBY=double.Parse(Console.ReadLine());Console.Write("請輸入B點的高程=");doubleHB=double.Parse(Console.ReadLine());//計算M到A的方位角doubleα0=導(dǎo)線.方位角(M*,MY,A*,AY);List<double>坐標(biāo)增量*集合=newList<double>();List<double>坐標(biāo)增量Y集合=newList<double/r