python gdal迭代写为ENVI的dat格式

def writeTiff(im_data, im_geotrans, im_proj, path):'''保存tif文件函数'''if 'int8' in im_data.dtype.name:datatype = gdal.GDT_Byteelif 'int16' in im_data.dtype.name:datatype = gdal.GDT_Int16else:datatype = gdal.GDT_Float32if len(im_data.shape) == 3:im_bands, im_height, im_width = im_data.shapeelif len(im_data.shape) == 2:im_data = np.array([im_data])im_bands, im_height, im_width = im_data.shape# 创建文件# dat格式需要将GTiff改成ENVI，保存文件名直接命名为datdriver = gdal.GetDriverByName("ENVI")dataset = driver.Create(path, int(im_width), int(im_height), int(1), datatype)  # 写入波段数：此时定义为1[W波段只用不保存]if(dataset != None):dataset.SetGeoTransform(im_geotrans)  # 写入仿射变换参数dataset.SetProjection(im_proj)  # 写入投影for i in range(im_bands):dataset.GetRasterBand(i+1).WriteArray(im_data[i])del dataset
if __name__=='__main__':tifNameList = os.listdir(tifFolder)tifPath = tifFolder + "/" + tifNameList[0]# 试读入第一个tif，为获取tif的基本信息dataset = readTif(tifPath)width = dataset.RasterXSize  # 栅格矩阵的列数height = dataset.RasterYSize  # 栅格矩阵的行数Tif_geotrans = dataset.GetGeoTransform()Tif_proj = dataset.GetProjection()# 数据波段for i in range(len(tifNameList)):tifPath = tifFolder + "/" + tifNameList[i]dataset = readTif(tifPath)ndvi_band = dataset.GetRasterBand(1)  # 数据波段ndvi = ndvi_band.ReadAsArray(0, 0, width, height)SavePath = "./data/modis_dat/" + \os.path.splitext(tifNameList[i])[0]+".dat"writeTiff(ndvi, Tif_geotrans, Tif_proj, SavePath)del tifPath, dataset, ndvi_band

python gdal迭代写为ENVI的dat格式相关推荐

python/gdal处理遥感影像（读取、投影转换、裁剪、建立图像金字塔等）
python/gdal处理遥感影像(读取.投影转换.裁剪.建立图像金字塔等) gdal库简单介绍 python使用gdal 一.安装python环境二.安装gdal库三.使用gdal处理遥感影像 ...
Python+GDAL几何校正任意自带经纬度数据的遥感影像
目录前言构建虚拟数据集几何校正完整代码前言李民录老师在他的博客中使用C++ GDAL的gdaltranslate.exe.gdalwarp.exe等工具对FY3A_MERSI数据进行了几何 ...
python基础教程：python+gdal+遥感图像拼接(mosaic)的实
这篇文章主要介绍了python+gdal+遥感图像拼接(mosaic)的实例,具有很好的参考价值,希望对大家有所帮助.一起跟随小编过来看看吧作为摄影测量与遥感的从业者,笔者最近开始深入研究gdal, ...
Python遥感图像处理应用篇(二十二)：Python+GDAL 批量等距离裁剪影像-续
之前写过一篇按照指定行列号数量来进行影像等距离裁剪的博客,链接如下: Python遥感图像处理应用篇(二十二):Python+GDAL 批量等距离裁剪影像_空中旋转篮球的博客-CSDN博客_pytho ...
基于Python+GDAL实现nc格式转geotiff格式
基于Python+GDAL实现nc格式转geotiff格式 1. 目的 2. 版本 3. 基础知识 3.1 什么是nc文件? 3.2 基于Python处理nc文件需要用到的库 4. 程序示例 5. ...
素数python函数怎么写_python如何判断一个数是否是素数
先来看下什么是质数: 质数(Prime number),又称素数,指在大于1的自然数中,除了1和该数自身外,无法被其他自然数整除的数(也可定义为只有1与该数本身两个因数的数). 简单来说就是,只能除以 ...
python机器学习手写算法系列——线性回归
本系列另一篇文章<决策树> https://blog.csdn.net/juwikuang/article/details/89333344 本文源代码: https://github.c ...
Python 可迭代的对象、迭代器和生成器
Python中可迭代对象和迭代器对象的理解我们在用for ... in ...语句循环时,in后面跟随的对象要求是可迭代对象,即可以直接作用于for循环的对象统称为可迭代对象(Iterable),如 ...
python实现迭代的快速排序（Iterative Quick Sort）
python实现迭代的快速排序(Iterative Quick Sort) 快速排序(Quicksort)是对冒泡排序算法的一种改进. 快速排序由C. A. R. Hoare在1960年提出.它的基本 ...
html文件怎么用Python做后端,利用python实现后端写网页（flask框架）.pdf
利利用用python实实现现后后端端写写网网页页 ((flask框框架架)) 如何用python做后端写网页-flask框架什么是Flask安装flask模块Hello World更深一步:数据绑 ...

python gdal迭代写为ENVI的dat格式

python gdal迭代写为ENVI的dat格式相关推荐

最新文章

热门文章