网站举报电话是多少,wordpress 自定义边栏,有源码后怎么做网站,营销战略咨询在Google Earth Engine的python版API更新后#xff0c;之前使用folium动态加载地图的代码就不能在正常运行#xff0c;因为整个Google Earth Engine的地图加载服务的URL发生了更新#xff0c;所以我们也需要更新相关绘制方法。下面我会讲解一种新的绘制方法#xff0c;大家…在Google Earth Engine的python版API更新后之前使用folium动态加载地图的代码就不能在正常运行因为整个Google Earth Engine的地图加载服务的URL发生了更新所以我们也需要更新相关绘制方法。下面我会讲解一种新的绘制方法大家可以根据自己的需求来扩充相关绘制方法直接上相关具体代码。 需要说明一下我所有的代码都是在Jupyter Lab中运行这些代码无法在pycharm中直接运行。
第一步、注册GEE。 第二步、导入外部库folium。 第三步、编写添加GEE图层的方法并且将其作为自定义方法加入到folium的Map中。 具体代码 # 定义加载GEE影像数或者矢量数据方法 def add_ee_layer(self, ee_object, vis_params, name): try: if isinstance(ee_object, ee.Image): map_id_dict ee.Image(ee_object).getMapId(vis_params) folium.raster_layers.TileLayer( tiles map_id_dict[tile_fetcher].url_format, attr Google Earth Engine, name name, overlay True, control True ).add_to(self) elif isinstance(ee_object, ee.FeatureCollection): ee_object_new ee.Image().paint(ee_object, 0, 2) map_id_dict ee.Image(ee_object_new).getMapId(vis_params) folium.raster_layers.TileLayer( tiles map_id_dict[tile_fetcher].url_format, attr Google Earth Engine, name name, overlay True, control True ).add_to(self) except: print(Could not display {}.format(name)) folium.Map.add_ee_layer add_ee_layer
代码解析
加载图层方法为add_ee_layer这个方法的具体参数参数有三个分别为 ee_object显示对象可以为image或者featureCollection vis_params显示规则其中image和JavaScript版配置一致featureCollection配置使用palette配置颜色即可 name图层的名称
第四步、调用上面定义的方法绘制影像Image和矢量数据FeatureCollection。 具体代码 dem ee.Image(USGS/SRTMGL1_003) vis { min: 0, max: 4000, palette: [006633, E5FFCC, 662A00, D8D8D8, F5F5F5] } roi ee.FeatureCollection(users/wangweihappy0/training03/bj_shp) center roi.geometry().centroid().getInfo()[coordinates] center.reverse() print(center is: {}.format(center)) my_map folium.Map(locationcenter, zoom_start8, height500) my_map.add_ee_layer(dem.updateMask(dem.gt(0)), vis, DEM) my_map.add_ee_layer(roi, {palette: FF0000}, roi) my_map.add_child(folium.LayerControl())
代码解析 vis配置的是DEM显示规则 center是地图要居中显示的位置 my_map是通过folium.Map配置显示的地图 my_map通过add_ee_layer添加定义的图层dem和roi
运行结果 第五步、总结第四步调用方法可以将其写为一个具体的方法然后直接调用这个方法。 具体代码 # 使用Folium展示地图 def Mapdisplay(center, layers, TilesOpensTreetMap,zoom_start10): mapViz folium.Map(locationcenter,tilesTiles, zoom_startzoom_start) for data in layers: mapViz.add_ee_layer(data.get(layer), data.get(vis), data.get(name)) mapViz.add_child(folium.LayerControl()) return mapViz
使用这个方法示例如下 具体代码 dem ee.Image(USGS/SRTMGL1_003) vis { min: 0, max: 4000, palette: [006633, E5FFCC, 662A00, D8D8D8, F5F5F5] } roi ee.FeatureCollection(users/wangweihappy0/training03/bj_shp) center roi.geometry().centroid().getInfo()[coordinates] center.reverse() print(center is: {}.format(center)) layers [ { layer: dem, vis: vis, name: DEM }, { layer: roi, vis: {palette: FF0000}, name: roi } ] Mapdisplay(center, layers, zoom_start8)
代码解析 这个程序运行结果和前面例子运行结果一样区别就是在这里使用的我们自定义的方法Mapdisplay()方法这个方法的参数分别为 center中心点坐标 layers要加载的所有图层信息列表配置规则是 [ { layer: xxx, //图层的数据可以是image或者featureCollection vis: xxx, //显示规则image配置和在JavaScript中配置一致 //featureCoolection配置用{palette: 颜色}就可以 name: xxx, //显示的图层名称 } ] Tiles这个使用默认就可以不用修改 zoom_start最初缩放比例比如这里使用的是8 来源请引用地理遥感生态网科学数据注册与出版系统。
