CN103514457B - 生态约束条件下土地利用结构优化的多目标改进差分进化模型 - Google Patents
生态约束条件下土地利用结构优化的多目标改进差分进化模型 Download PDFInfo
- Publication number
- CN103514457B CN103514457B CN201310322556.8A CN201310322556A CN103514457B CN 103514457 B CN103514457 B CN 103514457B CN 201310322556 A CN201310322556 A CN 201310322556A CN 103514457 B CN103514457 B CN 103514457B
- Authority
- CN
- China
- Prior art keywords
- land
- remote sensing
- formula
- vector
- land resource
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000006872 improvement Effects 0.000 title abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 82
- 238000005457 optimization Methods 0.000 claims abstract description 64
- 239000013598 vector Substances 0.000 claims abstract description 45
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 19
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 238000011160 research Methods 0.000 claims description 46
- 230000008901 benefit Effects 0.000 claims description 22
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 8
- 238000005520 cutting process Methods 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 6
- 230000006978 adaptation Effects 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 238000012805 post-processing Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 238000011161 development Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000008239 natural water Substances 0.000 description 2
- 230000000452 restraining effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000577395 Thenus Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000008713 feedback mechanism Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
一种生态约束条件下土地利用结构优化的多目标改进差分进化模型,其包括:1)遥感影像预处理;2)遥感影像分类,得到土地利用现状数据;3)利用差分进化算法,在输入种群规模NP、个体的维数D、差分矢量缩放因子F和交叉概率CR后,经过初始化种群、变异操作、交叉操作、选择操作,对土地利用现状数据进行土地资源优化处理,得到土地资源优化配置结果。
Description
技术领域
本发明涉及卫星遥感,尤其涉及基于卫星遥感的土地利用结构优化方法。
背景技术
以往土地资源优化配置大多采用线性规划法、目标规划法和***动力学等数学模型及在此基础上形成的集成优化方法,尚未解决多目标和多约束条件下的复杂优化问题。
线性规划法由于其自身的线性特征,在解决具有多目标多约束非线性的土地资源优化配置中显得力不从心,往往得不到理想的最优解。
目标规划模型虽然能够处理多目标土地利用优化配置问题,但它容易受用户已有知识水平和思维方式等的影响,在处理具有多个不同优先等级目标的多目标优化问题上难以把握。
***动力学模型虽然可以根据经济、社会和政策等多种因素的影响,在土地资源优化配置过程中不断调整优化目标和优化方案,但此类模型的建立需要对优化***有充分的认识,对***内各种反馈机制要非常了解,在因果关系不明确的情况下,不宜使用。
由著名学者StornR和PriceK于1995年提出的差分进化算法是一种随机的并行直接搜索算法,适于求解一些利用传统和常规的数学优化方法所无法解决的复杂优化问题,以其易用性、鲁棒性和强大的全局优化能力在多个领域成功应用(StornR,PriceK.DifferentialEvolution–asimpleandefficientadaptiveschemeforglobaloptimizationovercontinuousspaces[J].TechnicalReport,InternationalComputerScienceInstitute,Berkley.1995(8):22-25.)。土地利用***是一个涉及多目标多约束的复杂***,利用常规的数学优化方法来解决土地资源优化配置这一复杂优化问题的效果并不令人满意。由于差分进化算法具有算法简单、收敛速度快、所需领域知识少、比较适合于解决复杂的优化问题等特点,因此,利用差分进化算法进行土地资源优化配置具有很好的潜力。最近几年,国外一些学者已开始尝试利用差分进化算法进行农作物种植规划,并取得了一定的成果,但是针对区域土地资源优化配置的研究还非常少见;在国内,差分进化算法主要应用在人工智能、化工、生物、数据挖掘、影像处理等领域,利用差分进化算法进行土地资源优化配置的相关研究还未见报道。
发明内容
本发明旨在建立一种基于遥感影像的土地资源优化配置的方法,以解决区域生态约束条件下的土地资源优化配置问题。
本发明的基于遥感影像的土地资源优化配置方法,其包括:
1)遥感影像预处理;
2)遥感影像分类,得到土地利用现状数据;
3)利用差分进化算法,在输入种群规模NP、个体的维数D、差分矢量缩放因子F和交叉概率CR后,经过初始化种群、变异操作、交叉操作、选择操作,对土地利用现状数据进行土地资源优化处理,得到土地资源优化配置结果。
优选地,所述选择操作中,通过目标函数F(X)=f(X)+δ(t)*h(X)(式6)来作出选择;
其中,(式4);k为目标函数的个数;wm为第m个目标函数fm(X)的权重系数;
ps(X)=max{0,gs(X)}(式8);
式7为惩罚因子的计算;式8为对约束条件的违背程度;式9为惩罚函数的强度;式10为多段映射函数;式11为随迭代次数变化而不断调整的惩罚力度,其中,X表示群体,其由同一代的所有个体构成;s表示进化到第s代;t表示进化的代数;gs(X)表示不等式约束。
优选地,其中k=2;f1(X)为生态效益函数;f2(X)为经济效益函数。
优选地,所述差分矢量缩放因子F为自适应缩放因子Fi。
优选地,所述自适应缩放因子Fi的适应策略为:
式中,fa、fb和fc分别是个体矢量和的适应度,其中,将三个随机选择的用于变异的个体矢量按照适应度大小排列顺序,适应度最优的个体矢量记为适应度次优的个体矢量记为适应度最差的个体矢量记为Fl和Fu分别为缩放因子F的上限和下限;Fi在[Fl,Fu]之间随着和在矢量空间中的相对位置来自适应地变化;个体矢量和的原始值来自于土地利用现状数据。
优选地,所述交叉概率CR为自适应交叉概率CRi。
优选地,所述自适应交叉概率CRi的适应策略为:
CRu和CRl分别为交叉概率的上限和下限;g为进化过程中的迭代代数,G为进化的最大代数。
优选地,所述变异操作中的变异策略为:
为生成的变异矢量;r1,r2∈{1,2,…,NP},为互不相同且与当前目标矢量序号i也不相同的整数了;为当代种群中适应度最优的个体。
优选地,所述影像预处理包括:
几何精确校正与配准;
遥感影像镶嵌与裁剪;
最佳波段组合。
优选地,所述影像分类包括利用ISODATA非监督分类方法进行分类,得到初步分类结果;然后对初步分类结果进行分类后处理,根据野外调查数据、GPS定位数据、土地利用现状图和其他相关资料,采用目视解译法修改混分与错分像元,得到研究区土地利用分类结果,最后对分类结果进行精度检验。
本发明针对当前土地资源优化配置中广泛采用的目标规划法、***动力学和线性规划法等数学模型和方法存在的问题与不足,提出了一种基于遥感影像的利用差分进化算法进行土地资源优化配置的***及其方法,来解决区域生态约束条件下的土地资源优化配置问题。本发明提出的方法更为简洁,运行效率更高,也具有更好的鲁棒性和实用性。
附图说明
图1是研究区2000、2005、2010年TM影像(R5、G4、B3);
图2是遥感影像分类流程图;
图3是研究区2000、2005、2010年遥感影像分类图;
图4是研究区土地利用空间格局优化结果图;
图5是本发明的基于遥感影像的土地资源优化配置方法的流程图。
具体实施方式
下面结合附图对本发明进行详细说明。
如图2、5所示,基于遥感影像的利用差分进化算法进行土地资源优化配置***及其方法,包含以下模块或步骤:
(1)遥感影像预处理
由于遥感影像在数据采集过程中会受到空间、时间、波谱、辐射分辨率等多方面限制,不可避免的存在误差;这些误差将会降低遥感影像的数据质量和影响影像分析的精度;因此,在遥感影像分析之前,需要对原始的遥感影像进行预处理。遥感影像预处理包括:几何精校正、配准、影像镶嵌与裁剪、波段组合等几个环节。
①几何精校正与配准
几何精校正是在***校正的基础上,利用若干个地面控制点来建立一个模拟影像几何崎变的数学模型,在此基础上建立待校正影像空间与参考影像空间的校正变换函数。再根据这种变换函数把待校正影像空间中的像元变换到标准空间中去,进而实现遥感影像的几何精校正。
②遥感影像镶嵌与裁剪
镶嵌是指当研究区大于单幅遥感影像所覆盖的范围时,一般需要将两幅或多幅影像拼接在一起形成一幅或多幅覆盖全区的较大影像。首先,选择其中面积最大的一幅影像作为参考影像,用来决定镶嵌影像的对比度匹配、输出影像的像元大小、数据类型等;在遥感软件环境下,最终得到覆盖研究区全域的镶嵌影像。在镶嵌影像的基础上,以研究区行政界线文件在遥感软件中实现影像的矢量栅格的转化,运用掩模运算实现影像的不规则裁剪,从而得到研究区的遥感影像。
③最佳波段组合
即指对遥感影像的各波段进行合理的组合与分析,以期得到信息量最丰富、能够最有效地区分各土地利用类型的彩色影像。最佳波段组合能够降低波段维数,去除冗余信息波段,从而提高影像处理的效率和速度,是有效提取土地利用类别的重要措施。
遥感影像波段组合方案一般为3个波段合成,波段的选择对于最大程度地表现地物信息非常重要。最佳波段选择的原则有两点:所选择的波段组合信息量最大和对需要识别的地物类别之间最容易区分。一般来说,标准差越大,所包含的信息量就越大;波段相关系数越小,各波段的独立性越高,信息冗余度越小。
(2)遥感影像分类
遥感影像分类,采取以计算机分类为主,目视解译为辅的方法。经过对得到的研究区遥感影像进行分类实验分析,发现利用非监督分类方法比利用监督分类方法效果要好;因此,本研究利用ISODATA非监督分类方法进行分类,得到初步分类结果。然后,对初步分类结果进行分类后处理,根据野外调查数据、GPS定位数据、土地利用现状图和其他相关资料,采用目视解译法修改混分与错分像元,得到研究区土地利用分类结果,最后对分类结果进行精度检验。在分类结果的基础上,通过计算与统计手段,得到研究区土地利用现状数据,在此基础上,进行土地资源优化配置。
(3)控制参数改进
差分进化算法主要涉及种群规模NP、个体(变量)的维数D、差分矢量缩放因子F和交叉概率CR共四个控制参数,与NP和D比较起来,F和CR对算法的优化性能有着更重要的影响。其中,变量的维数D一般根据实际问题可以确定;种群规模NP一般介于5D到10D之间;缩放因子F取值范围为[0,2],交叉概率CR取值范围为[0,1]。
一般来说,合适参数的确定与实际问题有关,不同的优化问题对应不同的控制参数。本发明在总结现有研究的基础上提出缩放因子F和交叉概率CR两个控制参数的改进策略,以期使得到的优化结果更加科学合理。
①缩放因子F的自适应
差分矢量缩放因子F是差分进化算法的一个非常重要的控制参数,F控制扰动量的缩放幅度,对优化性能有着重要的影响。以前差分进化算法中缩放因子F大多数人为选取固定常数,即从优化开始到结束缩放因子F固定不变,这种做法人为因素影响较大,较难确定合适的参数,因此,需要制定缩放因子F的自适应策略,这种策略能够提供相对于每个解的最佳F值。
该策略的基本思路是根据变异操作中随机选取的三个个体在搜索空间的相对位置来自适应地调整缩放因子F的大小。如果随机选择的三个个体的位置在搜索空间相邻较近,则生成的差分矢量就会比较小,同时各维决策变量的扰动量也会较小,此时整缩放因子F就应该选择较大的值,否则扰动量太小,起不到变异的作用,这样就会导致在进化初期过早陷入局部搜索。如果随机选取的三个个体在搜索空间中的位置比较分散,则生成的差分矢量值就会太大,从而对各维决策变量的扰动幅度也相对较大,此时缩放因子F应选取相对较小的值,来控制扰动量太大,从而有利于算法进行局部搜索。因此,缩放因子F的取值应该是自适应的,能依据生成差分矢量的两个个体在搜索空间中的相对位置来调整自身变化和动态调整差分矢量的大小,从而达到平衡局部搜索和全局搜索之间的矛盾以及避免人为确定F的值带来的不确定性影响。具体的自适应策略如下:
式中,fa、fb和fc分别是个体矢量和的适应度(目标函数值),其中,将三个随机选择的用于变异的个体按照适应度大小排列顺序,适应度最优的个体矢量记为适应度次优的个体矢量记为适应度最差的个体矢量记为Fl和Fu分别为缩放因子F的上限和下限,一般Fl可取0.1,Fu可取0.9。Fi为自适应缩放因子,在[0.1,0.9]之间随着和在矢量空间中的相对位置来自适应地变化。
②交叉概率CR的自适应
交叉概率CR是通过控制种群中个体矢量的选择来影响多样性,从而影响算法的搜索速度与搜索成功率。在搜索的开始阶段,CR取值应稍小些,这有利于防止早熟和保持种群的多样性;但在搜索后期,交叉概率CR取值应大些,可以提高算法的局部搜索能力。本发明提出了一种交叉概率自适应策略,具体自适应策略如下:
式中,CRi为自适应交叉概率。CRu和CRl分别为交叉概率的上限和下限,一般CRl取值0.3,CRu取值0.9,当然,CRu和CRl可根据实际问题来确定。g为进化过程中的迭代代数,G为进化的最大代数。
(4)变异策略改进
差分进化算法中标准变异策略是DE/rand/1/bin,该变异操作中的和是从种群中随机选取的3个互不相同的个体矢量,这种策略有利于差分进化算法降低收敛速度和进行全局搜索。本发明利用改进的变异策略DE/best/1/bin,具体公式如下:
式中,为生成的变异矢量;r1,r2∈{1,2,…,NP},为互不相同且与当前目标矢量序号i也不相同的整数。为当代种群中适应度最优的个体。这种改进变异策略,以为变异的基向量,不再是完全的随机搜索,有利于加快算法的收敛速度。
(5)多目标函数处理方法
本发明中差分进化算法的多目标函数处理采用加权和法,具体公式为:
min(或max): 式4
式中,F(X)为新的效用函数(优化函数);k为目标函数的个数;wm为第m个目标函数的权重值,一般满足由于每个目标的权重值不同,通过调整权重系数的大小就能得到不同情景的多目标优化结果。
(6)约束条件的处理
大多数优化问题都包含不等式约束和等式约束两种类型,而且这些约束条件很多情况下都是非线性的,这就导致约束优化问题的求解比无约束优化问题的求解复杂和困难得多。一般来说,罚函数法是求解约束优化问题最常用的方法,它的基本思路是利用一个罚因子来对约束函数进行惩罚,再把得到的惩罚项加到优化目标函数中,从而构成一个无约束的广义目标函数,然后采用优化算法对新生成的广义目标函数进行求解,最终得到在惩罚项作用下的最优解。
约束优化问题一般表述为:
min(或max):f(X)
式中,f(X)表示目标函数,gq(X)表示不等式约束,hr(X)表示等式约束;通常在处理过程中要首先把等式约束转变为不等式约束。具体方法为:首先引入松弛因子ε,对于等式约束:h(X)=0,然后利用松弛因子把该等式约束转变为不等式约束:|h(X)|<ε;一般来说,ε是一个很小的正数,习惯上取其为l0-3、l0-4等实数。在此基础上,构造广义目标函数:
F(X)=f(X)+δ(t)*h(X)式6
式中,F(X)表示广义目标函数,f(X)表示原目标函数;δ(t)*h(X)表示惩罚项,δ(t)表示惩罚力度,h(X)表示惩罚因子。如果惩罚力度δ(t)为一个常数,在整个优化过程中固定不变,则称为固定罚函数法,如果δ(t)在优化过程中是变化的,则称为非固定罚函数法。非固定罚函数法比固定罚函数法有更多的优点,因为它在计算过程中,惩罚力度δ(t)可以随迭代次数t不断变化,给整个优化过程带来一定的灵活性,使求解过程变得更加快捷。本发明提出改进型非固定多段映射罚函数法进行约束条件的处理,数学表达为:
ps(X)=max{0,gs(X)}式8
式中,公式7式为惩罚因子的计算;公式8式为对约束条件的违背程度;公式9式为惩罚函数的强度;公式10式为多段映射函数;公式11式为随迭代次数变化而不断调整的惩罚力度;其中,X表示群体,也就是某一代中的所有个体;s表示进化到第s代;ps(X)表示对约束条件的违背程度;t表示进化的代数;gs(X)表示不等式约束,gs(X)为根据实际问题来确定的约束条件。
(7)土地资源优化配置***的编码框架
下面是基于遥感影像的利用差分进化算法进行土地资源优化配置***及其方法实现的编码框架。
*初始化参数值:D=6,NP=60,Fl=0.1,Fu=0.9,CRl=0.3,CRu=0.9,maxgen=500,xmin=0,xmax=138700
(初始化种群)
Fori=1toNP
Forj=1toD
x(i,j)=xmin+randomnumber[0,1]*[xmax-xmin]
Nextj
Nexti
*给迭代变量赋初值gen=1,然后在本次迭代中进行各种进化操作。
Forgen=1tomaxgen
Fori=1toNP
*首先将最大化函数转换为最小化函数,然后计算目标函数值,并记录
当代种群中目标函数值最优的个体
*检查各个变量是否违反约束条件,并采用改进型非固定多段映射罚
函数法进行约束条件处理;如果违反,则惩罚函数值。
f(xgen,i)=f(xgen,i)+δ(t)*h(X)。
Forj=1toD
(变异操作)
*从当前种群中任意选择三个互不相同而且与当前个体也不同的个体,得到每个个体相应的决策变量:x(gen,r1,j),x(gen,r2,j)andx(gen,r3,j),即要满足i≠r1≠r2≠r3,然后进行变异操作。
r1=randomnumber[0,1]*NP
r2=randomnumber[0,1]*NP
r3=randomnumber[0,1]*NP
*求出r1,r2,r3的适应度(目标函数值),然后根据其大小排列顺序。其中,xa表示适应度最优的个体矢量,xb表示适应度次优的个体矢量,xc表示适应度最差的个体矢量。然后,自适应计算缩放因子F。
*采用本发明提出的改进变异策略计算变异矢量。
*越界处理,即检查每个决策变量是否超出其取值范围,如果越界做如下处理。
Ifv(gen,i,j)<xminthenv(gen,i,j)=xmin+randomnumber[0,1]*(x(gen,i,j)-xmin)
Ifv(gen,i,j)>xmaxthenv(gen,i,j)=x(gen,i,j)+randomnumber[0,1]*(xmax-x(gen,i,j))
(交叉操作)
*采用本发明提出的自适应策略计算交叉概率。
*进行交叉操作。
p=randomnumber[0,1]
n=randomnumber[0,1]*Dn∈{1,2,…,D}(n为随机数)
Ifp<CRorj=n,Thenu(gen,i,j)=v(gen,i,j)
Elseu(gen,i,j)=x(gen,i,j)
Nextj
(选择操作)
*将经过变异和交叉操作后得到的个体与原个体进行比较,如果新个
体优于原个
体,就用新个体代替原个体并保留到下一代,否则,将原个体保留到下一代。
Iff(u(gen,i))<f(x(gen,i))thenx(gen+1,i)=u(gen,i)
Elsex(gen+1,i)=x(gen,i)
Nexti
Nextgen
*输出结果(Xbest,f(Xbest)),即输出适应度最优的目标函数值以及对应的决策变量。
End
实施例:
本发明实施案例选择具有典型滨海生态特色的辽宁省大洼县为研究区域。
(1)遥感影像预处理
本研究采用的原始数据主要有:2000年、2005年和2010年三期LandsatTM和LandsatETM+遥感影像,经过几何精校正、配准、影像镶嵌与裁剪、波段组合等几个环节的遥感影像预处理步骤,得到研究区2000年、2005年和2010年三期TM影像,见图1。
(2)遥感影像分类
①分类过程
为提高分类的精度,本文对遥感影像分类,采取以计算机分类为主,目视解译为辅的方法。经过对得到的研究区三期TM影像进行分类实验分析,发现利用非监督分类方法比利用监督分类方法效果要好;因此,本研究利用ENVI软件的ISODATA非监督分类方法进行分类。在分类的过程中,首先用ISODATA非监督分类方法把影像自动分为15类,然后再合并成需要的六种地类(耕地、林地、建设用地、湿地、水域、滩涂),得到研究区初步分类结果。然后,对初步分类结果进行分类后处理,根据野外调查数据、GPS定位数据、土地利用现状图和其他相关资料,采用目视解译法修改混分与错分像元,得到研究区三期土地利用分类结果,最后对分类结果进行精度检验。具体分类过程如图2所示。
②精度检验
根据与遥感影像同期的地面资料,通过随机选取若干样本区,利用ENVI软件计算其分类混淆矩阵和Kappa系数,从而对研究区三期TM遥感影像的分类结果进行精度检验。结果表明:2000年、2005年和2010年三期影像分类精度分别为82.43%、87.38%和84.45%;Kappa系数分别为0.81、0.85和0.83,即均达到最低允许判别精度0.70的要求。研究区2000年、2005年和2010年三期TM影像分类结果见图3。
(3)变量设置
根据研究区遥感影像分类结果,土地利用类型包括:耕地、林地、湿地、水域、滩涂和建设用地。因此,基于遥感影像的多目标改进差分进化算法的变量设置为:x1(耕地)、x2(林地)、x3(湿地)、x4(水域)、x5(滩涂)、x6(建设用地)。优化结果的求解就是找到6个变量当处于某一数量结构时,目标函数值达到最优。
(4)多目标函数的建立
目标函数主要有两个,即生态效益的最大化和经济效益的最大化。最后采用加权和法将多目标优化问题转变为单目标优化问题。
①目标函数一:生态效益最大化
土地资源优化配置的一个目标就是达到生态效益的最大化,生态效益最大化的计算过程就是每种土地利用类型面积乘以相应的单位面积生态服务价值,然后相加,最后得到区域总的生态效益价值。具体计算公式如下:
max: 式12
式中,xu为第u种土地利用类型,pu为第u种土地利用类型的单位面积生态服务价值。根据实地调查和文献研究,p1=0.692万元/hm2,p2=1.194万元/hm2,p3=1.082万元/hm2,p4=0.764万元/hm2,p5=0.761万元/hm2,p6=0。
②目标函数二:经济效益最大化
土地资源优化配置的另一个目标是达到经济效益的最大化,经济效益最大化的计算过程就是每种土地利用类型面积乘以相应的单位面积经济效益,然后相加,最后得到区域总的经济效益价值。具体计算公式如下:
max: 式13
式中,xv为第v种土地利用类型,qv为第v种土地利用类型的单位面积经济效益。根据实地调查和相关文献研究结果,综合确定每种土地类型的单位面积经济效益为:q1=3.0万元/hm2,q2=0.4万元/hm2,q3=50.0万元/hm2,q4=2.7万元/hm2,q5=0.2万元/hm2,q6=10.0万元/hm2。
③多目标函数的建立
优化目标是在体现生态效益优先并兼顾经济效益的条件下,进行土地资源优化配置。因此,在确定目标函数的权重时,应优先考虑生态效益,在前面生态服务价值测算基础上,采用专家咨询法来确定权重系数,最终得到目标函数一(生态效益最大化)的权重w1=0.65,目标函数二(经济效益最大化)的权重w2=0.35。最终多目标优化函数表达式为:
max: 式14
式中,f(X)为最终优化目标函数,在优化问题求解过程中,只需考虑f(X)的目标函数值。
(5)生态约束条件的构建
根据研究区自然条件、土地利用现状、研究区土地利用总体规划(2006-2020年)、研究区生态县建设规划(2008-2020年)和研究区国民经济和社会发展十二五规划,从生态环境保护、粮食安全、经济发展等多个方面来构建土地资源优化配置中的生态约束条件。以2020年为优化目标年为例,主要包括以下八个方面:
①土地总面积约束。各类土地面积之和应等于研究区土地总面积,即:
②耕地面积约束。随着研究区社会经济的快速发展,耕地面积仍会呈现下降趋势。但是,为了保障研究区的粮食生产总量,确保粮食安全,需要在一定时期内保留一定数量的耕地。根据研究区土地利用总体规划,截至2020年在规划期内的耕地保有量应不小于61675.00hm2。因此,耕地面积约束条件为:
61675.00hm2≤x1≤62768.29hm2式16
③林地面积约束。根据《大洼县生态县建设规划(2008-2020年)》,在规划期内,研究区将会大规模开展植树造林活动,增加林地面积。即:
x2≥2279.78hm2式17
④湿地面积约束。根据对研究区2000-2010年期间的土地利用变化分析,随着自然环境的变化以及人为因素的干扰,湿地面积呈现逐渐减少的趋势。但是,根据研究区生态县建设规划和生态用地重要性评价结果,研究区在优化期内将会采取严格的生态保护措施来保护湿地。因此,在2010-2020年的优化期内,湿地面积将有可能维持现状或小于现状面积。即:
x3≤14695.97hm2式18
⑤水域面积约束。研究区水域主要包括自然水面和人工养殖水面,根据研究区土地利用总体规划、生态县建设规划以及生态用地重要性评价结果,在保护自然水面的同时,应当适当控制或减少人工养殖水面,因此,在2010-2020年的优化期内,水域面积将小于现状面积。即:
x4≤14281.84hm2式19
⑥滩涂面积约束。根据研究区土地利用总体规划和生态用地重要性评价结果,研究区沿海滩涂在生态保护前提下,将被适当开发;但是,为了保护研究区生态环境,滩涂不能无限制开发,根据研究区生态县建设规划,到2020年滩涂面积要保持在22318.39hm2以上;因此,滩涂面积约束条件为:
22318.39hm2≤x5≤22905.71hm2式20
⑦建设用地面积约束。根据研究区社会经济发展和城乡建设用地增减挂钩项目的开展,研究区建设用地面积将会继续增加,但根据研究区土地利用总体规划,到2020年研究区建设用地面积不得超过22608hm2。因此,建设用地面积将大于现状面积而小于控制指标。即:
21768.41hm2≤x6≤22608.00hm2式21
⑧非负约束。即各类土地面积均不得为负值:
xi≥0,i=1,2,...,6式22
(6)土地资源优化配置***的建立与实现
根据前面建立的该***的编码框架,在VC++6.0环境下编程实现,其中,变量维数D设为6,种群规模NP设为60,缩放因子F和交叉概率因子CR通过提出的自适应策略确定,最大迭代次数为500,以2010年为基期,2020年为优化目标年,经过多次独立运行试验,最后得到研究区2020年土地利用数量结构优化结果(表1)和空间格局优化结果,见图4。
表1土地利用数量结构优化结果
说明:由于差分进化算法优化过程带有一定的随机性,优化结果将会在很小的范围内浮动,表中结果是取的几次优化结果的平均值。
Claims (9)
1.一种基于遥感影像的土地资源优化配置方法,其包括:
1)遥感影像预处理;
2)遥感影像分类,得到土地利用现状数据;
3)利用差分进化算法,在输入种群规模NP、个体的维数D、差分矢量缩放因子F和交叉概率CR后,经过初始化种群、变异操作、交叉操作、选择操作,对土地利用现状数据进行土地资源优化处理,得到土地资源优化配置结果;
所述选择操作中,通过目标函数F(X)=f(X)+δ(t)*h(X)(式6)来作出选择;
其中,(式4);k为目标函数的个数;wm为第m个目标函数fm(X)的权重系数;
ps(X)=max{0,gs(X)}(式8);
式7为惩罚因子的计算;式8为对约束条件的违背程度;式9为惩罚函数的强度;式10为多段映射函数;式11为随迭代次数变化而不断调整的惩罚力度,其中,X表示群体,其由同一代的所有个体构成;s表示进化到第s代;t表示进化的代数;gs(X)表示不等式约束。
2.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
其中k=2;f1(X)为生态效益函数;f2(X)为经济效益函数。
3.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述差分矢量缩放因子F为自适应缩放因子Fi。
4.如权利要求3所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述自适应缩放因子Fi的适应策略为:
式中,fa、fb和fc分别是个体矢量和的适应度,其中,将三个随机选择的用于变异的个体矢量按照适应度大小排列顺序,适应度最优的个体矢量记为适应度次优的个体矢量记为适应度最差的个体矢量记为Fl和Fu分别为缩放因子F的下限和上限;Fi在[Fl,Fu]之间随着和在矢量空间中的相对位置来自适应地变化;个体矢量和的原始值来自于土地利用现状数据。
5.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述交叉概率CR为自适应交叉概率CRi。
6.如权利要求5所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述自适应交叉概率CRi的适应策略为:
CRu和CRl分别为交叉概率的上限和下限;g为进化过程中的迭代代数,G为进化的最大代数。
7.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述变异操作中的变异策略为:
为生成的变异矢量;r1,r2∈{1,2,…,NP},为互不相同且与当前目标矢量序号i也不相同的整数;为当代种群中适应度最优的个体。
8.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述影像预处理包括:
几何精确校正与配准;
遥感影像镶嵌与裁剪;
最佳波段组合。
9.如权利要求1所述的基于遥感影像的土地资源优化配置方法,其特征在于:
所述影像分类包括利用ISODATA非监督分类方法进行分类,得到初步分类结果;然后对初步分类结果进行分类后处理,根据野外调查数据、GPS定位数据、土地利用现状图,采用目视解译法修改混分与错分像元,得到研究区土地利用分类结果,最后对分类结果进行精度检验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310322556.8A CN103514457B (zh) | 2013-07-29 | 2013-07-29 | 生态约束条件下土地利用结构优化的多目标改进差分进化模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310322556.8A CN103514457B (zh) | 2013-07-29 | 2013-07-29 | 生态约束条件下土地利用结构优化的多目标改进差分进化模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103514457A CN103514457A (zh) | 2014-01-15 |
CN103514457B true CN103514457B (zh) | 2016-06-29 |
Family
ID=49897149
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310322556.8A Expired - Fee Related CN103514457B (zh) | 2013-07-29 | 2013-07-29 | 生态约束条件下土地利用结构优化的多目标改进差分进化模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103514457B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2983574C (en) * | 2014-04-22 | 2021-11-16 | Blast Motion Inc. | Initializing an inertial sensor using soft constraints and penalty functions |
CN106991155A (zh) * | 2017-03-29 | 2017-07-28 | 北京图畅科技有限公司 | 基于3s技术进行地籍调查的方法 |
CN108564592B (zh) * | 2018-03-05 | 2021-05-11 | 华侨大学 | 基于动态多种群集成差分进化算法的图像分割方法 |
CN109508881B (zh) * | 2018-11-12 | 2021-04-16 | 自然资源部第二海洋研究所 | 海岛区域分类及生态资源价值评估方法 |
CN109472320B (zh) * | 2018-11-30 | 2021-09-07 | 南京农业大学 | 不确定条件下作物生育期模型品种参数自动校正方法 |
CN110990156A (zh) * | 2019-12-09 | 2020-04-10 | 深圳集智数字科技有限公司 | 一种资源自动配置方法及装置 |
CN111724271A (zh) * | 2020-04-27 | 2020-09-29 | 齐实 | 一种基于生长营养空间的高度盐碱化地区绿化植生方法 |
CN112001085B (zh) * | 2020-08-25 | 2022-05-17 | 南京林业大学 | 一种小流域水土保持土地利用优化方法及*** |
CN112149872A (zh) * | 2020-08-25 | 2020-12-29 | 南京森林警察学院 | 一种小流域土地利用结构优化方法及*** |
CN112183459B (zh) * | 2020-10-20 | 2022-10-04 | 安徽大学 | 一种基于进化多目标优化的遥感水质图像分类方法 |
CN112729130A (zh) * | 2020-12-29 | 2021-04-30 | 四川天奥空天信息技术有限公司 | 卫星遥感测量树木冠层高度的方法 |
CN114021913A (zh) * | 2021-10-22 | 2022-02-08 | 中冶南方工程技术有限公司 | 一种基于差分进化算法的高炉配料优化方法、电子设备以及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102325089A (zh) * | 2011-09-19 | 2012-01-18 | 西安电子科技大学 | 基于差分进化和捕食搜索策略的胖树型片上网络映射方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8249981B2 (en) * | 2008-12-16 | 2012-08-21 | Ge Corporate Financial Services, Inc. | Methods and systems for generating transition probability matrices through an optimization framework |
-
2013
- 2013-07-29 CN CN201310322556.8A patent/CN103514457B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102325089A (zh) * | 2011-09-19 | 2012-01-18 | 西安电子科技大学 | 基于差分进化和捕食搜索策略的胖树型片上网络映射方法 |
Non-Patent Citations (3)
Title |
---|
基于TM遥感影像的森林资源线性规划与优化配置研究;董斌等;《生态学报》;20120331;第32卷(第6期);1943-1950 * |
基于自适应交叉概率因子的差分进化算法及其应用;杨卫东等;《信息与控制》;20100430;第39卷(第2期);187-193 * |
多目标优化问题的差分进化算法研究;吴亮红等;《湖南大学学报(自然科学版)》;20090228;第36卷(第2期);53-57 * |
Also Published As
Publication number | Publication date |
---|---|
CN103514457A (zh) | 2014-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103514457B (zh) | 生态约束条件下土地利用结构优化的多目标改进差分进化模型 | |
Liu et al. | Random forest regression evaluation model of regional flood disaster resilience based on the whale optimization algorithm | |
CN107230113A (zh) | 一种多模型融合的房产评估方法 | |
CN106716479A (zh) | 从局部模型导出性能的计算机实现方法 | |
CN106529814A (zh) | 基于Adaboost聚类和马尔科夫链的分布式光伏超短期预测方法 | |
CN103063202A (zh) | 基于遥感影像的蓝藻生物量时空变化监测与可视化方法 | |
CN106845544A (zh) | 一种基于粒子群与支持向量机的小麦条锈病预测方法 | |
CN111119282A (zh) | 一种针对供水管网的压力监测点优化布置方法 | |
CN109711401A (zh) | 一种基于Faster Rcnn的自然场景图像中的文本检测方法 | |
CN104050514A (zh) | 一种基于再分析数据的海浪有效波高的长期趋势预测方法 | |
CN107358294A (zh) | 一种基于Elman神经网络的需水预测方法 | |
CN113837499A (zh) | 一种超短期风电功率预测方法及*** | |
CN104021308A (zh) | 基于ERA-Interim和ERA40预测海浪有效波高的方法 | |
CN108667069A (zh) | 一种基于偏最小二乘法回归的短期风电功率预测方法 | |
CN110069032A (zh) | 一种基于小波神经网络的茄子温室环境智能检测*** | |
CN111831966B (zh) | 一种基于高维概率分布函数的组合河道水位预报方法 | |
CN110490228A (zh) | 一种基于cpp的水文模型参数动态率定方法 | |
CN108549961A (zh) | 一种基于cmip5预估海浪有效波高的方法 | |
CN106228235B (zh) | 一种顾及格局与过程耦合效应的土地利用空间优化配置方法 | |
CN114723142A (zh) | 一种基于非支配排序遗传算法和flus模型的多目标土地利用模拟***与方法 | |
Zhao et al. | Regulation factors driving vegetation changes in China during the past 20 years | |
Yang et al. | Predictor selection for CNN-based statistical downscaling of monthly precipitation | |
Wang et al. | Simulating land use structure optimization based on an improved multi-objective differential evolution algorithm | |
Huang et al. | Agricultural water optimization coupling with a distributed ecohydrological model in a mountain-plain basin | |
CN109472320B (zh) | 不确定条件下作物生育期模型品种参数自动校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160629 Termination date: 20170729 |
|
CF01 | Termination of patent right due to non-payment of annual fee |