补种苗木变更申请（通用9篇）
篇1：补种苗木变更申请
新建铁路海南西环铁路海口至凤凰机场段XHZQ-6标一分部补种乔木变更申请
尊敬领导您们好！
武汉迪园园艺工程有限公司在建“新建铁路海南西环铁路海口至凤凰机场段XHZQ-6标一分部绿化工程”。应业主要求，对全线主线路堤边沟平台外侧超过2.5M的平台补种乔木，根据设计图纸要求我标段补种杜英。因本地区天气干燥并少雨，日照时间长，土壤贫瘠且以砂石为主，且路基全线封闭灌溉难度极大。杜英喜湿润环境，且光照时间不宜过长，土质要求偏高，不适合在路基内生长存活，特我司申请对补种苗木变更为更适合在本地区耐高温耐旱且生长存活较高的火焰木、发财树、盆架子，且规格与原设计规格一致。
特此汇报 此敬！
武汉迪园园艺工程有限公司
2015年11月25日
篇2：补种苗木变更申请
提要：绿化工程（主要指种植土厚度、地形堆坡、回填土质）概况，苗木死亡种类及死亡情况介绍、补种时间，种植穴大小不合格的方案
基层土不透水情况处理方案
种植土厚度不够处理方案
篇3：补种苗木变更申请
有观点认为:本案中,法院于2012年判决:“该房扩建部分未进行产权登记,不判决房屋归属,但属于家庭成员共同财产”。据此, B某理当是该房扩建部分的共有人。而当A某提出房屋所有权变更登记申请后,登记机构在未对该房现所有权状况审核的情况下为A某单独颁证,造成了房屋登记薄记载的情况与事实不符。故而认为登记机构登记有误, 理当严格审查申请主体的变化状况。
另一种观点认为:《房屋登记技术规程》中明确规定了房屋所有权转移登记申请材料表中含“婚姻关系证明”,而房屋所有权变更登记申请材料表中没有这项要求。本案中,登记机构在接受A某提出的房屋面积变更申请后,审核了A某的房屋产权证、建房执照和测绘部门出具的测绘报告,并在A某提交申请时就有关该房权利状况是否单独所有对A某进行了询问。A某填写的询问笔录为:“该房为其单独所有。”故登记机构将房屋所有权变更登记在A某一人名下无误。再因登记机构无权对实体民事法律关系加以审查,故而,在涉及对权利是否存在的审核时,登记机构无需审核实际的权利状态,而仅就登记簿记载的权利状态加以审查即可。总之,登记机构无权对申请人进行物权的确认。本案中,登记簿记载的该房原权利状态为A某单独所有,且建房执照登记在A某一人名下。登记机构在为A某办理房屋所有权变更登记时没有义务审核申请主体的变化状况。
笔者赞同第二种观点。
首先,如果申请人在登记申请时刻意隐瞒权利状况改变事实,如本案中,A某隐瞒了法院对该房扩建部分判决的事实,登记机构无法查证申请主体是否发生变化。《房屋登记办法》实施后,如果房屋在事实上存在共有人,但因当事人没有据实申请登记, 而导致登记名义人与实际权利人有差异时,登记机构不应对此承担责任。对于实际共有人因此受到的损失,只能通过其他法律途径要求处分人承担相应的赔偿责任。
篇4：补种苗木变更申请
有观点认为：本案中，法院于2012年判决：“该房扩建部分未进行产权登记，不判决房屋归属，但属于家庭成员共同财产”。据此， B某理当是该房扩建部分的共有人。而当A某提出房屋所有权变更登记申请后，登记机构在未对该房现所有权状况审核的情况下为A某单独颁证，造成了房屋登记薄记载的情况与事实不符。故而认为登记机构登记有误，理当严格审查申请主体的变化状况。
另一种观点认为：《房屋登记技术规程》中明确规定了房屋所有权转移登记申请材料表中含“婚姻关系证明”，而房屋所有权变更登记申请材料表中没有这项要求。本案中，登记机构在接受A某提出的房屋面积变更申请后，审核了A某的房屋产权证、建房执照和测绘部门出具的测绘报告，并在A某提交申请时就有关该房权利状况是否单独所有对A某进行了询问。A某填写的询问笔录为：“该房为其单独所有。”故登记机构将房屋所有权变更登记在A某一人名下无误。再因登记机构无权对实体民事法律关系加以审查，故而，在涉及对权利是否存在的审核时，登记机构无需审核实际的权利状态，而仅就登记簿记载的权利状态加以审查即可。总之，登记机构无权对申请人进行物权的确认。本案中，登记簿记载的该房原权利状态为A某单独所有，且建房执照登记在A某一人名下。登记机构在为A某办理房屋所有权变更登记时没有义务审核申请主体的变化状况。
笔者赞同第二种观点。
首先，如果申请人在登记申请时刻意隐瞒权利状况改变事实，如本案中，A某隐瞒了法院对该房扩建部分判决的事实，登记机构无法查证申请主体是否发生变化。《房屋登记办法》实施后，如果房屋在事实上存在共有人，但因当事人没有据实申请登记，而导致登记名义人与实际权利人有差异时，登记机构不应对此承担责任。对于实际共有人因此受到的损失，只能通过其他法律途径要求处分人承担相应的赔偿责任。
其次，依照相关法规，登记机关既没有权利也没有义务对当事人申请登记的涉及实质权利的内容进行调查，也无权对当事人的财产法律关系进行变更，因为这些做法违背权利人对自己财产的意思自治的原则。如果当事人之间发生物权归属及其内容的争议，依据 《物权法》第32条规定，物权的确认应当通过法院或仲裁机构进行。登记机构只要审核申请人递交的材料、查阅登记簿、在必要的时候进行实地查看就履行了登记机构审核的职责。本案中，申请人A某递交的材料与登记簿中记载的主体、客体一致，申请房屋面积变更的材料即登记簿中曾记载的建房执照资料和测绘部门重新勘察出具的测绘报告，符合登记程序。基于此，笔者认为，本案中登记机构无需在申请人申请房屋所有权变更登记时对其进行申请主体变化状况的审核。而针对B某提出的异议，建议其作为利害关系人申请异议登记，并向法院申请民事诉讼。
篇5：林中空地补植补种探讨
林中空地泛指森林里树木中间及周围的无林地。根据国家第八次森林资源清查的数据来看,山西省森林有林地面积282.41万hm2,其中人工林面积131.81万hm2,按照1.0%的比例估算,那么林中空地面积为2.82万hm2,要将这些林中空地全部改造成有林地,森林面积和蓄积量就会得到增加,森林质量就会得到提高。但同时,如此大的改造面积也给林业工作者带来巨大的考验,怎么全面完成这项任务呢?习近平总书记的讲话给我们指明了思路,他说:“我们说一张蓝图抓到底,不仅需要科学决策,需要思想境界。什么思想境界?就是功成不必在我”。这段论述深刻揭示了长远利益与眼前利益、革命利益与个人利益的辩证统一关系,对于解决在工作中存在的急功近利倾向,具有重要指导意义和现实针对性。尤其是针对林中空地补植补种工作,不能急于求成,要把这项工作作为长期攻坚战来完成,不必在林局、林场的一两任领导的任期内完成,要科学规划、因地制宜、分步分片实施。只有一批批林业人的坚持不懈的努力才能完成这项艰巨的任务,任重道远。
2结合森林抚育开展林中空地补植补种
对于林中空地的改造,国家新的《森林抚育技术规程》中有了规定,在第6条各种抚育方式适用的条件,6.5补植,符合以下条件之一的,即可采用补植。(1)人工林郁闭成林后的第一个龄级,目的树种、辅助树种的幼苗幼树保存率小于80%;(2)郁闭成林后的第二个龄级及以后各龄级,郁闭度小于0.5;(3)卫生伐后,郁闭度小于0.5的;(4)含有大于25m2林中空地的;(5)立地条件良好,符合经营目标的目的树种株数稀少的有林地;符合条件(5)的,应结合生长伐进行补植。在第7条控制指标,7.6(1)补植树种要选择能与现有树种互利生长或相容生长,并且其幼树具备从林下生长到主林层的基本耐阴能力的树种。对于人工用材林纯林,要选择材质好、生长快、经济价值高的树种;对于天然用材林,要优先补植材质好、经济价值高、生长周期长的珍贵树种或乡土树种;对于防护林,应选择能在冠下生长、防护性能良好并能与主林层形成复层混交的树种。(2)用材林和防护林经过补植后,林分内的目的树种或目标树种株数不低于每公顷450株,分布均匀,并且整个林分中没有半径大于主林层平均高1/2的林窗;(3)不损害林分中原有的幼苗幼树;(4)尽量不破坏原有的林下植被,尽可能减少对土壤的扰动;(5)补植点应配置在林窗、林中空地、林隙等处;(6)成活率应达到85%以上,3年保存率应达80%以上。
根据新的森林抚育技术规程,2015年太行山林局石源林场在森林抚育工作中,开始尝试在林中空地进行补植补种。具体如下:对地块的选择上,在槐树沟进行试点,在进行作业设计时,针对性选择有林中空地的小班做出补植要求。补植点设置在小班下部的林边空地、小班中间面积大的林窗和林隙,坚持适地适树的原则,补植树种为乡土树种。这次选择了乔木,油松、侧柏容器苗和山杏裸根苗木,以后也可以尝试灌木;整地采用小鱼鳞坑,规格为40cm×30cm×20cm,为减少对原有植被的破坏,株行距确定为3.0m×4.0m,品字形排列,每公顷840株。如遇到特殊地形,也可以采用不规则整地,株距和行距不低于2.0m,不超过5.0m,裸露的岩石可空出,根据实地情况可垒石坑,从附近取土填坑,保证每公顷密度不少于840株。严把苗木质量,对苗木有严格的要求,达不到标准的不能上山;栽植上下功夫,挖坑培土,扶正踏实,保证苗木成活率。隔年踏查,根据情况再进行补植,保证保存率。
3林中空地补植补种中存在的问题及解决的办法
3.1存在问题
3.1.1森林抚育作业设计难度加大,一是外业调查难度加大,以前是下中下布设标准样地,现在需要全部踏查,时间上和设计人员的体力消耗更多;二是林中空地的面积测量测算也给设计带来了新问题,补植补种面积不准,很难调绘在1∶10000的地形图上;三是林中空地补植补种的要求,实地情况很复杂,很难简单地以一种造林模式来概括,因此设计样式多,投资测算难。
3.1.2投资不足,尤其是面积大的林中空地补植补种上,正常的造林投资在每公顷4 500元~7 500元,而森林抚育每亩投资只是100元,即使是林局或林场有些配套的资金,但分配到每亩上也就是几块钱,基本上不起作用。实际上由于地形复杂,道路更远,分散不集中,林中空地的补植补种反而要更多的用工量,投资不足,更是影响了基层实施单位对林中空地改造的热情,林局和林场的领导们对此项工作也重视不起来。
3.1.3实际执行难,验收标准难统一。森林抚育技术规程也是刚开始实行,林中空地补植补种没有完整详细的实施方案或规范,也没有什么经验和教训,缺乏制度、机制及实施细则;验收难在面积不好计算,实地检查也很困难,不可能每个补植点都能像整地造林那样去验收,也不能像普通的森林抚育那样以标准样地验收。抚育是当年验收,造林却是3年验收,时间也不统一。
3.2解决方法
3.2.1从思想上高度认识,做出长期规划。行动上重视起来,由省林业主管部门牵头成立一个专门的林中空地补植补种领导小组,由主要领导挂帅,各市和各省直国有林管理局也一样,上到省,下到林场,拧成一股绳,思想统一形成合力。
3.2.2兵马未动粮草先行,多方筹措资金。针对林中空地补植补种有专项资金,加大投资力度,保证基层单位在实施过程中不贴钱或少贴钱。
3.2.3开展试点,总结经验和教训,形成实施方案和细则再逐步推开,有统一性,也要有区域的差异性。可在北、中、南分区选择试点推开,各市选择3个林场,省直林局选择3个林场,通过试点去形成适合各自的方案,然后再扩大试点范围,不断总结完善,3年内完成试点再全省推开。在树种配置上,以乡土树种为主,可乔木也可灌木;在整地规格上,密度要求以及栽植上形成实施细则,再用10年左右的时间推广完成。
3.2.4推广新科技。林中空地补植面积测算用GPS绕测,并根据补植面积给外业调查人员增加补助,调动其工作积极性。
3.2.5完善验收机制。实行责任追究制,林场自查,林局或林业局验收,各市和省直林局互查,加上省级抽查,明确责任追究到底,面积不实追究设计人员和林场实施人员责任,质量、成活率和保存率不达标,追究场领导和局领导责任,不只是当年或次年验收,每隔3年,再验收一次保存率,实地验收和航片、卫星结合起来,Google地图每5年更新一次,这些都能作为很好的参考,狠抓细抓验收,把这项工作做实做到位。
3.2.6结合小流域治理。完成林中空地的补植补种需要森林抚育,也要跳出森林抚育。现在省林业厅和林局提出一条沟的治理理念,而且也在试点推广,能抚育则抚育,能造林则造林,能改造则改造,因地制宜,把林边空地和林隙的补植补种放到造林设计中去,有效提高资金投资成效,弥补不足。
上述只是本人在森林抚育林中空地补植补种上的一些粗略的见解,希望大家在今后的工作中,不断尝试,发现问题,探讨研究,形成好的经验,逐步推广,一定能努力完成好林中空地改造这项艰巨而伟大的任务。
摘要：林中空地泛指森林中间及周围无林木生长的空地。介绍了山西林中空地的现状及改造的必要性,分析了林中空地补植补种中存在的问题,提出了结合森林抚育开展林中空地补植补种的建议。
篇6：补种苗木变更申请
执行立案没多久，童兰就在一次极限运动中意外死亡。明芳发现，童兰的存款、股票、基金、房产均被其父母继承。于是在律师的建议下，明芳向法院提出执行异议申请，申请变更童兰的父母为被执行人，要求其父母在继承童兰遗产的范围内履行还款的义务。
执行异议申请立案后，执行法官经过审查，裁定变更“明芳申请执行童兰一案”的被执行人为童兰的父母。由于童兰的父母不配合履行执行义务，执行法官依法采取了相应的执行措施，通过冻结划拨童兰银行账户存款及查封拍卖童兰房产的方式为明芳成功执结了100万元的借款及相应的迟延履行金。
[评析]根究相关司法解释规定，若债务人死亡，债务人生前拥有合法财产且已由继承人继承的，则债权人可以向法院申请追加继承人为被告，要求继承人在继承遗产的范围内承担偿还责任。本案中，童兰死亡后，其父母继承了童兰的遗产，明芳即可以申请变更童兰的父母为被执行人。据此，明芳只需按照法律程序的规定再立一个执行异议案件即可助其在童兰死亡后实现债权。
篇7：补种苗木变更申请
“责令补种”这一行政处罚措施因其符合最朴实的补偿思维, 最大限度地恢复了被破坏的生态环境, 符合人类与自然和谐发展的要求而被立法者、执法者及社会大众所接受, 越来越多地适用在现实生活中。此种规定虽在法律规定的静态层面上很好地遵循了科学规律, 但仔细推敲却缺乏对现实情况的进一步细致理性思考, 不可避免地给现实中的林业执法实践带来诸多困难, 主要体现在:
(一) 无明确规定补种地点
法律条文中没有明确地规定补种地点使得林业部门的检查监督工作无据可依, 为执法工作带来一定的障碍, 也成为违法者拒不补种的借口之一。关于补种地点我国实践中的一般做法为:在国有或者集体所有的荒山荒地上进行补种, 补种者将来也不享有林木的所有权。现阶段我国大部分地区都已进行集体林权改革, 新的形势下, 此种将补种树木植入国有或集体所有林地的做法过于狭隘, 已经不能合理地匹配现有体制, 需要进一步的完善与修订。
(二) 无具体确定的补种树种要求
林业工作人员在处罚决定做出时, 受到植物生长温度、气候、土地条件、树种特征等多种生物因素干扰, 难以确定具体的补种要求, 补种处罚书中缺少补种行为的核心标准, 继而造成林业部门无依据对其行为进行检查验收, 因此在具体的执法实践中林业管理部门很难判断补种的执行状况, 形成监督工作不可突破的障碍。
(三) “代为补种”规定的尴尬
根据我国《森林法》第三十九条的规定, 由林业主管部门做出“责令补种”的处罚决定, 同时规定在违法者拒不补种或补种不符合国家相关规定的, 由林业行政主管部门“代为补种”, 所需费用由违法者支付。如此规定使得林业主管部门陷入尴尬局面, 一方面在做出责令补种的处罚决定后, 如果违法者不执行, 林业部门本身不但没有强制执行的权力, 而且需要其自己进行树木的补种, 不仅没有对违法者实行惩戒, 却加重了林业部门自身的责任;另一方面, 规定要求补种费用需要有林业行政管理部门先行垫付, 后向违法者追讨。立法者坚守自然资源恢复原则的善良初衷并没有考虑到我国特定的社会现实状况, 为后续林业部门的执法、落实工作带来新的困扰:
一是, 我国行政部门并没有特定的款项作为盗伐林木行为补种树木的预算, 因此筹措补种树木资金比较困难;
二是, 即使林木部门先行补种, 补种款项的回收成为一个新的问题。首先, 现实生活中违法主体确定的难度较大。由于山林的覆盖面积广, 有些林户的承包林、自留山等分布在距住所比较远的位置, 交通不便利;加之树木成材周期长, 人工作业少等生长特点, 林户不能做到及时、全面的看护, 一旦发生盗伐林木行为, 未能及时发现, 当场抓获违法者的概率微乎其微。加之上述盗伐林木的行为特点, 越来越多作案行为呈现出计划预谋、团伙作案、分工明确等特点, 盗伐者可以利用较为先进的技术在作案后迅速逃离, 加大了森林公安的破案难度。其次, 盗伐者多为低收入群体, 以盗伐林木应急、治穷的案件, 罚款执行困难重重。
二、“责令补种”法律条款的修订与完善
针对当前“责令补种”条款面临的重重困境, 以我国特有的林业现状为依据, 结合我国基层的林业工作经验, 提出以下改进措施:
完善立法, 明确“责令补种”树木地点、树种、数量的方法, 为林业执法部门的处罚决定、检查监督提供坚实的法律依据。
(一) 补种树木地点的确定方法
补种地点的确定十分必要, 它不仅是违法者执行处罚的重要依据, 同时也可以体现出法律责任的根本属性和科学内涵。
理想条件下补种的树木应该补种在盗伐林木的原处, 此思路可以最为直观地补偿受害人, 但树木的生长必须遵照客观的生物规律, 如果被盗伐地点因自然条件、生长环境等不适宜补种新林木的, 可以责令违法者在被害人所承包经营的林地范围内, 由被害人行使补种地点的决定权。这样既可以实现遵循科学规律补种树木, 提高林木成活率, 同时也可以发挥被害人在行政处罚中的作用, 最大限度地满足其对于补偿的要求。
(二) 补种树种确定方法
法律虽未对补种树种做出明确的要求, 但将此权力交与林业行政主管部门, 笔者建议将此权力下放于被害人, 由被害人自行决定。由于被害人是被盗伐树木的所有权人或使用权人, 是日常种植、维护林木生长的监护人, 其最了解案发地适宜种植的树种种类, 同时, 树种的选择直接关系到补种的效果和经济补偿的具体数额, 被害人趋利避害选择同保持水土、维护国家生态安全有着同样的价值取向。
(三) 取消“林业部门代为补种”的相关规定
正如上文提及, “林业部门代为补种”的规定无法执行有损执法机关的权威形象, 同时, 代为补种后, 当事人交不出补种费用或拒绝缴纳时, 林业主管部门无法收回, 唯一追索途径只有向法院申请强制执行数额较小的补种费, 此种方法不仅周期长、费用高, 且多数盗伐者都是贫困的农民, 本来就是因为贫穷而盗伐树木, 执行成功率较低。
(四) 确定明确的补种监督检查机制
现行法律将“责令补种”的监督检查权赋予了林业行政主管部门, 笔者在此建议将该权利部分赋予被害人。由于其是盗伐林木行为的直接受害者、财产补偿的直接受益人, 对于补种树木的成活、生长及成才情况最为关注与了解, 他们可以对违法者的补种行为进行持久而有效地监督, 一定程度上减轻了林业行政主管部门的工作负担。
随着盗伐林木之“责令补种”处罚方式法律条款的不断明确与完善、立法思路日趋清晰与成熟, 其会在未来的盗伐林木行政处罚方式中发挥越来越重要的作用。
摘要：“责令补种”是我国《森林法》法律责任规定的行政处罚措施之一, 由于其代表的科学内涵而得到林业行政法的普遍适用, 但其在现实的司法实践中却因其立法漏洞而面临严峻的考验, 本文通过以“责令补种”条款在法律适用中的具体问题分析为基础, 提出完善“责令补种”条款之构想。
篇8：补种苗木变更申请
摘 要：针对医学图像CT图像像素不均匀对图像局部分割算法影响较大的问题，提出一种基于Lagrangian粒子增强补种算法的混合水平集医学CT图像分割算法。首先，针对局部图像的非均匀性，通过在计算水平集公式前先计算Lagrangian标记粒子来重建内嵌交接界面，从而提高水平集算法的质量守恒特性；其次，针对传统粒子方法在处理界面奇异性和复杂几何相关问题上的不确定性，通过增加速度矢量和单位法向量来促进奇异点和拓扑变化点速度场的收敛；最后，通过在合成数据测试集和真实CT图像上的仿真测试表明，所提算法在边缘分割收敛精度及运算速度上均要优于对比算法。
关键词：补种算法；奇异性；水平集； CT图像；分割
中图分类号：TP391 文献标识码：A
1 引 言
图像处理方向学者Osher与Sethian等[1]联合研究并提出了水平集（Level Set， LS）分割算法，主要思路是将低维闭合二维曲线映射到高维的水平集上进行处理，从而实现了图像分割算法稳定性的增强。尽管如此，图像界面奇异性和复杂几何相关性对水平集算法影响仍然很大，如何处理好上述问题，对于提高水平集分割算法在CT、X-ray等医学成像中的应用效果，至关重要，有利于增强算法在处理对比度差、像素不均匀、边界模糊等问题并存的医学影像成像分割问题上的表现[2-3]。
针对锋利边缘和极端变形等情况，不同学者提出各种改进方式，如文献[4]基于欧拉-拉格朗日粒子的方法对水平集算法进行改进，该方法在保持分割界面质量上表现不错。此外，还有局部水平集方法[5-6]、间断Galerkin方法[7-8]、梯度增强水平集方法[9-10]等。上述方法在处理锋利边缘和极端变形等情况有一定效果，但是当像素不均匀或交接界面存在奇异点或拓扑变化时，上述算法效果不佳。
对此，本文在文献[4]基础上，采用Lagrangian粒子增强补种算法与水平集算法相结合，通过增加速度矢量和单位法向量来促进奇异点速度场的收敛，提出了一种基于Lagrangian粒子增强补种算法的混合水平集医学CT图像分割算法（LPRLS）。
2 基于粒子的水平集算法
2.1 水平集算法
2.2 混合粒子水平集方法
该方法融合了欧拉水平集算法[16]和拉格朗日粒子算法[17]，采用拉格朗日无质量粒子来纠正图像界面待求解区域的水平集函数。两组粒子被随机放置在窄带的界面附近，并吸引界面到正确的一边，其中阳性粒子趋向>0的一边，而阴性粒子趋向≤0的一边，格朗日粒子将在给定速度下产生对流：
式中，xp为粒子位置；uxp为粒子的速度向量。由于演化方程的耗散缺失，因此流动的特征信息被完全保存。然后，利用三阶TVD Runge Kutta方法来求解时间导数的演化方程。
每个粒子都有一个位置和范围半径，用来弥补因水平集函数的局限性所引起的的位置误差。每个粒子的半径由基于网格的最大值和最小值进行界定。根据公式（6）粒子半径的最大和最小界定值，可由下式进行计算：
在界面清晰像素均匀的较容易解决图像表面地区，水平集函数的解是足够准确的，并且粒子不漂移距离交接界面区域过远。然而在像素不均匀、界面模糊的不容易解决区域，水平集函数计算会产生质量损失，粒子会漂移距离交接界面过远。界面上逃脱的正（或负）粒子，基于网格点、粒子半径和粒子符号可定义局部水平集函数。水平集函数集的纠正可通过对比网格水平集函数值和逃脱粒子的局部水平集函数值实现。基于网格水平集函数的符号，则逃脱阴阳离子的局部水平集函数可分别定义如下：
通过将>0和≤0的区域分别与网格水平集函数值+和-进行比较，然后基于误差校正方法对水平集函数进行重建。上述水平集函数值+和-分别位于阳性和阴性区域，然后通过逃脱粒子计算定义的局部水平集函数值。x>0和x≤0区域的每个角落，可通过下式计算：
水平集函数值+和-可通过设置+和-的等价变量来合成为一个单一的水平集函数，该水平集都具有在每个网格点的最小值。变量可定义如下：
3 粒子补种方法及改进
3.1 粒子补种方法
在存在拉伸和撕裂的流动界面中，特别是不均匀区域，这样的区域在执行上述混合粒子水平集算法时会缺乏足够数量的粒子产生。为了保持界面的准确性，在极度变形的交接界面附近进行粒子补种以保持适当的粒子分布是很重要的步骤。对此提出一种粒子补种算法，阴阳粒子补种范围，如图2所示：
如果粒子在距离交接界面3个网格点以内的总数量大于最大值1.5Np，或者小于最小值0.5Np，式中Np取值可定义为：
Np=16， for 2D screen64， for 3D screen（15）
粒子被随机删除（或添加）直到颗粒的数目在规定的最大值和最小值之间，如图2所示。补种过程不仅加强了区域特征信息实现对运动区域的完全描述，而且消除了不良粒子诱导。
在处理相关性问题中，在交接界面附近补种适量的粒子数量是很重要的。当基于水平集函数的几何信息进行补种时，表面的数量（面积或体积）可由下式给出：
式中，Ω为操作区域；H为Heaviside辅助函数，可定义如下：
这种基于计算数或几何信息补种程序方法最为简单，但这种补种标准并不明确。如果特征被冲击波或可压缩流动等外界干扰破坏，则原有的粒子补种机制无法删除含有错误信息的污染颗粒，从而影响到水平集切面的准确性。频繁的粒子补种操作并不是解决问题的根本方案，补种操作并不能找出导致界面错误的非正常粒子，反而会删除含有正确特征的粒子。
3.2 改进粒子补种算法
原始粒子补种算法对于处理界面奇异问题和复杂的几何相关性问题时效果不佳，比如单独或联合区域划分问题，因此需要一种更为严格和可靠的粒子补种方法。endprint
速度场和几何流会引起界面的奇异性，比如角落或边缘。通过分析速度场和界面之间的相关性可知包含错误信息的特定区域。当速度场（或特征线）收敛于界面上的一点，由于拉格朗日粒子特性保持的完全性，那么逃逸粒子将会沿该奇异点堆积。这些逃逸的受污染粒子将会通过补种算法清除污染粒子前的校正阶段影响交接界面划分的准确性。由于粒子方法无法区分拓扑事件，例如，合并或分割，因此拓扑结构的变化是影响水平集交接界面准确性的另一个因素。
图3 拓扑变化示意图
本文求解思路是在处理奇异性和拓扑合并或破裂问题时（如图3所示），采用求解双曲型方程的柔性解决方案。采取的策略是：当检测到区域中含有奇异点或拓扑结构变化时，粒子的校正阶段停止执行，水平集函数值保持原始柔性计算值。奇异点是由速度场收敛而出现，并根据相应的交接界面和检测的速度向量和单位法向量的角度值进行评测，如下：
式中，θ表示交接界面附近节点特征聚集程度，如果所有的θ都接近于0，则说明此处存在奇异点，在奇异点处的特征被融合并且传统的粒子算法不能正确计算。在这里假定θ∈-10o，10o时，该处为奇异点。
在图3中，拓扑合并和分裂时发生在两界面彼此非常接近时，这种拓扑结构的变化可通过检查周围环境和水平集值进行检测。如果像素点的水平集值在0，0.1Δx范围内，且水平集值为负值，则该位置附近将发生合并或者分裂的拓扑结构变化。例如，在图3中，节点Ei，j为发生拓扑结构变化的可疑点。
综上所述，增强粒子补种算法的计算步骤如下：
Step1：（初始化）读取水平集待分割图像，并初始化算法参数；
Step2：（奇异点检测）根据3.3节检测当前区域是否属于奇异点或拓扑结构变化点，若是转Step4，否则转Step3；
Step3：（粒子补种）根据3.1节，执行粒子补种算法；
Step4：（粒子状态检测）检查环绕点Ai-1，j-1、Bi+1，j-1、Ci+1，j和Di-1，j的粒子符号状态；
Step5：（粒子删除）判断是否满足删除条件，环绕点符号是否为“-”，若是则删除节点Ei，j附近的粒子。否则转Step6；
Step6：（水平集分割）根据2.1～2.2节执行混合粒子水平集分割算法。判断是否满足终止条件，是则终止算法，并输出结果，否则转Step2，继续执行粒子的补种和删除操作。
4 仿真实验与分析
仿真设备参数：处理器为i7 2.4GHz，内存为 4G ddr1333。仿真平台选取：matlab2012a。实验对象：选取常用医学用CT图像分别为血管CT、肿瘤CT、乳腺CT、脑CT以及胸腔CT五幅图像。评价指标选取图像分割算法的识别率、运算时间及方差。
仿真参数设置：角度θ∈-10o，10o、公式（7）参数a=0.1，b=0.5，粒子补种参数Np=16。所选取的仿真对比算法为标准LS算法以及本文所提出的LPRLS算法。仿真结果如图4（a）～（e）所示，通过对matlab水平集工具箱函数进行研究，并且结合本文所提算法理论，经过编写接口程序实现LPRLS算法。除上述参数设置外，其他水平集参数设置同工具箱内水平集算法参数设置。水平集算法的分割精度、时间消耗以及方差等评测指标仿真结果如表1所示，表1中的仿真结果为上述两种算法分别独立运行20次所求取的评价指标均值。
图4（a）～（e）显示的分别为标准LS算法和LPRLS算法在血管CT、肿瘤CT、乳腺CT、脑CT以及胸腔CT五组图像上的分割效果。（a）图血管CT图像分割中，LPRLS算法比LS算法在主干血管的识别率上要明显更高，但是存在的问题是，由于灵敏度过高，导致存在一定的识别噪声点，如图中散落的小圆圈；（b）图肿瘤CT图像分割中，LPRLS算法比LS算法更全面的对肿瘤轮廓进行分割，并且肿瘤内部较精细的轮廓也能识别出来；（c）图乳腺CT图像分割中，LPRLS算法对整个乳腺的识别全面性要明显高于LS算法，并且乳腺内部也可较为精确的进行分割识别；（d）图脑CT图像分割中，脑CT图像的纹理结构更加复杂，LS算法的分割效果不是很理想，从分割的精细程度上要明显不如LPRLS算法；（e）图胸腔CT图像分割中，LS算法在左胸中部轮廓识别中存在杂乱现象，整体轮廓也较杂乱，而LPRLS算法对于胸腔轮廓识别要更圆润，识别效果更好。图4从视觉直观性上给出了标准LS算法和LPRLS算法的图像分割效果对比，直观上LPRLS算法要明显优于标准LS算法。表1分别给出标准LS算法和LPRLS算法在血管CT、肿瘤CT、乳腺CT、脑CT以及胸腔CT五组图像上的分割数值对比结果。仿真对比指标
（a）血管CT图像
（b）肿瘤CT图像
（c）乳腺CT图像
（d）脑CT图像
（e）胸腔CT图像
选取识别率、运算时间和识别率方差三种，算法各允许20次求均值。在算法分割识别率指标上，LPRLS算法要明显高于标准LS算法，并且这种优势非常明显，在该指标上LPRLS算法要平均高出20%～30%，如在血管CT、肿瘤CT、乳腺CT三组图像中，标准LS算法的识别率近50%～60%，而LPRLS算法在这三幅图像上的分割识别率均达到85%以上；在脑CT和胸腔CT两组图像分割中，标准LS算法略高达到70%以上，但是LPRLS算法在这两组分割图像的分割识别率均达到90%以上。在计算时间指标对比中，标准LS算法在这五组图像上的运算时间在20s～35s之间，而LPRLS算法在这五组图像上的分割时间则分布在16s～20s之间，可看出在计算时间方面LPRLS算法也要优于标准LS算法。在代表算法稳定性方面的算法方差指标中，两者相差不大，LPRLS算法略优于标准LS算法。
5 结束语endprint
利用Lagrangian粒子增强补种算法对水平集算法进行改进，提高算法应对医学图像CT图像像素不均匀所带来的分割精度降低问题。通过在计算水平集公式前先计算Lagrangian标记粒子来重建内嵌交接界面，从而提高水平集算法的质量守恒特性；通过增加速度矢量和单位法向量来促进奇异点和拓扑变化点速度场的收敛，来处理界面奇异性和复杂几何相关问题。算法仿真结果显示，所提算法比原始算法性能大为提升。从仿真结果可看出，存在的不足是，算法存在误识别率，下一步将主要针对此问题进行进一步分析和研究。
参考文献
[1] OSHER S，SETHIAN J A. Fronts propagating with curvaturedependent speed： algorithms based on HamiltonJacobi formulations[J]. Journal of Computational Physics，1988，79（1）：12-49.
[2] 王斌， 李洁， 高新波. 一种基于边缘与区域信息的先验水平集图像分割方法[J]. 计算机学报， 2012， 35（5）： 1067-1071.
[3] 薛维琴， 周志勇， 张涛. 灰度不均的弱边缘血管影像的水平集分割方法[J]. 软件学报， 2012， 23（9）： 2489-2499.
[4] GIBOU F，MIN C，FEDKIW R. High resolution sharp computational methods for elliptic and parabolic problems in complex geometries[J]. Journal of Science Computing， 2013， 54（2）： 369-413.
[5] BOUMEHED M，ALSHAQAQI B，OUAMRI A. Moving Objects Localization by Local Regions Based Level Set： Application on Urban Traffic[J]. Journal of Mathematical Imaging and Vision， 2013， 46（2）： 258-274.
[6] JEFFREY C L，ZACHARY M. Level sets of the Takagi function： local level sets[J]. Monatshefte für Mathematik， 2012， 166（2）： 201-238.
[7] MARCHANDISE E，REMACLE J F，CHEVAUGEON N. A quadraturefree discontinuous Galerkin method for the level set equation[J]. Journal of Computational Physics， 2006， 212（1）： 338–357.
[8] MU L，WANG J P，WANG Y P. A computational study of the weak Galerkin method for secondorder elliptic equations[J]. Numerical Algorithms， 2013， 63（4）： 753-777.
[9] NAVE J C，ROSALES R R，SEIBOLD B. A gradientaugmented level set method with an optimally local coherent advection scheme[J]. Journal of Computational Physics， 2010， 229（10）： 3802–3827.
[10]LEE C，DOLBOW J，MUCHA P J. A narrowband gradientaugmented level set method for multiphase incompressible flow[J]. Journal of Computational Physics， 2014， 273（5）： 12–37.
[11]杨谊， 喻德旷. 基于改进水平集图像分割方法的乳腺超声病灶提取[J]. 计算机应用与软件， 2014， 31（11）： 217-221.
[12]Li Changyang， Wang Xiuying， Eberl S. A Likelihood and Local Constraint Level Set Model for Liver Tumor Segmentation from CT Volumes[J]. IEEE Transactions on Biomedical Engineering， 2013， 60（10）： 2967-2977.
[13]KACHROO P，RATLIFF L，SASTRY S. Analysis of the GodunovBased Hybrid Model for Ramp Metering and Robust Feedback Control Design[J]. IEEE Transactions on Intelligent Transportation Systems， 2014， 15（5）： 2132-2142.
[14]Zhuang Chijie， Zeng Rong， Zhang Bo. A WENO Scheme for Simulating Streamer Discharge With Photoionizations[J]. IEEE Transactions on Magnetics， 2014， 50（2）： #7007904.
[15]刘世兴， 宋端， 贾林. 辛RungeKutta 方法在求解LagrangeMaxwell方程中的应用研究[J]. 物理学报， 2013， 62（3）： #034501.
[16]OSHERA S，Cheng Litien， KANG M J. Geometric Optics in a PhaseSpaceBased Level Set and Eulerian Framework[J]. Journal of Computational Physics， 2002， 179（2）： 622-648.
篇9：补种苗木变更申请
1 材料与方法
1.1 材料
供试品种为黑龙江省农业科学院黑河分院选育的不同熟期大豆品种:黑河35 (生育期90d) 、黑河25 (生育期98d) 、黑河45 (生育期107d) 、黑河43 (生育期115d) 、黑河38 (生育期116d) 和黑河48 (生育期116d) 。
1.2 方法
试验于2012年在黑龙江省农业科学院黑河分院试验地 (N50°15′、E127°27′) 进行, 该地属于黑龙江省第四积温带。试验分为4个播期:第Ⅰ播期为5月11日 (为当地正常播期) 、第Ⅱ播期为5月21日、第Ⅲ播期为5月31日、第Ⅳ播期为6月10日。行长5 m, 4行区, 3次重复, 垄距60cm, 株距8cm。人工开沟双粒点播, 播种深度3~4cm。生育期间化学除草, 田间正常管理。
2 结果与分析
2.1 不同播期对大豆农艺性状的影响
由表1可知, 各品种株高随着播期的延后逐渐增高, 播期对不同大豆品种株高的影响存在差异, 早熟品种株高变化幅度大于晚熟品种。不同播期, 不同大豆品种的分枝数差异较大。各品种在第Ⅱ播期和第Ⅲ播期的分枝数较多。黑河35分枝数最多, 黑河45分枝数最少, 这是品种基因型决定的。各品种在不同播期条件下的主茎节数差异不大, 但前3个播期各品种主茎节数都高于第Ⅳ播期。黑河48在不同播期主茎节数差异最小。
2.2 不同播期对大豆产量及产量构成因子的影响
由表2可知, 随着播期的推迟, 各品种单株粒数、百粒重 (除了黑河25和黑河48外) 和产量均呈先升高后下降的趋势。早熟品种黑河35在第Ⅲ播期下, 获得较高百粒重和最高单株粒数及产量;黑河45在第Ⅱ播期下获得最高百粒重和产量;晚熟品种黑河38和黑河48在第Ⅱ播期下获得最高单株粒数和产量, 在第Ⅳ播期下, 单株粒数、百粒重和产量均最低。在当地正常播期条件下, 早熟品种产量比晚熟品种低, 随着播期的推迟, 早熟品种的产量高于晚熟品种, 综合产量构成因子的变化可以看出, 各品种适当晚播可提高产量。
3 结论
许多研究表明[2,3,4,5], 播期对大豆株高等农艺性状影响显著。随着播期的推迟, 株高逐渐增高, 说明适宜的晚播有利于大豆植株的增高, 但株高生长过快, 易形成弱苗, 引起倒伏。随着播期的推迟, 分枝数有先增加后减少的趋势, 说明过于早播和晚播均不利于分枝数的增加。
随着播期的推迟, 各品种单株粒数、百粒重 (除黑河25和黑河48外) 和产量均先升高后降低, 黑河35和黑河25在第Ⅲ播期产量最高, 黑河45、黑河43、黑河38和黑河48在第Ⅱ播期产量均最高, 所以由于某些原因造成当地主栽品种不能正常播种时, 应选择熟期较短的大豆品种以得到较高的产量。
参考文献
[1]高福成.爱辉区农业气候灾害分析及预防[J].黑河科技, 2000 (4) :12-13.
[2]宁海龙, 孙培乐, 宋兆华, 等.不同播期对春大豆生态性状的影响[J].大豆科学, 2011, 30 (1) :73-78.
[3]董全中.迟播对早熟大豆产量、品质及农艺性状的影响[J].黑龙江农业科学, 2007 (3) :13-16.
[4]李树臣, 陈学珍, 谢皓, 等.不同播期对夏播大豆产量及农艺性状的影响[J].北京农学院学报, 2003, 18 (4) :266-269.