1. 有關(guān) Matlab 獲取代碼關(guān)注WZZHHH回復(fù)關(guān)鍵詞，或者咸魚關(guān)注：WZZHHH123

? ? ?懷俄明探空站數(shù)據(jù)解算PWV和Tm：懷俄明探空站數(shù)據(jù)解算PWV和Tm

? ? ?懷俄明多線程下載探空站數(shù)據(jù)（包括檢查和下載遺漏數(shù)據(jù)的代碼）：懷俄明多線程下載

? ? ?對IGRAv2進(jìn)行質(zhì)量控制得到PWV和Tm的 matlab 代碼：IGRAv2進(jìn)行質(zhì)量控制得到PWV和Tm

? ? ?算 IGRAv2 探空站的 Tm 和 PWV：IGRAv2計(jì)算Tm和PWV

? ? ?提取探空站 IGRAv2 全部數(shù)據(jù)：

? ? ?ERA5 解算合集（溫度、氣壓、PWV、水汽壓和 Tm）代碼獲取：ERA5合集

MATLAB畫散點(diǎn)密度圖

? ? ? 代碼、測試數(shù)據(jù)的壓縮包在文末，可自行下載測試出圖效果

散點(diǎn)誤差圖出圖效果

代碼如下

具體使用，請自行修改

function?scatter_plot(x,y,derror)
%?點(diǎn)狀密度圖
%?根據(jù)測試數(shù)據(jù)畫的，實(shí)際應(yīng)用自己修改
%?--------------------改變密度大小------------------
numbins?=?50;
%?--------------------------------------------------X?=x;
Y?=y;
%?找到最大值最小值
MaxSP=max(x);MaxV=max(y);Maxi=1.1*max(MaxSP,MaxV);[values,?centers]?=?hist3([X?Y],?[numbins?numbins]);
centers_X?=?centers{1,1};
centers_Y?=?centers{1,2};
binsize_X?=?abs(centers_X(2)?-?centers_X(1))?/?2;
binsize_Y?=?abs(centers_Y(2)?-?centers_Y(1))?/?2;
bins_X?=?zeros(numbins,?2);
bins_Y?=?zeros(numbins,?2);
for?i?=?1:numbinsbins_X(i,?1)?=?centers_X(i)?-?binsize_X;bins_X(i,?2)?=?centers_X(i)?+?binsize_X;bins_Y(i,?1)?=?centers_Y(i)?-?binsize_Y;bins_Y(i,?2)?=?centers_Y(i)?+?binsize_Y;
end
scatter_COL?=?zeros(length(X),?1);
onepercent?=?round(length(X)?/?100);
for?i?=?1:length(X)if?(mod(i,onepercent)?==?0)fprintf('.');endlast_higher_X?=?NaN;id_X?=?NaN;c_X?=?X(i);last_lower_X?=?find(c_X?>=?bins_X(:,1));if?(~isempty(last_lower_X))last_lower_X?=?last_lower_X(end);elselast_higher_X?=?find(c_X?<=?bins_X(:,2));if?(~isempty(last_higher_X))last_higher_X?=?last_higher_X(1);endendif?(~isnan(last_lower_X))id_X?=?last_lower_X;elseif?(~isnan(last_higher_X))id_X?=?last_higher_X;endendlast_higher_Y?=?NaN;id_Y?=?NaN;c_Y?=?Y(i);last_lower_Y?=?find(c_Y?>=?bins_Y(:,1));if?(~isempty(last_lower_Y))last_lower_Y?=?last_lower_Y(end);elselast_higher_Y?=?find(c_Y?<=?bins_Y(:,2));if?(~isempty(last_higher_Y))last_higher_Y?=?last_higher_Y(1);endendif?(~isnan(last_lower_Y))id_Y?=?last_lower_Y;elseif?(~isnan(last_higher_Y))id_Y?=?last_higher_Y;endendscatter_COL(i)?=?values(id_X,?id_Y);
end%?擬合曲線
a=polyfit(x,y,1);
ax=linspace(0,Maxi+10,2000);
ay=a(1)*ax+a(2);Str_yx=['Y?=?',num2str(round(a(1),2)),'*X?+',num2str(round(a(2),2))];scatter(x,?y,?20,?scatter_COL,?'.'?);
hold?oncolormap('jet');
colorbar;xlabel('GNSS-PWV/mm','FontSize',12,'FontName','Times?New?Roman');
ylabel('ERA5-PWV/mm','FontSize',12,'FontName','Times?New?Roman');bx?=?linspace(0,120,2000);%?設(shè)置xy軸的界限
%?xlim([0?120]);ylim([0?120]);%?把誤差表示在圖上，位置根據(jù)自己需要改，我這是根據(jù)測試數(shù)據(jù)填的
text(5,115,Str_yx,'FontSize',12,'FontName','Times?New?Roman')
text(5,107,['Bias?=?',num2str(round(100*derror(1))/100,'%.2f')],'FontSize',12,'FontName','Times?New?Roman')
text(5,99,['RMSE?=?',num2str(round(100*derror(2))/100,'%.2f')],'FontSize',12,'FontName','Times?New?Roman')
text(5,91,['Cor?=?',num2str(derror(3),'%.2f')],'FontSize',12,'FontName','Times?New?Roman')%?畫擬合曲線
plot(ax,ay,'r')
hold?on%?畫對角線
plot(bx,bx,'k');
hold?ontitle('測試數(shù)據(jù)','FontSize',12)
box?on
grid?onset(gca,'LineWidth',1.2)
end

壓縮包地址如下：

? ? ? 壓縮包里面的內(nèi)容：

? ? ? 壓縮包地址：

? ? ? 關(guān)注: WZZHHH，回復(fù)?MATLAB畫散點(diǎn)密度圖

>.< 有問題，請你問你的老師，師兄師姐。我也挺忙的。