PROSAGA码农传奇-时空索引-XML导入INSPIRE GIS数据

这是一个GML文件，因此它可以被嵌入到的OGR驱动程序读取
rgdal
和
sf
包。因此：


    
  > sf::st_layers("./DB-Netz_INSPIRE_20171116.xml")
Driver: GML
Available layers:
               layer_name geometry_type features fields
1                 Network            NA        1     12
2     ConditionOfFacility            NA     7072     15
3              MarkerPost         Point    34325     11
4    TrafficFlowDirection            NA     7072     15
5        VerticalPosition            NA     1313     15
[etc]
</code>


  
    可以读取其中的空间数据（即具有非NA几何类型的层）
     
      sf::st_read
    
    ：
  

   
    
  > nodes = sf::st_read("./DB-Netz_INSPIRE_20171116.xml","RailwayNode")
Reading layer RailwayNode' from data source/home/rowlings/Downloads/SO/train/DB-Netz_INSPIRE_20171116.xml’ using driver `GML’
Simple feature collection with 21457 features and 20 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 6.021325 ymin: 47.39779 xmax: 15.03196 ymax: 54.90462
epsg (SRID):    4258
proj4string:    +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
plot(nodes$geom)
</code>
  

  
    这产生了一组我可以看到德国很好的点。
  

  
    也可以使用读取非空间数据
     
      st_read
    
     并在可能的情况下返回数据框：
  

   
    
  > ds = sf::st_read("./DB-Netz_INSPIRE_20171116.xml","DesignSpeed")
Reading layer DesignSpeed' from data source/home/rowlings/Downloads/SO/train/DB-Netz_INSPIRE_20171116.xml’ using driver `GML’
Warning message:
no simple feature geometries present: returning a data.frame or tbl_df 
</code>
  

  
    我想这是各个轨道部分的速度限制 - 您必须查找元数据以查看ID如何在这样的表和地理数据之间匹配：
  

   
    
  > head(ds)
   gml_id                    identifier applicableDirection fromPosition
1 Spd-2046676 urnoid:Spd-2046676                            0
2 Spd-2046677 urnoid:Spd-2046677                            0
3 Spd-2046678 urnoid:Spd-2046678                            0
4 Spd-2046679 urnoid:Spd-2046679                            0
5 Spd-2046680 urnoid:Spd-2046680                            0
6 Spd-2046681 urnoid:Spd-2046681                            0
[etc etc etc etc]     
</code>