Google Earth Engine(GEE)——用10米分辨率土地分类数据进行土地面积计算出现错误

错误:

FeatureCollection (Error)

Error in map(ID=0): Feature, argument 'geometry': Invalid type. Expected type: Geometry. Actual type: Feature.

 

原始代码:

var boten = ee.FeatureCollection("projects/ee-suraj7-28/assets/Boten");
var boten_shp = boten.geometry()
var Lunagprabang = ee.FeatureCollection("projects/ee-suraj7-28/assets/Luangprabang");
var Lunagprabang_shp = Lunagprabang.geometry()
var Vientiana = ee.FeatureCollection("projects/ee-suraj7-28/assets/Vientiana");
var Vientiana_shp = Vientiana.geometry()

// Buffer//
var boten_buffer=boten_shp.buffer(20000)
var Lunagprabang_buffer=Lunagprabang_shp.buffer(20000)
var Vientiana_buffer=Vientiana_shp.buffer(20000)


var features = ee.FeatureCollection([
  ee.Feature(boten_buffer, {label: 'Boten'}),
  ee.Feature(Lunagprabang_buffer, {label: 'Lunagprabang'}),
  ee.Feature(Vientiana_buffer, {label: 'vientiana'}),
])
print(features)

var years = ee.List.sequence(2021,2022)


var calculateClassAreaByYear = features.map(function(feature){
  var dw2=years.map(function(year){
  var dw1= ee.ImageCollection('GOOGLE/DYNAMICWORLD/V1')
              .filterDate(ee.Date.fromYMD(year, 1, 1), ee.Date.fromYMD(year, 12, 1)).filter(ee.Filter.bounds(feature)).select('label').reduce(ee.Reducer.mode()).eq(6);
              return dw1
              })
              return ee.FeatureCollection(dw2)
    })

print(calculateClassAreaByYear.flatten())

我们这里可以看到提示期待的类型不对,需求的类型是几何图形,但是这里给的确实矢量,所以我们需要再后面加一个.geometry()

修改后的代码:


var years = ee.List.sequence(2021,2022)


var calculateClassAreaByYear = features.map(function(feature){
  var dw2=years.map(function(year){
  var dw1= ee.ImageCollection('GOOGLE/DYNAMICWORLD/V1')
              .filterDate(ee.Date.fromYMD(year, 1, 1), ee.Date.fromYMD(year, 12, 31)).filter(ee.Filter.bounds(feature.geometry())).select('label').reduce(ee.Reducer.mode()).eq(6);
  var areaImage = dw1.multiply(ee.Image.pixelArea())    
  var area = areaImage.reduceRegion({
  reducer: ee.Reducer.sum(),
  geometry: feature.geometry(),
  scale: 10,
  maxPixels: 1e10
  })

 这里有一个矢量面积的计算函数:

ee.Image.pixelArea()
Generate an image in which the value of each pixel is the area of that pixel in square meters. The returned image has a single band called "area."

ee.Image.pixelArea()
生成一个图像,其中每个像素的值是该像素的面积,单位是平方米。返回的图像有一个称为 "面积 "的单带。

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
使用GEE进行计算森林面积可以按照以下步骤进行: 1. 在GEE上打开代码编辑器,选择JavaScript语言。 2. 在代码导入Globeland30数据: ```javascript var globeland30 = ee.Image('users/your_username/Globeland30'); ``` 其,`your_username`需要替换为你自己的用户名,`Globeland30`是你上传的土地利用栅格数据的名称。 3. 定义森林的分类值。在Globeland30数据,森林的分类值通常是11、14和20。你可以使用以下代码定义森林的分类值: ```javascript var forest = globeland30.eq(11).or(globeland30.eq(14)).or(globeland30.eq(20)); ``` 4. 计算森林面积。你可以使用以下代码计算森林面积: ```javascript var forest_area = forest.multiply(ee.Image.pixelArea()).divide(10000).reduceRegion({ reducer: ee.Reducer.sum(), geometry: your_geometry, scale: 30, maxPixels: 1e12 }); ``` 其,`your_geometry`是你感兴趣的区域的几何形状。这个区域可以是一个点、线或多边形,也可以是一个从图像提取的特定区域。在这个代码,我们使用了`reduceRegion`方法,该方法将森林图像的像素面积乘以每个像素的面积,并将其除以10000,以计算森林面积(单位为公顷)。 5. 输出森林面积。你可以使用以下代码输出森林面积: ```javascript print('Forest area:', forest_area.get('classification').toFixed(2), 'ha'); ``` 这将输出森林面积的值,保留两位小数,并附带单位(公顷)。 以上是在GEE上使用自己上传的土地利用栅格数据(Globeland30)进行计算森林面积的大致步骤,你可以根据自己的需要进行调整和修改。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值