利用Python下的FastKML包读取KML文件

2018-05-08

在网上下载的台站信息一般会被XML标准包装，如果需要提取很多台站信息的话，不可能从XML中手动提取出来。而有的台站信息是被kml形式包装的，其实KML(介绍)完全符合XML标准。下面会讲到如何用FastKML提取KML格式的数据。

修改IRIS上下载的XML

这是从IRIS下载的澳大利亚的一个台阵，kml文件部分信息显示如下：

<Document xmlns:atom="http://www.w3.org/2005/Atom" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns="http://www.opengis.net/kml/2.2">
  <name>GMAP Stations</name>
  <Style id="33CC33">
    <IconStyle>
      <Icon>
        <href>http://ds.iris.edu/static/img/markers/circle-dot-10x10-33CC33.27295a36c313.png</href>
      </Icon>
    </IconStyle>
  </Style>
  <Placemark>
    <name>AU.AS01</name>
    <styleUrl>#33CC33</styleUrl>
    <Point>
      <coordinates>133.950806,-23.6647</coordinates>
    </Point>
    <Style>
      <BalloonStyle>
        <text>&lt;div style='width: 275px;'&gt;&lt;h4&gt;AU AS01&lt;/h4&gt;Alice Springs Array Element 01&lt;br /&gt;Geoscience Australia&lt;br /&gt;&lt;br /&gt;Data Center: &lt;a href='http://ds.iris.edu'&gt;IRISDMC&lt;/a&gt;&lt;br /&gt;Operating Range: 2003-08-13 - 2599-12-31&lt;br /&gt;&lt;br /&gt;Geographic Location: -23.6647, 133.950806&lt;br /&gt;&lt;a href='https://ds.iris.edu/mda/AU/AS01?timewindow=2003-2599'&gt;More Information&lt;/a&gt;&lt;/div&gt;</text>
      </BalloonStyle>
    </Style>
  </Placemark>
 </Document>

补充缺失的说明信息

文件开头加入的语句
<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2">

文件末尾加入的语句
</kml>

修改后的XML形式

<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2">
<Document xmlns:atom="http://www.w3.org/2005/Atom" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns="http://www.opengis.net/kml/2.2">
  <name>GMAP Stations</name>
  <Style id="33CC33">
    <IconStyle>
      <Icon>
        <href>http://ds.iris.edu/static/img/markers/circle-dot-10x10-33CC33.27295a36c313.png</href>
      </Icon>
    </IconStyle>
  </Style>
  <Placemark>
    <name>AU.AS01</name>
    <styleUrl>#33CC33</styleUrl>
    <Point>
      <coordinates>133.950806,-23.6647</coordinates>
    </Point>
    <Style>
      <BalloonStyle>
        <text>&lt;div style='width: 275px;'&gt;&lt;h4&gt;AU AS01&lt;/h4&gt;Alice Springs Array Element 01&lt;br /&gt;Geoscience Australia&lt;br /&gt;&lt;br /&gt;Data Center: &lt;a href='http://ds.iris.edu'&gt;IRISDMC&lt;/a&gt;&lt;br /&gt;Operating Range: 2003-08-13 - 2599-12-31&lt;br /&gt;&lt;br /&gt;Geographic Location: -23.6647, 133.950806&lt;br /&gt;&lt;a href='https://ds.iris.edu/mda/AU/AS01?timewindow=2003-2599'&gt;More Information&lt;/a&gt;&lt;/div&gt;</text>
      </BalloonStyle>
    </Style>
  </Placemark>
 </Document>
 </kml>

FastXML安装

1	pip install fastkml

提取台站信息的py2脚本

from fastkml import kml
with open('gmap-stations.kml','rt') as f:
    doc=f.read()
    
# Create the KML object to store the parsed result
k = kml.KML()

# Read in the KML string
k.from_string(doc)

#  Check that the number of features is correct
#  This corresponds to the single "Document" 
features = list(k.features())

# Check that we can access the features as a generator

f2 = list(features[0].features())

#output the position of stations
with open('station.loc','w') as f:
	for sta in f2:
		f.write(sta.name+'\t'+str(sta.geometry.x)+'\t'+str(sta.geometry.y)+'\n')
		print sta.name, sta.geometry.x, sta.geometry.y

输出信息

AU.AS01 133.950806 -23.6647
AU.AS02 133.936905 -23.679199
AU.AS03 133.919693 -23.677401
AU.AS04 133.929199 -23.6597
AU.AS05 133.947495 -23.6492
AU.AS06 133.971405 -23.647499
AU.AS07 133.969696 -23.6656
AU.AS08 133.960007 -23.6814
AU.AS09 133.941406 -23.6994
AU.AS10 133.913895 -23.696899
AU.AS11 133.897797 -23.678301
AU.AS12 133.904404 -23.666401
AU.AS13 133.894394 -23.651899
AU.AS14 133.910294 -23.651699
AU.AS15 133.912201 -23.635599
AU.AS16 133.929993 -23.6369
AU.AS17 133.991699 -23.6644
AU.AS18 133.980804 -23.690001
AU.AS19 133.962494 -23.704399
AU.AS31 133.904007 -23.666

台站信息来源网站

FDSN查看台站信息，这里列出台网代号等供查询，可以查询永久台站和流动台站信息。

在IRIS下的gmap下载台站信息。这里需要自己输入台网代号和台站名字。台站信息存储的格式为KML格式