这样的图是怎么做出来的呀

hainllj

目前在MeteoInfoLab脚本汇总贴里王老师有个帖子脚本做出的效果与这个图差不多，但代入我的数据后效果不太理想曲线拟合：http://bbs.06climate.com/forum.php?mod=viewthread&tid=51274&extra=page%3D1

yufengwang

这个东西用python来画，贼好画，画散点图即可。
以下是CSDN的链接，可以学一下，就等于至少是4维的图了。

matplotlib.pyplot.scatter(x, y, s=20, c=None, marker='o', 
cmap=None, norm=None, vmin=None, vmax=None, alpha=None, 
linewidths=None, verts=None, edgecolors=None, hold=None, 
data=None, **kwargs） 
这个函数接收的参数很多，有默认值的平时也不需要我们指定，是可选的，这次我们用到的除了基本的x ，y参数，还有c，s，alpha和marker，c就是为点指定的颜色数组，s是点的面积大小，alpha是点的颜色的透明度，marker是指定点标记的形状
————————————————
版权声明：本文为CSDN博主「tianguiyuyu」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/tianguiyuyu/article/details/80497919

yebugu

hainllj 发表于 2021-3-29 08:26
非常感谢，但我太小白了，没看懂你的代码。是MATLAB吗？没学过呢

1.使用的是R语言的ggplot2包
2.物理意义没做解释但改成了正值。色阶你可以自行调整，即将湿度升、降序排列。
3.右边界拟合的不好，可能数据少的原因。

================================================

1. library(ggplot2)
   
2. library(scales)# humidity以百分比的形式出现
   
3. 
   
4. T=sort(runif(300,min=24,max=38))
   
5. O3=8*10^(-11)*T^7.99
   
6. O3=jitter(O3,amount=25)
   
7. O3=O3-min(O3)+10# 保证O3没有负值
   
8. Humidity=sort(runif(300,0,1),decreasing = TRUE)
   
9. df=data.frame(T=T,O3=O3,Humidity=Humidity)
   
10. 
    
11. 
    
12. #你有自己的数据集，重新拟合成理想的形式，如指数模型
    
13. ## 拟合过程也可以利用熟悉的Excel或者MATLAB替代
    
14. 
    
15. 
    
16. fit_zhishu <- nls(log(O3) ~ a*log(T)+b,  start = list(a = 8, b = 0))#
    
17. fit_O3=predict(fit_zhishu, T)
    
18. df$fit_O3=exp(fit_O3) #ln反转
    
19. 
    
20. 
    
21. # 以上为我构造的数据集，你可自行替换；第一列$T$，第二列$O_3$，第三列湿度,最后一列为拟合的$O_3$
    
22. p<-ggplot(df, aes(T, O3,colour = Humidity)) +geom_point(size=3,alpha=0.8,position=position_jitter(width=0, height=5))+geom_line(aes(T,fit_O3),color="blue",size=1.2,alpha=0.9)+scale_colour_gradientn(colours = rainbow(10),labels = percent)+xlab(bquote('T/'^o*C))+ylab(bquote(O[3]/ 'mg·'*m^-3))+theme_bw()

复制代码

================================================

hainllj

yebugu 发表于 2021-3-30 11:50
1.使用的是R语言的ggplot2包
2.物理意义没做解释但改成了正值。色阶你可以自行调整，即将湿度升、降序排 ...

非常感谢，{:eb511:}{:eb511:}