|
本帖最后由 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
|
|