我有一个描述地理坐标的大型 XML 文档(准确地说是 KML)；下面的代码片段应该可以让您了解它的外观。这里的问题是坐标是 double 的(小数点后 16 位)，这会在进一步处理时造成很多问题(此外，小数点后一位实际上是十分之一纳米 - 我们的 GPS 不是那么精确)。

我一直在寻找将精度降低到给定值的任何方法，例如小数点后 5 位给了我们一米的精度。我尝试用 Python(使用 lxml)解析 XML，更改值并保存新文档，但在处理过程中文档的格式发生了很大变化，并以某种方式中断了进一步处理。

因此，我正在寻找一种就地降低精度的方法，以便更改原始文件中的值。我认为 AWK 应该可以解决问题，但遗憾的是我的尝试没有成功。

这里是 my XML 的例子.

<Document xmlns="http://www.opengis.net/kml/2.2">
    <Folder><name>Export_Output02</name>
        <Placemark>
            <Style><LineStyle><color>ff0000ff</color></LineStyle><PolyStyle><fill>0</fill></PolyStyle></Style>
            <ExtendedData><SchemaData schemaUrl="#Export_Output02">
                <SimpleData name="species">1312</SimpleData>
                <SimpleData name="area">7848012</SimpleData>
                <SimpleData name="irrep_area">0.00000012742</SimpleData>
                <SimpleData name="groupID">2</SimpleData>
            </SchemaData></ExtendedData>
            <MultiGeometry>
                <Polygon>
                    <outerBoundaryIs>
                        <LinearRing>
                            <coordinates>-57.843052746056827,-33.032934004012787 -57.825312079170494,-33.089724736921667 -57.888494029914156,-33.073777852969904 -57.843052746056827,-33.032934004012787</coordinates>
                        </LinearRing>
                    </outerBoundaryIs>
                </Polygon>
                <Polygon>
                    <outerBoundaryIs>
                        <LinearRing>
                            <coordinates>-57.635769389832561,-33.032934004012787 -57.618028722946228,-33.089724736921667 -57.681210673689904,-33.073777852969904 -57.635769389832561,-33.032934004012787</coordinates>
                        </LinearRing>
                    </outerBoundaryIs>
                </Polygon>
            </MultiGeometry>
        </Placemark>
    </Folder>
</Document>

[编辑]

我的 Python 代码:

import glob
import argparse
from pykml import parser
from pykml.helpers import set_max_decimal_places

arg_parser = argparse.ArgumentParser(description='Script for batch reduction of precision of KML files', prog='KML precision reducer')

arg_parser.add_argument('-p', '--precision', type=int, default=5, help='Desired precision')
arg_parser.add_argument('-d', '--directory', default='./', help='Path to KML files')

args = arg_parser.parse_args()

path_to_kml = glob.glob(args.directory + '*.kml')
precision = args.precision

for kml_file in path_to_kml:
    print 'Processing ' + kml_file
    with open(kml_file) as file_read:
        doc = parser.parse(file_read)

    max_decimals={'longitude': precision, 'latitude': precision,}

    for element in doc.iter("*"):
        set_max_decimal_places(element, max_decimals)

    out_filename = kml_file.replace('.kml', '_out.kml')

    with open(out_filename, 'w') as file_write:
       doc.write(file_write, pretty_print=True, with_tail=True)

最佳答案

您可以使用 XSLT .下面的样式表使用 XSLT 2.0 .这也可以使用 XSLT 1.0，但它没有我在这里使用的 tokenize() 函数:

<?xml version="1.0" encoding="UTF-8"?>
<xsl:stylesheet 
    xmlns:xsl="http://www.w3.org/1999/XSL/Transform" 
    xmlns:gis="http://www.opengis.net/kml/2.2"
    version="2.0">

    <!-- This is an identity transform template - it copies all the nodes -->
    <xsl:template match="@*|node()">
        <xsl:copy>
            <xsl:apply-templates select="@*|node()"/>
        </xsl:copy>
    </xsl:template>

    <!-- this template has precedence over the identity template for the `coordinates` nodes -->
    <xsl:template match="gis:coordinates">
        <xsl:copy> <!-- it copies the element --> 
        <xsl:variable name="coords" select="tokenize(.,' ')"/> <!-- saves coordinate pairs in variable -->
        <xsl:for-each select="$coords"> <!-- for each coordinate pair, formats the values before and after the comma -->
            <xsl:value-of select="round(number(substring-before(.,','))*100000) div 100000"/>
            <xsl:text>,</xsl:text> <!-- puts the comma back between the coords -->
            <xsl:value-of select="round(number(substring-after(.,','))*100000) div 100000"/>
            <xsl:if test="position() != last()">
                <xsl:text> </xsl:text> <!-- puts the space back if it's not the last coord -->
            </xsl:if>
        </xsl:for-each>
        </xsl:copy>
    </xsl:template>
</xsl:stylesheet>

我在上面添加了一些注释来解释它是如何工作的。

如果将其应用于示例文档，它会将坐标截断为五位小数。这是一个显示转换后的 coordinates 元素的示例:

<LinearRing>
    <coordinates>-57.84305,-33.03293 -57.82531,-33.08972 -57.88849,-33.07378 -57.84305,-33.03293</coordinates>
</LinearRing>

这是一个 XSLT Fiddle 有一个工作结果。

我在上面的 XML Playground 中粘贴了您的完整 XML，它起作用了。我只是无法将 fiddle 与您的文件一起保存以将其链接到此处，因为文件太大，但您可以尝试将其粘贴到此处。转换整个文件大约需要 40 秒。

我不知道 Python 中的 XSLT 2.0 支持，但您可以使用命令行工具(例如 Saxon)运行转换，或者调用 Java 或其他支持 XSLT 2.0 的语言的程序(或者，如果您正在寻找针对这个特定问题的解决方案，可以使用在线工具来解决)。

关于python - 降低 XML 文档中值的精度，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/24224318/