广州市城市路桥隧道车辆通行费年票制试行办法
广东省广州市人民政府
广州市城市路桥隧道车辆通行费年票制试行办法
广州市人民政府令第6号
《广州市城市路桥隧道车辆通行费年票制试行办法》已经２０００年１２月４日市政府常务会议讨论通过，现予发布，自２００１年１月１日起施行。
市长 林树森
二○○○年十二月十一日
广州市城市路桥隧道车辆通行费年票制试行办法
第一条 为充分发挥本市城市道路的交通功能，平衡市区交通流量，规范贷款建设的路桥、隧道的收费行为，利于贷款的偿还，省人民政府同意本市实施城市路桥隧道车辆通行费年票制，结合本市实际情况，制定本办法。
第二条 路桥隧道车辆通行费（以下简称通行费）年票制，是指对在本市登记上牌的机动车辆（不包括番禺区、花都区、增城市、从化市，以下简称本市车辆）一次性征收全年通行费的制度。
对外地进入本市的机动车辆，在市区各入城路口收费站按次征收通行费（次票）。
征收的通行费用于偿还城市路桥、隧道建设的贷款。
第三条 实行年票制收费的路桥、隧道包括：内环路及其放射线、广园快速路西段、鹤洞大桥、解放大桥、江湾大桥、海印大桥、广州大桥、珠江隧道。
第四条 广州市市政管理部门是通行费征收管理的行政主管部门，负责组织实施本办法。
第五条 公安交通管理部门协助市市政管理部门做好车辆缴费的检查工作。公安交通管理部门应当向市市政管理部门提供本市车辆户籍及相关资料，并在办理车辆年审、新车入户、车辆报废、外地车迁入等手续时，检查车辆缴纳年票通行费情况。对没有缴费的车辆，责令车主到市市政管理部门缴费后，方可办理有关手续。
建设、交通、物价、财政、税务等有关部门，依照各自职能协同实施本办法。
第六条 通行费由市市政管理部门在市内缴费点和各进城路口收费站进行征收。 经华南快速干线、环城高速公路等收费路桥进入市区的外地车辆，由市市政管理部门委托其收费站代行征收通行费。
第七条 通行费收费标准如下：
（一）一类车：摩托车，年票１８０元，次票２元。
（二）二类车：２０座以下（含２０座）客车、２吨以下（含２吨）货车、机动三轮车、简易机动车，年票９８０元，次票１０元。
（三）三类车：２１座（含２１座）至５０座（含５０座）客车、２吨至５吨（含５吨）货车，年票２０４０元，次票２０元。
（四）四类车：５１座以上（含５１座）客车、５吨至１５吨（含１５吨）货车，年票２８８０元，次票３０元。
（五）五类车：１５吨以上货车和各种集装箱车，年票２６４０元，次票５０元。
（六）出租汽车（的士），年票９８０元，次票１０元。
（七）公共汽车，年票６００元，次票２０元。
第八条 本市车辆必须在每年办理年审手续前缴纳年票通行费。
新车入户、外地车迁入或者车辆报废，在办理有关手续时应当从当月起缴纳或者停止缴纳年票通行费。
已缴纳年票通行费的本市车辆在改装、换牌之日起３０日内持有关证件到原年票缴费点办理变更手续。
第九条 市市政管理部门对缴纳年票通行费的车辆，应当开具发票并发给与车牌号码、类别相一致的年票正、副本。车主应当将年票正本粘贴在车前挡风玻璃的右上角，摩托车贴于车头显眼处，副本随车携带，以备查验。
第十条 长期在市区行使的外地车辆，应当缴纳年票通行费；经常出入市区的外地车辆，可以选择缴纳年票通行费或者次票通行费。
次票通行费实行单向按次征收，凭票进入市区；进入一次缴纳一次。在市内留驻的，三日内有效；超过有效期的，必须另行缴纳通行费。缴费发票必须保存至离开市区为止，以备查验。
第十一条 军车、警车、消防车、省政府治理公路“三乱”督察队标志车以及救护车、殡葬车、抢险车、外国领事馆用车、残疾人专用机动车等免征通行费。
第十二条 通行费票证售出后，不予退换。年票如有遗失，车主必须提出书面补领申请，持《机动车行驶证》和缴费发票到原年票缴费点办理挂失手续，经核准后，给予补票；次票遗失不补，必须自行到收费站点重新补缴。
第十三条 通行费年票正副本、次票，任何人不得转借、冒用和伪造。
第十四条 通行费收费站、点应悬挂《广东省行政事业性收费许可证》，并向社会公开收费站、点名称和审批机关、主管部门、收费用途、收费标准、收费员工证号及投诉电话，接受社会监督。
第十五条 征收的通行费按行政事业性收费管理规定，纳入财政专户管理。市市政管理部门以及受委托的收费单位必须建立健全财务、审计、统计、票据管理制度和报表制度，及时做好通行费的征收和上缴工作。
第十六条 市市政管理部门应当定期向社会公布通行费的征收、管理、使用及贷款偿还情况。还清贷款后，应当停止征收通行费。
第十七条 市市政管理部门有权对市区内停车场、车站、码头和经过收费站出入口的车辆进行通行费缴纳情况的稽查。任何单位和个人不得拒绝接受市市政管理部门的检查。
对没有缴纳通行费的车辆，市市政管理执法人员应当即时处理；不能即时处理或者影响交通畅通的，可以责令驾驶人员到指定地点接受处理，处理完毕后，方可驶离。
第十八条 市市政管理执法人员执行公务时，应当统一着装，佩戴《行政执法证》。
第十九条 违反本办法，有下列行为之一的，由市政管理部门按下列规定处理：
（一）不按规定缴纳年票通行费的，责令其补缴年票通行费，并从应缴日起按日加收应当缴纳年票费额２‰的滞纳金。
（二）不按规定缴纳次票通行费的，除责令其补缴规定费额外，可并处５０元罚款。
（三）转借、冒用、使用伪造年票正副本的，除责令其补缴规定费额外，对单位并处５００元罚款，对个人并处１００元罚款。
（四）转借、冒用、使用伪造次票的，除责令其补缴规定费额外，对单位并处１００元罚款，对个人并处５０元罚款。
（五）不按规定将年票正本粘贴在车前挡风玻璃右上角、摩托车贴于车头显眼处的，处以警告或者２０元罚款。
（六）未购票车辆不按指定路线通过收费站，进入有票车道的，处以警告，可并处２０元罚款。
第二十条 阻碍、围攻、谩骂、殴打依法执行公务的市政管理执法人员或者违反本办法，应当给予治安管理处罚的，由公安部门依照《中华人民共和国治安管理处罚条例》处罚；构成犯罪的，由司法机关依法追究其刑事责任。
第二十一条 市市政管理部门工作人员玩忽职守，滥用职权，徇私舞弊的，由所在单位或者上一级主管部门给予行政处分；构成犯罪的，依法追究其刑事责任。
第二十二条 本办法自２００１年１月１日起施行。本市过去颁发的有关机动车辆通行路桥隧道征收通行费的有关规定。如与本办法不一致的，以本办法为准。
关于统一图幅理论面积与图斑椭球面积计算要求的通知
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
（国土调查办发〔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
