“
地学、水文、气象领域的自然科学数据通常以netcdf、hdf、二进制等方式存储,比如温度、降水、蒸发数据等;学会这些数据格式的读取和可视化是进行地学统计分析计算的关键,python提供了解析nc、hdf等数据格式的库,作者整理了如何利用python进行*.nc、*.hdf、*.dat格式数据的读取并进行栅格转换,欢迎转载和关注!
”arcpy读取二进制文件---以GSMaP降水数据集为例
温度、降水、蒸发等科学数据为了达到压缩文件大小,易于存储和维护等目的,通常在下载后往往以二进制*.bin、*.dat等方式呈现。因此在进行数据分析和统计计算时需要读取二进制文件进行数据还原。
注意: 读取二进制文件要特别要注意到 二进制文件数据存储的端序问题(字节序),小端序需要进行端序转换!
本地文件夹下存在2018年6月份前10天的日尺度降水数据集,下面利用arcoy脚本将该数据批量转换为tiff栅格影像,并在arcmap中展示。
源码#!/usr/bin/env python# -*- coding:utf-8 -*-# Name : gsmap2TIFF.py# Author : zengsk in NanJing# Created: 2019/8/24 23:26'''说明:1.该脚本是读取gsmap小时尺度降水数据,输出为tiff 2.运行环境需要安装python2 需要arcpy模块 3.使用arcgis自带的python环境(有arcpy模块) 4.运行结果可以直接用arcgis打开'''# 导入moduleimport osimport globimport numpy as npimport arcpyimport warningswarnings.simplefilter("ignore") # 忽略警告# 原始降水数据文件夹,可根据自己本地情况修改sPath = r'G:\NRT'oDir = r"G:\NRT\result\tiff"for fileName in glob.glob(sPath+'\*.dat'): print("Processing... {0}".fo