关于统一图幅理论面积与图斑椭球面积计算要求的通知
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
(国土调查办发〔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
内蒙古自治区实施《中华人民共和国野生动物保护法》办法
内蒙古自治区人大常委会
内蒙古自治区实施《中华人民共和国野生动物保护法》办法
内蒙古自治区人大常委会
(1991年12月24日内蒙古自治区第七届人民代表大会常务委员会第二十四次会议通过 1991年12月24日公布施行)
目 录
第一章 总 则
第二章 野生动物保护
第三章 野生动物管理
第四章 法律责任
第五章 附 则
第一章 总 则
第一条 为保护、发展和合理利用自治区野生动物资源,根据《中华人民共和国野生动物保护法》,结合自治区实际,制定本办法。
第二条 在自治区境内从事野生动物的保护管理、驯养繁殖、开发利用活动,必须遵守本办法。
本办法规定保护的野生动物是指:
(一)国务院批准的国家重点保护的陆生、水生野生动物;
(二)自治区人民政府批准的重点保护的陆生、水生野生动物;
(三)国家和自治区野生动物行政主管部门公布保护的有益或者有重要经济、科学研究价值的陆生野生动物。
第三条 旗县级以上人民政府林业、渔业行政主管部门分别负责本行政区域内的陆生、水生野生动物管理工作。
自治区武装森林警察部队执行保护野生动物资源的任务。
第四条 旗县级以上人民政府应当把保护发展和合理利用野生动物资源纳入国民经济和社会发展计划,加强对野生动物资源的管理,鼓励和支持开展野生动物的科学研究和驯养繁殖。
第五条 自治区人民政府对少数民族专业猎民队的狩猎生产和狩猎枪支的使用给予特许。
第六条 保护野生动物资源是每个机关、团体、部队、企业事业单位和公民的义务。
自治区各级人民政府要加强对野生动物资源保护的宣传和教育。
第七条 自治区各级人民政府对执行野生动物保护法和本办法成绩显著的,给予表彰奖励。
第二章 野生动物保护
第八条 自治区人民政府对本办法规定的野生动物实行重点保护。
重点保护的野生动物名录及其调整,由自治区野生动物行政主管部门提出,报自治区人民政府批准公布。
第九条 自然保护区,风景游览区,文物保护区,城镇市区,距国境线10公里以内、铁路和公路两侧1公里以内的地区,为禁猎区。
旗县级以上人民政府要将重点保护野生动物的繁殖、迁徙停歇、隐蔽、越冬等生息繁衍地划为禁猎区、禁渔区,并予公布。
第十条 每年3月15日至11月15日为禁猎期,5月1日至7月31日为禁渔期(季节性河流、湖泊除外)。
第十一条 单位或者个人发现国家和自治区重点保护野生动物伤病、危困、迷途,要及时采取紧急保护措施,并报当地野生动物行政主管部门组织抢救。
第十二条 影响或者破坏野生动物生存环境的建设项目,环境保护主管部门在审批环境影响报告书时,要征求同级野生动物行政主管部门的意见。
第十三条 保护管理野生动物资源所需经费,纳入旗县级以上人民政府财政预算,在野生动物行政主管部门的年度经费中列支。
自治区设立野生动物保护基金,基金管理办法由自治区野生动物行政主管部门会同财政部门制定,报自治区人民政府批准。
第十四条 保护国家和自治区重点保护野生动物,要避免造成农作物或者其他损失,造成损失的,由当地人民政府予以补偿。补偿办法和补偿标准由自治区人民政府制定。
第十五条 禁止在国家和自治区重点保护野生动物生息繁衍地使用有毒药物。特殊情况确需使用的,须经当地野生动物行政主管部门批准,并采取防范措施。
第三章 野生动物管理
第十六条 旗县级以上人民政府野生动物行政主管部门,要加强狩猎管理,制定保护管理和开发利用方案;定期组织野生动物资源调查,建立资源档案。
第十七条 因特殊需要猎捕国家重点保护野生动物的,按野生动物保护法第十六条规定办理。
猎捕自治区重点保护野生动物的,须向旗县级野生动物行政主管部门提出申请,经盟市野生动物行政主管部门审核,报自治区野生动物行政主管部门批准,核发特许猎捕证。
为了控制野生动物造成的危害,需要捕杀国家和自治区重点保护野生动物的,应当进行种群数量和危害情况调查,按上款规定办理。
第十八条 猎捕非国家和自治区重点保护野生动物的,必须向旗县级野生动物行政主管部门提出申请,领取狩猎证。
第十九条 猎捕者要按特许猎捕证、狩猎证的规定进行猎捕。持猎枪猎捕时,必须同时持有旗县级公安机关核发的持枪证。
第二十条 严禁使用军用武器、体育运动枪支、毒药、炸药和非人工操作制动的装置狩猎;严禁采取照明围猎、猎套、掏蛋等灭绝性手段狩猎。
第二十一条 进行狩猎生产,由旗县级野生动物行政主管部门划定猎区,经盟市野生动物行政主管部门审核,报自治区野生动物行政主管部门批准。
第二十二条 在野外对国家和自治区重点保护野生动物进行考察、科学研究,拍摄电影、录像,须经自治区野生动物行政主管部门批准。
外国人进行前款规定活动的,按野生动物保护法第二十六条规定办理。
第二十三条 驯养繁殖国家和自治区重点保护野生动物的,须提出申请,经国家或者自治区野生动物行政主管部门批准,由自治区野生动物行政主管部门核发驯养繁殖许可证,工商行政管理部门凭驯养繁殖许可证核发营业执照。
出售、转让驯养繁殖的国家和自治区重点保护野生动物的,由国家或者自治区野生动物行政主管部门批准。
第二十四条 收购、加工本办法规定保护的野生动物及其产品的单位和个人,必须持有旗县级以上野生动物行政主管部门核发的经营许可证;携带、运输、邮寄本办法规定保护的野生动物及其产品的单位和个人,必须持有旗县级以上野生动物行政主管部门核发的准运证。
第二十五条 出口国家重点保护野生动物及其产品的,进出口中国参加的国际公约所限制进出口的野生动物及其产品的,须持有贸易合同、经营许可证和营业执照,由所在地的旗县级野生动物行政主管部门签署意见,经自治区野生动物行政主管部门审核,报国务院野生动物行政主管部
门批准,取得国家濒危物种进出口管理机构核发的允许进出口证明书。
第二十六条 猎捕、经营利用野生动物及其产品的单位和个人,要缴纳野生动物资源保护管理费。收费标准和管理办法由自治区野生动物行政主管部门会同物价、财政部门制定,报自治区人民政府批准后施行。
第二十七条 各级野生动物行政主管部门要会同工商行政管理、物价、经贸、公安、环境保护、草原监理、交通、邮电、商业、供销、医药、海关等部门,对狩猎活动和野生动物及其产品的收购、出售、加工、运输、邮寄、进出口进行监督检查。
检查人员在执行公务时,要出示自治区野生动物行政主管部门及有关管理部门制发的检查证件。
第四章 法律责任
第二十八条 在禁猎区、禁猎期或者使用禁猎工具、禁猎方式狩猎的,由野生动物行政主管部门没收猎具、猎物和违法所得,并处以猎获物价格2至5倍的罚款;构成犯罪的,依法追究刑事责任。
第二十九条 未取得特许猎捕证或者未按特许猎捕证规定猎捕国家重点保护野生动物的,按野生动物保护法第三十一条规定处理;违法猎捕自治区重点保护野生动物构成犯罪的,依法追究刑事责任。
第三十条 未取得狩猎证或者未按狩猎证规定猎捕野生动物的,由野生动物行政主管部门没收猎获物和违法所得,处以猎获物价格2至5倍的罚款,并没收猎捕工具,吊销狩猎证。
未取得驯养繁殖许可证或者未按驯养繁殖许可证规定驯养繁殖国家和自治区重点保护野生动物的,由野生动物行政主管部门没收驯养繁殖的野生动物,吊销驯养繁殖许可证。
未取得经营许可证或者未按营业执照、经营许可证经营野生动物的,由工商行政管理机关和野生动物行政主管部门没收实物和违法所得,吊销营业执照和经营许可证,并处以实物价格2至5倍的罚款。
违反本办法规定,出售、收购野生动物及其产品的,由工商行政管理部门没收实物和违法所得,吊销有关证件,可以并处实物价格1至3倍的罚款。
在禁猎区破坏野生动物主要生息繁衍场所的,由野生动物行政主管部门责令停止破坏行为,限期恢复原状,处以实际损失1至3倍的罚款。
第三十一条 违反本办法规定,运输、携带、邮寄野生动物及其产品的,由工商行政管理部门没收实物和违法所得,并处以相当实物价格的罚款。没收的实物交当地野生动物行政主管部门处理。
第三十二条 伪造、倒卖、出租、转让特许猎捕证、狩猎证、经营许可证、准运证、驯养繁殖许可证、允许进出口证明书等证件的,由旗县级以上野生动物行政主管部门或者工商行政管理部门收缴、吊销有关证件,没收违法所得,处以违法所得2至5倍的罚款;构成犯罪的,依法追究
刑事责任。
第三十三条 违反本办法第二十二条规定的,由野生动物行政主管部门没收所获资料和标本。
第三十四条 野生动物行政主管部门及有关管理部门行政执法人员索贿受贿、徇私舞弊、玩忽职守的,由其上级主管部门给予行政处分;构成犯罪的,依法追究刑事责任。
第三十五条 当事人对行政处罚决定不服的,按野生动物保护法第三十九条规定办理。
第五章 附 则
第三十六条 本办法执行中的具体问题,由自治区野生动物行政主管部门解释。
第三十七条 本办法自发布之日起施行。
1991年12月24日