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

举报
此星光明 发表于 2022/08/19 23:17:25 2022/08/19
【摘要】 错误: FeatureCollection (Error) Error in map(ID=0): Feature, argument 'geometry': Invalid type. Expected type: Geometry. Actual type: Feature.   原始代码: var bo...

错误:

FeatureCollection (Error)

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

 

原始代码:


  
  1. var boten = ee.FeatureCollection("projects/ee-suraj7-28/assets/Boten");
  2. var boten_shp = boten.geometry()
  3. var Lunagprabang = ee.FeatureCollection("projects/ee-suraj7-28/assets/Luangprabang");
  4. var Lunagprabang_shp = Lunagprabang.geometry()
  5. var Vientiana = ee.FeatureCollection("projects/ee-suraj7-28/assets/Vientiana");
  6. var Vientiana_shp = Vientiana.geometry()
  7. // Buffer//
  8. var boten_buffer=boten_shp.buffer(20000)
  9. var Lunagprabang_buffer=Lunagprabang_shp.buffer(20000)
  10. var Vientiana_buffer=Vientiana_shp.buffer(20000)
  11. var features = ee.FeatureCollection([
  12. ee.Feature(boten_buffer, {label: 'Boten'}),
  13. ee.Feature(Lunagprabang_buffer, {label: 'Lunagprabang'}),
  14. ee.Feature(Vientiana_buffer, {label: 'vientiana'}),
  15. ])
  16. print(features)
  17. var years = ee.List.sequence(2021,2022)
  18. var calculateClassAreaByYear = features.map(function(feature){
  19. var dw2=years.map(function(year){
  20. var dw1= ee.ImageCollection('GOOGLE/DYNAMICWORLD/V1')
  21. .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);
  22. return dw1
  23. })
  24. return ee.FeatureCollection(dw2)
  25. })
  26. print(calculateClassAreaByYear.flatten())

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

修改后的代码:


  
  1. var years = ee.List.sequence(2021,2022)
  2. var calculateClassAreaByYear = features.map(function(feature){
  3. var dw2=years.map(function(year){
  4. var dw1= ee.ImageCollection('GOOGLE/DYNAMICWORLD/V1')
  5. .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);
  6. var areaImage = dw1.multiply(ee.Image.pixelArea())
  7. var area = areaImage.reduceRegion({
  8. reducer: ee.Reducer.sum(),
  9. geometry: feature.geometry(),
  10. scale: 10,
  11. maxPixels: 1e10
  12. })

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

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()
生成一个图像,其中每个像素的值是该像素的面积,单位是平方米。返回的图像有一个称为 "面积 "的单带。

文章来源: blog.csdn.net,作者:此星光明,版权归原作者所有,如需转载,请联系作者。

原文链接:blog.csdn.net/qq_31988139/article/details/126406202

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。