ASTER产品数据 (AST_08) 几何校正

2024-04-12 17:58

本文主要是介绍ASTER产品数据 (AST_08) 几何校正,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

前言

国内有关ASTER的博客等资料真的很少。因个人需求,用过这个数据的不少产品,分享一下校正的方法。
数据可从EarthData进行下载,需要注册账号。
在这里插入图片描述
这个数据的温度和发射率产品的分辨率可达90 m,但是回访周期特别长。几年可能也就几十幅数据。
官方有数据预处理工具HEG,安装后,不能正常运行,尝试了两台电脑都不行,所以放弃了。这里提供一种用gdal校正的思路。

AST_08 几何校正

AST_08是90 m分辨率的地表温度数据,乘以0.1可转换为开尔文单位。
由于是NASA的HDF数据,首先想到的是geoloc校正,即根据元数据里的地理位置信息自动校正。但是经过实验并不行。
经过产看文件格式,发现里面存放了具有11个像元间隔的经纬度信息。所以尝试使用添加Gcps,即控制点校正。

import gdal, osr
import numpy as npfile = r'E:\AST_08_00307242016030826_20210513075055_5224.hdf'dataset = gdal.Open(file)
sub_dataset = dataset.GetSubDatasets()def readSubdatasetAsArray(sub_dataset_name):return gdal.Open(sub_dataset_name).ReadAsArray()lst_dataset = gdal.Open(sub_dataset[0][0]) # lst的子数据集对象lats = readSubdatasetAsArray(sub_dataset[6][0])
lons = readSubdatasetAsArray(sub_dataset[7][0])row = np.linspace(0, 700, 11, endpoint=True)
col = np.linspace(0, 830, 11, endpoint=True)
cols, rows = np.meshgrid(row, col)gcps = [gdal.GCP(lons[row, col], lats[row, col], 0, rows[row, col], cols[row, col]) \for row in range(11) for col in range(11)]srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)
lst_dataset.SetGCPs(gcps, srs.ExportToWkt())correctdata = gdal.Warp(r'E:\ast08_utm.tif',  lst_dataset, dstSRS="EPSG:32649" , xRes=90, yRes=90, format='GTiff', tps=True)
del correctdata

我使用了EPSG代码指定了UTM 50N投影,有关各种投影的代码可以通过QGIS查询得到。
在这里插入图片描述
校正后影像,已和Landsat影像对比,校正的没有问题。
在这里插入图片描述

这篇关于ASTER产品数据 (AST_08) 几何校正的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/897839

相关文章

Spring 请求之传递 JSON 数据的操作方法

《Spring请求之传递JSON数据的操作方法》JSON就是一种数据格式,有自己的格式和语法,使用文本表示一个对象或数组的信息,因此JSON本质是字符串,主要负责在不同的语言中数据传递和交换,这... 目录jsON 概念JSON 语法JSON 的语法JSON 的两种结构JSON 字符串和 Java 对象互转

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

SpringBoot使用GZIP压缩反回数据问题

《SpringBoot使用GZIP压缩反回数据问题》:本文主要介绍SpringBoot使用GZIP压缩反回数据问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录SpringBoot使用GZIP压缩反回数据1、初识gzip2、gzip是什么,可以干什么?3、Spr

SpringBoot集成Milvus实现数据增删改查功能

《SpringBoot集成Milvus实现数据增删改查功能》milvus支持的语言比较多,支持python,Java,Go,node等开发语言,本文主要介绍如何使用Java语言,采用springboo... 目录1、Milvus基本概念2、添加maven依赖3、配置yml文件4、创建MilvusClient

SpringValidation数据校验之约束注解与分组校验方式

《SpringValidation数据校验之约束注解与分组校验方式》本文将深入探讨SpringValidation的核心功能,帮助开发者掌握约束注解的使用技巧和分组校验的高级应用,从而构建更加健壮和可... 目录引言一、Spring Validation基础架构1.1 jsR-380标准与Spring整合1

MySQL 中查询 VARCHAR 类型 JSON 数据的问题记录

《MySQL中查询VARCHAR类型JSON数据的问题记录》在数据库设计中,有时我们会将JSON数据存储在VARCHAR或TEXT类型字段中,本文将详细介绍如何在MySQL中有效查询存储为V... 目录一、问题背景二、mysql jsON 函数2.1 常用 JSON 函数三、查询示例3.1 基本查询3.2

SpringBatch数据写入实现

《SpringBatch数据写入实现》SpringBatch通过ItemWriter接口及其丰富的实现,提供了强大的数据写入能力,本文主要介绍了SpringBatch数据写入实现,具有一定的参考价值,... 目录python引言一、ItemWriter核心概念二、数据库写入实现三、文件写入实现四、多目标写入

使用Python将JSON,XML和YAML数据写入Excel文件

《使用Python将JSON,XML和YAML数据写入Excel文件》JSON、XML和YAML作为主流结构化数据格式,因其层次化表达能力和跨平台兼容性,已成为系统间数据交换的通用载体,本文将介绍如何... 目录如何使用python写入数据到Excel工作表用Python导入jsON数据到Excel工作表用

Mysql如何将数据按照年月分组的统计

《Mysql如何将数据按照年月分组的统计》:本文主要介绍Mysql如何将数据按照年月分组的统计方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录mysql将数据按照年月分组的统计要的效果方案总结Mysql将数据按照年月分组的统计要的效果方案① 使用 DA

鸿蒙中Axios数据请求的封装和配置方法

《鸿蒙中Axios数据请求的封装和配置方法》:本文主要介绍鸿蒙中Axios数据请求的封装和配置方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录1.配置权限 应用级权限和系统级权限2.配置网络请求的代码3.下载在Entry中 下载AxIOS4.封装Htt