|
|
本帖最后由 ononsiiii197 于 2021-12-23 13:30 编辑 ' T# q! r8 {! ?! o0 U
- ]+ {8 T# t; Q. i) i
: ~5 }7 e8 f. A0 V9 O8 P
- 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;
/ h7 V$ c+ \% J. m9 Z
3 i2 r( @% p0 w. J$ H+ ~8 T
4 {; F* U& H L$ W; r$ Z' u: j1 T
|
|