利用python GDAL库读写geotiff格式的遥感影像方法
作者：彩虹弯弯 发布时间：2023-08-31 13:15:06
标签：python,GDAL,geotiff
如下所示：
from osgeo import gdal
import numpy as np
def read_tiff(inpath):
ds=gdal.Open(inpath)
row=ds.RasterXSize
col=ds.RasterYSize
band=ds.RasterCount
geoTransform=ds.GetTransform()
proj=ds.GetTransform()
data=np.zeros([row,col,band])
for i in range(band):
dt=ds.GetRasterBand(1)
data[:,:,i]=dt.ReadAsArray(0,0,col,row)
return data
def array2raster(outpath,array,geoTransform,proj):
cols=array.shape[1]
rows=array.shape[0]
driver=gdal.GetDriverByName('Gtiff')
outRaster=driver.Create(newRasterfn,cols,rows,1,gdal.GDT_Byte)
outRaster.SetGeoTransform(geoTransform)#参数2,6为水平垂直分辨率，参数3,5表示图片是指北的
outband=outRaster.GetRasterBand(1)
outband.WriteArray(array)
outRaster.SetProjection(proj)#将几何对象的数据导出为wkt格式
outRaster.FlushCache()
if _name=="_main_":
data,geoTransform,proj=read_tiff('d:/a.tif')
array2raster("d:/b.tif",np.zeros[2400,2400],geoTransform,proj)
利用python GDAL库读写geotiff格式的遥感影像，并生成与原影像具有相同地理坐标和投影坐标的geotiff格式图片。
来源：https://blog.csdn.net/qq_20340733/article/details/78316347
0
投稿猜你喜欢
- 首先安装解析的第三方包：go get gopkg.in/yaml.v2示例：package main import ( "os&q
- 我们知道，做web开发，在调试时需要反复启动整个工程，那么上一个工程占用的端口，在下一次工程启动时就不能用了，因为占用的端口没有释放，但是手
- python 循环while和for in简单实例#!/uer/bin/env python# _*_ coding: utf-8 _*_l
- 前言如果你以前没有接触过面向对象的编程语言，那你可能需要先了解一些面向对象语言的一些基本特征，在头脑里头形成一个基本的面向对象的概念，这样有
- 前言：前面的一章我们展示了选择结构的三种基本形式，并给出了实例演示，但是大家在编写代码过程中仍然要注意一些细节，尤其是代码块的缩进，这对if
如下所示：def softmax(logits, axis=None, name=None, dim=None): ""- Blackfriday是在Go中实现的Markdown处理器。您可以安全地输入用户提供的数据，速度快，支持通用扩展（表，智能标点符号替换等）
- 共同点: 1.它们都是python的核心类型，是python语言自身的一部分核心类型与非核心类型 多数核心类型可通过特定语法来生成其对象，比
- 不多说，我们直接上源码：# -*- coding:UTF-8 -*-'''实现文件打包、上传与校验Created o
- 内容摘要：统计在线人数的方法很多，可以使用Application来统计在线人数，也可以使用IP来统计在线人数。各有优点。本文介绍了通过判断S
- 如何用数据库制作一个多用户版的计数器？代码和说明如下：count.asp' 计数器的核心程序<％Set c
- 相信为数不少的系统管理员每天都在做着同一样的工作——对数据进行备份。一旦哪一天疏忽了，而这一天系统又恰恰发生了故障，需要进行数据
- 文件名： VS7.1sp1-KB918007-X86.exe版本： &n
- 实现效果UI组件依然是使用 Quasar Framework。先来看一下效果：加减.gifinput type="number&q
- 如下所示：Uploadfiles = request.FILES.get('参数'， '')for i in
- 处于兴趣，写了一个遍历指定城市五天内的天气预报，并转为华氏度显示。把城市名字写到一个列表里这样可以方便的添加城市。并附有详细注释import
- 前言在使用传统物理机或云服务器上部署项目都会存在一些痛点比如：项目部署速度慢、资源浪费、迁移难且扩展低而使用 Docker 部署项目的优势包
- assert断言声明,遇到错误则立即返回在使用python语言开发深度学习模型时，经常会遇到模型运行结束时才会发现的崩溃状态，或者得到的结果
- 一、问题描述在实习的时候，需要将两个表格的内容进行匹配分类，比如两个不同的工程项目针对的对象都是A，那么就需要将这两个工程项目归类到A当中，
- --查看指定表的外键约束 select * from sysobjects where parent_obj in( select id f
