国土资源部


关于统一图幅理论面积与图斑椭球面积计算要求的通知

（国土调查办发〔2008〕32号）


各省、自治区、直辖市第二次土地调查领导小组办公室，国土资源厅（国土环境资源厅、国土资源局、国土资源和房屋管理局、房屋土地资源管理局），解放军土地管理局、新疆生产建设兵团国土资源局：


面积计算是第二次土地调查的一项重要内容，国务院第二次全国土地调查领导小组办公室组织有关专家，依据《第二次全国土地调查技术规程》，对图幅理论面积与图斑椭球面积计算公式进行了细化，明确了面积计算方法，统一了公式中的有关参数，现将《图幅理论面积与图斑椭球面积计算公式及要求》予以印发，请各地严格遵照执行。



附：图幅理论面积与图斑椭球面积计算公式及要求



二〇〇八年三月二十七日

图幅理论面积与图斑椭球面积计算公式及要求

一、 图幅理论面积计算公式
（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—图幅东西图廓的经差（单位：弧度）。
(B2﹣B1)—图幅南北图廓的纬差（单位：弧度），Bm﹦（B1﹢B2）/2。

二、椭球面上任意梯形面积计算公式
（2）
其中：A,B,C,D,E 为常数，按下式计算:
е²﹦(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
式中：a—椭球长半轴(单位：米)，b—椭球短半轴(单位：米)；
ΔL—图块经差(单位：弧度)； (B2﹣B1)—图块纬差(单位：弧度)
Bm﹦（B1﹢B2）/2。

三、高斯投影反解变换（ ）模型
（若坐标不带带号，则不需减去带号×1000000；）



+中央子午线经度值（孤度） （3）
式中：

公式说明：若坐标为没有带号前缀格式，则不需减去带号×1000000；若坐标为有带号前缀格式，则需减去带号×1000000。

四、计算用到的常数、椭球参数
在计算图幅理论面积与任意图斑椭球面积时，有关常数及保留的位数按给定数值计算。
常数：
π﹦3.14159265358979
206264.8062471
80椭球常数：
= 6378140 = 1/ 298.257
= 6356755.29
= 6.69438499958795E-03
= 6.73950181947292E-03
= 6399596.65198801
相关常数:
k0 = 1.57048687472752E-07
k1 = 5.05250559291393E-03
k2 = 2.98473350966158E-05
k3 = 2.41627215981336E-07
k4 = 2.22241909461273E-09

五、计算中的取位及要求
① 高斯投影反解变换后的B,L以秒为单位，保留到小数点后6位，四舍五入。
② 采用计算机计算时，所有变量数据类型均要定义为双精度。
③ 面积计算结果以平方米为单位，保留一位小数，四舍五入。
④ 各种比例尺标准分幅图经差、纬差见表1。
⑤ 在用大地坐标生成标准分幅图框时，要求在每条边框线的整秒处插入加密点。
表1 各种比例尺标准分幅图经差、纬差表
比例尺 1:100万 1:50万 1:25万 1:10万 1:5万 1:2.5万 1:1万 1:5千
经差 6º 3º 1º30′ 30′ 15′ 7′30″ 3′45″ 1′52.5″
纬差 4º 2º 1º 20′ 10′ 5′ 2′30″ 1′15″

六、任意图斑椭球面积计算方法
任意封闭图斑椭球面积计算的原理：将任意封闭图斑高斯平面坐标利用高斯投影反解变换模型，将高斯平面坐标换算为相应椭球的大地坐标，再利用椭球面上任意梯形图块面积计算模型计算其椭球面积，从而得到任意封闭图斑的椭球面积。
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）面积的代数和。

图1 椭球面上任意多边形计算面积
则任意多边形ABCD的面积P为：
P=ABCD= BCC1B1+ CDD1C1+ DAA1D1- ABB1A1
2、计算要求
① 利用图形坐标点将高斯坐标系下的几何图形反算投影到大地坐标系，进行投影变换。
② 任意指定一条经线L0，从选定多边形几何形状的起始点开始，沿顺时针方向依次计算相邻两点构成的线段，以及两点到指定经线的平行线构成的梯形面积。将该梯形沿纬度变化方向（Y轴）进行切割，至少需切割为2个部分。
③ 计算过程中应顺同一方向依坐标点逐个计算相邻两点连线与任意经线构成的梯形面积，坐标点不得有遗漏。若多边形包含内多边形（洞），则该多边形面积为外多边形面积减去所有内多边形面积之和。
④ 计算所有梯形面积的代数和即为该多边形的面积。
七、算法伪代码描述
为了确保编程使用的参数、算法一致，保证不同软件计算的椭球面积一致，我们用算法伪代码描述的方法对编程进行统一，在利用计算机编制椭球面积计算软件时，计算参数与计算顺序应严格按照以下代码执行。
1、参数说明
双精度类型：
圆周率值：PI = 3.14159265358979
中央经线:CenterL
RHO = 206264.8062471

A:ParamA
B:ParamB
C:ParamC
D:ParamD
E:ParamE

Const ZERO As Double = 0.000000000001

80椭球常数

椭球长半轴：aRadius = 6378140
椭球短半轴：bRadius = 6356755.29
椭球扁率：ParaAF = 1/ 298.257
椭球第一偏心率：ParaE1 = 6.69438499958795E-03
椭球第二偏心率：ParaE2 = 6.73950181947292E-03
极点子午圈曲率半径：ParaC = 6399596.65198801

k0:Parak0 = 1.57048687472752E-07
k1:Parak1 = 5.05250559291393E-03
k2:Parak2 = 2.98473350966158E-05
k3:Parak3 = 2.41627215981336E-07
k4:Parak4 = 2.22241909461273E-09

2、算法描述

初始化参数

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 = 0

For 起始点 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 / RHO
L = L / RHO
B1 = B1 / RHO
L1 = 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;
Next

End Sub

3、高斯坐标反解算法

Public Sub ComputeXYGeo(x As Double, y As Double, B As Double, L As Double, center As Double)

Dim y1 As Double
Dim bf As Double

y1 = y - 500000

Dim e As Double

e = Parak0 * x

Dim se As Double

se = Sin(e)
bf = e + Cos(e) * (Parak1 * se - Parak2 * Power(se, 3) + Parak3 * Power(se, 5) - Parak4 * Power(se, 7))

Dim v As Double
Dim t As Double
Dim N As Double
Dim nl As Double
Dim vt As Double
Dim yn As Double
Dim t2 As Double
Dim g As Double

g = 1

t = Tan(bf)
nl = ParaE1 * Power(Cos(bf), 2)
v = Sqr(1 + nl)
N = ParaC / v
yn = y1 / N
vt = Power(v, 2) * t
t2 = 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) / 720

B = TransArcToDegree(B)

Dim cbf As Double

cbf = 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 + center
L = TransArcToDegree(L)
End Sub

弧度转换为度
Public Function TransArcToDegree(arc As Double) As Double
Dim degree As Double
Dim min As Double
Dim sec As Double
Dim ret As Double
Dim tmp As Double
ret = arc * 180 / PI
degree = FormatValue(ret, 100, 100)
tmp = (ret - degree) * 60
min = FormatValue(tmp, 100, 100)
sec = (tmp - min) * 60
//秒保留到小数点后6位，四舍五入
sec = Format(sec, "####.000000") 'FormatValue(sec, 10000000, 100)
TransArcToDegree = degree * 3600 + min * 60 + sec
End Function

Private 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
国务院办公厅


国务院办公厅转发发展改革委财政部卫生部关于清理化解基层医疗卫生机构债务意见的通知

国办发〔2011〕32号


各省、自治区、直辖市人民政府，国务院各部委、各直属机构：
　　发展改革委、财政部、卫生部《关于清理化解基层医疗卫生机构债务的意见》已经国务院同意，现转发给你们，请认真贯彻执行。
　　　　　　　　　　　　　　　　　　　　　　　　　　　 　国务院办公厅
　　　　　　　　　　　　　　　　　　　　　　　　　　 二○一一年七月五日



关于清理化解基层医疗卫生机构债务的意见
发展改革委　财政部　卫生部

　　随着基本药物制度在基层全面实施，“以药补医”机制逐步扭转，基层医疗卫生机构以往形成的债务问题进一步显现，影响到基层医疗卫生机构正常运行。为确保基本药物制度顺利实施，促进基层医疗卫生机构建立新的运行机制和持续健康发展，根据深化医药卫生体制改革和清理化解乡村债务工作的有关要求，现就清理化解基层医疗卫生机构债务提出以下意见：
　　一、工作目标
　　各省（区、市）要按照“制止新债、锁定旧债、明确责任、分类处理、逐步化解”的总体要求，在严格制止发生新债的基础上，用2年左右时间全面完成基层医疗卫生机构长期债务的清理化解工作。
　　二、基本原则
　　（一）谁举债谁负责。各省（区、市）人民政府对本省（区、市）基层医疗卫生机构债务清理化解工作负总责，县级人民政府具体实施。中央财政按照重点支持中西部地区、适当兼顾东部地区的原则，对地方债务清理化解予以适当补助。
　　（二）先清理后化解。在全面摸清基层医疗卫生机构债务底数的基础上锁定债务，由各地结合实际情况确定化债的先后顺序，逐步化解。
　　（三）先承诺后补助。由国务院医改办公室与各省（区、市）医改领导小组签订责任书，各省（区、市）承诺在2年左右时间内全部完成基层医疗卫生机构债务化解任务。中央财政对在规定期限内完成债务化解任务的地方给予补助，国务院医改办公室对各地债务化解工作统一组织考核，未在规定期限内完成债务化解任务的扣回中央财政相应补助资金。
　　三、债务化解范围
　　纳入本次债务化解范围的基层医疗卫生机构是指由政府举办的乡镇卫生院和社区卫生服务机构。纳入化解范围的债务是基层医疗卫生机构发展建设过程中形成的长期债务，主要包括发生于业务用房、辅助用房建设维修和医疗设备购置等与基层医疗卫生机构发展建设直接相关的债务。债务计算时间原则上截止到2009年12月31日。2010年1月1日至本意见印发之日形成的债务，各地参照本意见进行化解，不纳入中央财政补助计算范围。
　　四、偿债资金来源
　　各地区可从以下渠道筹集偿债资金：一是统筹安排地方一般预算收入、上级财力性转移支付资金；二是中央财政安排的基层医疗卫生机构实施基本药物制度“以奖代补”资金和用于支持化解债务的专项补助资金；三是基层医疗卫生机构实施基本药物制度前的收支结余资金；四是从基层医疗卫生机构在“核定任务、核定收支”后超收的资金中安排一定的比例用于偿债；五是通过统筹有关非税收入途径筹集资金；六是社会捐资赞助的偿债资金。鼓励有条件的乡镇积极筹资用于偿债。各地要积极调整财政支出结构，落实偿债资金，省、地（市）两级财政要加大对财政困难县（市、区）化解债务工作的资金支持力度。中央财政以全国卫生财务年报反映的截至2009年12月31日的基层医疗卫生机构长期负债为基数确定补助总额，根据人口、财力等因素对各地进行适当补助，具体补助办法由财政部另行制定。
　　五、工作要求
　　（一）摸清底数，锁定债务。地方政府要统一组织审计、财政、卫生、监察等部门对每一个基层医疗卫生机构的每项、每笔债务认真清理核实，剔除不实债务，锁定实际债务。要坚持程序公开、过程公开和结果公开，在一定范围内公示各项债务情况，接受监督。各省（区、市）有关部门要对上报的基层医疗卫生机构债务进行审核和认定，确保债务数据合法、真实、完整和准确。
　　（二）明确化债主体，分类化解债务。按照举债主体、债务来源、债务用途对债务进行分类，严格划分县乡政府、卫生行政部门和基层医疗卫生机构等的责任，将审核认定的全部债务从基层医疗卫生机构剥离出来，交给当地政府。各级财政要在预算中特设专户，单独列支用于化解基层医疗卫生机构债务的支出，相关资金不作为预算安排正常卫生支出的基数。要区分轻重缓急，明确偿债次序，分类逐步化解债务。要优先化解已实施基本药物制度的基层医疗卫生机构的债务，优先偿还医务人员集资等个人的债务。
　　（三）落实政府投入责任，坚决制止发生新债。各地区、各有关部门要严格按照《国务院办公厅关于建立健全基层医疗卫生机构补偿机制的意见》（国办发〔2010〕62号）的要求，落实政府办基层医疗卫生机构的基本建设、设备购置、经常性收支差额补助等各项经费，不留经费缺口。基层医疗卫生机构建设项目和设备购置要按程序申报，经批准后实施，所需资金由政府纳入财政预算足额安排。未经批准、资金未落实的项目一律不得实施，经批准的项目在实施过程中不得随意扩大建设规模、提高建设标准。所有政府办基层医疗卫生机构都要认真执行国家有关规定，不得举借新债。各地要加强源头控制，确保基层医疗卫生机构建设与政府财力水平相适应，不得将应由地方政府承担的资金转给基层医疗卫生机构承担。
　　（四）加强组织领导，制定实施方案。各地区、各有关部门要充分认识化解基层医疗卫生机构债务的重要性和紧迫性，加强组织领导，明确责任分工，精心组织实施，强化协调配合，及时研究解决化解债务过程中出现的问题，确保各项工作平稳有序推进。各级医改领导小组要充分发挥统筹协调作用，发展改革、财政、卫生部门要会同审计、监察等部门按照各自职责，做好债务清理、核实、锁定和资金筹集等工作。在本意见印发60个工作日内，各省（区、市）要按照本意见要求，结合实际情况制定本省（区、市）清理化解基层医疗卫生机构债务的实施方案，并报国务院医改办公室、财政部、卫生部备案。2011年10月31日前各地要完成债务清理核实工作，按照债务来源和用途逐笔登记造册，建立债务台账和债权债务数据库；2011年12月31日前要完成债务认定工作，并将认定的债务从基层医疗卫生机构剥离给当地政府；当地政府要按债务协议，如期偿还结清各项债务。
　　（五）严格执行有关规定和财经纪律，建立责任追究制度。各地不得借新债还旧债，不得向群众摊派，不得挤占挪用其他医改专项资金，不得影响基层医疗卫生机构实施基本药物制度和推进综合改革。要将基层医疗卫生机构债务清理化解和不得举借新债落实情况作为对领导干部任期经济责任审计的重要内容。发展改革、财政、卫生等部门要密切配合，建立事前、事中、事后相结合的监督体系和制约机制。对违反规定搞建设、上项目、借新债的基层医疗卫生机构和部门，视情节轻重和数额大小，依法依规追究有关单位负责人和相关人员的行政责任。对虚报冒领、截留挪用、套取补助资金及造成资金损失的，除追回补助资金外，还要按照有关规定追究直接责任人和有关负责人的责任；涉嫌违法犯罪的，移交司法机关处理。