蚌埠市城市古树名木管理办法
安徽省蚌埠市人民政府办公室
关于印发蚌埠市城市古树名木管理办法的通知
蚌政办〔2011〕24号
蚌埠市人民政府办公室关于印发蚌埠市城市古树名木管理办法的通知
各县、区人民政府，市政府各部门、各直属单位：
《蚌埠市城市古树名木管理办法》已经市十四届人民政府第27次常务会议研究通过，现印发给你们，请认真贯彻执行。
二○一一年四月十一日
蚌埠市城市古树名木管理办法
第一条 为加强古树名木的保护管理，根据国务院《城市绿化条例》、建设部《城市古树名木保护管理办法》等法规规章规定，结合本市实际，制定本办法。
第二条 本办法适用于本市城市规划区内古树名木的保护和管理。
第三条 本办法所称古树名木，是指百年以上树龄的树木，稀有、珍贵树木，具有历史价值或重要纪念意义的树木。
第四条 市城市绿化行政主管部门是古树名木主管部门，其所属的城市园林绿化管理部门负责具体管理工作。
市林业、城市管理行政执法、规划、农业、水利等其他有关部门，按照各自职责，协同做好城市古树名木管理工作。
在城市规划区内，有关法律、法规规定由林业、水利等行政主管部门管理的绿化，其所管辖范围内的古树名木管理保护工作，由林业、水利等管理部门负责。
第五条 古树名木主管部门应当加强对古树名木的保护宣传和技术推广，广泛普及保护知识，不断提高保护和管理水平。
第六条 古树名木主管部门必须对管辖范围内的古树名木进行调查登记、鉴定分级、建立档案、设立标志、划定保护范围、制定保护措施、确定保护管理责任。
古树名木主管部门应定期对古树名木生长和管护情况进行检查。对长势濒危的古树名木提出抢救措施，并监督实施。
第七条 任何单位和个人都有保护古树名木及其附属设施的义务。对损伤、破坏古树名木及其附属设施的行为，有权劝阻、检举和控告。
第八条 鼓励单位和个人通过多种方式参与古树名木的保护管理。对保护管理古树名木成绩显著的单位或个人，古树名木主管部门应当给予表彰和奖励。
第九条 编制城乡建设规划应当在古树群周围划定绿线，保护古树群的生长环境和风貌。
第十条 古树名木的养护责任，按下列规定承担：
（一）生长在城市公共绿地、道路、公园等范围内的古树名木，由城市园林绿化管理部门负责养护。
（二）生长在铁路、公路、水库和河道用地范围内的古树名木，分别由铁路、公路、水库和河道管理部门负责养护。
（三）生长在林地、风景名胜区内的古树名木，由林地、风景名胜区管理部门负责养护。
（四）生长在机关、部队、学校、团体、寺庙、教堂、企事业单位管界内的古树名木，由所在单位负责养护。
（五）生长在居住小区内或私人庭院内的古树名木，由业主或业主委托的物业管理部门负责养护。
（六）生长在其他范围内的古树名木，由所在地的街道办事处和镇（乡）人民政府负责组织养护。
变更古树名木养护责任或责任人，应当到古树名木主管部门登记备案。
第十一条 古树名木养护费用由养护责任单位或责任人负担。养护古树名木确有困难的，可以向古树名木主管部门申请养护补助。
各级人民政府设立专项经费用于古树名木的资源普查、建档挂牌、复壮、抢救、养护补助等。
第十二条 古树名木养护责任单位或责任人，应当按照古树名木主管部门规定的养护管理措施实施养护管理，保障古树名木正常生长，并接受古树名木主管部门的指导、监督、检查。
第十三条 古树名木受到损害或长势衰弱，养护责任单位或责任人应当及时报告古树名木主管部门，并由古树名木主管部门组织抢救、复壮。
对已经死亡的古树名木, 必须上报并经市级古树名木主管部门确认、查明原因及明确责任后，方可处理，予以注销登记。
第十四条 禁止下列损坏古树名木的行为：
（一）刻划钉钉、缠绕绳索、架设电线；
（二）攀树、折枝、挖根或剥损树枝、树干、树皮；
（三）依树搭棚或借用树干做支撑物；
（四）擅自采摘树叶、果实和种子；
（五）距树冠垂直投影5米范围内堆放物料、挖坑取土、兴建临时设施建筑物，倾倒有害污水、污物垃圾，动用明火或排放烟气；
（六）在距树干边缘2米范围内地面，用水泥或其他建筑材料进行硬覆盖；
（七）擅自砍伐、修剪、移植、转让买卖；
（八）其他损害古树名木的行为。
第十五条 生产、生活设施妨碍和危害古树名木生长的，由古树名木主管部门责令有关单位或个人限期采取措施，排除妨碍和危害。
第十六条 建设项目在进行规划、设计、施工时，应当对古树名木采取避让和保护措施。保护措施应当由建设单位报经古树名木主管部门批准；未经批准的，不得施工。
因特殊情况，确需迁移古树名木的，建设单位必须提出妥善的迁移方案，经古树名木主管部门审查同意后，按古树名木级别报同级或上级人民政府批准，按照古树名木移植的有关规定组织施工。迁移古树名木所需全部费用及树木移植后5年内的恢复、养护费用，由申请人承担。
第十七条 违反本办法有关规定对古树名木造成侵害的，由古树名木主管部门责令其改正，需要给予行政处罚的，由城市管理行政执法部门依法予以处罚；构成犯罪的，由司法机关依法追究刑事责任。
第十八条 古树名木主管部门工作人员应当认真履行职责，切实做好古树名木的保护管理工作。因玩忽职守，致使古树名木损伤或死亡的，由所在单位或上级主管部门给予行政处分；构成犯罪的，由司法机关依法追究刑事责任。
第十九条 市辖各县古树名木的保护管理，可参照本办法执行。
第二十条 本办法自发布之日起施行。
关于统一图幅理论面积与图斑椭球面积计算要求的通知
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
（国土调查办发〔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
