Php KML到CSV（谷歌地球到Agisoft Photoscan Pro）

我正在尝试编写一个脚本，将Google Earth KML文件转换为CSV文件，以便在Agisoft的Photoscan Pro中使用。该脚本将使用户能够快速创建用于优化模型的GCP（地面控制点）的CSV文件，以便在Google Earth中更准确地覆盖这些模型

如何将其写入CSV

<?php

    $inputfile = simplexml_load_file('Points.kml');
    $placemarks = $inputfile->Document->Folder->Placemark;

    print "<pre>";
    print_r($inputfile->Document->Folder->Placemark[0]);

    foreach ($placemarks as $placemark){
        $name = $placemark->name;
        $longitude = $placemark->LookAt->longitude;
        $latitude = $placemark->LookAt->latitude;
        $altitude = $placemark->LookAt->altitude;
        //print $name.",".$longitude.",".$latitude.",".$altitude;
        //print "</br>";
        $array = array("$name", "$longitude", "$latitude", "$altitude");
        //var_dump($array);
    }
?>

---

simplexml_load_文件将返回SimpleXMLElement类，所以您应该使用该类中的方法来访问xml元素

---

例如，children方法-

如果您只需要点的坐标，并且它们的标签无关紧要，我建议您可以为PhotoScan编写一个简单的Python脚本，该脚本将获得kml文件的路径，然后逐行搜索文本，如果它在给定的行中，那么您可以从字符串中提取坐标，甚至使用应用的坐标创建新的标记

下面是一个示例代码（应另存为.py文件，并通过PhotoScan Pro中的“运行脚本”对话框启动），该代码要求输入kml文件，并将提取点的x、y、z坐标写入控制台窗格：

导入PhotoScan
path=PhotoScan.app.getOpenFileName（“指定KML文件：”）
文件=打开（路径“rt”）
feof=False
计数器=0
虽然不包括：
str=file.readline（）
如果不是len（str）：
feof=True
打破
如果str中有“”：
计数器+=1
str=str.split（“，1）[1]
str=str.rsplit（“，1）[0]
（x，y，z）=str.split（“，”）
打印（x、y、z）
file.close（）文件

如果希望将坐标保存到.txt文件中（例如，使用计数器变量作为GCP标签），可以修改此代码。但我认为将这些值直接导入地面控制窗格并不困难

<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
<Document>
    <name>Points.kml</name>
    <Style id="sn_open-diamond">
        <IconStyle>
            <color>ff0e12b7</color>
            <scale>0.7</scale>
            <Icon>
                <href>http://maps.google.com/mapfiles/kml/shapes/open-diamond.png</href>
            </Icon>
        </IconStyle>
        <LabelStyle>
            <color>b30e12b7</color>
        </LabelStyle>
    </Style>
    <StyleMap id="msn_open-diamond">
        <Pair>
            <key>normal</key>
            <styleUrl>#sn_open-diamond</styleUrl>
        </Pair>
        <Pair>
            <key>highlight</key>
            <styleUrl>#sh_open-diamond</styleUrl>
        </Pair>
    </StyleMap>
    <Style id="sh_open-diamond">
        <IconStyle>
            <color>ff0e12b7</color>
            <scale>0.827273</scale>
            <Icon>
                <href>http://maps.google.com/mapfiles/kml/shapes/open-diamond.png</href>
            </Icon>
        </IconStyle>
        <LabelStyle>
            <color>b30e12b7</color>
        </LabelStyle>
    </Style>
    <Folder>
        <name>Points</name>
        <open>1</open>
        <Placemark>
            <name>1</name>
            <LookAt>
                <longitude>-81.11867631881405</longitude>
                <latitude>28.58871312454101</latitude>
                <altitude>0</altitude>
                <heading>1.04958304488296e-13</heading>
                <tilt>0</tilt>
                <range>16</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_open-diamond</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-81.11867631881405,28.58871312454101,0</coordinates>
            </Point>
        </Placemark>
        <Placemark>
            <name>2</name>
            <LookAt>
                <longitude>-81.11878160604603</longitude>
                <latitude>28.58868018607449</latitude>
                <altitude>0</altitude>
                <heading>-5.038190918931147e-05</heading>
                <tilt>0</tilt>
                <range>16.00000207857145</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_open-diamond</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-81.11878160604603,28.58868018607449,0</coordinates>
            </Point>
        </Placemark>
        <Placemark>
            <name>3</name>
            <LookAt>
                <longitude>-81.11889392441689</longitude>
                <latitude>28.58897376333365</latitude>
                <altitude>0</altitude>
                <heading>-5.676793149912727e-05</heading>
                <tilt>0</tilt>
                <range>16</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_open-diamond</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-81.11889392441689,28.58897376333365,0</coordinates>
            </Point>
        </Placemark>
        <Placemark>
            <name>4</name>
            <LookAt>
                <longitude>-81.11879074625233</longitude>
                <latitude>28.58900436190265</latitude>
                <altitude>0</altitude>
                <heading>-7.3948102396262e-06</heading>
                <tilt>0</tilt>
                <range>16.00000167331963</range>
                <gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
            </LookAt>
            <styleUrl>#msn_open-diamond</styleUrl>
            <Point>
                <gx:drawOrder>1</gx:drawOrder>
                <coordinates>-81.11879074625233,28.58900436190265,0</coordinates>
            </Point>
        </Placemark>
    </Folder>
</Document>
</kml>

# Label,X,Y,Z
1,-81.11867632,28.58871312,0
2,-81.11878161,28.58868019,0
3,-81.11889392,28.58897376,0

import PhotoScan

path = PhotoScan.app.getOpenFileName("Specify KML file:")
file = open(path, "rt")
feof = False
counter = 0

while not feof:

    str = file.readline()
    if not len(str):
        feof = True
        break

    if "<coordinates>" in str:
        counter += 1
        str = str.split("<coordinates>", 1)[1]
        str = str.rsplit("</coordinates>", 1)[0]
        (x, y, z) = str.split(",")
        print(x, y, z)

file.close()