图幅理论面积与图斑椭球面积计算公式及要求
净面积计算
净面积（图斑地类面积）计算净面积=图斑面积（椭球面积）—线状地物面积—零星地物面积—田坎面积田坎面积=（图斑面积（椭球面积）—线状地物面积—零星地物面积）*田坎系数（地类图斑里）1、JBNTBHPK、DLTB与XZ相交，得到TB，目的是小组分割片块,并得到每个小图斑的TKXS。
2、shape_area重算（属性表中计算几何体），把面积太小的小图斑手动删除了，不然到时候算的时候那些小图斑面积太小而里面有线状地物的时候会出现好多负值，这步主要是减小误差。
2、XZDW与TB相交，用小图斑切割线状地物，得到TB_XZDW然后修复几何（选中删除没有几何形状的要素），再多部分转单部分，处理那些不是单一部分的，然后把分割后的线状地物长度（CD）重新算一下接下来就是分析切割后的线状地物的位置，这样相交后线状地物在边界上的会有两个重复的图元，这种都是现状地物在小图斑边界上的3、把TB_XZDW转换成点point（选中内部）生成在原来线的中点位置，新建字段，计算点的X,Y坐标，然后删除相同的（简化--融合）。
这样就删除了重复的线状地物，避免计算重复。
4、把TB转换成线，得到TB_LINE,然后与point空间链接（目标要素图层为融合后的point图层，注意搜索距离的设置，设置过大则Join_Count出现较大的数字，设置过小则可能会出现没有记录），得到point_spatial，给point附上边界属性，然后在point_spatial 属性表中新建mj字段，mj=CD*KD然后看那个join_count字段，也就是一个点上有几条线，判断线状地物的位置，0的在图斑里面，1的在一个图斑的边界上，2的是在两个图斑相交的地方，其中0的不扣，1和2的半扣，手动计算改。
point_spatial与TB空间链接(目标要素图层为TB图层)，并对点中的mj求和求出小图斑内总的线状地物的面积6、同样把零星地物空间链接附到小图斑上，把小图斑椭球面积重算。
椭圆的面积问题
椭圆的面积问题
椭圆是一种常见的平面图形，其具有独特的几何性质。
在本文中，我们将探讨椭圆的面积问题。
首先，让我们回顾一下椭圆的定义。
椭圆是平面上到定点F1、F2距离之和等于定长2a的点P的轨迹。
这两个定点称为焦点，2a称为长轴，沿着长轴的中心称为圆心。
为了计算椭圆的面积，我们需要找到一个公式，其中包含椭圆的重要参数。
我们知道，长轴的长度为2a，短轴的长度为2b。
那么，我们可以使用下面的公式来计算椭圆的面积：
S = πab
其中，π是一个常数，约等于3.14159。
因此，我们可以得出结论：一个椭圆的面积等于其长轴和短轴长度的乘积，再乘以π。
让我们通过一个具体的例子来看看如何计算椭圆的面积。
假设一个椭圆的长轴长度为6，短轴长度为4。
我们可以使用上述公式计算出椭圆的面积：
S = π × 6 × 4/2
= 12π
因此，该椭圆的面积约为37.7。
最后，让我们看一下如何应用椭圆的面积问题。
在实际应用中，我们可能需要计算椭圆形状的物体的面积。
例如，一个篮球场的形状可
能类似于一个椭圆，所以我们可以使用上述公式来计算篮球场的面积。
同样地，椭圆也经常在建筑设计和航天工程中使用。
综上所述，椭圆是一个非常重要的几何形状，在许多不同领域发挥
着作用。
通过使用上述公式，我们可以简单地计算椭圆的面积，并应
用它在实际问题中。
excel中计算椭球面积的方法
excel中计算椭球面积的方法
在Excel中计算椭球面积可以通过数学公式和Excel函数来实现。
椭球的面积计算公式为：S=4πab，其中a和b分别为椭球的两个半轴长度。
在Excel中，可以通过以下步骤来计算椭球的面积：
1. 首先，假设椭球的两个半轴长度分别为a和b，可以在Excel中选择两个单元格分别输入这两个数值。
2. 接下来，在另一个单元格中使用Excel的乘法函数“=ab”来计算a和b的乘积。
3. 然后，在另一个单元格中使用Excel的π函数“=PI()”来获取π的数值。
4. 最后，在另一个单元格中使用Excel的乘法函数将π的数值与a和b的乘积相乘，“=4πab”，即可得到椭球的面积。
这样，你就可以在Excel中使用这些简单的数学公式和Excel 函数来计算椭球的面积了。
当然，你也可以将这些步骤合并到一个
公式中，以便更加简洁地计算椭球的面积。
希望这些信息能够帮助到你。
【学习课件】第07章地球椭球与椭球计算理论
ppt课件
5
地球椭球与椭球计算理论
安徽理工大学
2.地球椭球参数间的相互关系
由前面式子得：
1
2 3 4
e2
a2 b2 a2
e'2
a2 b2 b2
1 e2
b2 a2
1 e2
a2 b2
5 6
并得： (1e2)(1e'2)1
7 8 9
推得：
e2
e' 2 1 e'2
e' 2
e2 1 e2
10
同理可得： a b1 e '2 b a1 e 2
5 Z轴与地球自转轴重合，X轴与地球赤道
8 9
道面的夹角叫做地
10 心纬度。该点的大
地经度L与地心纬度
12/9 构成地心纬度坐标
0 系。
ppt课件
安徽理工大学
12
地球椭球与椭球计算理论
安徽理工大学
(4)地心空间直角坐标系
地心空间直角坐标系是在大地体内建立
1 2
的坐标系O-XYZ．，它的原点与地球质
3 4
心重合，坐标轴系的配置方法如图所示。
0
返回本章首页
ppt课件
3
地球椭球与椭球计算理论
安徽理工大学
7.1地球椭球的基本几何参数及相互关系
1 2
1.地球椭球的基本几何参数
3 4
五个基本几何参数
a、b称为长度元素
5
椭圆的长半轴： a
6
椭圆的短半轴： b
扁率反映了椭球体的
7
椭圆的扁率：
扁平程度
8 9
ab
a
10
(椭圆的第一偏心率：
e a2 b2 a
椭圆的面积和周长公式
椭圆的面积和周长公式公式是数学题目的解题关键，那么椭圆的面积公式和周长公式是什么呢？下面是由小编为大家整理的“椭圆的面积和周长公式”，仅供参考，欢迎大家阅读。
椭圆的面积公式和周长公式椭圆周长公式：L=2πb+4（a-b）椭圆周长定理：椭圆的周长等于该椭圆短半轴长为半径的圆周长（2πb）加上四倍的该椭圆长半轴长（a）与短半轴长（b）的差。
椭圆面积公式：S=π（圆周率）×a×b，其中a、b分别是椭圆的半长轴，半短轴的长。
椭圆面积公式属于几何数学领域。
拓展阅读：椭圆的定义平面内与两个定点F1、F2的距离之和等于常数（大于|F1F2|）的点的轨迹叫做椭圆，这两个定点叫做椭圆的焦点，两焦点间的距离叫做焦距.。
注意：定义中椭圆标准方程为（x*x）/（a*a）+（y*y）/（b*b）=1 （a>b>0）a*a=b*b+c*c 离心率e=c/a 椭圆顶点（-a，0）（a，0）（0，b）（0-b）2a为长轴长2b为短轴长准线方程x=（a*a）/c。
性质：有对称性。
椭圆的简单几何性质可以总结为以下几种：（一）、对性质的考查：1、范围.2、对称性.3、顶点.4、离心率。
（二）、课本例题的变形考查：1、近日点、远日点的概念：椭圆上任意一点P（x，y）到椭圆一焦点距离的最大值：a+c与最小值：a-c及取最值时点P的坐标。
2、椭圆的第二定义及其应用；椭圆的准线方程及两准线间的距离、焦准距：焦半径公式。
3、已知椭圆内一点M，在椭圆上求一点P，使点P到点M与到椭圆准线的距离的和最小的求法。
4、椭圆的参数方程及椭圆的离心角：椭圆的参数方程的简单应用。
5、直线与椭圆的位置关系，直线与椭圆相交时的弦长及弦中点问题。
椭圆的面积公式是什么怎么算的
椭圆的⾯积公式是什么怎么算的椭圆⾯积公式：S=π(圆周率)×a×b，其中a、b分别是椭圆的长半轴，短半轴的长。
椭圆⾯积公式属于⼏何数学领域。
椭圆的⾯积公式怎么算点与椭圆点M（x0，y0）椭圆x²/a²+y²/b²=1；点在圆内：x0²/a²+y0²/b²<1；点在圆上：x0²/a²+y0²/b²=1；点在圆外：x0²/a²+y0²/b²>1；跟圆与直线的位置关系⼀样的：相交、相离、相切。
直线与椭圆y=kx+m①x²/a+y²/b²=1②由①②可推出x²/a²+(kx+m)²/b²=1相切△=0相离△<0⽆交点相交△>0可利⽤弦长公式：设A(x1,y1)B(x2,y2)求中点坐标根据韦达定理x1+x2=-b/a,x1*x2=c/a带⼊直线⽅程可求出y+y/2=可求出中点坐标。
|AB|=d=√(1+k²)[(x1+x2)²-4x1*x2]=√(1+1/k²)[(y1+y2)²-4x1*x2]椭圆⾯积⽤定积分怎么算椭圆⾯积⽤定积分算为S=abπ。
解题思路：设椭圆x^2/a^2+y^2/b^2=1取第⼀象限内⾯积有 y^2=b^2-b^2/a^2*x^2即 y=√(b^2-b^2/a^2*x^2)=b/a*√(a^2-x^2)由于该式反导数为所求⾯积，观察到原式为圆⽅程公式*a/b，根据(af(x))'=a*f'(x),且x=a时圆⾯积为a^2π/4可得当x=a时，1/4S=b/a*1/4*a^2*π=abπ/4即S=abπ。
第三次全国土地调查调查界线和控制面积_201805013
（四）任意图斑椭球面积计算
2. 计算方法
任意封闭区域总是可以分割成有限个梯形图块，因此，任
意封闭区域的面积
式中Si为按多边形边线
上的坐标点分割的梯形图块面积（i=1,2,…n）用公式
（D.1）计算。
（四）任意图斑椭球面积计算
求封闭多边形 区域（右图） ABCD的面积， 具体方法为：
椭球面上任意多边形计算面积
《第三次全国土地调查技术规程》中附录D《图幅 理论面积与图斑椭球面积计算公式及要求》对图 幅理论面积与图斑椭球面积计算公式、方法及要 求作了明确规定。
（一）椭球面上任一梯形图块面积计算公式
（一）椭球面上任一梯形图块面积计算公式
（二）图幅理论面积计算公式
椭球面上不同比例尺标准图幅的图形是一系列经差、纬差 不同的梯形。因此标准图幅理论面积计算公式可由公式 （D.1）推导为以公式。
各省县级控制界线制作完成后应提前上报全国土 地调查办审查备案。
（二）调查界线的调整
对照影像如果有需要进行零米线更新的，零米线 更新成果应以省为单位单独上报，提交的材料包 括：
——零米线更新情况的省级报告（包括零米线变 更说明、并附相关批准文件）；
——零米线更新所涉及的标准分幅图幅矢量数据 （线层、面层数据）；
县级及县级以上行政区域界线采用最新的全国陆地行政 区域勘界成果确定的界线。
乡镇级行政区域界线，采用各县（市、区）最新确定的 界线。
民政部门行政勘界形成的成果主要有：勘界过程中形成 的具有法律效力的各级行政区域界线协议书及附图、界 桩成果表、界桩登记表、界桩照片和文字记录等。
（一）调查界线的来源及使用
原则上一个土地调查单位只有一个控制范围，不得以权 属界线代替行政界线，形成一县多区域。
三调计算椭球面积
三调计算椭球面积椭球是一种具有特定形状的球体，它在数学和几何学中具有广泛的应用。
计算椭球面积的方法主要有数学推导和数值计算两种，下面将介绍这两种方法。
数学推导方法：首先，我们需要了解椭球的基本概念。
椭球是由一个旋转椭圆绕其短轴旋转一周形成的。
它由两个半轴（长半轴a和短半轴b）以及两个焦点确定。
椭球的体积可以通过长半轴、短半轴和焦距之间的关系来计算。
在数学中，椭球的面积可以通过椭球的参数方程来计算。
椭球的参数方程如下：x = a*cosθ*sinφy = b*sinθ*sinφz = c*cosφ其中，a、b、c分别为三个轴的半径，θ和φ分别为球面上的两个参数，取值范围分别为[0,2π]和[0,π]。
在计算椭球面积时，我们可以采用积分的方法。
首先，我们可以将椭球面分成许多小的面约元，然后通过对这些小面积的积分求和来得到整个椭球的面积。
对于椭球上的一个小面积约元，可以表示为：dS = a*b*sinφ*dθ*dφ其中，dθ和dφ分别表示θ和φ的微小变化量。
对于整个椭球的面积，可以通过以下积分来计算：A = ∫∫dS = ∫∫a*b*sinφ*dθ*dφ其中，积分范围为θ从0到2π，φ从0到π。
将上述积分化简可以得到：A = 4πab这就是椭球的表面积公式。
根据这个公式，我们可以计算任意椭球的面积。
数值计算方法：对于数值计算方法，我们可以使用数值积分的方法来计算椭球的面积。
数值积分是一种近似计算积分的方法，根据给定的精度要求，将积分区间分成若干小区间，然后计算每个小区间上函数值的平均值，并将其乘以小区间的宽度得到近似的积分值。
在计算椭球的面积时，我们可以将其分成许多小的面积元，然后对每个小面积元进行计算，最后将所有小面积元的面积相加得到椭球的表面积。
具体的计算步骤如下：1.选取适当的精度要求和分割小区间的个数；2.将椭球的参数方程代入小面积元的面积表达式中，得到面积元的表达式；3.将小面积元的表达式代入数值积分公式中，计算每个小面积元的面积；4.将所有小面积元的面积相加，得到椭球的表面积。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
图幅理论面积与图斑椭球面积计算公式及要求一、 图幅理论面积计算公式⎢⎣⎡-+---⨯∆=m12m 12m 122cos5(25Csincos3(23Bsin cos (21Asin 603604P B B B B B B B B B L πb ）））⎥⎦⎤-+--m 12m 12cos9(29Esincos7(27DsinB B B B B B ）） （1）式中：a —椭球长半轴(单位：米)，α—椭球扁率，b —椭球短半轴(单位：米)。
е²﹦(a ²﹣b ²)/a ²。
A ﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8。
B ﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8。
C ﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8。
D ﹦ （1/112）е6﹢ （45/2304）е8。
E ﹦ （5/2304）е8。
ΔL —图幅东西图廓的经差（单位：弧度）。
(B 2﹣B 1)—图幅南北图廓的纬差（单位：弧度），Bm ﹦（B 1﹢B 2）/2。
二、椭球面上任意梯形面积计算公式⎢⎣⎡-+---∆=m12m 12m 122cos5(25Csin cos3(23Bsin cos (21Asin 2S B B B B B B B B B L b ）））⎥⎦⎤-+--m 12m 12cos9(29Esincos7(27DsinB B B B B B ）） （2）其中：A,B,C,D,E 为常数，按下式计算: е²﹦(a ²﹣b ²)/a ²A ﹦1﹢（3/6）е²﹢（30/80）е4﹢（35/112）е6﹢（630/2304）е8B ﹦ （1/6）е²﹢（15/80）е4﹢（21/112）е6﹢（420/2304）е8C ﹦ （3/80）е4﹢ （7/112）е6﹢（180/2304）е8D ﹦ （1/112）е6﹢（45/2304）е8E ﹦ （5/2304）е8式中：a —椭球长半轴(单位：米)，b —椭球短半轴(单位：米)；ΔL —图块经差(单位：弧度)； (B 2﹣B 1)—图块纬差(单位：弧度) Bm ﹦（B 1﹢B 2）/2。
三、高斯投影反解变换（x y B L →）模型5000001000000y y '=--⨯带号 （若坐标不带带号，则不需减去带号×1000000；） 0E K x =)sin sin sin sin (cos 7453321E K E K E K E K E E B f -+-+=()()()()()624242222222'4590617201'935241'21⎪⎭⎫⎝⎛++-⎪⎭⎫ ⎝⎛-+++⎪⎭⎫ ⎝⎛-=N y t V t t N y t V t t N y t V B B fηη()()522222322'cos 186242851201'cos 12161'cos 1⎪⎭⎫ ⎝⎛⎪⎪⎭⎫⎝⎛+++++⎪⎭⎫ ⎝⎛⎪⎪⎭⎫⎝⎛++-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫⎝⎛=N y Bt t t N y Bt N y BL f f f ηηη +中央子午线经度值（孤度） （3） 式中：f tgB t = f B e 222cos '=η V C N /= b a C /2= 21η+=V。
为与椭球常数有关的量43210,,,,K K K K K 公式说明：若坐标为没有带号前缀格式，则不需减去带号×1000000；若坐标为有带号前缀格式，则需减去带号×1000000。
四、计算用到的常数、椭球参数在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算。
常数：π﹦3.14159265358979 =ρ 206264.806247180椭球常数：a 椭球长半轴 = 6378140 α椭球扁率 = 1/ 298.257b 椭球短半轴 = 6356755.29 2e椭球第一偏心率= 6.69438499958795E-032e '椭球第二偏心率= 6.73950181947292E-03c 极点子午圈曲率半径= 6399596.65198801相关常数:k 0 = 1.57048687472752E-07 k 1 = 5.05250559291393E-03 k 2 = 2.98473350966158E-05 k 3 = 2.41627215981336E-07 k 4 = 2.22241909461273E-09五、计算中的取位及要求① 高斯投影反解变换后的B,L 以秒为单位，保留到小数点后6位，四舍五入。
② 采用计算机计算时，所有变量数据类型均要定义为双精度。
③ 面积计算结果以平方米为单位，保留一位小数，四舍五入。
④ 各种比例尺标准分幅图经差、纬差见表1。
⑤ 在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。
表1 各种比例尺标准分幅图经差、纬差表六、任意图斑椭球面积计算方法任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。
1、计算方法：任意封闭区域总是可以分割成有限个任意小的梯形图块，因此，任意封闭区域的面积∑==ni isP 1，式中Si 为分割的任意小的梯形图块面积（i=1,2,…n ）用公式(2)计算。
求封闭区域（多边形如图1）ABCD 的面积 ，其具体方法为：（1）对封闭区域（多边形）的界址点连续编号（顺时针或逆时针）ABCD，提取各界址点的高斯平面坐标A(X1,Y1)，B(X2,Y2)，C(X3,Y3)，D(X4,Y4)；（2）利用高斯投影反解变换模型公式（3），将高斯平面坐标换算为相应椭球的大地坐标A(B1,L1)，B(B2,L2)，C(B3,L3)，D(B4,L4)；（3）任意给定一经线L0（如L0＝60°），这样多边形ABCD的各边AB、BC、CD、DA 与L0就围成了4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）；（4）由于在椭球面上同一经差随着纬度升高，梯形图块的面积逐渐减小，而同一纬差上经差梯形图块的面积相等，所以，将梯形图块ABB1A1按纬差分割成许多个小梯形图块AEiFiA1，用公式（2）计算出各小梯形图块AEiFiA1的面积Si，然后累加Si就得到梯形图块ABB1A1的面积，同理，依次计算出梯形图块BCC1B1、CDD1C1、DAA1D1的面积（注：用公式（2）计算面积时，B1、B2分别取沿界址点编号方向的前一个、后一个界址点的大地纬度，ΔL为沿界址点编号方向的前一个、后一个界址点的大地经度的平均值与L0的差）；（5）多边形ABCD的面积就等于4个梯形图块（ABB1A1、BCC1B1、CDD1C1、DAA1D1）面积的代数和。
B图1 椭球面上任意多边形计算面积则任意多边形ABCD的面积P为：P=ABCD= BCC1B1+ CDD1C1+ DAA1D1- ABB1A12、计算要求①利用图形坐标点将高斯坐标系下的几何图形反算投影到大地坐标系，进行投影变换。
②任意指定一条经线L0，从选定多边形几何形状的起始点开始，沿顺时针方向依次计算相邻两点构成的线段，以及两点到指定经线的平行线构成的梯形面积。
将该梯形沿纬度变化方向（Y轴）进行切割，至少需切割为2个部分。
③计算过程中应顺同一方向依坐标点逐个计算相邻两点连线与任意经线构成的梯形面积，坐标点不得有遗漏。
若多边形包含内多边形（洞），则该多边形面积为外多边形面积减去所有内多边形面积之和。
④计算所有梯形面积的代数和即为该多边形的面积。
七、算法伪代码描述为了确保编程使用的参数、算法一致，保证不同软件计算的椭球面积一致，我们用算法伪代码描述的方法对编程进行统一，在利用计算机编制椭球面积计算软件时，计算参数与计算顺序应严格按照以下代码执行。
1、参数说明双精度类型：圆周率值：PI = 3.14159265358979中央经线:CenterLRHO = 206264.8062471A:ParamAB:ParamBC:ParamCD:ParamDE:ParamEConst ZERO As Double = 0.00000000000180椭球常数椭球长半轴：aRadius = 6378140椭球短半轴：bRadius = 6356755.29椭球扁率：ParaAF = 1/ 298.257椭球第一偏心率：ParaE1 = 6.69438499958795E-03椭球第二偏心率：ParaE2 = 6.73950181947292E-03极点子午圈曲率半径：ParaC = 6399596.65198801k0:Parak0 = 1.57048687472752E-07k1:Parak1 = 5.05250559291393E-03k2:Parak2 = 2.98473350966158E-05k3:Parak3 = 2.41627215981336E-07k4:Parak4 = 2.22241909461273E-092、算法描述初始化参数Double e;Double a;e = ParaE2;ParaC = aRadius / (1 - ParaAF);ParamA = 1 + (3 / 6) * e + (30 / 80) * Power(e, 2) + (35 / 112) * Power(e, 3) + (630 / 2304) * Power(e, 4);ParamB = (1 / 6) * e + (15 / 80) * Power(e, 2) + (21 / 112) * Power(e, 3) + (420 / 2304) * Power(e, 4);ParamC = (3 / 80) * Power(e, 2) + (7 / 112) * Power(e, 3) + (180 / 2304) * Power(e, 4);ParamD = (1 / 112) * Power(e, 3) + (45 / 2304) * Power(e, 4);ParamE = (5 / 2304) * Power(e, 4);参数初始化结束中央经线转换为弧度CenterL = TransDegreeToArc(CenterL)选定本初子午线为参考经线StandardLat = 0For 起始点 To 倒数第二点由高斯坐标反解计算经纬度值ComputeXYGeo (PntColl.Point(i).y, PntColl.Point(i).x, B, L, CenterL)ComputeXYGeo (PntColl.Point(i + 1).y, PntColl.Point(i + 1).x, B1, L1, CenterL) 将经纬度转换为弧度值B = B / RHOL = L / RHOB1 = B1 / RHOL1 = L1 / RHO计算梯形面积Double AreaVal;//梯形面积值Double lDiference ;//经差Double bDiference; //纬差Double bSum;//纬度和Double ItemValue(5);//计算变量bDiference = (B1 - B0)；bSum = (B1 + B0) / 2；lDiference = (L1 + L) / 2;ItemValue(0) = ParamA * Sin(bDiference / 2) * Cos(bSum);ItemValue(1) = ParamB * Sin(3 * bDiference / 2) * Cos(3 * bSum);ItemValue(2) = ParamC * Sin(5 * bDiference / 2) * Cos(5 * bSum);ItemValue(3) = ParamD * Sin(7 * bDiference / 2) * Cos(7 * bSum);ItemValue(4) = ParamE * Sin(9 * bDiference / 2) * Cos(9 * bSum);AreaVal = 2 * bRadius * lDiference * bRadius * (ItemValue(0) - ItemValue(1) + ItemValue(2) - ItemValue(3) + ItemValue(4));areaSum = areaSum + AreaVal;NextEnd Sub3、高斯坐标反解算法Public Sub ComputeXYGeo(x As Double, y As Double, B As Double, L As Double, center As Double)Dim y1 As DoubleDim bf As Doubley1 = y - 500000Dim e As Doublee = Parak0 * xDim se As Doublese = Sin(e)bf = e + Cos(e) * (Parak1 * se - Parak2 * Power(se, 3) + Parak3 * Power(se, 5) - Parak4 * Power(se, 7))Dim v As DoubleDim t As DoubleDim N As DoubleDim nl As DoubleDim vt As DoubleDim yn As DoubleDim t2 As DoubleDim g As Doubleg = 1t = Tan(bf)nl = ParaE1 * Power(Cos(bf), 2)v = Sqr(1 + nl)N = ParaC / vyn = y1 / Nvt = Power(v, 2) * tt2 = Power(t, 2)B = bf - vt * Power(yn, 2) / 2 + (5 + 3 * t2 + nl - 9 * nl * t2) * vt * Power(yn, 4) / 24 - (61 + 90 * t2 + 45 * Power(t2, 2)) * vt * Power(yn, 6) / 720B = TransArcToDegree(B)Dim cbf As Doublecbf = 1 / Cos(bf)L = cbf * yn - (1 + 2 * t2 + nl) * cbf * Power(yn, 3) / 6 + (5 + 28 * t2 + 24 * Power(t2, 2) + 6 * nl + 8 * nl * t2) * cbf * Power(yn, 5) / 120 + centerL = TransArcToDegree(L)End Sub弧度转换为度Public Function TransArcToDegree(arc As Double) As DoubleDim degree As DoubleDim min As DoubleDim sec As DoubleDim ret As DoubleDim tmp As Doubleret = arc * 180 / PIdegree = FormatValue(ret, 100, 100)tmp = (ret - degree) * 60min = FormatValue(tmp, 100, 100)sec = (tmp - min) * 60//秒保留到小数点后6位，四舍五入sec = Format(sec, "####.000000") 'FormatValue(sec, 10000000, 100)TransArcToDegree = degree * 3600 + min * 60 + secEnd FunctionPrivate Function FormatValue(inputVal As Double, precsion As Long, scaleNum As Long) As Double FormatValue = (Int(inputVal * precsion) - Int(inputVal * precsion) Mod scaleNum) / precsion End Function。