在这个世界上,无论任何事物,都需要我们发现规律,寻找规律,最后利用规律。我们之前与学员分享了很多关于GEE和ChatGPT的规律,包括GEE的组织方式以及算法结构、ChatGPT的万能公式等。今天,我们利用这个方式,快速实现全国不同分辨率土地利用分类的结果。
01
ChatGPT的万能公式
赋予身份 + 内容设定 + 科学扩展 + 润色表达 + 科学排版
Prompt:
你是一名Google earth engine 的算法工程师,请根据GEE自带的土地利用分类产品,裁剪出10m(ESA)、30m(ee.ImageCollection("users/xxc/GLC_2000_2015"))、100m(Copernicus Global Land Service (CGLS))、500m(MOD12Q1)等不同分辨率的中国范围的产品数据,并显示其空间分布。
注意:这里需要把台湾等地区加进去,所以我们采用审图号的研究区范围,得到以下代码:
// 导入中国行政边界
var china = ee.FeatureCollection("projects/ee-tilmacatanla/assets/boundry/china_provinces").union();
// 1. 10米分辨率 - ESA WorldCover
var esaworldcover = ee.ImageCollection("ESA/WorldCover/v100").first();
print('esaworldcover',esaworldcover)
var esaworldcover_china = esaworldcover.clip(china);
// 定义可视化参数
var esaworldcoverVis = {
min: 10,
max: 100,
palette: [
'006400', // 树木
'ffbb22', // 灌木
'ffff4c', // 草地
'f096ff', // 农田
'fa0000', // 建筑
'b4b4b4', // 裸地
'f0f0f0', // 积雪和冰川
'0064c8', // 永久水体
'0096a0', // 湿地
'00cf75', // 红树林
'fae6a0' // 苔原
],
};
// 显示 ESA WorldCover
Map.addLayer(esaworldcover_china, esaworldcoverVis, 'ESA WorldCover 10m');
// 2. 30米分辨率 - GLC_2000_2015
var glc_2000_2015 = ee.ImageCollection("users/xxc/GLC_2000_2015").mosaic();
var glc_2000_2015_china = glc_2000_2015.clip(china);
Map.addLayer(glc_2000_2015_china, {}, 'GLC_2000_2015 30m');
// 3. 100米分辨率 - Copernicus Global Land Service
var copernicus100m = ee.ImageCollection('COPERNICUS/Landcover/100m/Proba-V/Global')
.filterDate('2015-01-01', '2015-12-31')
.first()
.select('discrete_classification');
var copernicus100m_china = copernicus100m.clip(china);
// 定义可视化参数
var copernicusVis = {
min: 0,
max: 200,
palette: [
'aec3d4', '152106', '225129', '369b47', '30eb5b',
'387242', '6a2325', 'c3aa69', 'b76031', 'd9903d',
'91af40', '111149', 'cdb33b', 'cc0013', '33280d',
'd7cdcc', 'f7e174', '743411'
]
};
// 显示 Copernicus 100m
Map.addLayer(copernicus100m_china, copernicusVis, 'Copernicus 100m');
// 4. 500米分辨率 - MODIS MCD12Q1
var modis500m = ee.ImageCollection("MODIS/006/MCD12Q1")
.filterDate('2015-01-01', '2015-12-31')
.first()
.select('LC_Type1');
var modis500m_china = modis500m.clip(china);
// 定义可视化参数
var modisVis = {
min: 1,
max: 17,
palette: [
'05450a', '086a10', '54a708', '78d203', '009900',
'c6b044', 'dcd159', 'dade48', 'fbff13', 'b6ff05',
'27ff87', 'c24f44', 'a5a5a5', 'ff6d4c', '69fff8',
'f9ffa4', '1c0dff'
],
};
// 显示 MODIS MCD12Q1 500m
Map.addLayer(modis500m_china, modisVis, 'MODIS MCD12Q1 500m');
// 调整地图视图到中国范围
Map.centerObject(china, 4);
通过调色,显示等,最终可以得到以下结果: