找回密码
 注册
关于网站域名变更的通知
查看: 254|回复: 2
打印 上一主题 下一主题

求助各位大神,关于matlab克里金插值问题!

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-12-23 10:52 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
本人现在有本市地图数据,以及自动气象站站点经纬度及温度数值,想用克里金插值方法画等温线图,怎么画,请各位帮助,其中wlmq_dt为地图数据,wlmq_zdz第一列和第二列为站点经纬度,第五列为温度。
; v  `( g5 B" W3 @$ ?
$ Y1 z5 ~9 |2 e# A: i! W. P

该用户从未签到

2#
发表于 2021-12-23 13:28 | 只看该作者
本帖最后由 ononsiiii197 于 2021-12-23 13:30 编辑
, m* w( |$ a1 r* O# a: \( A- H9 Q% f2 j5 J
' m7 V# _% T: i2 R# [
  • data1 = load('wlmq_dt.txt');
  • latlim = minmax(data1(:,2)');
  • lonlim = minmax(data1(:,1)');
  • MLabelLocation = linspace(lonlim(1),lonlim(2),6);
  • PLabelLocation = linspace(latlim(1),latlim(2),6);
  • figure('color','w',...
  •     'name','等温线图',...
  •     'numbertitle','off',...
  •     'units','normalized',...
  •     'pos',[0.2738,0.0638,0.4451,0.8125]);
  • axesm('MapProjection','mercator', ...
  •     'MapLatLimit',latlim, ...
  •     'MapLonLimit',lonlim, ...
  •     'Frame','on', ...
  •     'Grid','on', ...
  •     'MeridianLabel','on', ...
  •     'ParallelLabel','on', ...
  •     'MLabelLocation',MLabelLocation, ...
  •     'MLineLoc',MLabelLocation, ...
  •     'PLabelLocation',PLabelLocation, ...
  •     'PLineLoc',PLabelLocation);
  • tightmap;
  • plotm(data1(:,2),data1(:,1),'k');
  • data2 = xlsread('wlmq_zdz.xls');
  • lon = data2(:,1);
  • lat = data2(:,2);
  • tem = data2(:,5);
  • theta = [10 10]; lob = [1e-1 1e-1]; upb = [20 20];
  • [dmodel, perf] = dacefit([lat,lon], tem, @regpoly0, @corrgauss, theta, lob, upb);
  • LonLat = gridsamp([min(latlim) min(lonlim);max(latlim) max(lonlim)], 60);
  • TemNew = predictor(LonLat, dmodel);
  • LatNew = reshape(LonLat(:,1),[60,60]);
  • LonNew = reshape(LonLat(:,2),[60,60]);
  • TemNew = reshape(TemNew, size(LonNew));
  • geoshow(LatNew,LonNew,TemNew,'DisplayType','surface');
  • hold on
  • plotm(lat,lon,'k.');
  • colorbar;
    2 e' B: T% m' q) S- R

9 M; a% B& U; C$ [7 C' O* J* }# B8 T1 F# N

该用户从未签到

3#
发表于 2021-12-23 13:45 | 只看该作者
{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}{:soso_e179:}
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-6-10 08:48 , Processed in 0.078125 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表