国家税务总局关于中国粮油食品进出口(集团)有限公司重组改制有关契税问题的通知
作者:法律资料网 时间:2024-07-22 06:52:17 浏览:9417
来源:法律资料网
下载地址: 点击此处下载
国家税务总局关于中国粮油食品进出口(集团)有限公司重组改制有关契税问题的通知
国家税务总局关于中国粮油食品进出口(集团)有限公司重组改制有关契税问题的通知
国税函[2001]843号
2001-11-16国家税务总局
河北、黑龙江、江苏、浙江、福建、江西、山东、河南、湖北、新疆、广东省(自治区)财政厅，北京、山西、辽宁、上海、海南、云南省(直辖市)地方税务局:
根据国务院领导指示，中国粮油食品进出口(集团)有限公司(以下简称中粮公司)拟对下属部分企业进行重组改制。现就该公司重组改制过程中所涉及的契税问题通知如下:
一、关于吸收合并。中粮公司拟将部分下属企业进行吸收合并，被吸收企业不再保留法人地位。如合并各方均为中粮公司的全资子公司，则不征契税，其余征收契税。
二、关于股权划转。中粮公司拟将部分境内下属企业的股权划转至该公司在香港注册的全资子公司，根据国家对外贸易经济合作部的批复，被划转企业性质将变更为中外合资或外商独资，企业名称相应变更。由于股权划转不影响企业法人的存续，其土地、房屋权属不发生转移，故不征契税。
国家税务总局
二○○一年十一月十六日
下载地址: 点击此处下载
民事司法调解实务调研报告
周鹏龙
日益纷繁复杂，层出不穷的纠纷决定了要有与之相适应的多元化纠纷解决机制，司法调解作为多元化纠纷解决机制之重要方式之一，在定纷止争，促进社会和谐，维护社会正义、公平方面日益发挥着不可代替的优越性，尤其在我国大力构建和谐社会的时代背景下，充分发挥民事司法调解制度的作用和功效显得尤为重要和契合时机。因此在和谐社会大力推进的过程中发挥民事司法调解制度的作用和功能显得尤为重要。尤其面对实务界，理论界对民事司法调解制度的各种质疑和挑战，对诉讼调解制度进行实证研究更具有现实意义。因此本次调研本着从实务实证调查研究出发，在实务调研中了解民事司法调解制度在司法实践中的基本运行情况及其以调解方式办结案件的法律效果和社会效果，并在此基础上提出调解制度的现实意义及其功效，通过调查访问当事人及法官，参加法院审判活动及旁听个案等方法，提出民事司法调解在司法实践中运作中存在的问题并提出改进和完善建议，为民事司法调解制度的发展和完善尽一点微薄之力。
一、本次调研的基本情况。
本调研于2008年7月20日全面开始，分为两个时间段进行。第一阶段为2008年7月20日——2008年8月30日，该阶段主要通过在调研选择试点统计2008年度以调解方式办结案件的数量，记录调解案件案由、当事人地址或联系方式等基本信息，访问各案件承办法官及相关人员，以统计信息为基础通过实地走访当事人进行调查了解调解案件的法律效果及相关情况，并同时在每个环节做好记录。第二阶段为2009年1月20日——2月20日，该阶段在第一阶段的基础上继续在同一调研试点统计尚未统计的以调解办结的案件，记录调解案件案由，当事人地址或联系方式等基本信息，访问各办案法官及相关人员，以统计信息实地走访当事人进行调查了解调解案件的质量及相关情况，并同时做好记录。
二、本次调研的试点。
本次调研课题为《少数民族地区基层人民法院民事司法调解制度运行机制及法律效果调查》，主要针对少数民族地区基层法院的诉讼调解基本情况及法律、社会效果的实证调研。基于此本次调研地点选择在宁夏回族自治区海原县人民法院及其各派出法庭，包括海城镇法庭，李旺法庭和西安法庭三个派出法庭。
三、本调研的基本内容。
本次调研属于民事司法调解实务调研，其内容主要包括：
（1）、海原县人民法院及其各派出法庭（海城镇法庭、李旺法庭、西安法庭）民事司法调解基本运行机制及其特色之处。
（2）、实地走访当事人，调查以调解方式办结案件的法律效果及其基本情况。
（3）、实地访问法官对民事司法调解的态度，经验以及调解感想等。
（4）据以上调研提出民事司法调解在司法实践中的作用和意义等并提出存在的问题，同时给予改进和完善建议。
四、本调研基本构思以及调研报告结构安排
本次调研课题：“少数民族地区民事司法调解制度运行机制及其法律效果调查报告”，主要有两部分组成：（一）少数民族地区民事司法调解制度运行机制;（二）少数民族地区民事司法调解制度法律效果调查。
五、海原县人民法院诉讼调解制度存在之可行性分析
在构建社会主义和谐社会过程中，诉讼调解作为人民法院定争止纷的一种结案方式，通过审判人员在征得当事人同意后就案件达成调解协议，达到案结事了的目的，对促进纠纷的柔性解决、实现社会稳定发展起着举足轻重的作用。
（一）、审判压力过大与审判人员有限的矛盾十分突出，调解制度可以有效缓解案多人少的矛盾。（以海城镇法庭为例）
民事案件相关情况统计表
年 件 收案 结案 办案人数 审判人数 人均收案 人均收案率 人均结案 人均结案率
2006年 190 186 3 3 64 33.3% 62 33.3%
2007年 230 229 4 3 57.5 25% 76.3 57.25%
2008年 281 280 4 3 93 35% 93.3 70%
纵观上表可看出，海原县人民法院海城镇法庭审判人数由2006年3人增加至2008年4人，年均增长率为25%，人均收案率由 2006年33.3%增至2008年35%，人均结案率由2006年33.3%增至2008年70%，据此可见该院审判压力过大与审判人员有限的矛盾突出，而民事司法调解制度可以有效缓解案多人少的矛盾。
（二）、本院受理案件有近六成以调解方式结案。
民事案件调解情况统计表
年度 件 收案 结案 判决 调解 调解率
2006年 190 186 80 106 56.99%
2007年 230 229 105 124 54.15%
2008年 281 280 110 170 60.71%
由本表可知该院2006年至2008年诉讼调解案件由106件增至170件，调解率由56.99%增至60.71%，基本达到六成调解率，符合本院的司法形势，有效以调解方式解决民事案件，取得定纷止争的和谐结果。
（三）、诉讼调解有助于实现社会的公平正义。公平正义是法律永恒追求的基本价值目标，是构建和谐社会的基础。法官作为司法人员，作为公平正义的法律使者，行使着国家公权力，是社会稳定和谐的调节器。而法院诉讼调解的灵魂就在于它的公正性，调解的前提是分清责任、明辨是非。法官在分清责任、明辨是非后，通过运用调解手段居中化解当事人矛盾纠纷，从而让全社会各方面的利益关系得到妥善协调，人民内部矛盾和其他社会矛盾得到正确处理，社会公平和正义得到切实维护和实现。
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
（国土调查办发〔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
周鹏龙
日益纷繁复杂，层出不穷的纠纷决定了要有与之相适应的多元化纠纷解决机制，司法调解作为多元化纠纷解决机制之重要方式之一，在定纷止争，促进社会和谐，维护社会正义、公平方面日益发挥着不可代替的优越性，尤其在我国大力构建和谐社会的时代背景下，充分发挥民事司法调解制度的作用和功效显得尤为重要和契合时机。因此在和谐社会大力推进的过程中发挥民事司法调解制度的作用和功能显得尤为重要。尤其面对实务界，理论界对民事司法调解制度的各种质疑和挑战，对诉讼调解制度进行实证研究更具有现实意义。因此本次调研本着从实务实证调查研究出发，在实务调研中了解民事司法调解制度在司法实践中的基本运行情况及其以调解方式办结案件的法律效果和社会效果，并在此基础上提出调解制度的现实意义及其功效，通过调查访问当事人及法官，参加法院审判活动及旁听个案等方法，提出民事司法调解在司法实践中运作中存在的问题并提出改进和完善建议，为民事司法调解制度的发展和完善尽一点微薄之力。
一、本次调研的基本情况。
本调研于2008年7月20日全面开始，分为两个时间段进行。第一阶段为2008年7月20日——2008年8月30日，该阶段主要通过在调研选择试点统计2008年度以调解方式办结案件的数量，记录调解案件案由、当事人地址或联系方式等基本信息，访问各案件承办法官及相关人员，以统计信息为基础通过实地走访当事人进行调查了解调解案件的法律效果及相关情况，并同时在每个环节做好记录。第二阶段为2009年1月20日——2月20日，该阶段在第一阶段的基础上继续在同一调研试点统计尚未统计的以调解办结的案件，记录调解案件案由，当事人地址或联系方式等基本信息，访问各办案法官及相关人员，以统计信息实地走访当事人进行调查了解调解案件的质量及相关情况，并同时做好记录。
二、本次调研的试点。
本次调研课题为《少数民族地区基层人民法院民事司法调解制度运行机制及法律效果调查》，主要针对少数民族地区基层法院的诉讼调解基本情况及法律、社会效果的实证调研。基于此本次调研地点选择在宁夏回族自治区海原县人民法院及其各派出法庭，包括海城镇法庭，李旺法庭和西安法庭三个派出法庭。
三、本调研的基本内容。
本次调研属于民事司法调解实务调研，其内容主要包括：
（1）、海原县人民法院及其各派出法庭（海城镇法庭、李旺法庭、西安法庭）民事司法调解基本运行机制及其特色之处。
（2）、实地走访当事人，调查以调解方式办结案件的法律效果及其基本情况。
（3）、实地访问法官对民事司法调解的态度，经验以及调解感想等。
（4）据以上调研提出民事司法调解在司法实践中的作用和意义等并提出存在的问题，同时给予改进和完善建议。
四、本调研基本构思以及调研报告结构安排
本次调研课题：“少数民族地区民事司法调解制度运行机制及其法律效果调查报告”，主要有两部分组成：（一）少数民族地区民事司法调解制度运行机制;（二）少数民族地区民事司法调解制度法律效果调查。
五、海原县人民法院诉讼调解制度存在之可行性分析
在构建社会主义和谐社会过程中，诉讼调解作为人民法院定争止纷的一种结案方式，通过审判人员在征得当事人同意后就案件达成调解协议，达到案结事了的目的，对促进纠纷的柔性解决、实现社会稳定发展起着举足轻重的作用。
（一）、审判压力过大与审判人员有限的矛盾十分突出，调解制度可以有效缓解案多人少的矛盾。（以海城镇法庭为例）
民事案件相关情况统计表
年 件 收案 结案 办案人数 审判人数 人均收案 人均收案率 人均结案 人均结案率
2006年 190 186 3 3 64 33.3% 62 33.3%
2007年 230 229 4 3 57.5 25% 76.3 57.25%
2008年 281 280 4 3 93 35% 93.3 70%
纵观上表可看出，海原县人民法院海城镇法庭审判人数由2006年3人增加至2008年4人，年均增长率为25%，人均收案率由 2006年33.3%增至2008年35%，人均结案率由2006年33.3%增至2008年70%，据此可见该院审判压力过大与审判人员有限的矛盾突出，而民事司法调解制度可以有效缓解案多人少的矛盾。
（二）、本院受理案件有近六成以调解方式结案。
民事案件调解情况统计表
年度 件 收案 结案 判决 调解 调解率
2006年 190 186 80 106 56.99%
2007年 230 229 105 124 54.15%
2008年 281 280 110 170 60.71%
由本表可知该院2006年至2008年诉讼调解案件由106件增至170件，调解率由56.99%增至60.71%，基本达到六成调解率，符合本院的司法形势，有效以调解方式解决民事案件，取得定纷止争的和谐结果。
（三）、诉讼调解有助于实现社会的公平正义。公平正义是法律永恒追求的基本价值目标，是构建和谐社会的基础。法官作为司法人员，作为公平正义的法律使者，行使着国家公权力，是社会稳定和谐的调节器。而法院诉讼调解的灵魂就在于它的公正性，调解的前提是分清责任、明辨是非。法官在分清责任、明辨是非后，通过运用调解手段居中化解当事人矛盾纠纷，从而让全社会各方面的利益关系得到妥善协调，人民内部矛盾和其他社会矛盾得到正确处理，社会公平和正义得到切实维护和实现。
关于统一图幅理论面积与图斑椭球面积计算要求的通知
关于统一图幅理论面积与图斑椭球面积计算要求的通知
（国土调查办发〔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
