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

关于罪犯减刑制度的新探索——“分阶段预减刑期”的思考/汤斌

时间:2024-07-12 14:25:17 来源: 法律资料网 作者:法律资料网 阅读:9951
下载地址: 点击此处下载
关于罪犯减刑制度的新探索
…………“分阶段预减刑期”的思考

汤斌


内容提要:公正与效率是一对矛盾体,在我国几十年的司法实践中,执法问题一直在这一对矛盾中前进、发展、变革再前进再发展再变革。一直以来,理论界就在实践过程中,就减刑制度的存在与否、现行减刑制度的完善与发展、实际操作的规范化与科学化等问题,作了深入探讨和研究。建立在以监狱法为基础,以矫正理论为导向,与监狱工作实际相结合的,实施科学管理,严格执法的“分阶段预减刑期”是新时期罪犯减刑制度的新探讨,是我们追求公正与效率的统一是我们工作的方向和目标。

关键字:罪犯 减刑制度 探索

随着我国加入 WTO,政治经济体制改革正如火如荼的进行,作为监狱机关执法的一个重要方面,罪犯减刑制度的改革势在必行。以监狱法为基础,以矫正理论为导向,结合监狱工作实际,实施科学管理,严格执法的“分阶段预减刑期”这一全新的减刑模式是值得研究的重要课题。
一“分阶段预减刑期”的概述
“分阶段预减刑期”这一概念主要有三部分内容组成,一是分阶段理论,主要是指监狱机关将罪犯的改造刑期,根据不同阶段的不同改造行为标准分成改造初期、改造中期、改造后期三个阶段 ,将罪犯的改造行为纳入一个循序的、递进的改造过程中,分阶段考核,量化罪犯的改造行为理论。二是预减刑期理论,指由法律法规预先作出规定,对所有被判处有期徒刑(无期徒刑罪犯从减为有期徒刑之日起计算)的罪犯预先按既判刑期的一定比例(有待专家考证)计算出预减刑期,罪犯刑期将由实际执行刑期和预减刑期两部分组成的理论。三就是分阶段与预减刑期有机地结合起来,使罪犯知道只要在服刑的每个阶段的每一个过程中,只要能够遵守法律和监规纪律,自觉参加劳动,并主动接受教育,在经考核合格后,就可以被核减掉每个阶段内的预减刑期。相反,如果在阶段改造过程中出现违反法律行为或未达到法定的条件和要求,经评定不合格的,就将被依法撤消所取得预减刑期资格的一种减刑机制。
这一减刑制度的科学性就在于,监狱机关根据比例核算出每名罪犯可以得到的预减刑期,只要罪犯在执行期内没有违规行为发生,就能够获得预减刑期的减免,相反,如果因为违规或未达到法定条件,就不能获得预减刑期的减免。通过不断的正反激励措施,罪犯的改造动力不断增强,能够有效地防止部分罪犯受功利思想的影响,入狱后为了减刑而阶段性的表现积极或者讨好管教人员,甚至采取一些违法手段来实现减刑的“功利改造”“伪改造”行为,使罪犯能够始终如一地自觉改造,避免监管失控情况的发生。
行为管理科学认为,人做什么事情都需要有动力,罪犯的改造也不例外。从罪犯改造的功力构成来看,主要是由内驱力(正确的思想意识),推力 (管理措施环境氛围)和引力(主要是减刑,假释等奖励)三个因素构成这三股力量,经过整合形成罪犯的改造动力,促进着罪犯改造。
八十年代中期以来,我国监狱系统普遍推行了规范化管理和计分考核制度,并使之与分级管理、分级处遇和公正文明执法连成一体,形成了比较科学系统的管理激励体系,调动了罪犯的改造积极性,但在具体实践中,发现还存在着一些不足,诸如激励缺乏层次性和持续性,激励手段单一等。从而在一定程度上影响了罪犯改造,针对上述问题,我认为构建全过程、递进式的长效激励机制是有必要的,是新形势下,探索监管工作的重要课题。全程递进式的罪犯改造激励机制,主要表现为二种形式:一是目标激励就是让罪犯个人的改造目标和分监区,监狱的整体目标相结合,使其在制定和实现个人改造目标的过程中,不仅把握自身改造的努力方向,而且要符合分监区、监狱的总体目标,明确自身在整体中的价值和责任。二是数据激励,把罪犯的改选表现科学地分析,量化成具有可比性和说服力的数字。通过这一体现罪犯改造最真实性的考核数据,既可增强他们的改造信心和成就感,又能使他们看到问题和不足,在改造过程中使他们既感到有压力,又感到有吸引力,
全过程递进式的罪犯改造激励机制的运行主要还体现了三个原则:1、时效性原则,注意最佳时效,切勿坐失良机;2、持续性原则,三种激励方式要结合进行,不可间断,要自始自终地贯穿于罪犯改造的全过程;3、实效性原则,要讲究实际方式、方法,要制定符合规律的运作模式。
“分阶段预减刑期”这一全新的罪犯减刑制度就是全程激励机制的最好诠释,它既能把有效的激励模式贯穿于罪犯改造的全过程,又能在实际操作中符合科学性、规律性。
二、“分阶段预减刑期”的实施依据。
1、矫正理论依据
随着社会的不断进步,法治文明也体现出不同的特点——刑罚的预防性和矫正性。诞生于19世纪的教育刑理论观点认为,刑罚逐渐从犯罪行为转到了犯罪人方面,集中表现为刑罚的个别化,注重对犯罪人的矫正,防止他们再犯罪为宗旨。随着这一思想在20世纪的迅速发展人们开始重新为监狱定位,认为监狱不再纯碎是惩罚犯罪的国家暴力机器,其存在有自身的价值基础,是法治理性与工具理性的综合设定。监狱纪律的最终目的是矫正罪犯,而不是对其施以报复性的惩罚,政府应将良好的改造表现减少刑期的做法取代以时间计算刑期的做法,矫正理论的合理内涵使其在世界范围内受到广泛关注。
现代矫正理论的特性,就在于它强调了教育刑论的核心价值在于扩大刑罚的教育矫正机能,压缩刑罚的强制威慑机能,积极引导人们走出报复性惩罚的误区,转而采取理智和宽容的态度,即对犯罪人以理解、教育,感化和同情,并给予必要的信任和安慰,在潜移默化中使其增加对政府和热爱生活的良知,最终摆脱消极因素的桎梏重获新生。
2、法律依据
减刑是司法机关对于在服刑改造期间具备法定条件的罪犯 ,依照法定程序减轻其原判刑罚的刑罚执行变更制度。根据《刑法》第二十八条和《监狱法》第二十九条规定“对于监狱执行刑罚的被判处无期徒刑和有期徒刑的罪犯,在执行期间,如果认真遵守法规,接受教育改造,确有悔改表现或有立功表现,可以减刑,有重大立功表现的应当减刑”从这一法律明文规定。《刑法》第四条又规定:“对任何人犯罪,在适用法律上一律平等的,不允许有任何人有超越法律的特权。”这一适用法律平等原则是法律面前人人平等的宪法原则在刑事领域的具体化,它具体体现在定罪上一律平等,量刑上一律平等和行刑上一律平等为实现这一原则,修行后的刑法协调了相关的法定刑罚。可以看出,罪犯只要在服刑期间只要达到法定条件,就有获得法律奖励的权利,且这一权利具有普遍性,也就是说所有的服刑人员都可能被减刑,前提是符合法定的条件,而不应再有一些比例数字来制约和控制,用数字比例来限制这一具有普遍性的法律规定,有悖于现行法律法规的立法精神,“分阶段预减刑期”就很好地解决了这一矛盾。“分阶段预减刑期”着重说明了的是所有服刑人员都具有获得减刑的权利,至于能否减刑,主要看,其自身的改造表现是否符合标准达到法定条件。
3,激励理论依据
能够充分发挥减刑措施的激励作用。人类学的研究成果表明,人不仅有思想而且有感情,人与人之间、人与社会之间存在着一定的感情联系或者信任关系,这种感情和信任可以因某些因素而削弱或者因某些条件而增强,而且感情、信任往往成为主体积极实现某一愿望的力量源泉。恩格斯说,人们“行动的一切动力,都一定要通过他的头脑,一定要转变为他的愿望和动机,才能使他们行动起来”,这一著名论断成为现代管理科学激励机制的理论基础。20世纪后期美国心理学家弗鲁姆提出了颇具影响的“期望理论”。其要点是:人们在行动中的激励力量或动力是一个人某一行动的预期价值和该人认为将会在达到预期目标的概率之乘积。公式为:激励力量或动力二期望值x效价(这里激励力量指的是调动个人积极性、激发内部潜力的强度。期望值是据个人的经验判断达到目标的概率即可能性。效价则是达到目标后满足个人需要的价值,即主观有用性。)
从这一方程式可以看出,一个人对他所追求的目标的价值看得越大,估计能实现目标的概率越高,他的动力也就越大。或者说,激励力量越大,其内部的潜力也就越充分地调动起来。当一个人对达到某一目标漠不关心时,那么效价为零,其动力为零。而当一个人宁可不要达到某一目标时,效价为负数,动力自然为负。同样,即使效价再大,期望值为零或负值时也毫无动力可言。只有期望值和效价同时为最大时,动力才为极大值。这一理论对指导监狱刑罚执行中的减刑、假释工作具有重要意义。监管人员应当充分利用这一理论,采用多种形式手段来激发罪犯的感情、动机,把他们的期望值和效价推向最大值,这样才能促使他们的内心焕发出巨大的改造动力。
经过多年的实践,我国现行的减刑制度由于其自身存在一定的局限性,使其预期的激励作用未能很好地发挥出来。具体表现在:部分罪犯并不把减刑看成是政府对自身的宽大和宽容,不愿做艰苦的世界观改造,而采取“弄虚作假、欺上瞒下,用假改造来骗取减刑”,或者通过不正当途径获取减刑。当减刑目的达到或刑期无法再减时,就会故态重萌,消极等待刑满出狱。据调查,减刑后余刑不长或者无减刑希望的罪犯 绝大部分不积极改造,甚至出现少数人宁肯禁闭、严管也不愿参加劳动的情况。监狱又无权将减去的刑期再重新恢复,更为严重的是这种减刑前后的悬殊表现,给新人监的罪犯造成一定的消极影响,损害了减刑制度对罪犯不健康心理的矫正作用,使监狱处于不减不行减又不好的尴尬局面。实行“分阶段预减刑期”可以较好地克服这些问题。监狱通过预先告知减刑期的方式,使每一名罪犯都切切实实地看到可得的减刑期并不因为自己的罪行深重而有所区别,随之便萌发出最大期望值和效价,推动改造动力显示为极大值。不难看出,“分阶段预减刑期”既能使罪犯感受到政府对自己的宽大而主动接受改造,最大限度地发挥聪明才智来创造业绩,争取获得减刑早出狱,又能使其不敢松懈麻痹,投机取巧,确保改造积极性不发生大的波动,将减刑机制的激励作用充分发挥出来,真正做到“不想跑、不愿跑”。
4、实践依据
矫正理论,已越来越多地被各国的司法实践运用,且已经显现出较为突出的改造和矫正效果。西班牙规定,只要罪犯工作劳动两日即可折换刑期一日,意大利表现好的罪犯每服刑6个月即可减刑20日,监狱作出决定通知法院。巴基斯坦和缅甸两国实行的刑期减免制度分为普通减免和特殊减免,普通减免是凡判刑4个月以上,表现良好的罪犯,根据工种由监狱长决定每月减免6—8天的刑期,特殊减免指凡抢险救灾,劳动出色,学习努力或考试合格的罪犯,根据所判刑期获得减免2—4个月的刑期,并分别由典狱长,总监和联邦政府作出决定,南非把减刑称为赦免,一般罪犯可被赦免1/3的刑期,有前科而且被判处2年以上刑期的囚犯,赦免时间不能超过刑期的1/4,监禁期间又犯罪并被加刑的囚犯不被赦免。
在其它一些发达国家,减刑适用的标准更为灵活在美国的马里兰州,刑期折抵奖励规定罪犯在服刑期间可以获得3种刑期折抵奖:一是良好的行为,有良好行为的罪犯,每月可获得5天至10天的刑期折抵奖励,有违规时予以撤销。二是工作情况相当于参加劳动每月可得到5天的折抵奖励。三是特别活动,相当于参加文化学习,比赛得奖获得证书或学历等。可额外获得5天至10天刑期折抵奖励。
在加拿大,挣得赦免制度规定罪犯每月可挣得最多15天的赦免(10天参与计划,5天给予好行为)有轻微狱内犯罪行为或受到不良表现警告的丧失1天至2天,严重违反监视则失去当月赦免数,也可能丧失以前挣得的赦免。泰国、新加坡及香港地区的监狱犯人的减刑都是由法院直接授权给监狱审批的,凡是表现好的犯人,监狱有权给予减刑,减刑制度为原则期的1/3(屡犯为1/6)。受到减刑的犯人如表现不好在1/3(或1/6)的减刑期中扣除,直到扣完为止。
三、“分阶段预减刑期”的具体运用。
1、刑期分阶段
罪犯改造期间往往要经过初期、中期和后期三个阶段,这三个阶段可从其对改造环境的不适应,基本适应和适应及获释前的阶段来划分。每一阶段,罪犯的心理呈现出不同的状态,而不同的心理状态势必表现出不同的改造行为。不同阶段的激励有不同的效果,如何把减刑的效果发挥最大化,有待于“分阶段预减刑期”的实践运用。
改造初期,此时的罪犯改造面临着诸多的取得解决的问题。身份上由一个普通公民到一个服刑人员间不能适应,行为上由一个相对自由的人到相对不自由的人变不适应,对监狱这个陌生环境具有朦胧恐惧等感觉。此外,由于身份变化而带来的思想,心理变化也困扰着罪犯,如何使其尽快完成由不适应到基本适应这个过渡阶段,不仅取决于罪犯自身的感知能力也取决于外部力量的正确疏导。他们第一次减刑期的幅度、力度,将直接影响其下一阶段的改造心理和观念,如果达到初期郊果将为以后的改造打下坚实的基础。
改造中期。通过一段时间的改造对周围环境已由不适应到基本适应。惶惑、焦虑、紧张心理基本消除。但是由于犯罪心理的惯性存在,罪犯之间不良心理品质的影响,常态改造心理时常会尝到冲击,自我改造的动机表现出起伏。这时,需要及时地施以正反激励控制,引导矫正罪犯继续踏实改造。这时的减刑运用也是十分关键的。
改造后期。由于刑期已过大半最困难时期已经过去,对自由显得越加渴望,同时对改造投入的热情也较更高期望,获得最后一次减刑机会。与此同时也解决了部分罪犯获释前出现的低落现象。
2、考核办法
考核是监狱依照一定的标准和程序对罪犯在一定时期内的改造表现进行的综合考查和评定,考核的结果是实施罪犯“分阶段预减刑期”的依据,所以考核的公正性和科学性分别是“分阶段预减刑期”的能否顺利实施,能否发挥出最佳效果。
(1)、考核的主体是监狱机关,《监狱法》第五十六条规定“监狱应当建立罪犯日常考核制度。”由此可见,监狱是考核的实施者,是考核主体,考核权应当属于监狱。
(2)、考核的客体是罪犯在服刑期间的改造行为表现,主要包括:生活规范劳动规范,学习规范,文明礼貌规范等五个方面考核这五个方面的综合表现,其中以基本规范考核法为主。
(3)、考核依据主要是国家的有关法律法规主要是1994年颁发的《监狱法》。2004年3月2日司法部88号令《监狱服刑人员行为规范》及《规范实施细则》。待编的《新计分考核加、扣分实施细则》
(4)、考核方法:A、刑期共分为三个阶段:每个阶段采用百分考核法由基本规范、生活规范、劳动规范、学习规范、文明礼貌规范组成五个部分共一百分,每个部分为二十分;B、在阶段过程考核中,对罪犯的改造表现对照行为规范五个部分方面,根据要求凡是违反规定或不达标的均予以扣分;C、在阶段考核结束评审中,对总分低于80分认罪服法基本规范分低于12分的,不予呈报减刑期。
3、“分阶段预减刑期”实施程序
(1)罪犯在每一阶段考核结束前一个月提出申请。这是罪犯本人对自己在过程改造中对自身考核的一个评价,是罪犯行使权利的一部分,是自愿的原则。
(2)分监区组织罪犯小组评议,是实施互相监督、互相帮助、互相促进的措施,罪犯在改造中有什么漏洞,有多少违规,经过大家评议,行与不行,一目了然,使罪犯明白,自己在改造的同时有多少双眼睛在看着。
(3)分监区根据阶段计分考核结果罪犯现实表现作出阶段评审鉴定报告。一名罪犯在某一阶段的计分考核结果说明,该罪犯在阶段考核中是否达标,对照标准的客观因素罪犯的现实表现是在一定程度上弥补了客观考核的不足,罪犯中意见大、民警中意见多的,说明存在一定的问题要重新调查核实,只有综合考虑了主客观因素,才能作出公正的鉴定评审。
(4)由监狱、法院、检察院及社会监督职能机构共同组成的评审委员会对分监区提交的鉴定报告在十日内作出裁定,由多个职能机构联系执法是体现法律公平、公开、公正的一个方面。
(5)罪犯对裁定表示不服的,在按定裁定后十日内提出申诉。是体现罪犯具有申诉权的一方面。

关于矿产资源开发管理问题的通知

国家经贸委


关于矿产资源开发管理问题的通知
国家经贸委



各省、自治区、直辖市经贸委(经委、计经委):
我委上报国务院的《关于重要矿产资源开发继续实行行业管理问题的请示》(国经贸运行[1999]10号),已经国务院原则同意。现就有关问题通知如下:
根据《矿产资源开采登记管理办法》(国务院第241号令)规定,采矿权申请人在办理采矿许可证时向采矿登记管理机关提交:申请登记书和矿区范围;采矿权申请人资质条件的证明;矿产资源开发利用方案;依法设立矿山企业的批准文件;开采矿产资源环境影响评价报告;国务院地? 士蟛鞴懿棵殴娑ㄓμ峤坏钠渌柿稀I昵肟晒夜婊笄蛘叨怨窬镁哂兄匾壑档目笄诘目蟛试春凸沂敌斜;ば钥傻奶囟笾值模褂Φ碧峤还揖澄陀泄匦幸抵鞴懿棵诺呐嘉募?


1999年11月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