34_geemap学习笔记 | 像素值的交互式提取

34 extract values
主要函数：extract_values_to_points()

import os
import ee 
import geemap
geemap.set_proxy(10809)

Map = geemap.Map()
Map

#添加数据至地图
landsat7 = ee.Image('LE7_TOA_5YEAR/1999_2003') \
                      .select(0, 1, 2, 3, 4, 6)
landsat_vis = {
    'bands': ['B4', 'B3', 'B2'],
    'gamma': 1.4
}
Map.addLayer(landsat7, landsat_vis, "LE7_TOA_5YEAR/1999_2003")
Map.set_plot_options(add_marker_cluster=True)

#激活绘图工具
#点击Plotting功能，在地图上选择感兴趣区，默认绘制的是点

#导出像素值至shp或csv文件
out_dir = os.path.join(os.path.expanduser('~'), 'Downloads')
out_shp = os.path.join(out_dir, 'points.shp')
out_csv = os.path.join(out_dir, 'points.csv')

Map.extract_values_to_points(out_shp)

The shapefile has been saved to: C:\Users\Administrator\Downloads\points.shp

Map.extract_values_to_points(out_csv)

Interactive Region Reduction

m =geemap.Map()

#add data to the map
Collection = ee.ImageCollection('MODIS/006/MOD13A2') \
                        .filterDate('2015-01-01', '2019-12-31') \
                        .select('NDVI')
#Convert the image collection to an image
image = Collection.toBands()

ndvi_vis = {
    'min': 0.0,
    'max': 9000.0,
    'palette': [
        'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
        '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
        '012E01', '011D01', '011301'
    ]
}

m.addLayer(image, {}, 'MODIS NDVI Time-series')
m.addLayer(image.select(0), ndvi_vis, 'MODIS NDVI VIS')

m

#set reducer
m.set_plot_options(add_marker_cluster=True, maker=None)
m.roi_reducer = ee.Reducer.mean()

#Export data
out_dir = os.path.join(os.path.expanduser('~'), 'Downloads')
out_shp = os.path.join(out_dir, 'ndvi.shp')
m.extract_values_to_points(out_shp)

The shapefile has been saved to: C:\Users\Administrator\Downloads\ndvi.shp