【GEE笔记】在线分类流程,标注样本点、分类和精度评价

GEE在线分类流程

介绍

GEE(Google Earth Engine)是一个强大的地理信息处理平台,可以实现在线的遥感影像分析和处理。本文将介绍如何使用GEE进行在线的分类流程,包括标注样本点、分类和精度评价。本文以2020年5月至8月的哨兵2影像为例,对区域内的土地覆盖类型进行分类。

标注样本点

首先,加载原始影像,进行在线的标注。

加载的影像:

加载影像的代码如下:

javascript 复制代码
var geometry = ee.Geometry.Polygon(
        [[[121.81940156260009, 40.92383488850036],
          [121.81940156260009, 40.73887826797227],
          [121.99998933115478, 40.73887826797227],
          [121.99998933115478, 40.92383488850036]]], null, false)
var year=2020 //设置年份
var bandlist=['B2','B3','B4','B8','B11','B12'] //设置波段列表
var start = ee.Date(year+'-5-1'); //设置开始日期
var finish = ee.Date(year+'-8-1'); //设置结束日期

var dataset = ee.ImageCollection('COPERNICUS/S2_SR') //加载哨兵2影像集合
                  .filterDate(start, finish) //按日期过滤
                  .filterBounds(geometry) //按范围过滤
                  .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20)) //按云量过滤
                  
dataset=dataset.select(bandlist); //按波段选择            
var rgbVis = {
  min: 0.0, 
  max: 3000, 
  bands: ['B4', 'B3', 'B2'], //设置RGB波段
};


Map.centerObject(geometry) //地图中心定位到范围
print(dataset) //打印影像集合信息

var image=dataset.median().clip(geometry) //计算影像集合的中值,并裁剪到范围
Map.addLayer(image,rgbVis) //添加影像图层

标注样本点过程如下:

1、创建样本点集合对象(集合名称可以设置为类别名称便于区分)

2、设置类别字段及其属性(属性为类别数字代码)

分类

接下来,进行分类。选择随机森林算法作为分类器(也可选择GEE官方其他分类器)。将样本点集合分为训练集和测试集,用训练集来训练分类器,用测试集来评估分类器的性能。

分类主要包括合并样本点集合、划分训练集和测试集、训练分类器、应用分类器、添加分类图层。分类的结果如下图所示:

分类的代码如下:

javascript 复制代码
var sample=danshui.merge(haishui).merge(jianpeng).merge(nongtian).merge(luwie).merge(jianzhu)
.merge(tantu).merge(caodi) //合并样本点集合

var withRandom = sample.randomColumn('random'); //给样本点集合添加随机数列
var split = 0.7;  //七成训练 三成测试
var trainingPartition2 = withRandom.filter(ee.Filter.lt('random', split)); //按随机数列划分训练集
var testingPartition = withRandom.filter(ee.Filter.gte('random', split)); //按随机数列划分测试集

var trainingPartition=image.sampleRegions({ //从影像中提取训练集的像素值
    collection: trainingPartition2, //输入训练集
    scale:10, //设置空间分辨率
    properties: ['Map'], //设置类别字段
    })
var classifier = ee.Classifier.smileRandomForest(100).train({ //训练随机森林分类器
  features: trainingPartition, //输入训练集
  classProperty:'Map', //设置类别字段
  // inputProperties :['B4', 'B3', 'B2']
});
var classied=image.classify(classifier,"smileRandomForest"); //应用分类器

Map.addLayer(classied.randomVisualizer()) //添加分类图层

精度评价

最后需要评估分类器的性能,计算分类精度。

可以通过以下代码来评估分类精度:

javascript 复制代码
var pixelValues2=classied.sampleRegions({ //从影像中提取测试集的像素值
    collection: testingPartition, //输入测试集
    scale:10, 
});
print("pixelValues2",pixelValues2); //打印测试集信息
var confusionMatrix2 = pixelValues2.errorMatrix("Map", 'smileRandomForest');//,[10,20,30,40,50,60,80,90] //计算混淆矩阵
print('confusionMatrix', confusionMatrix2); //打印混淆矩阵
print('accuracy()', confusionMatrix2.accuracy()); //打印总体精度
// print('consumersAccuracy', confusionMatrix2.consumersAccuracy());
print('kappa()', confusionMatrix2.kappa()); //打印Kappa系数
//print('order()', confusionMatrix.order());
// print('producersAccuracy()', confusionMatrix2.producersAccuracy());
相关推荐
MUTA️1 分钟前
《MAE: Masked Autoencoders Are Scalable Vision Learners》论文精读笔记
人工智能·笔记·深度学习·transformer
使一颗心免于哀伤42 分钟前
《设计模式之禅》笔记摘录 - 6.原型模式
笔记·设计模式
刘孬孬沉迷学习1 小时前
5G标准学习笔记15 --CSI-RS测量
网络·笔记·学习·5g·信息与通信·信号处理
给自己记录用2 小时前
qt笔记(1)——Qtablewidget使用
笔记·qt
飞升不如收破烂~2 小时前
英语学习笔记2.0
笔记·学习
金色光环2 小时前
概率论:理解区间估计【超详细笔记】
笔记·数学·概率论·数理统计·区间估计
大佐不会说日语~2 小时前
JVM垃圾回收机制面试笔记
jvm·笔记·面试
超龄超能程序猿3 小时前
使用 Python 对本地图片进行图像分类
开发语言·人工智能·python·机器学习·分类·数据挖掘·scipy
梁辰兴3 小时前
企业培训笔记:Vue3前端框架配置
笔记·前端框架·npm·vue·vue3·node
陈洪奇10 小时前
注册中心学习笔记整理
笔记·学习