错误32:使用getMapId()时管道损坏

时间:2019-11-01 12:01:13

标签: python geojson google-earth-engine

我正在读取一个大的geojson文件,其中包含数千个建筑物的足迹。该文件如下所示:

{
  "type":"FeatureCollection",
  "features":
  [
    {"type":"Feature","geometry":{"type":"Polygon","coordinates":[[[-97.208203,50.90649],[-97.208203,50.906568],[-97.208367,50.906568],[-97.208367,50.90649],[-97.208203,50.90649]]]},"properties":{}},
    {"type":"Feature","geometry":{"type":"Polygon","coordinates":[[[-97.138545,50.91915],[-97.1387,50.919146],[-97.138692,50.919018],[-97.138537,50.919022],[-97.138545,50.91915]]]},"properties":{}},
    {"type":"Feature","geometry":{"type":"Polygon","coordinates":[[[-97.106312,50.949822],[-97.106178,50.949845],[-97.106255,50.950025],[-97.106389,50.950002],[-97.106312,50.949822]]]},"properties":{}},
.
.
.
  ]
}

我正在读取该文件并将其加载到ee.FeatureCollection()中,

import json

# Open file to parse
with open('/content/drive/My Drive/Earthengine/buildings.geojson') as f:
  data = json.load(f)

# Parse each of the polygons to features
polygons = []
for feat in data['features']:
  polygons.append(ee.Feature(feat))

# Create the feature collection
buildings = ee.FeatureCollection(polygons)

但是当我尝试通过以下方式将这些建筑物显示为大叶草时:

polyImage = ee.Image(0).byte().paint(buildings, 1)
polyImage = polyImage.updateMask(polyImage)

mapid = polyImage.getMapId()
map = folium.Map(location=[38., -100.], zoom_start=5)
folium.TileLayer(
    tiles=EE_TILES.format(**mapid),
    attr='Google Earth Engine',
    overlay=True,
    name='training polygons',
  ).add_to(map)
map.add_child(folium.LayerControl())
map

我收到以下错误:

---------------------------------------------------------------------------
BrokenPipeError                           Traceback (most recent call last)
<ipython-input-7-3c769c73ff13> in <module>()
      2 polyImage = polyImage.updateMask(polyImage)
      3 
----> 4 mapid = polyImage.getMapId()
      5 map = folium.Map(location=[38., -100.], zoom_start=5)
      6 folium.TileLayer(

15 frames
/usr/lib/python3.6/ssl.py in write(self, data)
    640         The 'data' argument must support the buffer interface.
    641         """
--> 642         return self._sslobj.write(data)
    643 
    644     def getpeercert(self, binary_form=False):

BrokenPipeError: [Errno 32] Broken pipe

是否可能因为文件太大而发生?可以解决这个问题吗?

1 个答案:

答案 0 :(得分:0)

我很确定这是因为Earth Engine仍在计算,并且在未收到响应后在客户端关闭了连接。有关断管的更多信息,请参见this answer

假设您有很多多边形,Earth Engine必须构造多边形,然后对其进行栅格化,这可能不会花费很多计算,但是可能会花费一些时间,因此会折断管道。解决此问题的一种方法是将建筑物多边形另存为资产,然后将资产直接加载到地图上。

您可以使用以下语法导出featureCollection:

task = ee.batch.Export.table.toAsset(
  collection=buildings,
  description='my-buildings', 
  assetId='users/myusername/buildings'
)
task.start()