所以我使用了 Google Earth Engine 并研究了他们存储库中的一些示例代码。我正在使用 Python 3.6。看起来 Google 将不再通过他们的 ee.mapclient() 支持 Python 3 中的映射功能。我想知道是否有人找到合适的解决方法？让我概述一下问题。

我尝试加载 ee.mapclient 来绘制 map 。

import ee
import ee.mapclient
ee.Initialize()

但是我得到一个错误:

ModuleNotFoundError                       Traceback (most recent call last)
<ipython-input-13-6d4860410653> in <module>()
      1 import ee
----> 2 import ee.mapclient
      3 ee.Initialize()

/media/krishnab/lakshmi/anaconda3/envs/pMining/lib/python3.6/site-packages/ee/mapclient.py in <module>()
     29 
     30 import collections
---> 31 import cStringIO
     32 import functools
     33 import math

ModuleNotFoundError: No module named 'cStringIO'

cStringIO 问题很容易解决，如下所示: python 3.4.0 email package install: ImportError: No module named 'cStringIO'

所以我去 Google Earth Engine Repo 上发布了一个问题，但发现了一个预先存在的问题:

https://github.com/google/earthengine-api/issues/16

在这个问题中，开发人员承认了这个问题，但表示由于底层 Tk 包的限制，他们不会修复它。

这里引用了这个问题:

We have not been actively maintaining the mapclient object, because it depends on Tk, a graphical user interface toolkit, which behaves differently on different machines. Could you describe your use case that requires mapclient? We may be able to suggest an alternative approach.

Google 开发人员提出要提交解决方法，但到目前为止还没有发布任何解决方法。

因此，我想知道是否还有其他人在 Python3.6 中找到了解决此问题的合适方法？

作为真正的代码示例，我可以提供以下来自 Google 示例存储库的代码:

import datetime
import ee
import ee.mapclient

ee.Initialize()
ee.mapclient.centerMap(-95.738, 18.453, 9)

# Filter the LE7 collection to a single date.
collection = (ee.ImageCollection('LE7_L1T')
              .filterDate(datetime.datetime(2002, 11, 8),
                          datetime.datetime(2002, 11, 9)))
image = collection.mosaic().select('B3', 'B2', 'B1')

# Display the image normally.
ee.mapclient.addToMap(image, {'gain': '1.6, 1.4, 1.1'}, 'Land')

# Add and stretch the water.  Once where the elevation is masked,
# and again where the elevation is zero.
elev = ee.Image('srtm90_v4')
mask1 = elev.mask().eq(0).And(image.mask())
mask2 = elev.eq(0).And(image.mask())
ee.mapclient.addToMap(
    image.mask(mask1), {'gain': 6.0, 'bias': -200}, 'Water: Masked')
ee.mapclient.addToMap(
    image.mask(mask2), {'gain': 6.0, 'bias': -200}, 'Water: Elev 0')

最佳答案

不确定在 Python 3 中导入 ee.mapclient 是否仍然存在问题，但这里可以使用 folium package 解决问题可以具有与 ee.mapclient 类似的方法。但是，这将在您的 Web 浏览器中呈现结果，而不是作为弹出窗口(除非使用 Jupyter 笔记本，否则您可以调用 folium map 对象进行呈现)。

import ee
import folium
import datetime
import webbrowser

ee.Initialize()

class eeMapHack(object):
    def __init__(self,center=[0, 0],zoom=3):
        self._map = folium.Map(location=center,zoom_start=zoom)
        return

    def addToMap(self,img,vizParams,name):
         map_id = ee.Image(img.visualize(**vizParams)).getMapId()
         tile_url_template = "https://earthengine.googleapis.com/map/{mapid}/{{z}}/{{x}}/{{y}}?token={token}"
         mapurl = tile_url_template.format(**map_id)
         folium.WmsTileLayer(mapurl,name=name).add_to(self._map)

         return

    def addLayerControl(self):
         self._map.add_child(folium.map.LayerControl())
         return


# initialize map object
eeMap = eeMapHack(center=[18.453,-95.738],zoom=9)

# Filter the LE7 collection to a single date.
collection = (ee.ImageCollection('LE7_L1T')
          .filterDate(datetime.datetime(2002, 11, 8),
                      datetime.datetime(2002, 11, 9)))
image = collection.mosaic().select('B3', 'B2', 'B1')
eeMap.addToMap(image, {'gain': '1.6, 1.4, 1.1'}, 'Land')

# Add and stretch the water.  Once where the elevation is masked,
# and again where the elevation is zero.
elev = ee.Image('srtm90_v4')
mask1 = elev.mask().eq(0).And(image.mask())
mask2 = elev.eq(0).And(image.mask())

eeMap.addToMap(image.mask(mask1), {'gain': 6.0, 'bias': -200}, 'Water: Masked')
eeMap.addToMap(image.mask(mask2), {'gain': 6.0, 'bias': -200}, 'Water: Elev 0')

# add layer control to map
eeMap.addLayerControl()

outHtml = 'map.html' # temporary file path, change if needed
eeMap._map.save(outHtml)

webbrowser.open('file://'+outHtml) 

您应该在浏览器中获得一个交互式 web map ，其中包含如下所示的分析结果:

关于python - Google Earth Engine Python API 的解决方法和 Python 3 中不支持 `ee.mapclient`，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/46960328/