CN107944477B - 一种针对城市群的融合多源数据的遗传分类方法 - Google Patents
一种针对城市群的融合多源数据的遗传分类方法 Download PDFInfo
- Publication number
- CN107944477B CN107944477B CN201711127403.2A CN201711127403A CN107944477B CN 107944477 B CN107944477 B CN 107944477B CN 201711127403 A CN201711127403 A CN 201711127403A CN 107944477 B CN107944477 B CN 107944477B
- Authority
- CN
- China
- Prior art keywords
- genetic
- pixel
- value
- ndvi
- ndbi
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/211—Selection of the most significant subset of features
- G06F18/2111—Selection of the most significant subset of features by using evolutionary computational techniques, e.g. genetic algorithms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/176—Urban or other man-made structures
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Software Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Mathematical Physics (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Physiology (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种针对城市群的融合多源数据的遗传分类方法,通过编码、种群初始化、构建适应度函数和遗传过程等步骤获得遗传进化后的种群,本算法以城市群为研究对象,通过结合三种遥感数据:VIIRS DNB、NDVI和NDBI,重新构建遗传算法的适应度函数,通过改进的遗传算法计算出最优阈值,能客观、快速、准确的获得城乡分类结果。
Description
技术领域
本发明涉及一种针对城市群的融合多源数据的遗传分类方法。
背景技术
近年来,全球范围内,城市化进程高速推进,快速城市化影响着全球变化的方方面面,在带来一系列积极变化的同时也产生了许多的环境问题,如城市热岛、城市能耗剧增导致的污染排放问题以及城市生态破坏等。快速、准确、有效地获取城市化时空演变信息对规划城市发展以及解决城市环境问题有重要意义。
遥感技术因其具备进行大面积同步观测的能力、兼顾时效性和经济性等优点可满足城镇扩展研究的需求。高、中分辨率的可见光影像被用于城镇信息提取,但是由于影像成本较高、覆盖范围有限,光谱范围难以区分城市不透水层和乡村独立居民点,可见光遥感数据应用于城镇土地利用覆盖变化的研究受到了一定的限制。
夜间灯光数据能够弥补可见光数据覆盖范围较小、难以区分城乡不透水层的局限性,适用于大范围的城镇信息提取。基于夜间灯光数据提取城镇信息主要分为以下三种方式:分类法、阈值法以及指数法。Cao(2009)提出基于支持向量机分类法,应用DMSP/OLS夜间灯光数据和SPOT-VGT数据提取城市信息,分类法由于应用阈值选择训练样本,样本的选择要求一定的先验知识,训练样本的选择主观性强对结果的影响比较大;Zhang(2013)提出结合NDVI和NTL夜间灯光数据构建调整型夜间灯光指数VANUI用以提取城镇信息,指数构建通常并非以城镇结构作为依据,而是通过统计规律及常识进行构建,主观性强,适用范围有限;阈值法(Henderson et al.2003;Imhoff et al.1997;Liu et al.2012;Sutton 2003)是应用最广泛的一种方法,但是最优阈值的选择一直是很难客观解决的问题。
发明内容
本发明的目的是解决目前传统城乡分类遗传算法主观性强,适用范围窄,计算速度慢,精度低的技术问题。
为实现以上发明目的,本发明提供一种针对城市群的融合多源数据的遗传分类方法,包括如下步骤:
(1)编码
将遥感影像与种群基因相关联,将夜间灯光数据像元值转化为基因值;
(2)种群初始化
根据种群大小、变异概率、交叉概率、最大遗传代数产生随机种群,种群中的每个个体由一个基因序列表示;
(3)构建适应度函数
适应度函数的构建分为以下四步:
a.初始分割
设个体的基因值对应的灯光亮度值为Chi,以Chi对VIIRS DNB进行初始分割,得到初始城乡分布;
b.掩膜
应用上一步得到的初始分割结果对VIIRS DNB、NDVI以及NDBI数据进行掩膜,得到城市覆盖类型和非城市覆盖类型下的VIIRS DNB、NDVI及NDBI影像;
c.计算类间方差和类内方差
根据掩膜结果,将VIIRS DNB、NDVI以及NDBI数据分别划分为初始城区范围和初始乡村范围;以n0和n1分别表示城市部分和乡村部分的像元数,城市部分像元定义为C0,乡村部分像元定义为C1:
C0={u1,u2,…,un0};
C1={n1,n2,…,nn1};
城市部分的像元亮度均值为μ0,乡村部分的像元亮度均值为μ1:
um表示对城市区域像元值u从像元下标m=1至n0进行逐个遍历;
nm表示对农村区域像元值n从像元下标m=1至n1进行逐个遍历;
C0和C1的方差为:
C0和C1占整幅影像的比例为:
Nimg表示整幅影像的像元数量;
类间方差σb 2和类内方差σin 2分别定义为:
μ表示整幅影像的像元平均值;d.计算适应度函数
以类间方差与类内方差的商构建影像分割函数f:
通过分别计算VIIRS DNB、NDVI和NDBI的分割函数fVIIRS、fNDVI和fNDBI来构建适应度函数,可得基因值Chi对应的适应度函数为:
fitness(Chi)=2×fVIIRS+fNDVI+fNDBI;
(4)进行遗传过程
通过选择过程、交叉过程和变异过程进行遗传进化,产生末代种群;
(5)结束遗传过程
当适应度函数值在遗传过程中连续五代没有发生改变时,或者,当遗传代数达到步骤(2)中设置的最大遗传代数时,结束遗传过程。
进一步地,步骤(1)中,依据夜间灯光数据的实际值域区间确定基因编码长度为12bit。
进一步地,步骤(4)中,所述选择过程选用轮转赌盘法进行选择。
与现有技术相比,本发明的有益效果是:
本算法以城市群为研究对象,通过结合三种遥感数据:VIIRS DNB、NDVI和NDBI,重新构建遗传算法的适应度函数,通过改进的遗传算法计算出最优阈值,能客观、快速、准确的获得城乡分类结果。
附图说明
图1是本发明的流程图;
图2是本发明以珠江三角洲为例的算法结果展示图;
图3是本发明计算结果的精度定性评价图;
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
本算法采用了以下研究数据:
1、夜间灯光数据
The Visible-Infrared Imager-Radiometer Suite(VIIRS)是新一代的辐射计,VIIRS Day/Night Band(DNB)数据相比于前一代夜间灯光数据DMSP/OLS有着明显的进步,主要总结为以下三个方面:空间分辨率从1km提高至500m;星上定标使数据质量更高;16bit的数据量所包含的信息要远远多于DMSP的DN值范围(0-63),能够大大缓解夜间灯光数据城市区域灯光饱和的问题。
2、MODIS数据
自2000年起,MODIS数据称为全球尺度研究应用最广泛的数据,应用500m MODIS8-day surface reflectance product(MOD09A1)计算得到Normalized DifferenceVegetation Index(NDVI),Normalized Difference Build-up Index(NDBI)以及ModifiedNormalized Difference water Index(MNDWI)。应用MNDWI数据做水体掩膜,可以将水体的影响剔除;归一化植被指数NDVI用于表征植被覆盖;归一化建筑指数NDBI用于表征建筑用地的土地覆盖类型。
本发明提出基于Genetic Algorithm(GA)遗传算法的分类方法用于自动产生最优阈值以快速、准确地提取城乡信息,如技术流程图,本算法主要分为以下四步:
1、编码
在遗传算法中,编码的过程就是将研究对象与遗传学中基因相关联,所以确定编码规则,是遗传算法应用于实际问题的第一步。将遥感影像与种群基因相关联的方法就是将像元值转化为基因值,本项目依据夜间灯光数据的实际值域区间确定基因编码长度为12bit。
2、种群初始化
确定种群编码规则后,根据种群大小、变异概率、交叉概率、最大遗传代数产生随机种群,本项目中种群中的每个个体由一个基因序列表示,基因的遗传变异会影响个体的表现型。
3、适应性函数构建
“优胜劣汰,适者生存”,适应度函数反映了个体的生存能力,遗传算法的最优解,也是遗传算法中最重要的一部分。本项目适应度函数的构建分为四步:
初始分割
个体的基因值对应的灯光亮度值为Chi,以Chi对VIIRS DNB进行初始分割,得到初始城乡分布。
掩膜
应用上一步得到的初始分割结果对VIIRS DNB,NDVI以及NDBI数据进行掩膜,可以得到城市覆盖类型和非城市覆盖类型下的DNB,NDVI及NDBI影像。
计算类间方差和类内方差
根据掩膜结果,将VIIRS DNB,NDVI以及NDBI数据分别划分为初始城区范围和初始乡村范围。n0和n1分别是城市部分和乡村部分的像元数,城市部分定义为C0,乡村部分定义为C1:
C0={u1,u2,…,un0}
C1={n1,n2,…,nn1}
城市部分的像元亮度均值为μ0,乡村部分的均值为μ1:
um表示对城市区域像元值u从像元下标m=1至n0进行逐个遍历;nm表示对农村区域像元值n从像元下标m=1至n1进行逐个遍历;
C0和C1的方差为:
C0和C1占整幅影像的比例为:
Nimg表示整幅影像的像元数量;
类间方差σb2,类内方差σin2分别定义为:
μ表示整幅影像的像元平均值;
计算适应度函数
已知当图像分类结果类间方差最大类内方差最小,表明影像的分类结果较好,因此以类间方差与类内方差的商构建影像分割函数f:
分别计算VIIRS DNB,NDVI以及NDBI的分割函数,构建适应度函数,计算可得基因值Chi对应的适应度函数值:
fitness(Chi)=2×fVIIRS+fNDVI+fNDBI
遗传过程
遗传算法通过遗传过程进行“优胜劣汰”的选择,遗传过程中的变异与交叉使得该算法能够避免陷于局部最优,从而得到全局最优解。选择过程,本项目选用轮转赌盘法进行选择,适应度函数值高的种群个体更容易被选中,从而进行遗传进化;交叉过程,通过自然选择的个体,以一定的交叉概率进行基因间互换,是遗传变异的一种形式;变异过程,通过个体基因点位的变化,产生新的个体,从而使种群样本更加丰富。
判断遗传过程是否结束
改进的遗传算法是一个循环过程,如果该过程中出现以下两种情况中的任一种,则判断该过程结束:当适应度函数值在遗传过程中连续五代没有发生改变时,说明算法已经搜索到全局最优值,则遗传过程结束;当遗传代数达到在种群初始化中设置的最大遗传代数时,算法结束。
通过构建适用于城乡分类的遗传算法适应度函数以及融合多源数据:VIIRS DNB,NDVI和NDBI可对城市群进行城乡时空信息的提取。如图2所示,以珠江三角洲为例,展示算法的分类结果。图2.(a)(b)(c)分别为珠江三角洲的VIIRS DNB,NDVI和NDBI影像,可以看出明显的城乡区别:城市像元值在VIIS DNB和NDBI影像上高而在NDVI影像上低;相反的,农村像元值在在VIIS DNB和NDBI影像上低而在NDVI影像上高。图2.(d)中本算法的分类结果与三种数据分别展示的城乡分异有相似的空间分布,为了进一步对本算法的计算结果进行评价,着重算法精度评价:
以30米的Landsat8OLI影像的土地利用分类结果作为验证数据,对本项目改进的遗传算法城乡分割结果进行精度评价。本项目中,将城市群作为整体进行研究,但是城市群范围较大,需要多景Landsat影像进行拼接,很难收集到时相相近且受云雾影像小的多幅验证影像。所以精度评价时,将选取城市群中的典型的城市进行评价,即以城市群为研究对象进行遗传算法分割,再从整个城市群的分割结果中裁剪出典型城市对其进行精度评价。
选取珠江三角洲城市群的典型城市广州市对改进的遗传算法进行精度评价,本次精度评价分为定性评价和定量评价两部分。如图3所示,3.(a)是Landsat8OLI的分类结果,通过将遗传算法的分类结果与其进行对比可以发现,应用本文中改进的遗传算法提取的城市区域(3.(c)),与代表真实地表覆盖的Landsat影像分类结果在空间分布上具有较好的一致性;3.(b)是Landsat影像的假彩色显示,对比假彩色影像可以很明显看出本文遗传算法提取的城镇空间信息与实际情况相符。综上,从定性的角度评价遗传算法提取城乡信息具有较高的精度。
表1 定量精度评价
从定量的角度对本文提出的算法进精度评价,如表1,可以看出遗传算法的分类结果与Landsat 8 OLI影像的分类结果相比,差异并不明显,总体精度达到0.91,Kappa系数达到0.71。综上,从定性和定量两个角度对本文提出的结合VIIRS DNB,NDVI和NDBI数据基于遗传算法的城乡分割算法具有较高的精度,能够满足快速、准确、高效的提取城镇信息的要求。
除上述实施方式外,本发明还可以有其他实施方式,凡采用等同替换或等效变换形成的技术方案,均落在本发明的保护范围内。
Claims (3)
1.一种针对城市群的融合多源数据的遗传分类方法,其特征在于,包括如下步骤:
(1)编码
将遥感影像与种群基因相关联,将夜间灯光数据像元值转化为基因值;
(2)种群初始化
根据种群大小、变异概率、交叉概率、最大遗传代数产生随机种群,种群中的每个个体由一个基因序列表示;
(3)构建适应度函数
适应度函数的构建分为以下四步:
a.初始分割
设个体的基因值对应的灯光亮度值为Chi,以Chi对VIIRSDNB进行初始分割,得到初始城乡分布;
b.掩膜
应用上一步得到的初始分割结果对VIIRSDNB、NDVI以及NDBI数据进行掩膜,得到城市覆盖类型和非城市覆盖类型下的VIIRSDNB、NDVI及NDBI影像;
c.计算类间方差和类内方差
根据掩膜结果,将VIIRSDNB、NDVI以及NDBI数据分别划分为初始城区范围和初始乡村范围;以n0和n1分别表示城市部分和乡村部分的像元数,城市部分像元定义为C0,乡村部分像元定义为C1:
C0={u1,u2,…,un0};
C1={n1,n2,…,nn1};
城市部分的像元亮度均值为μ0,乡村部分的像元亮度均值为μ1:
um表示对城市区域像元值u从像元下标m=1至n0进行逐个遍历;
nm表示对农村区域像元值n从像元下标m=1至n1进行逐个遍历;
C0和C1的方差为:
C0和C1占整幅影像的比例为:
Nimg表示整幅影像的像元数量;
类间方差σb 2和类内方差σin 2分别定义为:
μ表示整幅影像的像元平均值;
d.计算适应度函数
以类间方差与类内方差的商构建影像分割函数f:
通过分别计算VIIRSDNB、NDVI和NDBI的分割函数fVIIRS、fNDVI和fNDBI来构建适应度函数,可得基因值Chi对应的适应度函数为:
fitness(Chi)=2×fVIIRS+fNDVI+fNDBI;
(4)进行遗传过程
通过选择过程、交叉过程和变异过程进行遗传进化,产生末代种群;
(5)结束遗传过程
当适应度函数值在遗传过程中连续五代没有发生改变时,或者,当遗传代数达到步骤(2)中设置的最大遗传代数时,结束遗传过程。
2.如权利要求1所述的一种针对城市群的融合多源数据的遗传分类方法,其特征在于,步骤(1)中,依据夜间灯光数据的实际值域区间确定基因编码长度为12bit。
3.如权利要求2所述的一种针对城市群的融合多源数据的遗传分类方法,其特征在于,步骤(4)中,所述选择过程选用轮转赌盘法进行选择。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711127403.2A CN107944477B (zh) | 2017-11-15 | 2017-11-15 | 一种针对城市群的融合多源数据的遗传分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711127403.2A CN107944477B (zh) | 2017-11-15 | 2017-11-15 | 一种针对城市群的融合多源数据的遗传分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107944477A CN107944477A (zh) | 2018-04-20 |
CN107944477B true CN107944477B (zh) | 2019-07-02 |
Family
ID=61932247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711127403.2A Active CN107944477B (zh) | 2017-11-15 | 2017-11-15 | 一种针对城市群的融合多源数据的遗传分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107944477B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3614268B1 (de) * | 2018-08-21 | 2021-07-07 | Siemens Aktiengesellschaft | Verfahren und vorrichtung zum verarbeiten von daten mittels codierter operationen |
CN111742329B (zh) * | 2020-05-15 | 2023-09-12 | 安徽中科智能感知科技股份有限公司 | 一种基于多源遥感数据融合和深度神经网络的矿山开采典型地物动态监测方法及平台 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100383805C (zh) * | 2005-11-03 | 2008-04-23 | 复旦大学 | 星载微波和红外遥感融合对城区地物特征进行分类的方法 |
CN102592181A (zh) * | 2011-12-30 | 2012-07-18 | 中国农业科学院农业资源与农业区划研究所 | 一种作物种植面积统计数据空间分布的优化方法 |
US20150012471A1 (en) * | 2013-07-02 | 2015-01-08 | SenseGon Technologies Ltd. | System for user psychosocial profiling |
CN105528596B (zh) * | 2016-02-03 | 2018-10-09 | 长江大学 | 利用阴影的高分辨率遥感影像建筑物自动提取方法及系统 |
-
2017
- 2017-11-15 CN CN201711127403.2A patent/CN107944477B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107944477A (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Attribution of climate and human activities to vegetation change in China using machine learning techniques | |
Xu et al. | Classification of local climate zones using ASTER and Landsat data for high-density cities | |
CN103942420B (zh) | 一种建筑物尺度的太阳直射辐射能快速估算方法 | |
Li et al. | Exploring the relationship between the 2D/3D architectural morphology and urban land surface temperature based on a boosted regression tree: A case study of Beijing, China | |
Huang et al. | Detection and prediction of land use change in Beijing based on remote sensing and GIS | |
Wang et al. | Impacts of urbanization on soil organic carbon stocks in the northeast coastal agricultural areas of China | |
Wang et al. | Classification of land-cover types in muddy tidal flat wetlands using remote sensing data | |
Zheng et al. | The desaturation method of DMSP/OLS nighttime light data based on vector data: Taking the rapidly urbanized China as an example | |
Liu et al. | High-resolution mapping of mainland China’s urban floor area | |
CN107944477B (zh) | 一种针对城市群的融合多源数据的遗传分类方法 | |
Han et al. | How do 2D/3D urban landscapes impact diurnal land surface temperature: Insights from block scale and machine learning algorithms | |
Ma et al. | A multisource remotely sensed data oriented method for “ghost city” phenomenon identification | |
Guo et al. | New morphological features for urban tree species identification using LiDAR point clouds | |
Qin et al. | Estimating aboveground carbon stocks of urban trees by synergizing ICESat-2 LiDAR with GF-2 data | |
Sun et al. | Population spatial distribution based on Luojia 1–01 nighttime light image: A case study of Beijing | |
Zhou et al. | Automatic responsive-generation of 3D urban morphology coupled with local climate zones using generative adversarial network | |
Yu et al. | Downscaling mapping method for local climate zones from the perspective of deep learning | |
Li et al. | Oasis functional stability evaluation based on multiple Indicators, Northwest China | |
Lu et al. | Mapping the terraces on the Loess Plateau based on a deep learning-based model at 1.89 m resolution | |
Li | Dynamic monitoring algorithm of natural resources in scenic spots based on MODIS Remote Sensing technology | |
Zhou et al. | Comparison of object-oriented and Maximum Likelihood Classification of land use in Karst area | |
Liu et al. | Exploring the seasonal effects of urban morphology on land surface temperature in urban functional zones | |
Le et al. | Case Study of GIS Application in Analysing Urban Heating Island Phenomena in Tropical Climate Country | |
Alphan | Incorporating visibility information into multi-criteria decision making (MCDM) for wind turbine deployment | |
Cong et al. | Temporal and spatial changes of the Yellow River Delta Wetland based on multi-source data during 30 years |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |