为了防控疫情，继续宅在家里～～。

国内某些在国外受过高等教育的人竟然认为公开承认上图红色注释会引起某些麻烦。

在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，这个方法的具体参数参数有三个分别为：

第四步、调用上面定义的方法绘制影像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(location=center, zoom_start=8, height=500)
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())

代码解析：

运行结果：

第五步、总结第四步调用方法可以将其写为一个具体的方法，然后直接调用这个方法。

具体代码：

# 使用Folium展示地图
def Mapdisplay(center, layers, Tiles="OpensTreetMap",zoom_start=10):  mapViz = folium.Map(location=center,tiles=Tiles, zoom_start=zoom_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_start=8)

代码解析：

这个程序运行结果和前面例子运行结果一样，区别就是在这里使用的我们自定义的方法Mapdisplay()方法，这个方法的参数分别为：

[

{

"layer": xxx, //图层的数据可以是image或者featureCollection

"vis": xxx, //显示规则，image配置和在JavaScript中配置一致

//featureCoolection配置用{"palette": 颜色}就可以

"name": xxx, //显示的图层名称

}

]

大家如果有问题需要交流或者有项目需要合作，可以微信联系我，加微信好友请留言加上“GEE”。

知乎专栏：https://zhuanlan.zhihu.com/c_123993183

CSDN：https://blog.csdn.net/shi_weihappy