用于读取 LANDSAT XML 图元文件的 Python 脚本

标签 python landsat

我需要获取多个 Landsat 图像的场景 ID。图像的元数据包含在 xml 文件中:

<?xml version="1.0"?>
<searchResponse xsi:schemaLocation="http://upe.ldcm.usgs.gov/schema/metadata
http://earthexplorer.usgs.gov/EE/metadata.xsd" 
xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"    
xmlns="http://upe.ldcm.usgs.gov/schema/metadata">
    <metaData>
        <browseAvailable>Y</browseAvailable
        <browseURL>http://earthexplorer.usgs.gov/browse/tm/146/41/1998/LT51460411998300XXX01.jpg</browseURL> 
        <sceneID>LT51460411998300XXX01</sceneID> 
        <sensor>LANDSAT_TM</sensor> 
        <acquisitionDate>1998-10-27</acquisitionDate> 
        <dateUpdated>2012-07-31</dateUpdated> 
        <path>146</path> 
        <row>41</row> 
        <...>
     </metaData>
     <metaData>
        <sceneID>LT51460411998300XXX01</sceneID> 
        <sensor>LANDSAT_TM</sensor> 
        <acquisitionDate>1998-10-27</acquisitionDate> 
        <dateUpdated>2012-07-31</dateUpdated> 
        <path>146</path> 
        <row>41</row> 
        <...>
    <etc etc>

下面的代码将列出元数据中的所有节点(browseAvailable、browseURL、sceneID),但我真正想做的只是获取 sceneID 列表。

#!/usr/bin/python

#import os
#import sys
from xml.dom.minidom import parse

import xml.etree.ElementTree as ET
#import grass.script as grass

tree = ET.parse('C:/Users/Simon/Documents/tif/Metadata/metadata_test1.xml')
root = tree.getroot()

for metadata in root:  
    for data in metadata:          
    sceneID = data.text
        grass.message('data -> %s' % (sceneID))

我尝试使用 Python xml.etree.ElementTree 上描述的工具,但没有成功。谁能帮我调整上面的代码以仅返回 sceneID?

最佳答案

试试这个:

#!/usr/bin/python

from xml.dom.minidom import parse

dom = parse('C:/Users/Simon/Documents/tif/Metadata/metadata_test1.xml')

ids  = dom.getElementsByTagName('sceneID')
for id in ids:
    print "%s" % id.childNodes[0].nodeValue

关于用于读取 LANDSAT XML 图元文件的 Python 脚本,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13256732/

相关文章:

python-3.x - 如何使用 Amazon S3 下载 Landsat 8 图像

r 栅格覆盖函数 - landsat 堆栈错误

time-series - 谷歌地球引擎: how to map a function over a collection of all Landsat sensors to create NDVI timeseries

javascript - Python GEE 从 Landsat 图像中提取特征集合时间序列

python - 在一行中构建 for 循环并收到错误的输出,这是一个 Nonetype

python - 为什么导入 python 模块不导入嵌套模块?

python - flask 日志记录 - 无法将其写入文件

Python Matplotlib – 在 x 轴上代表采样位置的条形图

javascript - 谷歌地球引擎: Region of Landsat Image

Python lambda自组合函数: How to return a function that is repeated n amount of times?