Python绘制卫星降雨校正散点密度图

最新推荐文章于 2025-03-20 21:20:11 发布

地狱道

最新推荐文章于 2025-03-20 21:20:11 发布

阅读量355

点赞数

文章标签： python 开发语言

本文链接：https://blog.youkuaiyun.com/2301_77925375/article/details/131389630

版权

这段Python代码展示了如何利用matplotlib和scipy库绘制卫星降雨数据的散点密度图，同时进行回归分析。通过gaussian_kde计算点密度并用颜色表示，还计算了相关性、BIAS、RMSE等统计指标，以及绘制了1:1线和回归线。最终保存为plot图像。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy.stats import gaussian_kde
from sklearn.metrics import mean_squared_error
from matplotlib.pyplot import MultipleLocator
from statistics import mean
from matplotlib import rcParams
config = {"font.family":'Times New Roman',"font.size": 16,"mathtext.fontset":'stix'}
rcParams.update(config)
# 读取数据
filename=r'./data613.xlsx'
df2=pd.read_excel(filename)#读取文件
p1=df2['p1']
p2=df2['p2']
def get_regression_line(real,pred,data_range=(0,120)):
    # 拟合
    def slope(xs,ys):
        m = (((mean(xs) * mean(ys)) - mean(xs * ys)) / ((mean(xs) * mean(xs)) - mean(xs * xs)))
        b = mean(ys) - m * mean(xs)
        return m, b
    k, b = slope(real,pred)
    regression_line = []