MATLAB绘制相关系数空间分布图
作者：第八星系-石头人邮箱：2205455617@qq.com
数据：era5数据
读取数据
%% 相关系数空间分布
clc;clear;close all
%% 读取数据
path = "E:\z\hhh.nc"; %设置数据文件路径
ncdisp(path) %显示数据文件的变量以及变量的精度等
mlat = double(ncread(path,'latitude')); %纬度
mlon = double(ncread(path,'longitude'));%经度
=meshgrid(mlat,mlon); %将向量转化为矩阵
t = double(ncread(path,'t2m')); %温度变量
r = double(ncread(path,'sro')); %径流量变量
sjyData = readmatrix('E:\z\hhh.txt'); %区域范围
VarName1 = sjyData(:,1)'; %区域经度
VarName2 = sjyData(:,2)'; %区域纬度
lat = reshape(nlat,57*29,1); %将矩阵转化为向量，用于画散点图
lon = reshape(nlon,57*29,1); %57和29分别代表原数据行数和列数
画图
%% 作图
figure('position',) %设置图片位置、大小
m_proj('miller','lon',,'lat',);%设置投影类型
m_contourf(nlon,nlat,cor2t,'linestyle','none');%相关系数填色
hold on
m_plot(VarName1,VarName2,'k','linewidth',1);%绘制区域边界
hold on;
shading interp;
hold on
m_scatter(lon,lat,cor5t,'k','linewi',1) %显著性检验打点
hold on
m_grid('linestyle','none'); %设置网格
hold off
colorbar
colormap(nclCM(382)) %nclCM为colormap配色补充包的函数
完整代码
%% 空间相关性作图
clc;clear;close all
%% 读取数据
path = "E:\z\hhh.nc"; %设置数据文件路径
ncdisp(path) %显示数据文件的变量以及变量的精度等
mlat = double(ncread(path,'latitude')); %纬度
mlon = double(ncread(path,'longitude'));%经度
=meshgrid(mlat,mlon); %将向量转化为矩阵
t = double(ncread(path,'t2m')); %温度变量
r = double(ncread(path,'sro')); %径流量变量
sjyData = readmatrix('E:\z\hhh.txt'); %区域范围
VarName1 = sjyData(:,1)'; %区域经度
VarName2 = sjyData(:,2)'; %区域纬度
lat = reshape(nlat,57*29,1); %将矩阵转化为向量，用于画散点图
lon = reshape(nlon,57*29,1); %57和29分别代表原数据行数和列数
%% 处理数据
% 计算相关性（corrcoef函数）
% 变量t与变量r的关系
cor1t = zeros(1,1); %定义新矩阵
cor2t = zeros(57,29);
cor3t = zeros(57,29);
for i = 1:57 %循环计算每一个个点的相关性，R表示相关性，P为显著性检验结果
for j = 1:29
for k = 1:252
cor1t(k,1) = t(i,j,k);
cor1t(k,2) = r(i,j,k);
end
= corrcoef(cor1t);
cor2t(i,j) = R(1,2);
cor3t(i,j) = P(1,2);
end
end
% 是否通过显著性检验（P<0.01）
cor4t = zeros(57,29);
for i = 1:57
for j = 1:29
if cor3t(i,j)>=0.01
cor4t(i,j) = NaN;
elseif cor3t(i,j)<0.01
cor4t(i,j) = cor3t(i,j);
end
end
end
% 筛选出需要的区域（inpolygon函数）
for i = 1:57
for j = 1:29
GZs = inpolygon(nlon(i,j),nlat(i,j),VarName1,VarName2);
if GZs == 0
cor2t(i,j) = NaN;
cor4t(i,j) = NaN;
end
end
end
cor4t(isnan(cor4t) ~= 1) = 1;
cor5t = reshape(cor4t,57*29,1);
%% 作图
figure('position',) %设置图片位置、大小
m_proj('miller','lon',,'lat',);%设置投影类型
m_contourf(nlon,nlat,cor2t,'linestyle','none');%相关系数填色
hold on
m_plot(VarName1,VarName2,'k','linewidth',1);%绘制区域边界
hold on;
shading interp;
hold on
m_scatter(lon,lat,cor5t,'k','linewi',1) %显著性检验打点
hold on
m_grid('linestyle','none'); %设置网格
hold off
colorbar
colormap(nclCM(382)) %nclCM为colormap配色补充包的函数
微信搜索“第八星系人造大气理论爱好者”公众号，关注获取文章数据
Awesome Daily Sport Predictions Tips
For the lady asking about baseball mlb picks, tenis pronosticos este semana, cricket predictions for tonight, tenis pronosticos hoy & manana, winning tennis picks strategy, american football nfl computer predictions, h2h stats picks, expert camel racing predictions against the spread, cricket computer picks, crypto casinos usa,I highly recommend this top rated useful sport prediction link or camel racing picks this weekend, best nfl betting promo codes, fanduel / betmgm nfl promo codes, nfl tips predictions for this week, basketball nba predictions against the spread, crypto casinos usa, soccer mls picks this weekend, best tennis betting promo codes, darts picks tomorrow, nfl & nba picks for today, alongside all this straight from the source about awesome sport prediction link on top of best nfl bets predictions platform, cricket predictions today straight up, tennis picks monday night, tennis predictions tomorrow, table tennis predictions today straight up, daily picks for basketball nba, best cricket betting picks site, espn / draftking soccer mls promo codes, darts predictions today, up to date camel racing news, on top of this updated best sport prediction url which is also great. Also, have a look at this look at this on new sport prediction blog as well as tennis picks each week, daily predictions news for nfl tips, nfl predictions tomorrow, tennis h2h prediction news today, camel racing predictions maker, winning ice hockey nhl predictions, winning nfl & nba picks strategy, best nfl & nba picks model, nfl & nba picks each week, h2h stats picks straight up, as well as this what do you think on new sport prediction site with daily predictions for american football nfl, basketball nba prediction news of the day, crypto casino promo codes, table tennis picks, tennis predictions maker, cricket picks monday night, basketball nba picks straight up, free darts picks today & tomorrow, tennis predictions for tonight, tennis h2h predictions maker,for good measure. Check more @ Awesome TAJIR4D Site 7236fa7Excellent Google Indexing Site
FrankJScott ??? 2024-9-6 00:37For the lady asking about baseball mlb picks, tenis pronosticos este semana, cricket predictions for ...
In response to the man asking about links google, google website crawl request, search working, request indexing search console, live url test, seo robots, crawl as googlebot, google re index page, indexing api wordpress, submit my site,I highly recommend this her explanation about google indexing blog or search crawl, google search index update, noindex website, google search console xml sitemap, auto index google, submit your url to google, my website not appearing on google, check if a site is indexed by google, index of google, google crawl my website, not to mention this recommended google indexing advice as well as google search my site, test noindex, google fast indexing, instant backlink indexer, meta robots noindex tag, url inspection google, free google index checker, tell google about your website, google seo crawler, add my site to google, and don't forget this top google indexing url which is also great. Also, have a look at this right here on google indexing details not to mention quick indexing, google check page seo, submit domain to google, bulk url submission, indexing crawling, google search console no index, page index check, remove a page from google index, website page indexing, submit new website to google, alongside all this a fantastic read about google indexing info with link search, google search console crawl rate, google console indexing, google check index, increase google crawl rate, best backlink indexing tool, submit page for indexing google, indexing wordpress site, new website not showing in google search, crawl page google,for good measure. Check more @ Updated 7RAJATOGEL Guide 6fa7e03
Recommended AMANSLOT88 Guide
FrankJScott ??? 2024-9-6 00:37For the lady asking about baseball mlb picks, tenis pronosticos este semana, cricket predictions for ...
For the lady talking about judi sgp, sdy slot 4d, sgp judi, sydney singapore hongkong prize, online toto 4d, togel online 2023, data sgp hk sydney, judi toto online, judi hk online, judi 4d,I highly suggest this helpful resources about AMAN88 blog or sydney singapore hongkong prize, singapore hk, hk sdy, toto singapore online, hk sgp sdy, sdy slot login, judi hk online, togel online 2023, data hk sdy sgp, sdy hk, not to mention this website about AMANSLOT88 forum which is worth considering with judi toto, nama 4d, hk sdy sgp, data sgp hk sydney, result sdy sgp hk, sdy sgp hk pools, hk singapore, 4d setiap hari apa, slot 4d singapore, 4d asia, not forgetting sites such as this here are the findings for AMAN88 info which is also great. Also, have a look at this visit website about AMANSLOT88 advice not to mention togel online 2023, bandar judi toto, singapore prize hk, bandar togel idn, sgp judi, singapore 4d slot, hk singapore, nama judi togel, online toto 4d, judi sgp online, alongside all this watch this video for AMAN88 forum not to mention bandar judi 4d, sgp sydney, hk sdy sgp, singapore 4d slot, toto singapore online,such a good point on as well as online toto 4d, sdy hk sgp, toto 4d new, judi toto online, result sdy sgp hk,for good measure. Check more @ Best DVL TOTO Guide decb1e7
页:
[1]