ChatGPT-o1+GEE|快速实现全国不同分辨率土地利用分类的结果

GEE遥感训练营 · 公众号 · · 2024-09-13 23:49

正文

请到「今天看啥」查看全文

在这个世界上，无论任何事物，都需要我们发现规律，寻找规律，最后利用规律。我们之前与学员分享了很多关于GEE和ChatGPT的规律，包括 GEE的组织方式以及算法结构、 ChatGPT的万能公式等。今天，我们利用这个方式，快速实现全国不同分辨率土地利用分类的结果。

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 WorldCovervar 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 WorldCoverMap.addLayer(esaworldcover_china, esaworldcoverVis, 'ESA WorldCover 10m');
// 2. 30米分辨率 - GLC_2000_2015var 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 Servicevar 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 100mMap.addLayer(copernicus100m_china, copernicusVis, 'Copernicus 100m');
// 4. 500米分辨率 - MODIS MCD12Q1var 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 500mMap.addLayer(modis500m_china, modisVis, 'MODIS MCD12Q1 500m');
// 调整地图视图到中国范围Map.centerObject(china, 4);

通过调色，显示等，最终可以得到以下结果：