阿拉善盟行政公署关于印发阿拉善盟行政公署工作规则的通知
内蒙古自治区阿拉善盟行政公署
阿拉善盟行政公署关于印发阿拉善盟行政公署工作规则的通知
阿署发〔2009〕1号
各旗人民政府，盟直各部门，开发区、示范区，各大企业：
《阿拉善盟行政公署工作规则》已经2009年3月22日行政公署第一次常务会议通过。现印发给你们，请认真遵照执行。
二○○九年三月二十六日
阿拉善盟行政公署工作规则
第一章 总则
一、根据《中华人民共和国地方各级人民代表大会和地方各级人民政府组织法》、《中华人民共和国民族区域自治法》和《内蒙古自治区人民政府工作规则》,制定本规则。
二、行政公署工作的指导思想是：以邓小平理论和"三个代表"重要思想为指导，深入贯彻落实科学发展观，认真执行党的方针、政策和自治区党委、政府及盟委各项部署，全面履行政府职能，努力建设服务政府、责任政府、法治政府和廉洁政府。
三、行政公署工作的准则是，实行科学民主决策，坚持依法行政，推进政务公开，健全监督制度，加强廉政建设。
第二章 组成人员职责
四、行政公署由下列人员组成:盟长、副盟长、行政公署秘书长、各组成部门主任、局长。
五、行政公署实行盟长负责制,盟长领导行政公署的工作。副盟长协助盟长工作。
六、行政公署工作中的重大事项,必须经行政公署全体会议或常务会议讨论决定。
七、 副盟长按分工负责处理分管工作,受盟长委托,负责其他方面的工作或专项任务,并可代表行政公署进行外事活动。盟长出访期间,由常务副盟长主持工作。
八、行政公署秘书长在盟长领导下,协助处理行政公署日常工作。
九、各局局长、各委员会主任在行政公署领导下负责本部门工作。各部门决定工作中的重要事项,应事先请示分管副盟长,并定期向分管副盟长报告工作,重大事项向盟长报告。
第三章 行政公署职能
十、 行政公署全面履行经济调节、市场监管、社会管理和公共服务职能。
十一、贯彻执行国家宏观调控政策,主要运用经济、法律手段和必要的行政手段,引导和调控经济运行,促进国民经济又好又快发展。
十二、严格市场监管，推进公平准入，完善监管体系，规范市场执法行为，形成统一开放竞争有序的现代市场体系。
十三、加强社会管理，强化政府促进就业和调节收入分配职能，完善社会保障体系，健全基层社会管理体制，妥善处理社会矛盾，维护社会公平正义和社会稳定，健全各种突发事件应急管理机制。
十四、强化公共服务，完善公共政策，健全公共服务体系，增强基本公共服务能力，促进基本公共服务均等化。
第四章 行政公署决策
十五、行政公署及各部门要健全重大事项决策的规则和程序，完善群众参与、专家咨询和政府决策相结合的决策机制。
十六、全盟国民经济和社会发展规划及计划、财政预算、改革开放和经济社会发展的政策措施、重要社会管理事务、行政公署发布的规范性文件、大型项目建设等重大决策,由行政公署全体会议或常务会议讨论和决定。
十七、行政公署各部门提请行政公署研究决定的重大事项，
都必须经过深入调查研究，并经专家或研究、咨询机构等进行必要性、可行性和合法性论证；涉及相关部门的，应充分协商；涉及地区的,应事先征求意见;涉及重大公共利益和人民群众切身利益的，要向社会公开征求意见，必要时应举行听证会。
十八、行政公署在做出重大决策前,根据需要通过多种形式听取民主党派、人大代表和政协委员、社会团体、专家学者等方面的意见和建议。
十九、行政公署各部门必须坚决贯彻落实行政公署的决定，及时跟踪和反馈执行情况。
第五章 依法行政
二十、行政公署及各部门要严格按照法定权限和程序履行职责，行使行政权力。
二十一、行政公署根据经济社会发展的需要，适时制定规范性文件，及时修改或废止不相适应的规范性文件、行政措施或决定。拟定和制定与群众利益密切相关的规范性文件，原则上都要公布草案，向社会征求意见。
规范性文件实施后要进行评估，发现问题，及时完善。
二十二、提请行政公署讨论的规范性文件草案，由有关部门按计划组织起草，经行政公署法制部门审查后按程序办理。规范性文件的解释工作由行政公署法制部门负责。
二十三、以行政公署和行政公署办公厅名义制发的规范性文件，须经行政公署法制部门进行合法性审查。各部门制定的规范性文件，必须符合法律、法规和行政公署的决定和命令，并征求相关部门的意见；涉及两个及以上部门职权范围的事项，应由行政公署制定规范性文件、发布决定和命令，或由有关部门联合制定规范性文件。其中，涉及群众切身利益、社会关注度高的事项及重要涉外、涉港澳台事项，应当事先请示行政公署；部门联合制定的规范性文件，发布前须经行政公署批准。各部门制定的规范性文件要及时报行政公署备案，由行政公署法制部门审查并定期向行政公署报告。
二十四、严格行政执法责任制和执法过错追究制，有法必依、违法必究，公正执法、文明执法。
第六章 政务公开
二十五、行政公署及各部门要大力推进政务公开，按照《中华人民共和国政府信息公开条例》的要求，健全政府信息发布制度，完善各类公开办事制度，提高政府工作的透明度。
二十六、行政公署全体会议和常务会议讨论决定的事项，行政公署各部门制定的政策性措施，除需要保密的外，应及时公布。
二十七、凡涉及群众利益，需要群众广泛知晓的事项以及法律法规和行政公署规定需要公开的其他事项，均应通过阿拉善盟政府门户网站、《阿拉善政报》、新闻发布会以及报刊、广播、电视等方式，依法、及时、准确地向社会公开。
第七章 监督制度
二十八、行政公署要自觉接受盟人大工委的监督，认真负责地报告工作，接受询问和质询，依法备案规范性文件，接受执法检查，办理议案和建议；接受盟政协的民主监督，及时通报情况，办理建议案和提案，认真听取意见和建议。
二十九、行政公署各部门要依照有关法律的规定接受司法机关实施的监督，同时要自觉接受监察、审计等部门的专门监督。对监督中发现的问题，要认真查处，限期整改，并向行政公署报告。
三十、加强行政系统内部监督，健全政府层级监督制度。行政公署及各部门要严格执行行政复议法和规范性文件备案审查制度，及时撤销或修改违反法律、法规的规范性文件，纠正违法或者不当的行政行为，并主动征询和认真听取各旗（开发区）及其部门的意见和建议。
三十一、行政公署及各部门要重视人民群众来信来访工作，进一步完善信访制度，确保信访渠道的畅通，维持信访秩序。坚持属地管理、分级负责，谁主管、谁负责，依法、及时、就地解决问题与疏导教育相结合的原则，严格实行信访责任制和责任追究制，有效防范和化解社会矛盾。
三十二、行政公署及各部门要接受新闻舆论和群众的监督。对新闻媒体报道和各方面反映的重大问题，各部门要积极主动地查处，限期整改，并向行政公署报告。
三十三、行政公署及各部门要推行行政问责制度和绩效管理制度，明确问责范围，规范问责程序，严格责任追究，提高政府执行力和公信力。
第八章 廉政建设
三十四、行政公署及各部门要从严治政。对职权范围内的事项要按程序和时限积极负责地办理，对不符合规定的事项要坚持原则不得办理；对因推诿、拖延等官僚作风及失职、渎职造成影响和损失的，要追究责任；对越权办事、以权谋私等违规、违纪、违法行为，要严肃查处。
三十五、行政公署及各部门要严格执行财经纪律，规范公务接待。要艰苦奋斗、勤俭节约，切实降低行政成本，建设节约型机关。
三十六、行政公署组成人员要廉洁从政，严格执行党中央、国务院、自治区党委、人民政府及盟委有关廉洁自律的规定，不得利用职权和职务为本人或特定关系人谋取不正当利益；要严格要求亲属和身边工作人员，不得利用特殊身份谋私利。
第九章 会议制度
三十七、行政公署实行全体会议、常务会议和盟长办公会议制度。
三十八、行政公署全体会议由盟长召集并主持,一般每半年召开一次,由行政公署组成人员组成,根据需要可安排行政公署副秘书长、旗长（开发区管委会主任）和有关部门、单位主要负责人列席会议。
全体会议的主要任务是:传达贯彻党中央、国务院, 自治区党委、人民政府和盟委、盟人大工作委员会的重要决定事项;通报全盟经济社会发展形势和有关重要情况,安排部署行政公署的重要工作;讨论决定行政公署工作中的重大事项等。
三十九、行政公署常务会议由盟长或盟长委托常务副盟长召集并主持,一般每月召开两次以上。行政公署常务会议由盟长、副盟长、行政公署秘书长组成。行政公署副秘书长、行政公署办公厅副主任列席会议,根据需要可安排有关旗（开发区）、部门和单位负责人列席会议。
常务会议的主要任务是:研究审议全盟国民经济和社会发展的长远规划、年度计划、盟本级财政预算;讨论通过以行政公署或行政公署办公厅名义制发的规范性文件;研究决定行政公署日常工作中涉及全面工作以及需要统筹安排的重要事项；讨论决定行政公署工作中涉及重大公共利益、人民群众切身利益等重要事项。
四十、 盟长办公会议由盟长或盟长委托常务副盟长召集并主持,不定期召开。盟长办公会议由盟长、副盟长、行政公署秘书长组成,根据需要可安排行政公署相关副秘书长及有关旗（开发区）、部门和单位负责人列席会议。
盟长办公会议的主要任务是:通报或听取全区性会议精神、重要访问和考察情况;听取旗（开发区）、部门及有关方面重要工作汇报; 研究提交行政公署全体会议、常务会议讨论的事项; 研究盟长或常务副盟长提出的其他重要事项。
四十一、根据工作需要，行政公署可召开专题会议，研究、协调、处理工作中的具体问题和专项事宜。行政公署专题会议一般由副盟长召集并主持，有关部门、单位负责人参加。
行政公署专题会议凡涉及资金、项目、机构编制安排的,严格按照规定的审批权限和程序办理。
四十二、提请行政公署全体会议讨论的议题,由行政公署常务会议或盟长确定;提请行政公署常务会议和盟长办公会议讨论的议题,由盟长确定,或由分管副盟长、行政公署秘书长协调审核后报常务副盟长确定。行政公署全体会议、常务会议和盟长办公会议纪要由行政公署秘书长审定，报请常务副盟长或盟长签发。行政公署专题会议纪要由分管副秘书长审定，报请副盟长签发。
未经上述程序的事项,不得列入会议议题,不得临时动议。
各类会议的会务组织工作由行政公署办公厅负责,会议文件于会前送达与会人员。对议题有意见或建议的，可在会前提出。行政公署全体会议、常务会议讨论决定的事项,宜于公开的,应及时报道。新闻稿须经行政公署相关副秘书长审核,行政公署秘书长审定签发,如有需要报请盟长或常务副盟长审定签发。
四十三、行政公署领导同志不能出席全体会议、常务会议或盟长办公会议,需向盟长或主持会议的常务副盟长请假;其他人员不能出席或列席会议,需向行政公署秘书长请假。行政公署办公厅定期通报各类会议的参加和请假情况。
四十四、行政公署及各部门召开的工作会议,要减少数量,控制规模,严格审批,并尽可能采用电视电话会议、视频会议等快捷、节俭的形式召开。应由各部门召开的全盟性专业工作会议,不得以行政公署或行政公署办公厅名义召开,不邀请各旗人民政府（开发区）主要负责人出席。确需邀请的,须报行政公署批准。
第十章 公文审批制度
四十五、各地区、各部门报送行政公署审批的公文,应当符合《国家行政机关公文处理办法》的规定。公文审批按照制度化、规范化、程序化和行政公署领导同志分工负责的原则办理,并根据需要由行政公署领导同志转请其他领导同志核批，重大事项报盟长审批。
四十六、各部门提请以行政公署或行政公署办公厅名义制发的规范性文件,经行政公署法制部门审查、分管副盟长审核、行政公署常务会议讨论通过后,报盟长或常务副盟长签发;行政公署发布的决定、命令、指导性意见,经分管副盟长审核后,报盟长签发;向自治区人民政府的请示、报告性公文,经分管副盟长审核后,报盟长签发，盟长外出期间由主持工作的常务副盟长签发。以行政公署名义就某一事项作出安排和部署,或就某一请示事项作出批复,以及致函不相隶属机关商洽工作、询问和答复问题、请求批准等,由副盟长或行政公署秘书长签发；涉及两个或两个以上副盟长分管工作的公文,经有关副盟长会审后签发;属于重大问题,经分管副盟长审核后,由盟长或常务副盟长签发。
四十七、以行政公署办公厅名义发文,经行政公署分管副秘书长审核并向分管副盟长报告后, 行政公署秘书长复审签发,也可由行政公署秘书长委托行政公署分管副秘书长向分管副盟长报告后审核签发。
四十八、各旗（开发区）、各部门报送行政公署的公文,由各旗（开发区）、各部门的主要负责同志或主持工作的负责同志签发,并严格按程序报送,不得直接报送行政公署领导同志个人。凡报送行政公署的公文,统一由行政公署办公厅签收、登记,按行政公署领导同志分工呈批。
四十九、各旗（开发区）、各部门向行政公署报送的请示性公文,需转有关职能部门研究办理时,凡涉及两个或两个以上部门的,由行政公署办公厅明确主办部门。主办部门主要负责人要主动与有关部门协商,属于部门职权范围内的事项,由主办部门牵头商请有关部门直接答复报文单位。需行政公署审批的事项, 主办部门要及时商请有关部门提出办理意见。
五十、各旗（开发区）、各部门代行政公署和行政公署办公厅拟制的文稿, 应将申请发文的理由、背景和有关说明材料一并报送。凡内容涉及其他部门职权范围的, 主办部门应主动与有关部门协商,取得一致意见后方可报送文稿。如有分歧,主办部门的主要负责人应出面协调; 仍不能取得一致意见的, 主办部门可列明各方理据, 提出办理建议, 并与有关部门会签后报行政公署研究决定。
五十一、提高公文处理工作效率,确保公文运转及时、准确。各旗人民政府（开发区）、各部门报送行政公署和行政公署办公厅的各类公文,原则上应在15个工作日内办结,其中转有关部门直接答复或征求意见的,各部门要在10个工作日内办结,逾期不能办结的,应书面说明情况和理由,否则，行政公署办公厅将予以通报。紧急公文随到随办。
五十二、行政公署及各部门要进一步精简公文。部门职权范围内的事项,由部门自行发文或联合发文,不得要求行政公署批转或行政公署办公厅转发;确需行政公署批准的事项并经批准后,可加注"经行政公署同意"字样,由部门自行发文。行政公署及行政公署办公厅一般不转发自治区委、办、厅、局和自治区非常设机构的公文。
各旗（开发区）、各部门报送行政公署或行政公署办公厅的公文,凡不符合《国家行政机关公文处理办法》及有关规定的,退回原报文单位重新办理。
第十一章 督查落实
五十三、行政公署对国家、自治区和盟重大决策和重要工作部署实行督查工作制度。各旗人民政府（开发区）、行政公署各部门对行政公署部署的工作负有落实责任。要进一步健全和完善运转协调、责任明确、反馈及时的督查工作机制,确保行政公署工作高效有序运行和政令畅通。
五十四、督查工作的重点是党中央、国务院、自治区党委、政府和盟委、行政公署重大决策和重要工作部署以及会议决定的重大事项的落实情况;领导批示交办事项的落实情况;全国和自治区人大代表建议、自治区和盟政协提案的办理情况和办复质量;群众普遍关心、社会关注的热点、难点问题的解决情况等。
五十五、督查工作要注重实效,采用网络督查、跟踪督查、现场督查、暗访督查等多种形式,全面深入准确地了解情况,报实情、重实效。督查工作实行责任制、考核制和通报制。以行政公署名义开展的督查工作,由行政公署督查室综合协调、组织实施。
第十二章 作风和纪律
五十六、行政公署组成人员要坚决贯彻执行党和国家的路线方针政策以及工作部署，严格遵守纪律，有令必行，有禁必止。
五十七、行政公署组成人员必须坚决执行行政公署的决定,如有不同意见可在内部提出,在没有重新做出决定前,不得有任何与行政公署决定相违背的言论和行为;代表行政公署发表讲话和文章,以及个人发表涉及未经行政公署研究决定的重大问题及事项的讲话和文章,事先须经行政公署同意。
五十八、行政公署各部门发布涉及行政公署主要工作部署、经济社会发展重要问题，与群众利益密切相关事项的信息，要经过严格审定，重大情况要及时向行政公署报告。
五十九、行政公署组成人员要严格遵守保密纪律和外事纪律，严禁泄漏国家秘密、工作秘密或者因履行职责掌握的商业秘密，维护国家安全、荣誉和利益。
六十、行政公署组成人员要做学习的表率，行政公署及各部门要努力创建学习型机关。
六十一、盟长、副盟长要深入基层,调查研究，指导工作,解决实际问题。下基层要轻车简从，减少陪同和随行人员,简化接待。
六十二、行政公署各部门要精简公务活动，严格执行领导同志出席各类会议和公务活动的有关规定，确需行政公署领导参加的，应履行报批手续。
六十三、行政公署组成人员要严格执行请销假制度。行政公署副盟长、行政公署秘书长离开巴彦浩特出访、出差、学习、因病休养或休假,须事先报告盟长，由行政公署办公厅通报行政公署其他领导同志。行政公署各部门主要负责人离开巴彦浩特,应事先向行政公署办公厅报告，由行政公署办公厅向行政公署领导同志报告。
行政公署直属特设机构、直属机构、办事机构、直属办事单位适用本规则。
关于统一图幅理论面积与图斑椭球面积计算要求的通知
国土资源部
关于统一图幅理论面积与图斑椭球面积计算要求的通知
（国土调查办发〔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
