这篇文章主要介绍“Python如何实现批量绘制遥感影像数据的直方图”，在日常操作中，相信很多人在Python如何实现批量绘制遥感影像数据的直方图问题上存在疑惑，小编查阅了各式资料，整理出简单好用的操作方法，希望对大家解答”Python如何实现批量绘制遥感影像数据的直方图”的疑惑有所帮助！接下来，请跟着小编一起来学习吧！

首先，明确一下本文需要实现的需求：现需对多幅栅格数据文件进行依据其像元数值的直方图绘制，具体绘制内容即各栅格图像像素数值的分布情况；所有栅格数据都保存在同一目标路径下，且均为

.tif

.tif

.tif

知道了需求，我们便开始代码的书写。具体代码如下：

# -*- coding: utf-8 -*-
"""
Created on Tue Jul 20 12:17:05 2021

@author: fkxxgis
"""

import os
import numpy as np
import matplotlib.pyplot as plt
from osgeo import gdal

lai_file_path="G:/Postgraduate/LAI_Glass_RTlab/LAI_Difference_Global/"
pic_save_path="G:/Postgraduate/LAI_Glass_RTlab/LAI_histogram.png"

file_list=os.listdir(lai_file_path)
tif_list=[]
for file in file_list:
    if os.path.splitext(file)[1]==".tif":
        if file[0:3]!="MCD":
            tif_list.append(file)

sub_plot_num=1
for tif in tif_list:
    lai_algorithm=tif[0:3]
    print(lai_algorithm)
    lai_raster=gdal.Open(lai_file_path+tif)
    lai_raster_array=lai_raster.ReadAsArray()
    plt.subplot(2,2,sub_plot_num)
    plt.hist(lai_raster_array)
    sub_plot_num=sub_plot_num+1
plt.show()
plt.savefig(pic_save_path, dpi=300)

其中，

lai_file_path

pic_save_path

首先，借助

os.listdir()

lai_file_path

.tif

os.path.splitext()

file

".tif"

lai_file_path

.tif

随后，还需要将我们不需要的

.tif

.tif

MCD

接下来，开始直方图的绘制。在这里我选择了将几幅直方图以子图的形式绘制在一个总图中，因此需要借助

sub_plot_num

Array

plt.hist()

最后，通过

plt.savefig(pic_save_path, dpi=300)