CN111666918A - 一种基于多因素的海岸线变化识别方法 - Google Patents
一种基于多因素的海岸线变化识别方法 Download PDFInfo
- Publication number
- CN111666918A CN111666918A CN202010572802.5A CN202010572802A CN111666918A CN 111666918 A CN111666918 A CN 111666918A CN 202010572802 A CN202010572802 A CN 202010572802A CN 111666918 A CN111666918 A CN 111666918A
- Authority
- CN
- China
- Prior art keywords
- data
- coastline
- clustering
- time
- time period
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000008859 change Effects 0.000 title claims abstract description 26
- 230000003628 erosive effect Effects 0.000 claims abstract description 13
- 230000006870 function Effects 0.000 claims description 57
- 238000012937 correction Methods 0.000 claims description 27
- 230000005855 radiation Effects 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 12
- 238000002310 reflectometry Methods 0.000 claims description 11
- 238000005192 partition Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000004062 sedimentation Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000001932 seasonal effect Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 239000000443 aerosol Substances 0.000 claims description 3
- 238000013145 classification model Methods 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000002834 transmittance Methods 0.000 claims description 3
- 230000006399 behavior Effects 0.000 abstract description 3
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 238000006467 substitution reaction Methods 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 15
- 238000012800 visualization Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000013499 data model Methods 0.000 description 2
- 238000013524 data verification Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 241000282412 Homo Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000004071 biological effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 238000003920 environmental process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 206010020718 hyperplasia Diseases 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/13—Satellite images
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Geometry (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种基于多因素的海岸线变化识别方法,将海岸线变化(侵蚀、增长、稳定性);NDVI(归一化差分植被指数);和沉降演变作为随时间的变量影响因素,通过结合多个变量的替代方法,从而对与复杂环境系统相关的问题进行建模,并消除不精确和主观的概念,充分挖掘了利用地理信息中的遥感数据的潜力,显示了局部动力学的复杂行为,从而为环境问题和海岸带综合管理增加了有用的、实质性的信息。
Description
技术领域
本发明归属于海岸线变化识别技术领域,具体涉及一种基于多因素的海岸线变化识别方法。
背景技术
海岸是人类社会经济活动最活跃、最集中的地区,与海运交通、领土主权、工程建设、资源开发、空间利用等诸多方面密切相关。海岸天然的地理优势,孕育了一个又一个自“大航海”时代以来的世界经济中心。进入21世纪,沿海城市在资源、环境、交通上优势更加明显,世界沿海国家经济重心向滨海地区转移;新时代我国加快推进海洋强国建设,海岸作为探索海洋的出发点与归宿点,正确理解、合理认知、科学探究海岸空间为实现海洋强国战略目标提供重要保障。
海岸线是海洋和陆地的分界线,是海岸乃至地球上最重要的边界线。全球海岸线总长44万千米,分布在100多个国家和南北极地区,全球已有超过一半的人口居住在距海岸线100千米内的沿海地区;我国拥有超过1.8万千米的大陆岸线和1.4万千米的岛屿岸线,沿岸滩涂超过2万平方千米,探明海滨矿藏储量15.3亿吨,沿海地区集中了40%的人口和60%的产值。海岸线的位置、走向和形态不仅体现了全球及海岸带环境过程,更是社会经济活动、人类综合作用的结果与反映。此外,海岸线作为陆海交界,除受到一般陆地和地表过程影响外,还受到潮汐、波浪、环流、信风等海洋动力作用以及生物活动影响,具有极大的复杂性和时空不确定性。
发明内容
本发明提供一种基于多因素的海岸线变化识别方法,提出了用于海岸区人类影响分类的模糊模型,其集成了海岸线变化、归一化差分植被指数和沉降影响,以通过地理信息系统(GIS)的图形可视化过程来增强数值语言模糊分类。显示了局部动力学的复杂行为,从而为环境问题和海岸带综合管理增加了有用的、实质性的信息。
一种基于多因素的海岸线变化识别方法,包括如下步骤:
S1.输入海岸线基线,采用遥感卫星图像覆盖所述海岸线区域的面积延伸;
S2.将来自每个所述遥感卫星场景的海岸线分割为多个长条多边形形状;
S3.采用模糊模型使用所述遥感卫星图像来提取关于时间变化的信息:(i)海岸线变化;(ii)归一化差分植被指数;和(iii)沉降演变;
S4.分别计算长条多边形的侵蚀、堆积和稳定的面积,所述侵蚀、堆积和稳定的面积被计算为相对于每个多边形区总面积的百分比,并将百分比用作第一输入(海岸线变化)的三个变量(X1、X2和X3);
S5.计算归一化差分植被指数(NDVI)作为第四变量X4,NDVI=(NIR-RED)/(NIR+RED),NIR和RED分别表示近红外和红色反射率;
S6.分别建立针对每个多边形区和时间图像分类区域的建筑物累积变量及沉降演变,作为海岸线模糊分类模型的第五变量X5和第六变量Y;
S7.采用模糊模型对每个变量X1、X2、X3、X4和Y进行计算,
S8.以交互形式对每个变量X1、X2、X3、X4和Y进行验证及微调,直到满足验证标准;
S9.对验证后的变量X1、X2、X3、X4和Y值采用特定的语言进行标记:低、中、高,根据标记结果确定海岸带的人为影响分类等级。
进一步地,定义快照质量函数(Snapshot Quality)和历史代价函数(HistoryCost),快照质量函数衡量的是当前时间段变量X1、X2、X3、X4和Y的聚类划分质量,表示聚类结果和当前时间段的数据的匹配程度;历史代价函数衡量的是相邻时间段变量X1、X2、X3、X4和Y的聚类划分是否具有较好的时序平滑特性,表示的是历史模型与当前时间段的聚类结果的差异;对于包含T个时间序列数据S={X1,X2,X3,X4,X5},且每个时间点的数据集X包含N个样本,xτ={x1,x2,x3,x4,x5}。定义如下的快照质量函数sq(Cτ,Mτ)以及历史代价函数hc(Cτ,Cτ-1):
其中,Cτ表示τ时间段数据的聚类结果,Mτ是时间段数据的相似度矩阵。
进一步地,步骤S8中对每个变量X1、X2、X3、X4和Y的验证方法包括如下步骤:
S11.对τ时间段的数据进行聚类时,将目标函数定义为:
将函数sq(Cτ,Mτ)及hc(Cτ,Cτ-1)带入目标函数J得到全局目标函数:
式中,c、N分别是聚类数和样本数,表示τ时间段的数据与τ时间段第i个分段(第i个聚类)的距离,表示τ时间段数据属于第i个聚类的隶属度,且满足每个样本属于所有类的隶属程度之和等于1,Cτ指τ时间段的隶属矩阵,Cτ∈{0,1}c*k,而Ci,k=1表示τ时间段数据属于第i类,各数据在各时间段内只能属于某个确定的类,即∑τCi,k=1;
S12.最小化全局目标函数以求得最优分割参数,通过拉格朗日定理求解,定义下式为拉格朗日目标函数:
S14.将数据的时间坐标作为一个额外变量并使其参与聚类过程,迭代计算方式如下:
计算τ时间段c个聚类中心vi:
计算聚类中心的模糊协方差矩阵Fi:
计算距离函数值D2(xk,vi):
利用距离函数值D2(xk,vi)更新模糊划分矩阵U:
当满足条件||U(l+1)-U(l)<ε||时,中止聚类算法,否则增加迭代次数,使l=l+1,转到步骤S11,重复上述步骤,直至满足该条件为止。
进一步地,将所述遥感卫星图像进行图像预处理,进行大气校正,进而获得每个图像元素校正后的全光谱信息,且所述遥感卫星图像在所述海岸线区域上具有小于10%的云量覆盖。
进一步地,利用大气校正模型去除天空光和大气散射效应,使遥感图像更准确地反映特征的光谱值,大气校正模型由下式获得:
ρ为地表反射率,ρΤ(θs,θv,φs,φv)为大气上界反射率,θs是太阳天顶角,θs是太阳方位角,θv是传感器方位角,T(θs)是大气透射率,Tg(θs,θv)是太阳-目标大气路径透射率,T(θv)是目标-传感器大气路径透射率,ρR+α是分子散射和气溶胶散射所构成的路径辐射反射率,s是大气半球反射率。
进一步地,所述图像预处理采用暗对象减影方法,将由于太阳照度、传感器观察几何形状和季节性变化而引起的卫星图像辐射度差异进行校正和归一化。
进一步地,将所述遥感卫星图像进行图像预处理,进行几何校正,具体的操作步骤是地面控制点选择,像素坐标变换和像素亮度值的重新采样,像元坐标变换采用的是基于坐标变换的的方法,变换的公式为
x’=a0+a1x+a2y+a3x2+a4xy+a5y2
y'=b0+b1x+b2y+b3x2+b4xy+b5y2
其中,x’,y’为校正后输出的控制点坐标,x,y为地面控制点在原图像中的坐标,RMSerror为每个控制点的均方根误差。
进一步地,模糊逻辑可以根据模型设计分类为低、低/中、中/高或高。
由于受潮汐、海岸地形等因素的影响较大,卫星过境时水边线刚好位于平均大潮高潮线的影像是很难获取的,海岸线自动解译算法的研宄大多是基于遥感影像中水边线的提取,即卫星过境时刻所记录的海陆分界线。本发明提出了用于海岸区人类影响分类的模糊模型,其集成了海岸线变化、归一化差分植被指数和沉降影响,以通过地理信息系统(GIS)的图形可视化过程来增强数值语言模糊分类。诊断全世界沿海地区的人为影响(即与人类活动有关的影响),例如沿海开发和规划、过度捕捞等。沿海环境的分类、风险和脆弱性评估是必不可少的环节,因此,本发明的算法最终显示了局部动力学的复杂行为,从而为环境问题和海岸带综合管理增加了有用的、实质性的信息。
本发明首次提出将海岸线变化(侵蚀、增长、稳定性);NDVI(归一化差分植被指数);和沉降演变作为随时间的变量影响因素,通过结合多个变量的替代方法,从而对与复杂环境系统相关的问题进行建模,并消除不精确和主观的概念,充分挖掘了利用地理信息中的遥感数据的潜力,并通过包括社会经济数据在内的其他影响因素以增强数值模糊分类的准确率。通过生成的空间地图实现图形可视化,这有利于地理特征的映射增强对每个部门演化模式识别的区分。
本发明的模型输出表示分数是从0到1的数字,其可转换为模糊语言分类变量;即低、中和高。同时,通过NDVI(归一化差分植被指数)使用GIS通过图形增强了可视化程度。
本发明在对当前时间段的数据进行聚类时会考虑上一时间段的聚类结果。通过本发明的算法,可以有效地处理分析具有非线性及强耦合特性的海岸带数据,使得改进算法能够实时分析线性时间序列数据,加入了快照质量函数和历史代价函数,增强了聚类结果和当前时间段数据的匹配程度,且使得相邻时间段的聚类划分具有良好的时序平滑特性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。
图1为本发明海岸线变化识别方法的流程图。
具体实施方式
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
沿海岸带探测人类影响的方法包括海岸线侵蚀/增生模式评估,具体地,通过(a)地形剖面分析(考虑到跨岸形态以及作用在海滩上的破坏力和建设力之间的平衡);(b)海岸线变化率;(c)土地利用/土地覆盖监测结合利用地理信息系统可视化的海岸线变化率。
本发明结合多个变量的替代方法,从而对与复杂环境系统相关的问题进行建模,并消除不精确和主观的概念,充分挖掘了利用地理信息中的遥感数据的潜力,通过包括社会经济数据在内的其他影响因素以增强数值模糊分类的准确率。通过生成的空间地图实现图形可视化,这有利于地理特征的映射增强对每个部门演化模式识别的区分。
一方面,模型输出表示分数是从0到1的数字,其可转换为模糊语言分类变量;即低、中和高。另一方面,通过NDVI(归一化差分植被指数)使用GIS通过图形增强了可视化程度。
从遥感数据、海岸线变化、NDVI计算和结算影响中提取海岸线位置。模糊模型被设计为具有五个变量,即侵蚀、增长、稳定性、NDVI和累积。所有的语言学标签(模糊集)、隶属函数、模糊规则和去模糊化提供了输出,该输出是代表海岸带人为影响分类的明确数字。
第一步:数据处理
输入基线,使用遥感卫星图像来覆盖研究区域的面积延伸。卫星数据集是考虑到每年的相同/最近月份(例如八月和九月)而选择的,以寻求通过最小化季节变化来增加土地利用类别的分离。此外,所有选择的图像在研究区域上应当具有小于10%的云量覆盖。
由于遥感数据受到大气影响等诸多因素的影响,因此这些数据集无法用于进一步分析。卫星图像只能在执行若干图像预处理步骤(包括大气校正、几何校正)之后才能使用,以消除或最小化大气影响,进而获得每个图像元素(像素)的校正后的全光谱信息。本发明的大气校正方法可将由于太阳照度、传感器观察几何形状和季节性变化而引起的卫星图像辐射度差异进行校正和归一化。
作为输入,考虑三个方面,模糊模型使用前面描述的卫星图像来提取关于时间变化的信息:(i)海岸线变化;(ii)NDVI(归一化差分植被指数);和(iii)沉降演变。
以海岸岛为例,由于研究区域在所有侧面(岛)都被水包围,因此在类似的缩放水平(1:5000的均匀比例)下,使用屏幕上手动数字化技术,将来自每个卫星场景的海岸线分割为多边形形状。这样可以有效描绘海岸线。
计算侵蚀和堆积面积(长条多边形),对于每两个连续的多边形,使用GIS环境中的空间联合工具,侵蚀、堆积和稳定性的面积被计算为相对于每个多边形区的总面积的百分比,并且那些值(%)被用作第一输入(海岸线变化)的三个变量(X1、X2和X3)。
NDVI(归一化差分植被指数)模糊模型中使用的第二个输入数据集是NDVI,该指数由图像中每个像素的新计算值组成,范围从-1到+1。NDVI由方程(1)和两个所需的输入波段,即近红外(NIR)和红色(red)反射率计算。
NDVI=(NIR-RED)/(NIR+RED)
在执行图像预处理(包括大气校正、几何校正)后,计算更具代表性的植被覆盖指数时需要反射值,得到每个多边形区域的NDVI平均值,并用作表示第四个变量(X4)的第二输入。
卫星传感器在获取地物光谱的过程中,来自太阳的辐射能量第一次穿越大气层会被吸收一部分,剩余部分到达地面后,会以透射和反射的方式与地物相互作用,这两种方式的辐射能量之和等于上述剩余部分的辐射能量,在上升通过大气层的过程中,地面物体反射的辐射能量将被大气吸收,最终传感器获得的辐射能量,一部分来自地面物体的反射辐射,另一部分来自大气的反射和散射的辐射能量。遥感影像记录的大气辐射能量,给影像造成一定的失真,为了正确评估遥感影像地物辐射值,就需要去除这些辐射噪声,从遥感图像中去除辐射噪声的过程称为大气校正,目前存在较多的大气校正理论模型,针对不同的遥感影像类型和状况选取合适的模型才可以获得较好的处理结果。
本发明遥感卫星图像在海岸线区域上具有小于10%的云量覆盖。本发明设计具有较好的辐射校正精度的大气校正模型,该模型考虑了地面非朗伯体情况。该模型使用的地面目标反射率与传感器入口处反射率关系式由式(1)所示:
在卫星飞行的过程中,会出现侧滚和距离地面高度的变化,这将会使获取的遥感影像发生几何形变。地面数据接收站会对获取的遥感影像做几何粗校正处理,在实际应用中,一般需要对遥感数据再进行几何精校正,在几何精细校正过程中,需要地面控制点GCPS(地面控制点)。具体的操作步骤是地面控制点选择,像素坐标变换和像素亮度值的重新采样。像元坐标变换采用的是基于坐标变换的的方法,变换的公式为
x’=a0+a1x+a2y+a3x2+a4xy+a5y2 (2)
y’=b0+b1x+b2y+b3x2+b4xy+b5y2 (3)
其中,x’,y’--校正后输出的控制点坐标,x,y--地面控制点在原图像中的坐标,RMSerror--每个控制点的均方根误差。
采用上述变换公式,将原始影像中的像元值坐标(x’,y’)转换为输出影像中的像元坐标值(x,y),通过计算每个地面控制点的均方根误差(RMSerror),可以验证校正模型的有效性。
其中,ρ为地表反射率,ρΤ(θs,θv,φs,φv)为大气上界反射率,θs是太阳天顶角,θs是太阳方位角,θv是传感器方位角,T(θs)是大气透射率,Tg(θs,θv)是太阳-目标大气路径透射率,T(θv)是目标-传感器大气路径透射率,ρR+α是分子散射和气溶胶散射所构成的路径辐射反射率,s是大气半球反射率。
利用大气校正模型能够去除天空光和大气散射效应,实现遥感图像的大气校正。使遥感图像更准确地反映特征的光谱值。
沉降演变所使用的第三输入数据(Y)是沉降影响(构建区域)。海岸线附近的基础设施和建筑物可直接影响海岸侵蚀以及洪水。在海岸线附近保持地形学方面的最小距离对于减少海岸带对人类的影响是非常重要的,然而,在本发明中,可以观察到,即随着时间的推移,在海岸线附近的沉降增加。采用数学函数(三角形或L-函数)根据特定范围和变量单位(X1、X2、X3、X4和Y)对每个特定语言标记(低、中、高)。
而且,本发明基于对象的算法,通过分析集成邻域信息,将任何图像视为对象,这将增强分析并增加分类图像即LULC的准确性。因此,对于此,使用分割方法才用特征提取工具从每个卫星图像提取LULC类别(建筑物、植物和其它)。在该过程期间,测试缩放和合并级以获得三个类别的最佳结果,包括在所有卫星图像中的构建区域。由于模糊输入需要精确的结果,因此,将分割的光栅转换为矢量数据集,以更精确地描绘三个LULC类别在编辑会话期间GIS环境中的每年的多边形区域。
建立针对每个扇区和时间图像分类区域的变量,然后用作海岸线模糊分类模型的第五变量(X5)。
模糊模型设计综合三个输入:海岸线变化、NDVI和沉降影响(建成区),开发模糊模型设计。这三个输入由五个变量(X1、X2、X3、X4和X5)组成。从卫星图像中提取基线信息,提取所检测到的所有输入变量,根据具体的变量特征,输入不同的范围和单位。在这种情况下,考虑到多边形海岸线变化的总量(%),X1、X2、X3(海岸线变化)在0到100(%)的范围内。NDVI(变量X4),范围从-1.0到1.0,X5(建筑物)范围从0到100,然后通过时间变化进行评估。该模糊模型的输出是范围从0到1的数字,对海岸带人为影响进行分类分级。当输出数接近1时,其表明高的人为影响分类,而接近0是指低的人为影响分类,在该范围之间,模糊逻辑可以根据模型设计分类为低、低/中、中/高或高。所提出的模糊模型中使用的推理方法是基于模糊规则概念和表示输出的模型种类,通过由每个推理规则的聚合产生模糊集,在模糊模型中,第一输入(海岸线变化)被分成三个变量(基于海岸线的状态)、侵蚀(X1)、堆积(X2)和稳定(X3),考虑到比较连续几年检测到的变化。为该变量考虑的语言标签被命名为低、中和高。根据表中给出的参数,所选择的隶属函数的类型是三角形或L函数。
最后,使用三个输入(海岸线变化、NDVI和沉降演变影响)、(建立区域)与五个变量(X1、X2、X3、X4和X5)构成模糊规则。规则由五个变量(X1,X2,X3,X4和X5)和它们各自的语言学标记(A1,A2,A3),(B1,B2,B3),(C1,C2,C3),(D1,D2,D3),(E1,E2,E3)组成。通过对五个变量进行积分来定义最终模糊规则输出Y(F1,F2,F3),其语言学标签使用"如果-则"规则格式。
例如,If X1εA1 And X2εB1 And X3εC3 And X4εD3 And X5εE1 Then YεF1,使用例如语言变量表达相同规则的方式:“如果侵蚀低且增长低且稳定高且NDVI高且累积低,则输出海岸带人类影响分类低”。所有这些变量和规则的集合需要被验证,否则其可能被分类为任意估计。在这种情况下,进行一些初步测试以便以交互形式微调这些功能的规则和参数,直到满足验证标准。验证步骤用于确定所实现的最终输出(模糊海岸带人为影响分类)的准确性和质量。
传统的数据验证问题中,数据集往往是静态且数量固定的。但在海岸线变化问题中,数据会具有两个重要特性,即演化特性和数量特性。演化特性是指数据的统计性质(例如均值、方差等指标)随着时间动态变化,具有潜在的内部变化规律以及时变性。这种内部变化规律既可能是由于该数据所描述的海岸线系统自身的自然且合理的长期变化,也有可能是由于数据噪声干扰导致的短期变动。数量特性则是指数据的规模庞大。由于在每一个固定的时间单位都会产生一批新的数据,系统的数据常常处于不断更新中,从而造成了总数目不断增加。因此,具有演化特性和数量特性的数据被称为演化数据。该类数据给数据验证带来了较大的难度。数据的演化特性要求算法能够处理数据随时间变化的趋势,合理的分析相邻时间段之间数据所发生的转变。数据的数量特性则要求算法能够在有限次遍历内完成分析,而非将整体数据全部存储于内存再运行算法聚类。
为了使聚类结果能够体现数据在时间域上的关联性并实时获得聚类的结果,本发明在对当前时间段的数据进行聚类时会考虑上一时间段的聚类结果。具体来说,每个新的时间段都有一批新的数据等待着被聚类算法归入某些聚类簇中。如果他们没有背离历史期望,而是与历史聚类结果相似,那么整个数据的聚类结果应该是与历史聚类结果相近的,聚类算法生成的整体聚类模型应该也是类似的。但如果他们的内部结构变化剧烈,那么算法应该修改聚类结果以反映出这种变化。演化聚类算法一方面要求在分析当前时间段的数据的同时能够结合之前时间段的模型,得出多个时间段之间数据统计信息的联系,另一方面要求聚类算法能够应对大规模的数据进行实时聚类得到结果。因此,为了刻画这两方面因素对聚类结果的影响程度,定义如下两个度量标准函数:快照质量函数(SnapshotQuality)和历史代价函数(HistoryCost)。快照质量函数衡量的是当前时间段变量X1、X2、X3、X4和Y的聚类划分质量,表示聚类结果和当前时间段的数据的匹配程度。历史代价函数衡量的是相邻时间段变量X1、X2、X3、X4和Y的聚类划分是否具有较好的时序平滑特性,表示的是历史模型与当前时间段的聚类结果的差异。
所以,全局目标函数包含两个部分。一部分是对当前样本数据聚类质量的度量,通过快照质量函数sq(Cτ,Mτ)衡量。Cτ表示τ时间段数据的聚类结果,Mτ是时间段数据的相似度矩阵。该函数值越小,表示聚类后的模型和该时间段数据的潜在真实模型越符合。另一部分则是历史数据聚类结果与当前时间段的聚类结果之间的差异度,通过历史代价函数hc(Cτ,Cτ-1)衡量。该函数值反映了当前时间段的数据模型和上一时间段的数据模型之间的差异度。该值越小,表示相邻时间段的聚类结果越相似。故在对τ时间段的数据进行聚类时,为了同时考虑上述两个度量函数,将目标函数定义为:
对于包含T个时间序列数据S={X1,X2,X3,X4,X5},且每个时间点的数据集X包含N个样本,xτ={x1,x2,x3,x4,x5}。定义如下的快照质量函数sq(Cτ,Mτ)以及历史代价函数hc(Cτ,Cτ-1):
将上述函数(3)、(4)带入式(2),得到全局目标函数:
式中,c、N分别是聚类数和样本数,表示τ时间段的数据与τ时间段第i个分段(第i个聚类)的距离,表示τ时间段数据属于第i个聚类的隶属度,且满足每个样本属于所有类的隶属程度之和等于1,Cτ指τ时间段的隶属矩阵,Cτ∈{0,1}c*k,而Ci,k=1表示τ时间段数据属于第i类,各数据在各时间段内只能属于某个确定的类,即∑τCi,k=1。
最小化全局目标函数以求得最优分割参数,通过拉格朗日定理求解,定义函数(6)为拉格朗日目标函数:
在时序分割问题中,每个聚类的数据点都必须符合自身固有的时间顺序,因此将数据的时间坐标作为一个额外变量并使其参与聚类过程,所得到的聚类结果便是动态聚类结果的计算方法。迭代计算方式如下:
计算τ时间段c个聚类中心vi:
计算聚类中心的模糊协方差矩阵Fi:
计算距离函数值D2(xk,vi):
利用距离函数值D2(xk,vi)更新模糊划分矩阵U:
当满足条件||U(l+1)-U(l)<ε||时,中止聚类算法,否则增加迭代次数,使l=l+1,转到步骤(2),重复上述步骤,直至满足该条件为止。
通过本发明的算法,可以有效地处理分析具有非线性及强耦合特性的海岸带数据,使得改进算法能够实时分析线性时间序列数据,加入了快照质量函数和历史代价函数,增强了聚类结果和当前时间段数据的匹配程度,且使得相邻时间段的聚类划分具有良好的时序平滑特性。
以上所述仅是对本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (8)
1.一种基于多因素的海岸线变化识别方法,其特征在于,包括如下步骤:
S1.输入海岸线基线,采用遥感卫星图像覆盖所述海岸线区域的面积延伸;
S2.将来自每个所述遥感卫星场景的海岸线分割为多个长条多边形形状;
S3.采用模糊模型使用所述遥感卫星图像来提取关于时间变化的信息:(i)海岸线变化;(ii)归一化差分植被指数;和(iii)沉降演变;
S4.分别计算长条多边形的侵蚀、堆积和稳定的面积,所述侵蚀、堆积和稳定的面积被计算为相对于每个多边形区总面积的百分比,并将百分比用作第一输入(海岸线变化)的三个变量(X1、X2和X3);
S5.计算归一化差分植被指数(NDVI)作为第四变量X4,NDVI=(NIR-RED)/(NIR+RED),NIR和RED分别表示近红外和红色反射率;
S6.分别建立针对每个多边形区和时间图像分类区域的建筑物累积变量及沉降演变变量,作为海岸线模糊分类模型的第五变量X5和第六变量Y;
S7.采用模糊模型对每个变量X1、X2、X3、X4和Y进行计算,
S8.以交互形式对每个变量X1、X2、X3、X4和Y进行验证及微调,直到满足验证标准;
S9.对验证后的变量X1、X2、X3、X4和Y值采用特定的语言进行标记:低、中、高,根据标记结果确定海岸带的人为影响分类等级。
2.根据权利要求1所述的一种基于多因素的海岸线变化识别方法,其特征在于,定义快照质量函数(Snapshot Quality)和历史代价函数(History Cost),快照质量函数衡量的是当前时间段变量X1、X2、X3、X4和Y的聚类划分质量,表示聚类结果和当前时间段的数据的匹配程度;历史代价函数衡量的是相邻时间段变量X1、X2、X3、X4和Y的聚类划分是否具有较好的时序平滑特性,表示的是历史模型与当前时间段的聚类结果的差异;对于包含T个时间序列数据S={X1,X2,X3,X4,X5},且每个时间点的数据集X包含N个样本,xτ={x1,x2,x3,x4,x5}。定义如下的快照质量函数sq(Cτ,Mτ)以及历史代价函数hc(Cτ,Cτ-1):
其中,Cτ表示τ时间段数据的聚类结果,Mτ是时间段数据的相似度矩阵。
3.根据权利要求2所述的一种基于多因素的海岸线变化识别方法,其特征在于,步骤S8中对每个变量X1、X2、X3、X4和Y的验证方法包括如下步骤:
S11.对τ时间段的数据进行聚类时,将目标函数定义为:
将函数sq(Cτ,Mτ)及hc(Cτ,Cτ-1)带入目标函数J得到全局目标函数:
式中,c、N分别是聚类数和样本数,表示τ时间段的数据与τ时间段第i个分段(第i个聚类)的距离,表示τ时间段数据属于第i个聚类的隶属度,且满足每个样本属于所有类的隶属程度之和等于1,Cτ指τ时间段的隶属矩阵,Cτ∈{0,1}c*k,而Ci,k=1表示τ时间段数据属于第i类,各数据在各时间段内只能属于某个确定的类,即∑τCi,k=1;
S12.最小化全局目标函数以求得最优分割参数,通过拉格朗日定理求解,定义下式为拉格朗日目标函数:
S14.将数据的时间坐标作为一个额外变量并使其参与聚类过程,迭代计算方式如下:
计算τ时间段c个聚类中心vi:
计算聚类中心的模糊协方差矩阵Fi:
计算距离函数值D2(xk,vi):
利用距离函数值D2(xk,vi)更新模糊划分矩阵U:
当满足条件||U(l+1)-U(l)<ε||时,中止聚类算法,否则增加迭代次数,使l=l+1,转到步骤S11,重复上述步骤,直至满足该条件为止。
4.根据权利要求2所述的一种基于多因素的海岸线变化识别方法,其特征在于,将所述遥感卫星图像进行图像预处理,进行大气校正,进而获得每个图像元素校正后的全光谱信息,且所述遥感卫星图像在所述海岸线区域上具有小于10%的云量覆盖。
6.根据权利要求1-3所述的一种基于多因素的海岸线变化识别方法,其特征在于,所述图像预处理采用暗对象减影方法,将由于太阳照度、传感器观察几何形状和季节性变化而引起的卫星图像辐射度差异进行校正和归一化。
8.根据权利要求1-7所述的一种基于多因素的海岸线变化识别方法,其特征在于,模糊逻辑可以根据模型设计分类为低、低/中、中/高或高。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010572802.5A CN111666918B (zh) | 2020-06-22 | 2020-06-22 | 一种基于多因素的海岸线变化识别方法 |
KR1020217020898A KR20220000898A (ko) | 2020-06-22 | 2021-02-24 | 다인자를 기반으로 한 해안선 변화 식별방법 |
PCT/CN2021/077565 WO2021258758A1 (zh) | 2020-06-22 | 2021-02-24 | 一种基于多因素的海岸线变化识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010572802.5A CN111666918B (zh) | 2020-06-22 | 2020-06-22 | 一种基于多因素的海岸线变化识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111666918A true CN111666918A (zh) | 2020-09-15 |
CN111666918B CN111666918B (zh) | 2023-08-15 |
Family
ID=72389255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010572802.5A Active CN111666918B (zh) | 2020-06-22 | 2020-06-22 | 一种基于多因素的海岸线变化识别方法 |
Country Status (3)
Country | Link |
---|---|
KR (1) | KR20220000898A (zh) |
CN (1) | CN111666918B (zh) |
WO (1) | WO2021258758A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112926465A (zh) * | 2021-03-02 | 2021-06-08 | 中国人民解放军战略支援部队信息工程大学 | 基于点云类型的海岸线性质识别方法及装置 |
CN113673582A (zh) * | 2021-07-30 | 2021-11-19 | 西南交通大学 | 基于系统聚类分析的铁路动态基准点多层级分群方法 |
WO2021258758A1 (zh) * | 2020-06-22 | 2021-12-30 | 大连海洋大学 | 一种基于多因素的海岸线变化识别方法 |
CN115439748A (zh) * | 2022-09-13 | 2022-12-06 | 中山大学 | 一种海岸线侵蚀程度的监测方法、装置以及电子设备 |
CN116124181A (zh) * | 2023-04-14 | 2023-05-16 | 国家海洋技术中心 | 一种潮汐观测设备的现场校准方法及系统 |
CN116797798A (zh) * | 2023-07-03 | 2023-09-22 | 生态环境部卫星环境应用中心 | 基于尺度变换的海岸线快速提取方法和装置 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109598310B (zh) * | 2018-12-25 | 2023-04-14 | 核工业北京地质研究院 | 一种多因子敏感设施识别方法 |
CN114462247B (zh) * | 2022-02-14 | 2022-10-21 | 中国人民解放军61540部队 | 一种北太平洋海表盐度年代际模态识别方法及系统 |
KR102526638B1 (ko) | 2022-04-21 | 2023-04-28 | (주) 지오씨엔아이 | 다중시기 위성영상과 인공지능을 활용한 시계열 해안선 탐지 장치 및 그 장치의 구동방법 |
CN114781537B (zh) * | 2022-05-07 | 2024-04-16 | 自然资源部第二海洋研究所 | 一种基于高分辨率卫星影像的入海排口疑似排污识别方法 |
CN115909044B (zh) * | 2022-07-13 | 2023-07-07 | 中国科学院地理科学与资源研究所 | 一种国土空间结构时空演变模式挖掘方法 |
CN115144411B (zh) * | 2022-09-05 | 2022-12-06 | 国家卫星海洋应用中心 | 基于卫星散射计的海冰检测的方法、装置、设备和介质 |
CN115511390B (zh) * | 2022-11-14 | 2023-06-27 | 南方科技大学 | 一种沿海沿江脆弱性评估方法、系统、终端及存储介质 |
CN115661617B (zh) * | 2022-12-28 | 2023-03-14 | 成都中轨轨道设备有限公司 | 一种面向遥感大数据的动态自适应分布式协同工作方法 |
CN116258961B (zh) * | 2023-01-18 | 2023-12-01 | 广州市绿之城园林绿化工程有限公司 | 一种林业图斑变化快速识别方法及系统 |
CN116721346B (zh) * | 2023-06-14 | 2024-05-07 | 山东省煤田地质规划勘察研究院 | 一种基于深度学习算法的岸线智能识别方法 |
CN116521904B (zh) * | 2023-06-29 | 2023-09-22 | 湖南大学 | 一种基于5g边缘计算的船舶制造数据云融合方法及系统 |
CN116778303B (zh) * | 2023-08-25 | 2023-10-31 | 山东省国土测绘院 | 一种基于无人机遥感的植被覆盖度测量方法 |
CN117668958B (zh) * | 2023-12-06 | 2024-05-17 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种海岸线分形维数自动计算方法、系统及设备 |
CN117541626B (zh) * | 2024-01-08 | 2024-05-07 | 广东泰一高新技术发展有限公司 | 基于深度学习的遥感影像变化检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105740794A (zh) * | 2016-01-27 | 2016-07-06 | 中国人民解放军92859部队 | 一种基于卫星影像的海岸线自动提取与分类方法 |
WO2018028191A1 (zh) * | 2016-08-10 | 2018-02-15 | 福州大学 | 一种基于波段比模型和太阳高度角的tavi计算方法 |
CN109190538A (zh) * | 2018-08-24 | 2019-01-11 | 华北水利水电大学 | 一种基于遥感技术的多泥沙河流三角洲海岸带演化分析方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256419B (zh) * | 2017-12-05 | 2018-11-23 | 交通运输部规划研究院 | 一种利用多光谱解译提取港口码头图像的方法 |
CN110991393A (zh) * | 2019-12-17 | 2020-04-10 | 北京航天泰坦科技股份有限公司 | 一种海岸线变迁遥感监测与分析的方法及装置 |
CN111666918B (zh) * | 2020-06-22 | 2023-08-15 | 大连海洋大学 | 一种基于多因素的海岸线变化识别方法 |
-
2020
- 2020-06-22 CN CN202010572802.5A patent/CN111666918B/zh active Active
-
2021
- 2021-02-24 KR KR1020217020898A patent/KR20220000898A/ko not_active Application Discontinuation
- 2021-02-24 WO PCT/CN2021/077565 patent/WO2021258758A1/zh active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105740794A (zh) * | 2016-01-27 | 2016-07-06 | 中国人民解放军92859部队 | 一种基于卫星影像的海岸线自动提取与分类方法 |
WO2018028191A1 (zh) * | 2016-08-10 | 2018-02-15 | 福州大学 | 一种基于波段比模型和太阳高度角的tavi计算方法 |
CN109190538A (zh) * | 2018-08-24 | 2019-01-11 | 华北水利水电大学 | 一种基于遥感技术的多泥沙河流三角洲海岸带演化分析方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021258758A1 (zh) * | 2020-06-22 | 2021-12-30 | 大连海洋大学 | 一种基于多因素的海岸线变化识别方法 |
CN112926465A (zh) * | 2021-03-02 | 2021-06-08 | 中国人民解放军战略支援部队信息工程大学 | 基于点云类型的海岸线性质识别方法及装置 |
CN113673582A (zh) * | 2021-07-30 | 2021-11-19 | 西南交通大学 | 基于系统聚类分析的铁路动态基准点多层级分群方法 |
CN113673582B (zh) * | 2021-07-30 | 2023-05-09 | 西南交通大学 | 基于系统聚类分析的铁路动态基准点多层级分群方法 |
CN115439748A (zh) * | 2022-09-13 | 2022-12-06 | 中山大学 | 一种海岸线侵蚀程度的监测方法、装置以及电子设备 |
CN115439748B (zh) * | 2022-09-13 | 2023-09-26 | 中山大学 | 一种海岸线侵蚀程度的监测方法、装置以及电子设备 |
CN116124181A (zh) * | 2023-04-14 | 2023-05-16 | 国家海洋技术中心 | 一种潮汐观测设备的现场校准方法及系统 |
CN116797798A (zh) * | 2023-07-03 | 2023-09-22 | 生态环境部卫星环境应用中心 | 基于尺度变换的海岸线快速提取方法和装置 |
CN116797798B (zh) * | 2023-07-03 | 2024-01-09 | 生态环境部卫星环境应用中心 | 基于尺度变换的海岸线快速提取方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
KR20220000898A (ko) | 2022-01-04 |
CN111666918B (zh) | 2023-08-15 |
WO2021258758A1 (zh) | 2021-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111666918B (zh) | 一种基于多因素的海岸线变化识别方法 | |
CN110287932B (zh) | 基于深度学习图像语义分割的道路阻断信息提取方法 | |
Mesev | The use of census data in urban image classification | |
Davis et al. | The method for object-based diagnostic evaluation (MODE) applied to numerical forecasts from the 2005 NSSL/SPC Spring Program | |
CN113642849B (zh) | 考虑空间分布特征的地质灾害危险性综合评价方法及装置 | |
CN111028255A (zh) | 基于先验信息与深度学习的农田区域预筛选方法及装置 | |
CN111985389A (zh) | 一种基于流域属性距离的流域相似判别方法 | |
CN115965812B (zh) | 无人机影像对湿地植被物种和地物分类的评估方法 | |
Pawłuszek et al. | Landslides identification using airborne laser scanning data derived topographic terrain attributes and support vector machine classification | |
CN116310853A (zh) | 一种基于多源数据的中小城市边缘区提取方法 | |
CN115019163A (zh) | 基于多源大数据的城市要素识别方法 | |
Aahlaad et al. | An object-based image analysis of worldview-3 image for urban flood vulnerability assessment and dissemination through ESRI story maps | |
Zhou et al. | Automatic responsive-generation of 3D urban morphology coupled with local climate zones using generative adversarial network | |
Ganjirad et al. | Google Earth Engine-based mapping of land use and land cover for weather forecast models using Landsat 8 imagery | |
CN117171533A (zh) | 一种地理测绘作业数据实时采集处理方法及系统 | |
CN116205522B (zh) | 一种多维cnn耦合的滑坡易发性评价方法及系统 | |
Uehara et al. | Time-series metrics applied to land use and land cover mapping with focus on landslide detection | |
Kong et al. | A graph-based neural network approach to integrate multi-source data for urban building function classification | |
Kato | Integrating openstreetmap data in object based landcover and landuse classification for disaster recovery | |
Owusu et al. | Towards a scalable and transferable approach to map deprived areas using Sentinel-2 images and machine learning | |
Awuah | Effects of spatial resolution, land-cover heterogeneity and different classification methods on accuracy of land-cover mapping | |
Chen et al. | Evaluation of landslide potential due to land use in the slope | |
Deng et al. | Automated and refined wetland mapping of Dongting Lake using migrated training samples based on temporally dense Sentinel 1/2 imagery | |
Qi et al. | High-Resolution Land-Cover Mapping Based on a Cross-Resolution Deep Learning Framework and Available Low-Resolution Labels | |
Bello et al. | Large scale mapping: an empirical comparison of pixel-based and object-based classifications of remotely sensed data |
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 |