文章预览
前面给大家分享过使用rgee获取IGBP土地覆被数据的代码 rgee土地覆被变化检测——2001~2019中国土地覆被变化 ,但是由于GEE API更新,rgee无法正常使用了,推荐使用吴秋生老师开发的Python版本GEE,非常方便,而且Python是GEE官方提供的接口语言,不会像R那样由于更新无法使用,功能相比JS版上限更高,推荐大家学习。 将上述R中的rgee代码转换为Python中的geemap代码,可以使用以下步骤和示例代码: import ee import geemap import os # Initialize Earth Engine ee.Initialize() # Define a region of interest with a shapefile using geemap ee_roi = geemap.shp_to_ee( './SHP/ChinaSizhi.SHP' ) # Call MODIS land cover dataset (MCD12Q1) cmr_lucc = ee.ImageCollection( "MODIS/006/MCD12Q1" ) \ .select( 'LC_Type1' ) \ .filterBounds(ee_roi) # Mosaic the image collection and clip to the region of interest cmr_lucc1 = cmr_lucc.mosai
………………………………