我正在读取一个大的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
是否可能因为文件太大而发生?可以解决这个问题吗?
答案 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()