Блог

Главная - Вопросы по программированию - Как извлечь данные из файла GML

Как извлечь данные из файла GML

#python #gml

#python #gml

Вопрос:

У меня есть текстовый файл, и я хотел бы извлечь <gml:pos>73664.300 836542.700</gml:pos> его из него. Точнее, я хотел бы получить систему координат GPS [73664.300 836542.700] из pos тега. Файл содержит несколько <wfs:member> , и каждый из них имеет <gml:pos> (самый глубокий слой).

 <?xml version='1.0' encoding='UTF-8'?>
<wfs:FeatureCollection xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wfs/2.0 http://schemas.opengis.net/wfs/2.0/wfs.xsd http://www.opengis.net/gml/3.2 http://schemas.opengis.net/gml/3.2.1/gml.xsd http://www.deegree.org/app https://web.de/feature_descr?SERVICE=WFSamp;amp;VERSION=2.0.0amp;amp;REQUEST=DescribeFeatureTypeamp;amp;OUTPUTFORMAT=application/gml+xml; version=3.2amp;amp;TYPENAME=app:lsa_dataamp;amp;NAMESPACES=xmlns(app,http://www.deegree.org/app)" xmlns:wfs="http://www.opengis.net/wfs/2.0" timeStamp="2020-11-18T15:01:17Z" xmlns:gml="http://www.opengis.net/gml/3.2" numberMatched="unknown" numberReturned="0">
  <!--NOTE: numberReturned attribute should be 'unknown' as well, but this would not validate against the current version of the WFS 2.0 schema (change upcoming). See change request (CR 144): https://portal.opengeospatial.org/files?fact_id=6798.-->
  <wfs:member>
    <app:dat_set xmlns:app="http://www.deegree.org/app" gml:id="app:dat_set_1">
      <app:point>2</app:point>
      <app:art>K         </app:art>
      <app:L_Name>westt / woustest             </app:L_Name>
      <app:geom>
        <!--Inlined geometry 'data_1_APP_GEOM'-->
        <gml:MultiPoint gml:id="data_1_APP_GEOM" srsName="EPSG:25832">
          <gml:pointMember>
            <gml:Point gml:id="GEOMETRY_ad608059-f297-4554-8464-cdde248cb531" srsName="EPSG:25832">
              <gml:pos>73664.300 836542.700</gml:pos>
            </gml:Point>
          </gml:pointMember>
        </gml:MultiPoint>
      </app:geom>
    </app:lsa_pointdata>
  </wfs:member>
  <wfs:member>
    <app:dat_set xmlns:app="http://www.deegree.org/app" gml:id="app:dat_set_2">
      <app:point>3</app:point>
      <app:art>K         </app:art>
      <app:L_Name>route / riztr        </app:L_Name>
      <app:geom>
        <!--Inlined geometry 'data_2_APP_GEOM'-->
        <gml:MultiPoint gml:id="data_2_APP_GEOM" srsName="EPSG:25832">
          <gml:pointMember>
            <gml:Point gml:id="GEOMETRY_440d8630-b674-4768-a5b7-3fab46d9ac8c" srsName="EPSG:25832">
              <gml:pos>74354.900 837456.300</gml:pos>
            </gml:Point>
          </gml:pointMember>
        </gml:MultiPoint>
      </app:geom>
    </app:lsa_pointdata>
  </wfs:member>
  <wfs:member>
    ...
...

Как я могу получить эти координаты gps?
Заранее благодарю вас.

Ответ №1:

Вы можете использовать lxml и XPATH.

 data = b'''
<?xml version='1.0' encoding='UTF-8'?>
<wfs:FeatureCollection xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wfs/2.0 http://schemas.opengis.net/wfs/2.0/wfs.xsd http://www.opengis.net/gml/3.2 http://schemas.opengis.net/gml/3.2.1/gml.xsd http://www.deegree.org/app https://web.de/feature_descr?SERVICE=WFSamp;amp;VERSION=2.0.0amp;amp;REQUEST=DescribeFeatureTypeamp;amp;OUTPUTFORMAT=application/gml+xml; version=3.2amp;amp;TYPENAME=app:lsa_dataamp;amp;NAMESPACES=xmlns(app,http://www.deegree.org/app)" xmlns:wfs="http://www.opengis.net/wfs/2.0" timeStamp="2020-11-18T15:01:17Z" xmlns:gml="http://www.opengis.net/gml/3.2" numberMatched="unknown" numberReturned="0">
  <!--NOTE: numberReturned attribute should be 'unknown' as well, but this would not validate against the current version of the WFS 2.0 schema (change upcoming). See change request (CR 144): https://portal.opengeospatial.org/files?fact_id=6798.-->
  <wfs:member>
    <app:dat_set xmlns:app="http://www.deegree.org/app" gml:id="app:dat_set_1">
      <app:point>2</app:point>
      <app:art>K         </app:art>
      <app:L_Name>westt / woustest             </app:L_Name>
      <app:geom>
        <!--Inlined geometry 'data_1_APP_GEOM'-->
        <gml:MultiPoint gml:id="data_1_APP_GEOM" srsName="EPSG:25832">
          <gml:pointMember>
            <gml:Point gml:id="GEOMETRY_ad608059-f297-4554-8464-cdde248cb531" srsName="EPSG:25832">
              <gml:pos>73664.300 836542.700</gml:pos>
            </gml:Point>
          </gml:pointMember>
        </gml:MultiPoint>
      </app:geom>
    </app:dat_set>
  </wfs:member>
  <wfs:member>
    <app:dat_set xmlns:app="http://www.deegree.org/app" gml:id="app:dat_set_2">
      <app:point>3</app:point>
      <app:art>K         </app:art>
      <app:L_Name>route / riztr        </app:L_Name>
      <app:geom>
        <!--Inlined geometry 'data_2_APP_GEOM'-->
        <gml:MultiPoint gml:id="data_2_APP_GEOM" srsName="EPSG:25832">
          <gml:pointMember>
            <gml:Point gml:id="GEOMETRY_440d8630-b674-4768-a5b7-3fab46d9ac8c" srsName="EPSG:25832">
              <gml:pos>74354.900 837456.300</gml:pos>
            </gml:Point>
          </gml:pointMember>
        </gml:MultiPoint>
      </app:geom>
    </app:dat_set>
  </wfs:member>
</wfs:FeatureCollection>
'''

from lxml import etree
from io import BytesIO
f = BytesIO(data)

ns = {"gml": "http://www.opengis.net/gml/3.2"}
tree = etree.parse(f)
for e in tree.findall("//gml:pos", ns):
    print(e.text)

Комментарии:

1. Большое спасибо за ваш ответ @Justin. Это работает нормально. Однако существует огромный текстовый файл, поэтому я не могу скопировать и вставить его в скрипт python. Итак, я добавил data = b'''` at the beginning and «‘ » в конце данных, как и вы, затем я попытался открыть его : data_set = open("path_data/col_data.txt","a") , и я получаю сообщение об ошибке: TypeError: a bytes-like object is required, not '_io.TextIOWrapper' at f = BytesIO(data_set) . Есть ли способ управлять этим? Спасибо

2. Данные не обязательно включать в ваш скрипт на python. Я просто не хотел создавать отдельный файл только для того, чтобы дать вам ответ. Это было для демонстрационных целей.

3. Хорошо, я понимаю. Однако я получаю ошибку и при чтении файла. data_set = open("path_data/col_data.txt","r") . data_set выходы <_io.TextIOWrapper name='path_data/col_data.txt' mode='r' encoding='UTF-8'> . Не могли бы вы, пожалуйста, поддержать меня и там?

4. Пробовали tree = etree.parse(data_set) ли вы после этого data_set = open("path_data/col_data.txt","r") ?

5. Смотрите lxml.de/tutorial.html#parsing-from-strings-and-files , в частности, раздел о «Функции parse()»