Open Geospatial Engine教程:坡度、坡向山阴的计算和可视化和导出

张开发
2026/6/5 6:43:39 15 分钟阅读
Open Geospatial Engine教程:坡度、坡向山阴的计算和可视化和导出
OGE的界面不同于GEE这里的控制台被放在了左下角代码的以python为主没有JavaScript API因此我们熟悉python的结果就非常容易上手不用再学习JavaScript来二次学习。下满是OGE在线开发的主要界面算子代码# 导入开放地球引擎Open Geospatial Engine的Python SDKimportoge# 初始化OGE环境建立与后端服务的连接oge.initialize()# 创建服务对象用于调用OGE的数据获取、处理等功能serviceoge.Service()# 从服务中获取指定覆盖数据Coverage# coverageID: 数据集的唯一标识此处应为ASTGTM_N28E056注意原代码中的空格是笔误正确写法无空格# productID: 产品类型此处为ASTER GDEM数字高程模型分辨率30米demservice.getCoverage(coverageIDASTGTM_N28E056,productIDASTER_GDEM_DEM30)# 调用坡向计算处理过程Coverage.terrAspect# .execute(dem, 1) 执行计算参数1可能表示坡度计算时的邻域大小或平滑因子具体取决于算子定义aspectservice.getProcess(Coverage.terrAspect).execute(dem,1)# 定义坡向图层的可视化参数# min/max: 坡向值范围坡向通常为0°~360°但此处设置为-1~1可能是归一化后的范围或特定算法输出# palette: 颜色渐变列表从深灰到浅灰用于渲染不同坡向vis_params{min:-1,max:1,palette:[#808080,#949494,#a9a9a9,#bdbebd,#d3d3d3,#e9e9e9]}# 将可视化参数应用到坡向数据生成地图瓦片并命名为terrAspect# .getMap() 触发地图渲染返回可用于显示或保存的地图对象aspect.styles(vis_params).getMap(terrAspect)# 使用内置地图客户端将视图中心定位到指定经纬度及缩放级别# 参数经度56.25°纬度28.40°缩放级别11级别越大显示越精细oge.mapclient.centerMap(56.25,28.40,11)可视化结果控制台展示的结果拓展下面是我执行slope代码的结果这里建议自己先自己写代码然后再看所有代码另外建议尝试使用hillshadow来计算一下和可视化。全部代码importoge oge.initialize()serviceoge.Service()demservice.getCoverage(coverageIDASTGTM_N28E056,productIDASTER_GDEM_DEM30)aspectservice.getProcess(Coverage.terrAspect).execute(dem,1)slopeservice.getProcess(Coverage.terrSlope).execute(dem,1)vis_params{min:-1,max:1,palette:[#808080,#949494,#a9a9a9,#bdbebd,#d3d3d3,#e9e9e9]}aspect.styles(vis_params).getMap(terrAspect)slope.styles(vis_params).getMap(terrSlope)oge.mapclient.centerMap(56.25,28.40,11)导出处理后的结果很多时候我们初学者来说如果你一直傻傻等待那么就错了你需要点击任务面板中查看任务点击运行就可以到出你要导出的影像importoge oge.initialize()serviceoge.Service()demservice.getCoverage(coverageIDASTGTM_N28E056,productIDASTER_GDEM_DEM30)aspectservice.getProcess(Coverage.terrAspect).execute(dem,1)vis_params{min:-1,max:1,palette:[#808080,#949494,#a9a9a9,#bdbebd,#d3d3d3,#e9e9e9]}aspect.styles(vis_params).export(terrAspect)oge.mapclient.centerMap(56.25,28.40,11)导出的结果可以看到任务名称坐标系和分辨率的设定。因此这里我们可以根据自己的需求来设定。在这里点击下载就可以将tif下载到本地。具体含义差异.getMap(“terrAspect”)表示获取地图瓦片或地图对象通常用于在 OGE 内置地图客户端中直接可视化渲染坡向图层。调用后地图会显示在界面上配合 oge.mapclient.centerMap() 可实现交互式浏览。这里批处理的就是用export.export(“terrAspect”)表示导出数据例如导出为 GeoTIFF、NetCDF 等格式将计算结果保存到本地或云端存储不主动显示地图。后续可能需要手动加载导出的文件才能查看。地形的算法

更多文章