灯湖矿区控制网设计书,在卫星遥感影像图上进行控制网设计的有关设计方法和过程设计书.doc

应用大地测量学基础课程设计报告书灯湖矿区控制网设计书一 灯湖矿区控制网设计任务通知书测绘队：灯湖矿区位于省煤田东北域。本矿已经列入国家经济建设计划，准备进行重点开发。根据建设计划，首先需要施测80平方公里的1:2000地形图。为了加快建井速度，需在竖井间进行对向掘进，最大贯通距离为3km，矿井最大单翼长度为5km，井深一般不超过800m,。为了满足矿山设计、建井、生产三阶段测绘各种比例尺地形图、井巷贯通以及工业场地施工测量的需要，应在全矿区统一的具有足够精度密度的平面控制网与高程控制网。要求你队在任务通知书下达后一周内完成控制网技术设计，一月内完成建网任务。 省煤炭管理局 年 月 日二 测区概况简述1. 地理概况本测区中心位置为东经11713，北纬3412。测区地面高程为+30+244m。位于省县境内，属于华东丘陵地带，南部及中部为山地。其上林木繁茂，不理通视，其余为平原，密布村庄和厂矿，建筑群遍地，影响观测。灯湖区位于测区东部，湖面开阔，中部有大堤。2. 交通情况测区内有铁路支线通过，公路四通八达，村间大道可行汽车，交通方便。3. 气候情况测区地处华东近海区，气候宜人。全年平均降雨量为1040mm，雨量集中于6、7、8三个月。全年平均气温+15度，夏季气温较高，一般为30度左右，冬季有雪，但不寒冷，最低温度为-5度，冻土线深度为0.1米。全年平均风力为23级，夏季略受台风影响。宜于野外作业时间为311月份，年平均作业利于率为21天/月。4. 居民及居民点测区内地少人多，人口稠密，多为汉族，少数民族为回族。测量作业所需人力、物力、财料及食宿均可就地解决。5. 作业目的及任务范围灯湖矿区位于广州市灯湖煤田东北域，煤藏量为10000亿吨。要求于2005年XX月之前，完成建网任务和提交阳春一矿的54 的1：2000和23 11000地形图。在全矿区建立统一的具有足够精度密度的平面控制网与高程控制网。另外，为贯通工程需布设一个专用控制网，要求两近井点横向相对中误差小于30mm。灯湖矿区的测图总面积为77km2三 已有测绘成果资料及其分析利用1. 三角网成果及其精度测区内及附近有国家二等网点三个：玉山、太山、广具山。系1959年由测绘队施测，作业所依据的规范为一、二、三、四等三角测量细则（1958年）。三点标石保存完好。坐标系统为1954年北京坐标系，三度分带，中央子午线为117。该二等网的主要情况如下：（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)玉山二379340020506000195.4广具山二379072520515900135.0太山二378680020510250244.0II2二35.000II4二46.000用已知数据经过计算以下结果：（1）.三角形平均边长为12km.（2）.最小求距角38。（3）.三角形最大闭合差+2.82，闭合差正负号的分布符合偶然误差的特性。（4）.按三角形闭合差计算所得的测角精度中误差为0.82” ,平差后为0.91”.（5）.最弱边长相对中误差为1:170000。故以上点均可作为起算点进行整体网平差，解算未知点坐标。1平面控制部分测区内有国家二等网中的控制点三个：玉山，太山，广具山。系1959年由西安科技大学测绘队实测。采用的是1954年北京坐标系。平差后网的精度：测角中误差为0.91秒，最弱边边长的相对的中误差为1：17000。造标的情况良好。经勘测标石的质量良好，可以应用。 2高程控制部分本测区内有三个二等的水准点，但依点之记只找到其中的两个：II2，II-4，且标石的质量良好。系由1979年西安科技大学测绘队实测。测绘的精度埋石的情况与质量均符合现行的规范要求。为1956年黄海高程系统。本次高程控制网采用了这两个点。 3.坐标系统的确定因为该矿区位于3带第四十带内，其中央子午线为120，矿区中心处的横坐标值以 ,地面平均高程值 。现算出在地面上的长度投影变形比为1:63300,可见，在国家的统一的坐标的系统里，矿区的投影长度变形并未超过1/4万的要求，故采用国家统一的坐标系。五 平面控制网布设方案1. 首级网等级和布网方案、加密网设计首级网：采用E级GPS网，由三个已知点和八个未知点组成，平均边长25千米。采用GPS静态相对定位测量，观测时段45分钟，GPS接收机用3或4台。设计出几个同步时间段，构成若干异步环。估算最弱点中误差。布点位置见图。加密网：采用一级符合导线网和闭合导线网。每条导线平均长度为500m，全网在GPS1号点和已知点太山之间布设一条附和导线，在GPS6号点处布设一条闭合导线，导线网未知点总共为23个。详细情况见图。观测方法采用全站仪，按规范要求测角、测边。估算最弱点中误差。对于1:2000比例尺测图，首级控制网或加密控制网精度估算后，相邻点最弱点误差以5cm为限。布设图如下:2. 各级控制网精度估算及其结果GPS网：从每个观测时段选取两条独立基线，最好各时段选取的独立基线不重复，将未知点坐标设为未知数，列出用未知数表示观测值的方程组，从而可以得出未知数的系数矩阵，即B矩阵。间接平差误差方程 V=AX+L; 法方程： NX+B=0;单位权中误差估值：m0=VTPV(n-k-2) n为误差方程数，k为卫星数 平差值的精度估算：mi=1pi pi=1Qii根据布网结果所得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.00.00.01.00.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.00.01.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.0-1.00.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.00.0-1.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.00.00.01.00.00.00.0;0.00.00.00.00.00.00.00.00.00.00.00.00.01.00.00.0;0.00.00.00.00.00.00.00.01.00.00.00.00.00.00.00.0;0.00.00.00.00.00.00.00.00.01.00.00.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.00.00.00.00.01.00.0;0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.01.0;0.00.00.00.0-1.00.01.00.00.00.00.00.00.00.00.00.0;0.00.00.00.00.0-1.00.01.00.00.00.00.00.00.00.00.0;0.00.00.00.0-1.00.00.00.01.00.00.00.00.00.00.00.0;0.00.00.00.00.0-1.00.00.00.01.00.00.00.00.00.00.0;0.00.00.00.00.00.0-1.00.01.00.00.00.00.00.00.00.0;0.00.00.00.00.00.00.0-1.00.01.00.00.00.00.00.00.0;0.00.00.00.00.00.00.00.0-1.00.01.00.00.00.00.00.0;0.00.00.00.00.00.00.00.00.0-1.00.01.00.00.00.00.0;0.00.00.00.00.00.00.00.00.00.0-1.00.00.00.01.00.0;0.00.00.00.00.00.00.00.00.00.00.0-1.00.00.00.01.0; 0.00.00.00.00.00.00.00.00.00.01.00.00.00.0-1.00.0;0.00.00.00.00.00.00.00.00.00.00.01.00.00.00.0-1.0;0.00.00.00.00.00.00.00.00.00.00.00.0-1.00.01.00.0;0.00.00.00.00.00.00.00.00.00.00.00.00.0-1.00.01.0; GPS网以3km为单位权，用公式=10mm+10PPMD可计算出单位权方差为3.14 cm，计算出每条独立基线的方差，之后根据公式p=0*0/i*i可以计算出每条独立基线的权值，如下（取对角线）：0.93,0.93,0.57,0.57,1.33,1.33,0.64,0.64,0.64,0.64,0.75,0.75,0.86,0.86,0.86,0.86,0.81,0.81,0.81,0.81,0.58,0.58,0.78,0.78,1.52,1.52,0.62,0.62,0.83,0.83,0.75,0.75,1.23,1.23,1.08,1.08,1.08,1.08,1.60由MATLAB软件可解算出各未知点的协因数阵， 点号坐标GPS1GPS2GPS3GPS4GPS5GPS6GPS7GPS8X0.44220.66160.61740.31080.38830.29730.63420.3259Y0.44220.66160.61740.31080.38830.29730.63420.3259通过比较得最弱点为GPS2号点，其精度为25.72mm,符合限差要求（限差为5cm）；符合导线：起始点GPS1,终点太山，中间布设13个点，相邻点之间距离平均为500m,且该导线经过主井石红。详细情况见图。设未知点坐标为未知数，根据附和导线的关系列出误差方程，求得B系数矩阵:以1km为单位权，根据相应等级的每千米中误差为10mm,用公式得各段的权值；由MATLAB软件可解算出各未知点的协因数阵:点号坐标D1S23D2D3D4D5S5D6D7S7S8D8D9X0.43890.55600.86951.06190.95030.99500.86770.79900.66980.66620.48720.30530.2644Y0.38400.44290.70590.66070.62910.71510.63100.58640.46260.28780.23010.21320.1295通过比较可得最弱点Qx=1.0619,Qy=0.6607, 故可求的最弱点x精度为10.304mm，y精度为8.128mm，均符合限差要求（限差为50mm）； 闭合导线：以G6为起始点，经过10个点后到达起始点，构成闭合环。风井上村在此闭合环中。详细情况见图。点号坐标D10D11S26D11S27S28D12D13S30D14D15X0.36440.52680.59580.64160.56850.32190.32850.29460.18230.12760.1168Y0.23040.44990.53510.54450.46610.25100.21490.18650.15350.13250.1079 D11号点为最弱点，QX11=0.6416,QY11=0.5445, 所以x精度为8.0009mm,y精度为7.3790mm. 3. 精度统计表GPS网精度统计表：等级点数（含已知点）平均边长最弱点中误差GPS网E级113.2318km25.72mm附和导线统计表：等级点数（含已知点）平均边长最弱点中误差一级15494m10.304mm闭合导线统计表：等级点数（含已知点）平均边长最弱点中误差一级11485m8.009mm六 水准高程控制网布设方案1. 水准网采取四等符合或闭合水准网，平均12km一个水准点。从已知水准点II2到II4布设三条符合水准路线，以II4为起始点布设两条闭合水准路线。水准点大多布设于公路道路旁，详细情况见图。2. 网中最弱点精度估算及结果（1）符合水准路线1：从II2开始，经过17个点，包括主井石红到达已知点II4。以1.2km定为单位权，各水准测段权值如下（取对角线）：0.942,1.35,2.69,3.77,2.69,1.57,1,1.256,2.355,1.18,0.9,0.96,1.35,1.884,0.75,1,1.34,1.884精度估算结果（协因数阵取对角线）：点号1234G25678协因数0.97841.56261.82532.0022.22792.57072.98803.21473.2974点号9101112131415G5协因数3.38283.33423.12312.87562.64851.89441.15670.5100从中可以看出9号点为最弱点，Q9=3.3828，故最弱点中误差MF=18.392mm.（2）符合水准路线2：以II2为起始点，进过17个未知点到达已知点II4。路线见图。各测段权值如下：1.05,3.76,2.7,1.29,1.256,1.256,1.57,1.71,1.45,1.05,1.884,1.45,1.71,1.256,1.45,1.35,2.7,2.09计算得协因数阵：点号252423因数0.87491.09151.37301.88632.30662.61862.79012.88652.9251点号47G3464544434142协因数2.84482.73282.51522.26711.83531.37360.78730.4589从中可以看出，48号点处最弱，Q48=2.9251，故此水准路线最弱点中误差为MF=18.735mm.（3）符合水准路线3：以II2为起始点，经过10个未知点到达已知点II4。路线见图。各测段权值如下：2.355,1.884,1.57,1.05,1,2.7,0.9,1.35,1.35,2.7,2.09;计算所得协因数阵：点号25242322212019181716协因数0.40010.83131.24771.66451.83671.83151.59211.24610.75090.4473最弱点为21号点，Q21=1.8367，最弱点中误差为MF=15.37mm.(4)闭合水准路线1：以II4为起点，经过20个点，包括风井上村，到达起始点II4形成闭合环。各段权值如下：1.884,1.34,1,0.75,1.884,1.35,1.256,1.45,1.57,2.09,1.64,1.884,1.11,2.09,1.45,1.177,3.77,1.345,1.11,1.05,3.76,1.05;计算后协因数阵：点号G51514131211G6协因数0.51281.17271.94522.77603.04353.35643.6144点号26272829303253协因数3.77233.84603.89883.90063.86343.71783.5982点号5455561415G5协因数3.37433.01462.88352.46771.86911.1233从中可以看出最弱点为29号点，Q29=3.9006, 最弱点中误差MF=19.749mm.(5).闭合水准路线2：以II4为起点，经过20个点，到达起始点II4形成闭合环。各段权值如下：1.884,1.34,1,1.256,1.345,1.177,1.45,2.09,1.57,2.9,1.884,1.256,0.942,1.256,1.05,1.05,1.18,0.9,1.05,3.76,1.05;计算后协因数：点号G5151456555453协因数0.51321.17531.95372.48422.90833.30843.5670点号3132G833343536协因数3.71153.85973.91883.98073.99783.92023.7634点号373839404142协因数3.47123.06592.61011.87691.12580.8958最弱点为34号点，Q34=3.9978，最弱点中误差为MF=19.994mm七 技术依据及作业方法1. 执行规范为： 一、二、三、四等三角测量细则（1958年），全球定位系统（GPS）测量规范 GB/T 18314-2009，工程测量规范GB50026-2007，国家三、四等水准测量规范GB12898-91,1:500 1:1000 1:2000地形图图式GB/T20257.1-2007.2. 标石图结构为寻常标，材料为钢标，中心标石一般用混凝土灌制而成，也可用坚硬石料凿成。标石的顶部中央均镶嵌一金属磁质标志，建标时应注意：标形端正，基础坚实，标架稳固，标石的标志中心与标石的标志中心应居于同一条铅垂线上。标的顶端安装供观测照准用的标心柱或照准圆筒。3. 仪器的选择及检验项目要求GPS单频接收机，DS3水准仪，全站仪检验项目：水准仪标尺零点差、基辅分化值检验，指标差检验，全站仪检验等。 4观测方法及各项限差GPS网为首级网，观测时分10个观测时段进行观测，每测段时间不少于45分钟，计算后最弱点中误差小于5cm；导线网为加密网，可用全站仪观测，可从一个起始点开始观测，到另一个已知点结束，GPS网点可作为导线网的已知点，导线网最弱点中误差小于5cm；水准网用水准仪进行往返观测，每站前后视距差小于1米，整条路线小于3米，每条路线的长度小于16km，最弱点中误差小于2cm.6. 概算内容和平差方法GPS网概算内容：数据预处理及其基线解算，观测结果的外业检核导线网概算内容：外业观测成果资料的检查，观测成果归算到标石中心的计算，地面观测值归算到椭球面的改正，椭球面上的观测值归算到高斯平面上的计算，主要包括角度归心改正、高差计算、边长斜距改平距、归心改正、化算到高斯平面的改正、计算导线方位角条件和环形条件闭合差、计算测角中误差、计算测距中误差、计算导线相对闭合差。GPS网、导线网、水准网均可以用间接平差法进行精度估算，平差大致过程为：绘制控制网平差略图、编制起算数据表和观测值表、选择参数、将观测量的平差值表达成所选参数的函数，并转化程误差方程、解算法方程求的未知数、计算改正数和平差值及其精度。7. 提供成果资料图上设计用的地形图、控制网设计图、控制点标石图、平面和高程控制网网形图八 工作量综合计算及工作进程计划表外业计算：用全站仪进行导线网观测时，应及时计算出每测回的2c值、2C互差值（一级导线小于13），多测回方向角（方向角互差小于9）；水准仪观测时应及时计算出前后视距差（小于1米）、累计视距差（小于3米）；内业计算：概算和平差（见前述）。工作进程计划表：时间任务备注1天踏勘找点（已知点）查看标石情况7天选位和造标埋石GPS点、水准点、导线点需待标石凝固符合要求后方可观测，导线点需要通视1天GPS网观测每时段观测时间不少于45分钟1天GPS网数据解算及结果精度分析确保精度达标7天水准仪检校水准网施测水准路线需要往返观测，且路线要至少经过4个GPS点3天水准网高程计算和GPS网高程拟合（包括错误测站重测时间和休息时间）确保精度达标3天全站仪检校导线网施测及时计算发现改正错误1天导线网计算及错误测站复测确保精度1天休息1天GPS网、水准网、导线网全面检查查看标石情况2天整理资料及结果、上交图纸1天单位检验验收九 装备、仪器、器材及经费预算装备：运送人员车辆一辆、运送材料卡车一辆；仪器：选择3台GPS单频接收机（配带脚架、卷尺），2台D3水准仪（配带脚架、标尺、尺垫）、全站仪一台（配带脚架、棱镜、气温气压计、卷尺）；器材：草帽、对讲机、标石、混凝土、重力计等经费预算：工 资工 程 师3000元人月*1人技 术 员2000元人月*1人测 工1200元人月*3人临 时 工800元人月*2人单频接收机每台15000元*3台+全站仪每台20000元*1台+DS3水准仪每台2000元*1+标石每座100元*73座+生活补贴费每人/每天100元*30天*7人+交通住宿费每人/每天100元*30天*7人+总工资=122500元十. 上交资料清单 控制网设计图纸和控制网技术设计说明书、仪器检验资料、 外业观测手簿、记簿及点之记、电算资料、 最后成果表、工作总结。附件：计算过程中所编写的C+程序#include#includeclass angledouble d,f,m;public:angle()d=0;f=0;m=0;angle(double xd,double xf,double xm)d=xd;f=xf;m=xm;double getd()return d;double getf()return f;double getm()return m;void put()coutgetd()度getf()分getm()秒endl;void input();void operator=(angle &a)d=a.d;f=a.f;m=a.m;friendangle add(angle &a,angle &b);friendangle cha(angle &a,angle &b);angle jian();angle chu(int r)m=m+f*60+d*60*60; d=0; f=0; m=m/r; return *this;angle set(int r) d=d+r*180; return *this; ;/void angle:input() double a,b,c; cout请输入： 度： 分： 秒：abc; this-d=a;this-f=b;this-m=c;/angle add(angle &a,angle &b) angle c; c.d=a.d+b.d; c.f=a.f+b.f; c.m=a.m+b.m; if(c.m60)c.m=c.m-60;c.f=c.f+1; if(c.f60)c.f=c.f-60; c.d+=1; return c;angle cha(angle &a,angle &b) angle c;if(a.mb.m) c.m=a.m-b.m;else a.f=a.f-1; a.m+=60; c.m=a.m-b.m;if(a.fb.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; return c;/angle angle:jian() while(this-m=60)this-m-=60;this-f+=1; while(this-mm+=60; this-f-=1; while(this-f=60)this-f-=60;this-d+=1; while(this-ff+=60; this-d-=1; while(this-d=360)this-d-=360; while(this-dd+=360; return *this;/void main() int i=0,j,k=1;double pi=3.1415925,p=2062.65;angle bian99,jiao100; cout请输入导线的转角个数：j;while(ij) cout请输入第i+1个转角的角度：endl;jiaoi.input();i+; coutendl;/angle qsb,zb;cout请输入起始边的方位角：endl;qsb.input();cout请输入终结边的方位角：endl;zb.input();/coutendl;/angle bhc(0, 0, 0),fj(360,0,0); /角度闭合差bhc.operator =(add(bhc,qsb);for(int t=0;t0)bhc.operator =(cha(bhc,fj); if(bhc.getd()0) bhc.chu(1);cout输出闭合差：endl;bhc.put();coutendl;/ 角度闭合差分配bhc.chu(-j); for(int m=0;mj;m+) jiaom.operator =(add(jiaom,bhc); for( m=0;mj;m+) cout第m+1个角度平差后其值为：; jiaom.put(); coutendl;bian0.operator =(add(qsb,jiao0);bian0.set(-1);bian0.jian();while(kj-1) biank.operator =(add(biank-1,jiaok); biank.set(-1);biank.jian(); k+;for( m=0;mj-1;m+)cout第m+1条边的方位角为：;bianm.put();coutendl;double qsdx,qsdy,zdx,zdy,chax,chay,zc=0,qsdx1,qsdy1;cout请输入起始点坐标（先横坐标，后纵坐标）：qsdxqsdy;qsdx1=qsdx; qsdy1=qsdy;cout请输入终点坐标（先横坐标，后纵坐标）：zdxzdy;double s99,x99,y99;for( m=0;mj-1;m+)cout请输入第m+1sm;coutendl;zc+=sm;/坐标差计算 for( m=0;mj-1;m+) double w=bianm.getd()*pi/180+bianm.getf()/60*pi/180+bianm.getm()/60/60*pi/180; xm=sm*cos(w); ym=sm*sin(w);/ 坐标闭合差 for( m=0;mj-1;m+)qsdx1=qsdx1+xm;qsdy1=qsdy1+ym;chax=qsdx1-zdx; chay=qsdy1-zdy;cout全长横坐标闭合差总和：chaxendl; cout全长纵坐标闭合差总和：chayendl;for( m=0;mj-1;m+) /坐标闭合差分配xm-=(sm*chax/zc);ym-=(sm*chay/zc);coutendl;/ for( m=0;mj-1;m+)cout第m+1条边的坐标差：endl; cout横坐标差：xm 纵坐标差：ymendlendl;coutendl;double zbx99,zby99;zbx0=qsdx+x0; zby0=qsdy+y0;for(m=1;mj-2;m+)zbxm=zbxm-1+xm;zbym=zbym-1+ym;for(m=0;mj-2;m+)cout输出第m+1个点的近似坐标：;cout 横坐标：zbxm 纵坐标：zbymendl;/double a99,b99,c99,d99,e99,f99;for(m=0;mj-1;m+)am=-xm/sm; bm=-ym/sm;cm=xm/sm; dm=ym/sm;cout注意啦：在边长及角度系数中，凡大于j-1，或小于等于0的系数都为已知点系数，其值均为零！endl;a0=0; b0=0;cj-2=0; dj-2=0;for(m=0;mj-1;m+) coutendl; cout第m+1条边的系数：endl; coutxm: am ym: bm xm+1: cm ym+1: dmendl;/a0=0; b0=0;c0=0;d0=0;e0=-p*y0/(s0*s0);f0=p*x0/(s0*s0);for(m=1;mj-2;m+)am=p*ym-1/(sm-1*sm-1);bm=-p*xm-1/(sm-1*sm-1);cm=p*(ym/(sm*sm)-ym-1/(sm-1*sm-1);dm=-p*(xm/(sm*sm)-xm-1/(sm-1*sm-1);em=-p*ym/(sm*sm);fm=p*xm/(sm*sm);a1=0; b1=0;/aj-2=p*yj-3/(sj-3*sj-3);bj-2=-p*xj-3/(sj-3*sj-3);cj-2=p*(yj-2/(sj-2*sj-2)-yj-3/(sj-3*sj-3);dj-2=-p*(xj-2/(sj-2*sj-2)-xj-3/(sj-3*sj-3);ej-2=0;fj-2=0;/aj-1=p*yj-2/(sj-2*sj-2);bj-1=-p*xj-2/(sj-2*sj-2);cj-1=0;dj-1=0;ej-1=0;fj-1=0;/cout注意啦：在边长及角度系数中，凡大于j-1，或小于等于0的系数都为已知点系数，其值均为零！endl;for(m=0;mj;m+) coutendl第m+1个转角的系数为：endl; coutxm-1: am ym-1: bm xm:cm ym: dm xm+1: em ym+1: fmendl;coutendl欢迎使用，敬请回顾！endl;在卫星遥感影像图上进行控制网设计任务书：对选定的一幅遥感影像图，面积约为100200平方千米，标定已知指点，设计测绘地形图（1:5000）工程所需控制点，按课程设计大纲应要求进行控制网设计。选图设计如下：一已有成果资料图中已知点有D018,D025,D026,其三维坐标如下：点号X（米）Y（米）H（米）备注D01837926703951822240.446平面坐标为D级GPS点，高程为国家二等水准点D02537897333951101538.663D026378900739515904182.514 注： 1980年西安坐标系、中央子午线经度117 三个已知点保存情况完好，可作为首级网的起算点。二 布设方案 GPS网：如图所示，根据云龙湖测区面积大小，共选定12个点为GPS布设点，其均匀分布于整个测区内，布设时点与点之间距离应适中。点号为E01、E02、E03.、E11、E12.GPS网可用三台GPS单频接收机按城市GPS网E级进行观测，分为15个时段，各时段见图。每个时段观测时间不少于45分钟。每个时段选取2条独立基线，全网共