某省某某煤田东北域灯湖矿区控制网计划书Word下载.docx
《某省某某煤田东北域灯湖矿区控制网计划书Word下载.docx》由会员分享,可在线阅读,更多相关《某省某某煤田东北域灯湖矿区控制网计划书Word下载.docx(39页珍藏版)》请在冰点文库上搜索。
该二等网的主要情况如下:
(1).三角形平均边长为12km.
(2).最小求距角38°
(3).三角形最大闭合差+2.82″,闭合差正负号的分布符合偶然误差的特性。
(4).按三角形闭合差计算所得的测角精度中误差为±
0.82”,平差后为±
0.91”.
(5).最弱边长相对中误差为1:
170000。
(6).仪器检验项目符合规范要求,归心元素的测定正确,观测成果的取舍合理。
(7).造表埋石质量良好。
经现场踏勘,三点标石完好。
所以玉山、太山、广具山三点的坐标可作为起算点坐标,进行不同坐标系之间的转换。
2.高程成果及其精度
国家二等水准路线由西向东横穿测区北部。
根据二等水准路线图,本测区内及其附近应有三个二等水准点,依点之记只找到其中两个:
II2和II4,标石保存完好。
该二等水准路线系国家测绘队于1979年施测。
施测精度及埋石质量均符合规范要求。
高程为1956年黄海高程系。
3.地形图资料
1984年※※航测队航测本地区1:
50000地形图,成图质量良好,可供技术设计图上选点之用。
4.成果的分析利用
收集测区原有资料,了解作业单位、作业时间、执行规范、平差方法等。
注意要点:
起算边长的选择问题、起算方位角的选择问题、考查原有控制点点位的可靠性、针对本次控制网服务的目的作出原有成果的利用方案、测区内有没有国家高级点可以利用。
四.坐标系统的选择及处理方法的论证、起始数据的配置和处理
坐标系统选择1954年北京坐标系,高程系统为1956年黄海高程系统。
如需其它系统的坐标可求的坐标转换参数进行求解。
二等三角点、水准点坐标和高程:
点名
等级
X(m)
Y(m)
H(m)
玉山
二
3793400
20506000
195.4
广具山
3790725
20515900
135.0
太山
3786800
20510250
244.0
II2
35.000
II4
46.000
用已知数据经过计算以下结果:
故以上点均可作为起算点进行整体网平差,解算未知点坐标。
1.平面控制部分
测区内有国家二等网中的控制点三个:
玉山,太山,广具山。
系1959年由西安科技大学测绘队实测。
采用的是1954年北京坐标系。
平差后网的精度:
测角中误差为0.91秒,最弱边边长的相对的中误差为1:
17000。
造标的情况良好。
经勘测标石的质量良好,可以应用。
2.高程控制部分
本测区内有三个二等的水准点,但依点之记只找到其中的两个:
II—2,II-4,且标石的质量良好。
系由1979年西安科技大学测绘队实测。
测绘的精度埋石的情况与质量均符合现行的规范要求。
为1956年黄海高程系统。
本次高程控制网采用了这两个点。
3.坐标系统的确定
因为该矿区位于3º
带第四十带内,其中央子午线为120º
,矿区中心处的横坐标值以,地面平均高程值。
现算出在地面上的长度投影变形比为1:
63300,可见,在国家的统一的坐标的系统里,矿区的投影长度变形并未超过1/4万的要求,故采用国家统一的坐标系。
五.平面控制网布设方案
1.首级网等级和布网方案、加密网设计
首级网:
采用E级GPS网,由三个已知点和八个未知点组成,平均边长2~5千米。
采用GPS静态相对定位测量,观测时段45分钟,GPS接收机用3或4台。
设计出几个同步时间段,构成若干异步环。
估算最弱点中误差。
布点位置见图。
加密网:
采用一级符合导线网和闭合导线网。
每条导线平均长度为500m,全网在GPS1号点和已知点太山之间布设一条附和导线,在GPS6号点处布设一条闭合导线,导线网未知点总共为23个。
详细情况见图。
观测方法采用全站仪,按规范要求测角、测边。
对于1:
2000比例尺测图,首级控制网或加密控制网精度估算后,相邻点最弱点误差以5cm为限。
布设图如下:
2.计算过程中所编写的C++程序
#include<
math.h>
classangle
{
doubled,f,m;
public:
angle(){d=0;
f=0;
m=0;
}
angle(doublexd,doublexf,doublexm){d=xd;
f=xf;
m=xm;
doublegetd(){returnd;
doublegetf(){returnf;
doublegetm(){returnm;
voidput(){cout<
<
this->
getd()<
"
度"
getf()<
分"
getm()<
秒"
endl;
voidinput();
voidoperator=(angle&
a){d=a.d;
f=a.f;
m=a.m;
friendangleadd(angle&
a,angle&
b);
friendanglecha(angle&
anglejian();
anglechu(intr){m=m+f*60+d*60*60;
d=0;
f=0;
m=m/r;
return*this;
angleset(intr){d=d+r*180;
};
////////////////////////
voidangle:
:
input()
{doublea,b,c;
cout<
请输入:
度:
分:
秒:
cin>
>
a>
b>
c;
this->
d=a;
f=b;
m=c;
};
//////////////////////////////////
angleadd(angle&
b)
{anglec;
c.d=a.d+b.d;
c.f=a.f+b.f;
c.m=a.m+b.m;
if(c.m>
60)
{c.m=c.m-60;
c.f=c.f+1;
if(c.f>
{c.f=c.f-60;
c.d+=1;
returnc;
anglecha(angle&
if(a.m>
b.m){c.m=a.m-b.m;
else{a.f=a.f-1;
a.m+=60;
c.m=a.m-b.m;
if(a.f>
b.f){c.f=a.f-b.f;
else{a.d-=1;
a.f+=60;
c.f=a.f-b.f;
c.d=a.d-b.d;
/////////////////////////////////////
angleangle:
jian()
{
while(this->
m>
=60)
{this->
m-=60;
f+=1;
m<
0)
m+=60;
f-=1;
f>
f-=60;
d+=1;
f<
f+=60;
d-=1;
d>
=360)
d-=360;
d<
d+=360;
////////////////
voidmain()
inti=0,j,k=1;
doublepi=3.1415925,p=2062.65;
anglebian[99],jiao[100];
请输入导线的转角个数:
j;
while(i<
j)
{
cout<
请输入第"
i+1<
个转角的角度:
jiao[i].input();
i++;
}
///////////////////////////////
angleqsb,zb;
请输入起始边的方位角:
qsb.input();
请输入终结边的方位角:
zb.input();
//////////////////////////////
////////////////////////////////
anglebhc(0,0,0),fj(360,0,0);
///////////////角度闭合差
bhc.operator=(add(bhc,qsb));
for(intt=0;
t<
t++)
{bhc.operator=(add(bhc,jiao[t]));
bhc.set(j);
bhc.operator=(cha(bhc,zb));
bhc.jian();
if(bhc.getd()>
{bhc.operator=(cha(bhc,fj));
}
if(bhc.getd()<
{bhc.chu
(1);
输出闭合差:
bhc.put();
/////////////////////////角度闭合差分配
bhc.chu(-j);
for(intm=0;
m++)
{jiao[m].operator=(add(jiao[m],bhc));
for(m=0;
{cout<
第"
m+1<
个角度平差后其值为:
;
jiao[m].put();
bian[0].operator=(add(qsb,jiao[0]));
bian[0].set(-1);
bian[0].jian();
while(k<
j-1)
{bian[k].operator=(add(bian[k-1],jiao[k]));
bian[k].set(-1);
bian[k].jian();
k++;
j-1;
{cout<
条边的方位角为:
bian[m].put();
doubleqsdx,qsdy,zdx,zdy,chax,chay,zc=0,qsdx1,qsdy1;
请输入起始点坐标(先横坐标,后纵坐标):
qsdx>
qsdy;
qsdx1=qsdx;
qsdy1=qsdy;
请输入终点坐标(先横坐标,后纵坐标):
zdx>
zdy;
doubles[99],x[99],y[99];
{
条边的长度(单位为米):
cin>
s[m];
zc+=s[m];
//////////////////////////坐标差计算
for(m=0;
{doublew=bian[m].getd()*pi/180+bian[m].getf()/60*pi/180+bian[m].getm()/60/60*pi/180;
x[m]=s[m]*cos(w);
y[m]=s[m]*sin(w);
//////////////////////////////////////////坐标闭合差
{qsdx1=qsdx1+x[m];
qsdy1=qsdy1+y[m];
chax=qsdx1-zdx;
chay=qsdy1-zdy;
全长横坐标闭合差总和:
chax<
全长纵坐标闭合差总和:
chay<
m++)///////////坐标闭合差分配
{x[m]-=(s[m]*chax/zc);
y[m]-=(s[m]*chay/zc);
///////////////////////////////////////
条边的坐标差:
横坐标差:
x[m]<
纵坐标差:
y[m]<
endl<
doublezbx[99],zby[99];
zbx[0]=qsdx+x[0];
zby[0]=qsdy+y[0];
for(m=1;
j-2;
{zbx[m]=zbx[m-1]+x[m];
zby[m]=zby[m-1]+y[m];
for(m=0;
输出第"
个点的近似坐标:
横坐标:
zbx[m]<
纵坐标:
zby[m]<
//////////////////////////////////////////////////
doublea[99],b[99],c[99],d[99],e[99],f[99];
{a[m]=-x[m]/s[m];
b[m]=-y[m]/s[m];
c[m]=x[m]/s[m];
d[m]=y[m]/s[m];
注意啦:
在边长及角度系数中,凡大于"
j-1<
,或小于等于0的系数都为已知点系数,其值均为零!
!
a[0]=0;
b[0]=0;
c[j-2]=0;
d[j-2]=0;
条边的系数:
x"
"
a[m]<
y"
b[m]<
x"
c[m]<
d[m]<
c[0]=0;
d[0]=0;
e[0]=-p*y[0]/(s[0]*s[0]);
f[0]=p*x[0]/(s[0]*s[0]);
a[m]=p*y[m-1]/(s[m-1]*s[m-1]);
b[m]=-p*x[m-1]/(s[m-1]*s[m-1]);
c[m]=p*(y[m]/(s[m]*s[m])-y[m-1]/(s[m-1]*s[m-1]));
d[m]=-p*(x[m]/(s[m]*s[m])-x[m-1]/(s[m-1]*s[m-1]));
e[m]=-p*y[m]/(s[m]*s[m]);
f[m]=p*x[m]/(s[m]*s[m]);
a[1]=0;
b[1]=0;
a[j-2]=p*y[j-3]/(s[j-3]*s[j-3]);
b[j-2]=-p*x[j-3]/(s[j-3]*s[j-3]);
c[j-2]=p*(y[j-2]/(s[j-2]*s[j-2])-y[j-3]/(s[j-3]*s[j-3]));
d[j-2]=-p*(x[j-2]/(s[j-2]*s[j-2])-x[j-3]/(s[j-3]*s[j-3]));
e[j-2]=0;
f[j-2]=0;
a[j-1]=p*y[j-2]/(s[j-2]*s[j-2]);
b[j-1]=-p*x[j-2]/(s[j-2]*s[j-2]);
c[j-1]=0;
d[j-1]=0;
e[j-1]=0;
f[j-1]=0;
////////////////////////////////////////////
在边长及角度系数中,凡大于
个转角的系数为:
m-1<
e[m]<
f[m]<
欢迎使用,敬请回顾!
3.各级控制网精度估算及其结果
GPS网:
从每个观测时段选取两条独立基线,最好各时段选取的独立基线不重复,将未知点坐标设为未知数,列出用未知数表示观测值的方程组,从而可以得出未知数的系数矩阵,即B矩阵。
间接平差误差方程V=AX+L;
法方程:
NX+B=0;
单位权中误差估值:
n为误差方程数,k为卫星数
平差值的精度估算:
根据布网结果所得B矩阵为:
1.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0;
0.01.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0;
0.00.00.00.01.00.00.00.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.01.00.00.00.00.00.00.00.00.00.00.0;
-1.00.01.00.00.00.00.00.00.00.00.00.00.00.00.00.0;
0.0-1.00.01.00.00.00.00.00.00.00.00.00.00.00.00.0;
-1.00.00.00.00.00.01.00.00.00.00.00.00.00.00.00.0;
0.0-1.00.00.00.00.00.01.00.00.00.00.00.00.00.00.0;
1.00.00.00.00.00.0-1.00.00.00.00.00.00.00.00.00.0;
0.01.00.00.00.00.00.0-1.00.00.00.00.00.00.00.00.0;
0.00.01.00.00.00.00.00.00.00.00.00.00.00.00.00.0;
0.00.00.01.00.00.00.00.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.00.01.00.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.00.00.01.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.00.0-1.00.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.00.00.0-1.00.00.00.00.00.00.00.00.0;
0.00.00.00.00.00.00.00.0