zonal_stats分区统计值均为nan，导致计算pearsonr时报错Constant Input Warning: An input array is constant.

计算pearson相关系数时报错：Constant Input Warning: An input array is constant. ....warning.warn (stats.constantinputwarning (msg))

原因：两个变量的矩阵均为nan，或者说均为一样的值，没有变异性，导致分母为零，计算报错

就要找两个变量的矩阵均为nan的原因

zonal_stats分区统计值均为nan：因为是两块区域构成的一个矢量图层，在分区统计时计算结果均为Nan

解决方式： 合并零散的区域成一个整体union_region = gdf.unary_union

gdf = gpd.read_file(vectors_file)
# 合并零散得区域成一个整体
union_region = gdf.unary_union
# 使用zonal_stats计算栅格数据在矢量区域内的统计量
stats = zonal_stats(union_region, file_path, stats='mean')