python投影Himawari-8数据

IDL是可以构建GLT对影像进行校正的，在python里面应该是用gdal.warp()实现的，据传可以编写VRT文件写入经纬度文件的路径就可以实现和GLT校正同样的功能，可惜我太菜了不会搞，就弄个笨办法来用了。
输入文件为HSD数据提取的影像或者NC文件中得到的影像数据，具体长这样：

import gdal
inputfile='*\\HS_H08_20190916_0300_FLDK_B3_R20.dat'
outtif='*********.tif'
memDs = gdal.Open(inputfile)
cols = memDs.RasterXSize
rows = memDs.RasterYSize
if rows == 11000 and cols == 11000:res = 0.01
else:res = 0.02
# 几何校正
# 定义空间参考
srs = osr.SpatialReference()
# 定义地球长半轴a=6378137.0m，地球短半轴b=6356752.3m，卫星星下点所在经度140.7，目标空间参考
srs.ImportFromProj4('+proj=geos +h=35785863 +a=6378137.0 +b=6356752.3 +lon_0=140.7 +no_defs')
memDs.SetProjection(srs.ExportToWkt())
memDs.SetGeoTransform([-cols*res*50000, int(res*100000), 0, cols*res*50000, 0, int(-res*100000)])
dstFilePath = os.path.join(outfolder,'H8_'+str(cols)+'.tif')
if os.path.exists(outtif):os.remove(outtif)
warpDs = gdal.Warp(outtif, memDs, dstSRS='EPSG:4326', outputBounds=(60.0, -90.0, 222.0, 90.0), xRes=res, yRes=res)
del warpDs

(60.0, -90.0, 222.0, 90.0)的范围是利用IDL计算影像像元坐标时算出来的经纬度范围估算的外边界范围，用起来误差还是不太大的。校正完影像长这样（跟矢量很贴合了）：

放大看这里：

tips:
srs.ImportFromProj4(’+proj=geos +h=35785863 +a=6378137.0 +b=6356752.3 +lon_0=140.7 +no_defs’)
代码中的数据换成FY4A的就可以对FY4A影像进行投影的，亲测可用。
outputBounds=(60.0, -90.0, 222.0, 90.0)替换为outputBounds=(80.0, -60.0, 200.0, 60.0)，生成的影像即为NC文件对应的经纬度范围内的影像（省去裁剪了），不过NC影像大小为6001X6001，使用gdal.Warp生成的为6000X6000，问题不大，6000X6000的影像如下:

python投影Himawari-8数据相关推荐

arcgis怎么运行python_在arcgis上用python脚本(arcpy)做数据批处理
arcgis软件是测绘地理信息行业空间数据处理与数据分析比较普及的一套软件系统,具有良好的扩展功能,其提供的基于python的开发包arcpy为数据的批处理提供了很好的基础.本文就是介绍在arcgis ...
基于MaxCompute分布式Python能力的大规模数据科学分析
简介: 如何利用云上分布式 Python 加速数据科学. 如果你熟悉 numpy.pandas 或者 sklearn 这样的数据科学技术栈,同时又受限于平台的计算性能无法处理,本文介绍的 MaxCom ...
python读取nc文件数据生成专业气象图
python读取nc文件数据生成专业气象图如色斑图.等值线.风向标,如下图所示: 代码如下: # -*- coding: utf-8 -*- """ Created on ...
Python数据挖掘：数据转换-数据规范化
来源:天善智能韦玮老师课堂笔记作者:Dust ---------- 数据转换: 简单变换 1.数据变换的目的是降数据转化为更方便分析的数据. 2.简单变换通常使用函数变换的方式进行,常见的函数变换包 ...
【python】使用python脚本将LFW数据中1672组同一个人多张照片拷贝出来
使用python脚本将LFW数据中1672组同一个人多张照片拷贝出来 dataCleaning4multiple.py 源码如下: import os, random, shutil import s ...
python游戏中调整箭头下落速度_入门 | 三行Python代码，让数据预处理速度提高2到6倍...
原标题:入门 | 三行Python代码,让数据预处理速度提高2到6倍选自TowardsDataScience 作者:George Seif,机器之心编译在 Python 中,我们可以找到原生的并行 ...
Python使用pandas设置数据列中float数据类型的有效小数位数、抑制科学计数法
Python使用pandas设置数据列中float数据类型的有效小数位数.抑制科学计数法目录
Python使用matplotlib绘制数据去重前后的柱状图对比图（在同一个图中显示去重操作之后同一数据集的变化情况）
Python使用matplotlib绘制数据去重前后的柱状图对比图(在同一个图中显示去重操作之后同一数据集的变化情况) #仿真数据 import pandas as pd import numpy a ...
python使用fpdf生成数据报告pdf文件
python使用fpdf生成数据报告pdf文件目录 python使用fpdf生成数据报告pdf文件 # 报告生成整体代码
python与R行列数据过滤（row column filtering）:dplyr、 isnull、isna、drop、select、iloc、loc、isin、filter
python与R行列数据过滤(row column filtering):dplyr. isnull.isna.drop.select.iloc.loc.isin.filter 很多工程师可能刚开始的 ...

python投影Himawari-8数据

python投影Himawari-8数据相关推荐

最新文章

热门文章