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

绍兴市城市供水管理办法

时间:2024-07-03 21:50:43 来源: 法律资料网 作者:法律资料网 阅读:9525
下载地址: 点击此处下载

绍兴市城市供水管理办法

浙江省绍兴市人民政府


绍兴市城市供水管理办法

市政府令〔2003〕60号


  现发布《绍兴市城市供水管理办法》,自2003年5月1日起施行。


市 长
二OO三年四月七日

绍兴市城市供水管理办法

  第一条 为了加强城市供水管理,发展城市供水事业,保障城市生活、生产用水和其他各项建设用水,保障供水企业和用户的合法权益,根据《城市供水条例》、《浙江省城市供水管理办法》等有关法规、规章的规定,结合本市实际,制定本办法。
  第二条 在本市市区范围内从事城市供水工作和使用城市供水,应当遵守本办法。
  第三条 城市供水工作坚持开发保护水源和计划用水、节约用水相结合,实行统一规划、优化资源配置。
  第四条 优先发展城市公共供水,加强城市供水科学技术研究,推广先进技术,提高城市供水的现代化水平。
  第五条 市供水行政主管部门负责市区城市供水管理工作。
环保、卫生、价格、质量技术监督、水行政等主管部门应按照各自职责,协助市供水行政主管部门共同做好城市供水管理工作。
  第六条 市城市规划行政主管部门应当会同有关部门根据城市建设和社会发展需要编制城市供水发展规划,并将其纳入城市总体规划 。
  第七条 市水行政、环保、卫生和供水等行政主管部门按照国家有关规定,划定饮用水水源保护区,经市人民政府批准后公布,并对饮用水水源实施保护和管理。对跨市的饮用水水源保护区,按照国家有关规定进行保护和管理。
市人民政府在饮用水水源保护区范围内设置保护标志,任何单位和个人未经批准不得毁坏和移动。
  在饮用水水源保护区内,禁止一切污染水质的活动。
  第八条 城市供水工程的建设,应当按照城市供水发展规划及其年度建设计划执行。城市公共供水管道及附属设施的建设应当与城市规划区范围内的道路等各项基础设施和公共设施同步建设。
  第九条 城市供水工程的设计、施工和监理应当由持有相应资质证书的单位承担,并遵守国家有关技术标准和规范。
禁止无证或超越资质证书规定的经营范围承担城市供水工程的设计、施工和监理任务。
  第十条 城市供水工程竣工后,应当按照国家有关规定组织验收;未经验收或验收不合格的,不得投入使用。
  第十一条 城市新建、改建、扩建工程项目需要增加用水的,其工程项目总概算应当包括供水工程建设投资。
  第十二条 实际用水地点在城市公共供水管网覆盖范围以外,无法使用城市公共供水,确需自建设施供水的,应当经有关行政主管部门批准后方可建设。
  在城市公共供水管网覆盖范围内的自建设施供水,要逐步核减许可取水量,直至完全取消。
  第十三条 供水企业从事城市供水经营活动,必须依法经有关行政主管部门批准。国家对供水企业从事供水活动另有规定要求的,从其规定。
  第十四条 供水企业应当建立、健全水质检测制度,设立符合国家规定的水质检验采样点,定期检验水源水、出厂水和管网水水质,确保供水水质符合国家规定的饮用水卫生标准。
用户对水质有特殊要求的,应当自行进行特殊处理。
第十五条 二次供水设施选址、设计、施工及所用材料,应保证不使饮用水水质受到污染,并有利于清洗和清毒,各类蓄水设施应加强卫生防护,定期清洗和清毒。
  第十六条 供水企业应当按照国家有关规定设置管网测压点,做好水压监测工作,确保供水管网的压力符合国家和省规定的标准。
  第十七条 在供水系统正常运行的情况下,供水企业应当保持不间断供水。因工程施工、设施检修等原因确需临时停止供水或者降低供水水压的,应当提前二十四小时通知用户;因发生灾害或者紧急事故,不能提前通知的,应当在抢修的同时通知用户,尽快恢复正常供水。
供水设施抢修时,有关单位和个人应当给予支持和配合,不得阻扰或者干扰抢修工作进行,有关手续可在抢修结束后补办。
  第十八条 使用公共供水不能间断供水的用户,应当自建贮水设备或采取其他保障措施;对水压有特殊要求而确需使用加压设备的用户,应通过贮水设备或采取其他隔离保障措施进行加压。
  第十九条 供水企业与用户应当签订《供用水合同》,明确双方的权利和义务。用户要求改变用水性质,进行更名、过户与销户等变更的,应向供水企业办理有关手续。
  第二十条 供水企业应按约定期限抄表。由于用户原因造成不能抄见水表的,供水企业可根据该用户前六个月中的最高月用水量估算当月用水量;如用户在三个月内不能解决妨碍抄表问题的,供水企业不退还多估水费;用户实际用水量超过估算用水量的,用户应按实际用水量支付水费。若遇水表故障的,供水企业可根据用户前三个月平均用水量计收水费。
  第二十一条 用户应当按照规定的计量标准和水价标准按时交纳水费。用户接到水费缴纳通知单之日起十五日内不交纳水费的,按应交纳水费额每日加收千分之五的滞纳金。用户无正当理由连续两次不缴纳水费的,供水企业可依据合同规定对其停止供水,但应当提前五天通知用户。
  第二十二条 水费结算以结算水表数据为准,结算水表由法定检验机构按有关规定进行检定。用户对结算水表计量准确性有异议的,可依法申请检定。经检定,误差超过规定标准的,当月按验表结果核收水费,并免交其他费用;误差在规定标准以内的,用户除缴纳正常水费外,还应承担验表及水表复装所必需的费用。
  第二十三条 供水企业可根据用户实际用水情况调整结算水表口径。用户全月平均小时用水量低于其结算水表最小流量时,供水企业可根据供用水合同规定改小结算水表口径,并由用户承担相应的工料费。
  第二十四条 用户无正当理由连续六个月以上停止用水的,供水企业可拆除结算水表并对其销户。用户要求复装的,应当重新办理用水手续。
  第二十五条 因城市房屋拆迁要求停止供水的,由房屋拆迁人负责办理用水销户手续,并与供水企业结清被拆迁人应付水费。房屋拆迁人在结清水费后可向被拆迁人追偿。
房屋产权或企业产权发生变更时,变更双方应明确水费支付责任并办妥用水过户手续。
  第二十六条 城市供水价格实行政府定价。城市供水价格的制定和调整,由供水企业向市价格行政主管部门提出书面申请,经价格行政主管部门按规定权限和程序依法审核批准后实施。
  城市供水服务项目的价格标准,由市价格行政主管部门依法核定。
  第二十七条 市政、园林与环卫等部门因作业确需用水的,应向供水企业办理用水手续,由供水企业合理设置专用供水设施。各使用部门应按规定向供水企业缴纳水费。
  第二十八条 禁止擅自改变用水性质、偷盗或者向其他单位和个人转供城市公共供水。
  第二十九条 供水企业对其取水口、泵站、净水厂和配水管网等设施应定期检查维修,确保安全运行。
  第三十条 城市供水设施管理、维护责任以结算水表为分界点进行划分:
  (一)结算水表以前(含结算水表)靠近水源侧的供水设施由城市供水企业负责;
  (二)算水表以后(不含结算水表)靠近用户侧的供水设施由用户负责,也可由用户有偿委托供水企业承担。
  第三十一条在城市公共供水设施安全保护范围内,禁止从事下列活动:
  (一)建造永久或临时建筑物、构筑物;
  (二)开挖沟渠或者挖坑取土;
  (三)打桩或者顶进作业;
  (四)其他可能损坏供水设施或危及供水设施安全的活动。
  城市公共供水设施的安全保护范围确定应当遵守国家有关给水设计规范的要求。接入市区的公共供水主管道及其附属设施的地面和两侧安全保护范围为:口径1000毫米以上、不满1400毫米管道两侧各5米,口径1400毫米以上、不满1800毫米管道两侧各6米,口径1800毫米以上管道两侧各8米。
  第三十二条 涉及城市供水设施的建设工程,建设单位或者施工单位应当在开工前向城市供水企业查明地下供水管网情况,可能影响城市供水设施安全的,建设单位或者施工单位应当与城市供水企业商定相应的保护措施,并组织实施。
  第三十三条 任何单位和个人不得擅自改装、迁移或者拆除城市公共供水设施。工程建设确需改装、迁移或者拆除城市公共供水设施的,建设单位应当在申请工程规划许可证前,与供水企业商定相应补救措施,报经市供水行政主管部门批准。
  第三十四条 市区新建住宅的供水应实行“一户一表、计量出户”,对管网压力条件允许的新建住宅,可利用城市供水管网直接供水。未实行“一户一表、计量出户”的住宅要按有关规定进行“一户一表、计量出户”改造。
  第三十五条 由用户自行建设的结算水表前供水设施,其设计方案应征得供水企业同意;施工过程中,供水企业应加强检查、监督。供水设施经供水企业验收合格并移交其统一管理后方可投入使用。
建设单位建设结算水表后用水设施的,应通知供水企业参加其设计图纸会审,用水设施必须按有关规范进行冲洗、消毒与水压试验。
  第三十六条  自建供水系统的用户不得擅自将内部自建供水设施与城市公共供水管道相连接。
  禁止生产或者使用有毒、有害物质单位的生产用水管网系统及设备与城市公共供水管网系统直接连接。
  第三十七条 禁止在城市公共供水管网上直接装泵抽水或者在与城市公共供水管网直接相联的用户用水设施上装泵抽水。
  第三十八条 公共消火栓为消防专用,禁止任何单位或个人损坏、擅自启用公共消火栓。
  第三十九条 使用城市公共供水的单位,应根据国家有关规定设置内部消防自救系统,供水企业对其内部消防自救系统采用间接供水。
确需利用城市公共供水直接作为单位内部消防系统水源的,应征得供水企业同意,订立消防用水合同并支付相应费用。消防用水管线应与其正常生产、生活用水管线相分离,无火警不得启用消火栓。
  第四十条 供水管材、设备、净水剂、用水器具及其他与饮用水接触的材料应当符合国家标准,国家尚未制定统一标准的,应当符合地方或行业标准。
禁止生产、销售和使用不符合标准或明令淘汰的供水管材、设备、净水剂、用水器具及其他与饮用水接触的材料。
  第四十一条 违反本办法规定的,由有关行政主管部门依法进行处罚,构成犯罪的,依法追究刑事责任。
  第四十二条 城市供水管理人员玩忽职守、滥用职权、徇私舞弊的,由其所在单位或上级机关给予行政处分;构成犯罪的,依法追究刑事责任。
  第四十三条 本办法自2003年5月 1日起施行,1996年8月30日绍兴市人民政府颁布的《绍兴市城市供水管理办法》同时废止,各县(市)的城镇供水管理可参照本办法执行。


重庆市查处无照经营行为条例

重庆市人大常委会


重庆市查处无照经营行为条例
重庆市人民代表大会常务委员会
重庆市人民代表大会常务委员会公告第158号


(2001年5月25日经重庆市第一届人民代表大会常务委员会第三十二次会议通过,自2001年7月1日起施行)


第一条 为查处无照经营行为,维护社会主义市场经济秩序,根据有关法律、法规的规定,结合本市实际,制定本条例。
第二条 在本市行政区域内从事商品生产、经营或营利性服务,法律、法规规定应登记的,经营者应依法向工商行政管理部门申请登记,领取营业执照。
临时占道经营,符合城市容貌管理规定的,应持临时占道许可证及相关文件,到工商行政管理部门申请登记,领取营业执照后,在批准的地点和期限内营业。
本条例所称无照经营行为,是指违反前两款规定,擅自从事商品生产、经营或营利性服务的行为。
第三条 市和区、县(自治县、市)工商行政管理部门负责依法查处无照经营行为。
市和区、县(自治县、市)人民政府应加强对查处无照经营行为工作的领导,市政、规划、公安等有关部门、乡镇人民政府和街道办事处应按各自职责,积极协助工商行政管理部门做好查处无照经营行为的工作。
第四条 查处无照经营行为,应坚持引导与查处,教育与处罚相结合的原则。
第五条 经营者应将营业执照置放在生产场所或营业场所的醒目位置。
经营者应依照法律、法规的规定,办理验照、年检、变更登记或其他手续。
任何组织和个人不得为无照经营者非法提供营业执照、合同文本、银行帐户、发票等经营条件。
第六条 工商行政管理部门查处无照经营行为时,可行使下列职权:
(一)调查涉嫌无照经营的经营者及相关的组织和个人;
(二)查阅、复制与涉嫌无照经营有关的合同、发票、帐册等经营资料和财务资料;
(三)检查与涉嫌无照经营相关的场所和物品;
(四)扣留、封存与涉嫌无照经营有关的资料和财物;
(五)法律、法规规定的其它职权。
第七条 工商行政管理部门实施扣留、封存措施、应遵循下列规定:
(一)须经市或区、县(自治县、市)工商行政管理部门主要负责人批准,交通不便地区或不及时采取扣留、封存措施可能影响查处时,可先予扣留、封存后,在三个工作日内补办批准手续;
(二)须由两名以上执法人员执行,并应着装和出示有效执法证件;
(三)出具扣留、封存财物通知书,并告知当事人申请行政复议和提起诉讼等权利。对扣留、封存的资料和财物,应当场开列清单一式二份,由执法人员和当事人签名或盖章,一份交当事人。当事人拒绝签名或盖章或当事人不在现场的,应有两名以上见证人签名,并由执法人员注明情况。
(四)对扣留、封存的资料和财物,应妥善保管,严禁动用、调换或损毁。被扣留、封存物品易腐烂变质的,工商行政管理部门在征得当事人同意后,可先行处理。
(五)扣留、封存的期限不得超过三十日,因案情复杂不能在规定期限内作出处理决定需延长期限的,应经市或区、县(自治县、市)工商行政管理部门批准,但延长期限不得超过三十日。
第八条 未依法登记为有限责任公司或者股份有限公司,而冒用有限责任公司或者股份有限公司名义,或未经登记,以合伙企业、个人独资企业名义,从事生产、经营或营利性服务的,由工商行政管理部门依照《中华人民共和国公司法》、《中华人民共和国合伙企业法》、《中华人民共和国个人独资企业法》等法律和行政法规的规定予以处罚。
从事前款规定以外其它无照经营行为的,由工商行政管理部门视其情节轻重,责令停止,并可没收相关生产经营工具、设备、物资及违法所得,处以五千元以下罚款。
第九条 违反本条例第五条规定,法律、法规没有规定的,按以下规定处罚:
(一)未将营业执照置于生产经营场所醒目位置的,由工商行政管理部门责令改正;拒不改正的,处以二百元以下罚款。
(二)不依法办理验照、年检、变更登记或其他手续的,由工商行政管理部门责令限期办理,可并处二千元以下罚款。
(三)为无照经营者非法提供营业执照、合同文本、银行帐户、发票等经营条件的,由工商行政管理部门责令改正,没收违法所得,并处五千元以下罚款。对无照经营者非法持有的营业执照、合同文本、发票等,由工商行政管理部门或有关部门予以收缴或依法处理。
第十条 当事人非法动用、调换或隐匿、转移被依法封存的财物的,由工商行政管理部门处以被隐匿、转移财物价值一倍以上三倍以下罚款;被隐匿、转移财物价值难以确定的,处一千元以上一万元以下的罚款。
第十一条 当事人拒不执行工商行政管理部门作出的罚款、没收违法所得处罚决定的,作出处罚决定的部门可将依法扣留的款项用于抵缴罚款或没收款。有扣留、封存的物品的,依照《中华人民共和国拍卖法》的规定予以拍卖或按国家有关规定处理,拍卖所得价款或变价款用于抵缴罚款或没收款。
第十二条 工商行政管理所可对个人违反本条例规定的行为作出处罚。但罚款额超过一千元或者没收违法所得额超过二千元的处罚,应由市或区、县(自治县、市)工商行政管理部门作出。
第十三条 工商行政管理部门实施处罚时,对符合登记条件的经营者,应告知其办理营业执照。
第十四条 当事人对工商行政管理部门作出的具体行政行为不服的,可依法申请行政复议或者提起行政诉讼。
当事人逾期不申请复议也不向人民法院起诉,又不履行行政处罚决定的,工商行政管理部门可依法申请人民法院强制执行。
第十五条 当事人使用暴力或以暴力相威胁,阻碍、抗拒执法行为的,由公安机关依照《治安管理处罚条例》处理;构成犯罪的,依法追究刑事责任。
第十六条 工商行政管理部门工作人员在查处无照经营行为工作中应当公正执法、文明执法。滥用职权、徇私舞弊、玩忽职守的,由有关部门给予行政处分;造成人身、财物损害的,依法予以赔偿;构成犯罪的,依法追究刑事责任。
第十七条 本条例自2001年7月1日起施行。


2001年5月25日

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

国土资源部


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

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