download_by_canvas.py 777 B

1234567891011121314151617181920212223242526272829
  1. import json
  2. import ee
  3. from ee_plugin import Map
  4. from ee_plugin.contrib import palettes
  5. dem = ee.Image('JAXA/ALOS/AW3D30/V2_2').select('AVE_DSM')
  6. Map.addLayer(dem, { 'min': 0, 'max': 3000 }, 'DEM', True)
  7. # MANUAL STEP: use "Create layer from extent" tool and activate the resulting layer
  8. # get first feature geometry from active layer
  9. layer = iface.activeLayer()
  10. feature = next(layer.getFeatures())
  11. geom = feature.geometry()
  12. geom_json = json.loads(geom.asJson())
  13. # show geometry (double-check)
  14. geom_ee = ee.Geometry.Polygon(geom_json['coordinates'], 'EPSG:3857', False)
  15. Map.addLayer(geom_ee, {}, 'geom')
  16. # download dem using given geometry as region
  17. url = dem.getDownloadURL({
  18. 'name': 'dem',
  19. 'scale': 30,
  20. 'region': json.dumps(geom_ee.getInfo())
  21. })
  22. print(url)