上传一份克里金差值的源程序，希望有帮助

jiajian0512 · 发表于 2014-10-19 15:23:27

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

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 onplotm(lat,lon,'k.');colorbar;

hcy5877 · 发表于 2014-10-20 08:36:02

沙发么，呵呵，学习了

hcy5877 · 发表于 2014-10-20 08:36:32

谢谢，谢谢楼主

香草拿铁 · 发表于 2015-2-9 22:50:21

学习了
最近在找这个
谢谢分享

navy.blue · 发表于 2015-2-9 23:35:36

楼主有协同克里格插值的源程序吗？

xyj-20012318 · 发表于 2015-3-1 18:03:58

多谢了，学习一下

ilworld · 发表于 2015-3-1 19:06:24

学习学习看看~

jxboy99 · 发表于 2015-3-2 15:03:52

牛，学习，好东西。

鳄梨与梨 · 发表于 2017-9-24 17:17:08

请问压缩包里那么多的程序哪个是克里金插值的程序啊?

		自动登录	找回密码
密码			立即注册

[源程序] 上传一份克里金差值的源程序，希望有帮助

登录后查看更多精彩内容~

评分

相关帖子