您的位置: 首页 » 法律资料网 » 法律法规 »

财政部、国家计委、国家粮食储备局、中国农业发展银行关于中央对地方粮食风险基金补助实行包干办法的通知

时间:2024-07-22 16:26:51 来源: 法律资料网 作者:法律资料网 阅读:8082
下载地址: 点击此处下载

财政部、国家计委、国家粮食储备局、中国农业发展银行关于中央对地方粮食风险基金补助实行包干办法的通知

财政部 国家计委 等


财政部、国家计委、国家粮食储备局、中国农业发展银行关于中央对地方粮食风险基金补助实行包干办法的通知


各省、自治区、直辖市财政厅(局)、计委(计经委)、粮食局(厅)、农业发展银
行分行:
为了进一步深化粮食流通体制改革,充分调动地方调整粮食种植结构、促进粮食购销企业顺价销售和加强粮食风险基金管理的积极性,根据《国务院关于进一步完善粮食流通体制改革政策措施的通知》(国发〔1999〕11号)文件精神,从1999年起,中央对地方粮食风险基金补助实行包干办法。
一、粮食风险基金包干的基本原则
(一)坚持粮食流通体制改革的“三项政策、一项改革”不变。实行包干办法后,国有粮食购销企业要坚持按保护价敞开收购农民余粮;除经批准处理陈化粮外,粮食购销企业要坚持顺价销售;粮食收购资金坚持封闭运行。
(二)对国有粮食购销企业的补贴标准整体水平不变。实行包干办法后,地方政府可以结合本地实际情况,制定鼓励国有粮食购销企业售粮的措施,适当降低超储费用补贴标准,拿出一部分资金在企业售粮后再给予适当奖励,以防止粮食购销企业坐享超储补贴,但费用补贴标准的整体水平不能减少。财政对国有粮食购销企业的超储利息补贴要及时足额用于支付农业发展银行利息。财政对国有粮食购销企业的超储库存所需利息要据实结算。
(三)实行包干后,地方财政要确保对国有粮食购销企业的各项补贴及时足额到位。包干如有超支由地方政府负责补足;包干如有结余,可用于补助处理国有粮食购销企业商品周转粮库存中陈化粮的差价损失或用于消化粮食企业新老财务挂账等粮食方面的开支,但不得挪作他用。
二、粮食风险基金包干方案
从1999年开始,中央财政对各省、自治区、直辖市粮食风险基金补助款实行总额包干。包干总额的计算依据为:
(一)超储库存基数。以各省、自治区、直辖市1998粮食年度超储库存最高3个月的平均数为包干基数。
(二)超储补贴标准。利息根据超储库存包干基数和当期利率据实计算,费用按原粮每年每斤3分钱计算。
(三)省级储备粮油利息、费用按省级人民政府批准的1998年实际建立的省级储备粮油数量和利息、费用补贴标准计算。
(四)对部分地区1998年受灾减产,影响库存等特殊因素予以适当照顾。
(五)中央对各省、自治区、直辖市的包干的补助资金,按国务院批准的1999年粮食风险基金中央财政筹资比例计算。
三、国有粮食购销企业商品周转粮中陈化粮的处理。由于陈化粮数量和处理数量、价差损失目前均难以核定,因此,处理的差价损失补助办法,由有关部门另行制定。
四、按现行政策规定,对地方自营出口粮食亏损用超储库存1年所需利息费用补贴的,已包含在中央对地方粮食风险基金补助包干总额内。对黑龙江、吉林、辽宁、内蒙古四省(区)地方自营玉米出口亏损补贴标准超过超储库存1年所需利息费用的部分,由地方财政和中央财政按现行政策规定的比例另外筹集资金解决。
五、中央对地方的包干补助款按季均衡拨付各地在农业发展银行开设的粮食风险基金专户,地方自筹配套资金也必须按季足额拨付到专户。财政部门会同粮食、农发行等有关部门,将本季度应拨付的粮食超储补贴和省级储备粮油利息、费用补贴及时拨补到购销企业。
六、对1999年1月1日以前发生的粮食风险基金支出,中央财政对各省、自治区、直辖市仍按原政策规定据实清算。财政部门对粮食购销企业仍按原政策规定的补贴标准和实际超储库存据实补贴,欠拨的补贴资金要尽快拨补到位。
七、粮食风险基金预算管理办法、专户管理办法,财务、会计处理办法不变。实行包干办法后,中央财政会同有关部门仍要继续加强粮食风险基金管理,中央对地方粮食风险基金包干补助和地方配套资金继续在农业发展银行专户管理。地方财政仍需按现行比例配套资金,地方配套不足的,中央财政通过扣款强制到位。各省、自治区、直辖市财政厅(局)必须按时上报粮食风险基金专户月报。针对目前某些地方粮食风险基金滞留、欠拨过多的情况,财政部将会同农业发展银行制定具体的监管办法。
八、实行包干办法后,各级财政部门要会同有关部门制定促进粮食购销企业扩大销售的办法,进一步加强粮食购销企业财务监管,监督企业如实反映粮食收购、销售和库存情况,防止弄虚作假骗取财政补贴;监督企业严格按现行财会制度规定,如实核算成本、费用,防止多转或少转成本、费用等短期行为,扩大亏损或潜亏挂账。
九、本办法由财政部负责解释。各省、自治区、直辖市财政厅(局)可依据本办法会同有关部门制定具体实施办法。



国家版权局废止的有关著作权管理的规章、规范性文件目录(1984年-1999年)

国家版权局


国家版权局废止的有关著作权管理的规章、规范性文件目录(1984年-1999年)



中华人民共和国国家版权局令
 (第4号)

  为促进政府职能转变,提高依法行政水平,保障著作权行政管理秩序,经2003年7月16日国家版权局第2次局务会议审议通过,决定废止一批规章及规范性文件,现予公布,自公布之日起生效。

                        国家版权局局长 石宗源
                           2003年12月4日
    

┏━━┯━━━━━━━━┯━━━━━━━━━━━━━━━━━━━┯━━━━━━━━━━┓
┃序号│发文机关    │文件名称               │文号及发文时间   ┃
┠──┼────────┼───────────────────┼──────────┨
┃1  │文化部出版局  │关于颁发《图书、期刊版权保护试行条例》│文出字(84)第849号 ┃
┃  │        │的通知                │84-6-15      ┃
┃  │        │附:图书、期刊版权保护试行条例    │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃2  │文化部出版局  │文化部关于转发《书籍稿酬试行规定》的通│文出字(84)第1791号┃
┃  │        │知                  │84-10-19     ┃
┃  │        │附:文化部出版局关于试行《书籍稿酬试行│          ┃
┃  │        │规定》的报告             │          ┃
┃  │        │书籍稿酬试行规定           │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃3  │文化部出版局  │关于颁发《图书、期刊版权保护试行条例实│文出字(85)第37号 ┃
┃  │        │施细则》和《图书约稿合同》、《图书出版│85-1-1      ┃
┃  │        │合同》的通知             │          ┃
┃  │        │附:图书、期刊版权保护试行条例实施细则│          ┃
┃  │        │图书约稿合同             │          ┃
┃  │        │图书出版合同             │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃4  │国家版权局   │关于期刊一年专有出版权的说明     │(86)权字第20号  ┃
┃  │        │                   │86-4-11      ┃
┠──┼────────┼───────────────────┼──────────┨
┃5  │国家版权局   │关于内地出版港澳同胞作品版权问题的暂行│(86)权字第30号  ┃
┃  │        │规定                 │86-5-24      ┃
┠──┼────────┼───────────────────┼──────────┨
┃6  │国家版权局   │关于清理港、澳、台作者稿酬的通知   │(87)权字第52号  ┃
┃  │        │                   │87-11-16     ┃
┠──┼────────┼───────────────────┼──────────┨
┃7  │国家版权局   │关于汇总使用港澳台同胞作品付酬情况的通│(88)权字第22号  ┃
┃  │        │知                  │88-5-9      ┃
┠──┼────────┼───────────────────┼──────────┨
┃8  │国家版权局   │关于台、港、澳作者来大陆领取稿酬的有关│(88)权字第33号  ┃
┃  │        │问题的通知              │88-8-1      ┃
┠──┼────────┼───────────────────┼──────────┨
┃9  │国家版权局   │关于《图书、期刊版权保护试行条例》第十│(88)权字第56号  ┃
┃  │        │五条第四款的解释           │88-12-22     ┃
┠──┼────────┼───────────────────┼──────────┨
┃10 │国家版权局   │关于版权纠纷处理决定书规范化的通知  │(88)权字第53号  ┃
┃  │        │                   │88-12-27     ┃
┠──┼────────┼───────────────────┼──────────┨
┃11 │文化部、国家版权│关于今后解释《图书、期刊版权保护试行条│(89)权字第2号   ┃
┃  │局       │例》的联合通知            │89-1-6      ┃
┠──┼────────┼───────────────────┼──────────┨
┃12 │新闻出版署、国家│关于征订台、港、澳作者的图书应出示合同│(89)新出联字第2号 ┃
┃  │版权局     │审核登记号的通知           │89-2-17      ┃
┠──┼────────┼───────────────────┼──────────┨
┃13 │国家版权局   │关于认真执行对台、港、澳版权贸易有关规│(90)权字第3号   ┃
┃  │        │定的通知               │90-2-2      ┃
┠──┼────────┼───────────────────┼──────────┨
┃14 │国家版权局   │关于版权贸易合同审核登记问题的补充规定│(90)权字第7号   ┃
┃  │        │                   │90-5-14      ┃
┠──┼────────┼───────────────────┼──────────┨
┃15 │国家版权局   │关于适当提高书籍稿酬的通知      │(90)权字第11号  ┃
┃  │        │附:书籍稿酬暂行规定         │90-6-15      ┃
┠──┼────────┼───────────────────┼──────────┨
┃16 │国家版权局办公室│关于下发对台、港、澳版权贸易示范《出版│(90)权字15号   ┃
┃  │        │合同》的通知             │90-8-2      ┃
┃  │        │附:《出版合同》           │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃17 │国家版权局   │关于加强音像版权管理的通知      │(91)权字16号   ┃
┃  │        │附:国家版权局“版权贸易合同审核登记证│91-7-2      ┃
┃  │        │书”样本               │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃18 │国家版权局   │关于维护出版社出版外国作品专有出版权的│(91)权字第34号  ┃
┃  │        │通知                 │91-11-29     ┃
┠──┼────────┼───────────────────┼──────────┨
┃19 │国家版权局   │关于颁发著作权许可使用合同标准样式的通│(92)权字第4号   ┃
┃  │        │知                  │92-1-24      ┃
┃  │        │附:《图书出版合同》(标准样式)   │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃20 │国家版权局   │关于《图书出版合同》(标准样式)   │国权(1992)16号  ┃
┃  │        │基本稿酬加印数稿酬计算公式的订正   │92-4-4      ┃
┠──┼────────┼───────────────────┼──────────┨
┃21 │新闻出版署、国家│关于认真做好报刊转摘作品付酬及收转工作│新出联(1993)8号  ┃
┃  │版权局     │的通知                │93-5-14      ┃
┠──┼────────┼───────────────────┼──────────┨
┃22 │国家版权局   │关于《对侵犯著作权行为行政处罚的实施办│国权(1994)47号  ┃
┃  │        │法》的通知              │94-6-23      ┃
┃  │        │附:对侵犯著作权行为行政处罚的实施办法│          ┃
┠──┼────────┼───────────────────┼──────────┨
┃23 │国家版权局   │关于印发《对侵犯著作权行为行政处罚的实│国权(1995)28号  ┃
┃  │        │施办法》的通知            │95-7-25      ┃
┃  │        │附:对侵犯著作权行为行政处罚的实施办法│          ┃
┠──┼────────┼───────────────────┼──────────┨
┃24 │国家版权局   │著作权行政处罚实施办法        │局令第2号      ┃
┃  │        │                   │97-1-28      ┃
┠──┼────────┼───────────────────┼──────────┨
┃25 │国家版权局   │关于发布著作权法实施条例第五条中“表演│国权(1999)43号  ┃
┃  │        │”的具体应用问题的解释的通知     │99-12-9      ┃
┃  │        │附:著作权法实施条例第五条中“表演”的│          ┃
┃  │        │具体应用问题的解释          │          ┃
┠──┼────────┼───────────────────┼──────────┨
┃26 │国家版权局   │关于制作数字化制品的著作权规定    │国权(1999)45号  ┃
┃  │        │                   │99-12-9      ┃
┗━━┷━━━━━━━━┷━━━━━━━━━━━━━━━━━━━┷━━━━━━━━━━┛




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

国土资源部


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

(国土调查办发〔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