Python-GDAL读取遥感影像直方图统计

通过GDAL+Python对遥感影像进行直方图统计

from osgeo import gdal
import numpy as np
import matplotlib.pyplot as plt# 通过GDAL读取栅格影像
filename = "/Volumes/王金阳的仓库/植被固碳计算/analyse/result/2006result/old/arable2/2006arable2.tif"
dataset = gdal.Open(filename)
im_width = dataset.RasterXSize
im_height = dataset.RasterYSize
im_data = dataset.ReadAsArray(0, 0, im_width, im_height)
print(im_data.shape)###########显示灰度直方图########
# 遍历影像中的每一个像元的像元值
data = []
for i in range(im_data.shape[0]):for j in range(im_data.shape[1]):# print(j)data.append(im_data[i][j])
data.sort()# 统计最大最小值
data = np.array(data)
print(data.min(), data.max())# 由于空值通常用-3.4028235e+38表达，避免错误统计需要提前剔除
a = sum(data == 0)
print(a)
print(data.shape)
data = data[a:]
# 根据影像中最大最小值设定坐标轴
bins = np.linspace(10, 100, 100)
# 绘制直方图，设定直方图颜色
plt.hist(data, bins, facecolor="blue")
# 横坐标轴名称
plt.xlabel('像元值')
# 纵坐标轴名称
plt.ylabel('频数')
# 图表头名称
plt.title('灰度分布直方图')
# 显示中文字体
plt.rcParams['font.sans-serif'] = ['Source Han Sans CN']
# 导出绘制得到的图片
plt.savefig('./test2.jpg')
plt.show()
band1_data = im_data

直方图统计结果

Python-GDAL读取遥感影像直方图统计相关推荐

python/gdal处理遥感影像（读取、投影转换、裁剪、建立图像金字塔等）
python/gdal处理遥感影像(读取.投影转换.裁剪.建立图像金字塔等) gdal库简单介绍 python使用gdal 一.安装python环境二.安装gdal库三.使用gdal处理遥感影像 ...
8影像计算ndvi landsat_使用GDAL读取遥感影像的信息
读取影像数据集的元数据 GDAL已经提供了足够方便的函数,可以读取影像的一些元数据信息, 从而方便对数据进行处理.GDAL一般是以字典的形式对元数据进行组织的, 但是对于不同的栅格数据类型,元数据的类 ...
3_读取遥感影像信息
读取遥感影像信息本文中的所有方法可以读取遥感影像信息,但是仅限于概括地获取元数据信息,无法具体到对像元进行处理. 打开已有的GeoTiff GDAL不能通过with打开影像,会报错 GDAL打开影像 ...
gdal进行遥感影像读写_如何使用遥感影像进行矿物勘探
gdal进行遥感影像读写 Meet Jose Manuel Lattus, a geologist from Chile. In the latest Soar Cast, he discusses ...
python文件读取操作练习题（统计单词）
python文件读取操作练习题(统计单词) 学习路线:python的文件读取基础入门(read(),readlines(),with.open()) ->python文件读取操作练习题(统计单词 ...
Python地学分析 — GDAL对遥感影像重投影
欢迎关注博主的微信公众号:"智能遥感". 该公众号将为您奉上Python地学分析.爬虫.数据分析.Web开发.机器学习.深度学习等热门源代码. 本人的GitHub代码资料主页(持续 ...
原创程序|基于GDAL的遥感影像批量处理工具介绍（三）
本文主要介绍基于C#+GDAL-Python实用工具开发的遥感影像批量处理工具,该工具目前主要包括影像批量切片生成KML文件和影像批量切片生成Tiles文件.该工具.Net框架版本为4.0,GDAL版 ...
基于c语言的遥感图像处理,基于GDAL的遥感影像显示（C#版）
在matlab中实现遥感影像和shp文件的结合显示 clc;close all;clear; road=shaperead('boston_roads.shp'); %读取shape文件 figure ...
[Python图像处理] 三十七.OpenCV直方图统计两万字详解（掩膜直方图、灰度直方图对比、黑夜白天预测）
该系列文章是讲解Python OpenCV图像处理知识,前期主要讲解图像入门.OpenCV基础用法,中期讲解图像处理的各种算法,包括图像锐化算子.图像增强技术.图像分割等,后期结合深度学习研究图像识别 ...
Python+pandas读取Excel文件并统计演员参演电影数量
Excel样本数据请参考Python读取Excel文件并统计演员参演电影 >>> import pandas as pd >>> df = pd.read_exce ...

Python-GDAL读取遥感影像直方图统计

通过GDAL+Python对遥感影像进行直方图统计

直方图统计结果

Python-GDAL读取遥感影像直方图统计相关推荐

最新文章

热门文章