将GeoDataFrame多边形转换为kml文件

时间:2016-03-25 15:28:29

标签: python-2.7 kml polygons geopandas

我有一个地缘数据GeoDataFrame,其中包含各种多边形和颜色,用于绘制气象数据(我问here的另一个问题):

        color   geometry
0   #fbfdd1 (POLYGON ((-97.12191717810094 32.569, -97.1194...
1   #f3fabf (POLYGON ((-97.12442748846019 32.569, -97.1219...
2   #ebf7b1 (POLYGON ((-97.12944810917861 32.569, -97.1269...
3   #daf0b2 (POLYGON ((-97.18969555780023 32.569, -97.1879...
4   #cbeab3 (POLYGON ((-97.18969555780023 32.5710632999095...
5   #afdfb6 (POLYGON ((-97.18467493708175 32.569, -97.1821...
6   #92d4b9 (POLYGON ((-97.17463369564484 32.5730575804109...
7   #74c9bc (POLYGON ((-97.17714400600408 32.5764063816167...
8   #5bbfc0 (POLYGON ((-97.17714400600408 32.5790959050363...
9   #40b5c3 (POLYGON ((-97.17463369564484 32.5814268890055...
10  #31a6c2 (POLYGON ((-97.17714400600408 32.5852716913413...
11  #2397c0 (POLYGON ((-97.17714400600408 32.5878055733984...
12  #1e83b9 (POLYGON ((-97.17714400600408 32.5895482376014...
13  #206eaf (POLYGON ((-97.17714400600408 32.5911487379959...
14  #2259a5 (POLYGON ((-97.17714400600408 32.5927834911588...
15  #23479d POLYGON ((-97.17463369564484 32.59421434681196...
16  #243594 POLYGON ((-97.17463369564484 32.5962866795434,...
17  #1a2b7d POLYGON ((-97.1721233852856 32.59996829071199,...

我想将其转换为kml / kmz文件,但之前我从未使用过该文件类型,因此我不确定如何继续。我尝试过使用this脚本,但它需要一些我没有的高度字段。在python中有一个好/简单的方法吗?如果可能的话,我想避免使用在线转换工具。

2 个答案:

答案 0 :(得分:1)

所以我可能找到了解决方案......

我安装了Geospatial Data Abstraction Library并且一直在使用ogr2ogr功能。

正如我在我的问题中解释的那样,我有一个带有多边形和相关颜色的geopandas GeoDataFrame,我将其写入json文件:

ogr2ogr -f KML /Users/Me/Documents/mydata.kml /Users/Me/Documents/mydata.json

在终端/命令行中,我输入:

import subprocess
subprocess.call("ogr2ogr -f KML /Users/Me/Documents/mydata.kml /Users/Me/Documents/mydata.json",shell=True)

从技术上讲,您可以使用库'subprocess'在python脚本中调用此命令:

get

这会生成一个带有基于lat / lon的多边形的kml文件。但是,它会自动将所有线条颜色设置为红色而没有填充颜色(即使我的json文件中有颜色)。我还没有找到一个很好的解决方案,所以我一直在手工编辑KML文件以获得我想要的样式。

答案 1 :(得分:1)

fiona包装的

geopandas库非官方地支持您必须手动启用的KML驱动程序。

import fiona

fiona.supported_drivers['KML'] = 'rw'
gdf.to_file('test.kml', driver='KML')

请注意,它也可以读取KML文件,但在“嵌套的” kml文件上不能很好地工作