EASE-GRID投影的产品文件(nc/hdf)转Tiff
由于最近用到的很多产品都是EASE-GIRD投影,转tif废了很大的功夫,写个文章来记录一下。
一、产品
AMSR2 产品、GlobSnow 产品、FY-3B等
二、EASE-GRID投影
关于这个投影的官方具体介绍可以看这里!
gEASE Grids Map Projection & Grid Definitions | National Snow and Ice Data Center (nsidc.org)
三、转Tiff
先看代码【python】
cols = 721
rows = 721
# 设置影像的显示范围
geo = (-9036842.76, 25067.53, 0, 9036842.76, 0, -25067.53)
# 构造projection
src_srs = osr.SpatialReference()
src_srs.ImportFromEPSG(3408) # 定义输出的坐标系为NSIDC_EASE_Grid_North"AUTHORITY["EPSG","3408"]
src_srs_wkt = src_srs.ExportToWkt() # 给新建图层赋予投影信息
#需要写入TIFF的数据
dayarray = fp['SWE_NorthernDaily_D'][:]
#输出文件名
out_file =os.path.join(output_dir, hdf_file[42:-18] + "fy.tif")ap.env.overwriteOutput = 1 # 输出文件夹里面已经有内容的,就覆盖掉
ap.CheckOutExtension("Spatial")
#转tif
write_tiff(out_file, geo, src_srs_wkt, rows, cols, dayarray)
由上段代码可知,要想转Tiff,需要行列号、geo(地理范围)、投影、数据这些信息。
(1)以GlobSnow为例
先用matlab查看数据信息
得到行列号、地理范围、投影信息
之后就可以输入代码转tiff啦
(2)如果用Matlab查看数据信息,没有地理范围,也没有写确定的投影坐标,那么,可以查看投影官方信息,根据行列号和分辨率大小确定投影和地理范围,具体如下:
例如确定分辨率为25km,根据行列号721×721,可以确定投影为NSIDC_EASE_Grid_North["EPSG","3408"],还有地理范围。
就可以输入代码转Tiff啦!
注意:这里只能是半球产品才能参考该网站所给的地理范围,如果是某个小于半球尺度的区域,需要想其他办法确定地理范围。
参考文献:
https://blog.csdn.net/weixin_38953602/article/details/101158084
EASE-GRID投影的产品文件(nc/hdf)转Tiff相关推荐
- 使用MATLAB 将EASE-Grid 2.0投影坐标系下 的NC文件转换为相同坐标系下的geotiff文件
目录 以SMOS L3 土壤水分产品数据为例 EASE-Grid2.0 不同分辨率网格参数 使用maprefcells创建对应网格参照系 转换后的tiff文件检查 以SMOS L3 土壤水分产品数据为 ...
- IDL实现MODIS Grid(正弦投影)产品的重投影及拼接处理
IDL实现MODIS Grid(正弦投影)产品的重投影及拼接处理 前言 map_proj_image函数使用关键 单个文件的重投影示例 多个文件的重投影+拼接 后记 前言 关于MODIS正弦投影产品的 ...
- ENVI_IDL:批量重投影ModisSwath产品(调用二次开发接口)+解析
目录 1. 课堂内容 1. 获取Modis Swath数据(这里只获取Lat.Lon.Aod(气溶胶厚度)三个数据集以及aod数据集的两个属性),并对aod数据进行简单的处理 2. 调用二次开发接口以 ...
- nc\hdf\h5数据格式批量提取为tif格式
需求:目前存在有nc,hdf,h5等数据格式,每一个文件按时间序列存在多个波段,需要将所有的文件提取为单波段的tif格式,可使用如下插件.
- linux互传文件nc命令
使用nc命令可以很快的在两台主机传递文件,且不需要在同一网段,只要设置好端口即可. 一.安装(CentOS下) yum install -y nc (需要root权限,可以用加上sudo) 二.使用 ...
- 【j2ee spring】38、巴巴运动网的产品文件的上传
巴巴运动网的产品文件的上传 1.项目图解 2.我们开始做我们的相应的功能模块 页面的素材我会上传的,链接是:http://download.csdn.net/detail/cutter_point/8 ...
- Python读取hdf文件并转化为tiff格式输出
目录 1.Python3读取hdf文件 2.使用Python2读取hdf文件 3. 如何使Python2和3兼容 3.1 修改Python2.7的名字 3.2 pip更新 3.2 使用pip2下载py ...
- 【Python】读取数据(nc,hdf,二进制)
目录 读取nc文件相应变量: 读取hdf文件: 解析二进制文件 : 读取nc文件相应变量: # H8data = H8data['latitude','longitude','CLTYPE'] 这样提 ...
- 代码分享 | tiff,nc,hdf,mat,grib,raw遥感数据python读取
1. tiff数据 read TIFF image # 读取tiff数据需要GDAL库 from osgeo import gdal def Read_tiff(file_path):# img_fi ...
最新文章
- codevs1137 计算系数
- MATLAB语法基础
- 东方通 中间件_东方通:中间件国产替代进程中的艰难领军者
- 虚拟示波器OSC802介绍、拆机
- 【剑指Offer学习】【全部面试题汇总】
- 【CyberSecurityLearning 70】DC系列之DC-1渗透测试(Drupal)
- boost::fusion::filter_view用法的测试程序
- springmvc rest风格化案例
- oracle中的存储过程教程,oracle 存储过程
- react直接使用bootstrap失效的原因
- H5端input标签通过手机软键盘进行搜索,并关闭软键盘
- 关于定量控制污水流量计的安装、维护与定期比对及计量
- VUE中数组赋值push与=的区别
- 云弹性:它如何影响云计算?
- In-class Initialization
- java打星星_java初学打印星星以及九九乘法表
- Hyperledger Fabric从源码分析交易
- PDManer[元数建模]-v4.1.0发布:一款简单好用的数据库建模平台
- 英伟达最新雷达工作:NVRadar:一种实时的雷达障碍检测和占位栅格预测方法
- java找不到符号_java找不到符号,请帮忙看一下,程序在附件中