用GDAL绘制8、16位tif图像的像素分布图

易错点

1. tif格式还是用GDAL读取最适配
2. 16位的图像不是所有库都能读取，比如用PIL就会报错
3. 横坐标显示所有像素值会太密集，只选取几个具有代表性的像素值进行显示

代码实现

from osgeo import gdal
from matplotlib import pyplot as plt
import numpy as np
# 打开图像文件
ds = gdal.Open("1.tif")

# 读取图像数据
data = ds