matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

这两行清空MATLAB的工作区，并关闭所有图形窗口。

1. fphn=fopen('hunan.txt','r');
   6 q5 T8 m) q$ v: \) o6 m4 \' j
2. 
   & q' O* Y0 n\" o4 K8 @) e
3. hnb=fgetl(fphn);
   ; a\" u+ j9 E/ `) `, z! h
4. ) z0 o+ V2 @6 N
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);/ b& J% u5 Z+ a
   
6. 
   % g6 B! `# k) P# F- y) Y$ `5 e
7. fclose(fphn);

复制代码

这部分打开名为 'hunan.txt' 的文件，读取文件中的湖南省的经纬度数据。fgetl 用于读取文件的第一行，而 fscanf 用于读取文件的剩余部分。hnmap 存储了湖南省59个位置的经纬度，每一列对应一个位置。文件的每一行应包含两个浮点数。

1. hnmap=hnmap';
   5 `8 }1 z% Z* K0 Y9 j6 a
2. 7 C2 k& h5 E\" t: R: D% d: M: e
   
3. xa=hnmap(:,[1]);
   0 i4 Y$ F9 g. m2 s
4. 
   ! f+ ^' n# {' a3 n
5. ya=hnmap(:,[2]);
   $ g# r% f8 q0 X1 y( J5 m$ C
6. : ^/ ~: x7 v2 |\" h
   

复制代码

这里对 hnmap 进行了转置，然后将经度和纬度分别存储在 xa 和 ya 中。

1. fp=fopen('LATLON57.txt','r');
   % N3 ~0 V9 g* P
2. 
   : A\" N! `& S- h. B+ e
3. LL57=fscanf(fp,'%d %f %f',[3,97]);\" w. z  @- O% G- L
   
4.   u# f) S: r* n7 R* P
   
5. fclose(fp);
   7 y\" V8 W, C& o; t
6. 8 R1 S6 G! i3 h5 {/ G- B0 [7 ~6 Y4 A; A
   

复制代码

这部分打开名为 'LATLON57.txt' 的文件，读取文件中的湖南省97个县的编号和经纬度数据。LL57 存储了97个县的相关信息，每一列对应一个县。文件的每一行应包含一个整数和两个浮点数。

1. LL57=LL57';
   & d  f: r& s8 M\" n5 i
2. 
   - [* m  i- Z\" L. k4 Z8 d
3. x=LL57(:,[3])/10;, O7 r\" i3 Y( E
   
4. 
   & P; q* `% K& H1 u
5. y=LL57(:,[2])/10;

复制代码

这里对 LL57 进行了转置，然后将经度和纬度分别存储在 x 和 y 中。在这里，经度和纬度都除以10，可能是由于数据单位不同而进行的缩放。

1. fpy=fopen('etw00100.txt','r');) `0 G7 q$ y$ k0 b, ~
   
2. 
     R( s, w. A; @$ |% h
3. ymd57=fscanf(fpy,'%d',[3,1]);
   ( {% f& B' O! k* K
4. 
   9 V\" `0 q6 r8 D  N
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);2 |( Q- ^  o; c7 C
   
6. & P, p. O7 E8 m3 Y
   
7. fclose(fpy);

复制代码

这部分打开名为 'etw00100.txt' 的文件，读取文件中的湖南省97个县的温度数据。yu97 存储了97个县的温度信息，每一列对应一个县。文件的每一行应包含一个整数和两个浮点数。

1. yu97=yu97';. M6 r/ x- s2 I, {% d
   
2. 
   3 S2 B, d. ]0 T4 i
3. z=yu97(:,[2]);

复制代码

这里对 yu97 进行了转置，然后将温度数据存储在 z 中。

1. hold on;; j0 U: M) O0 ^* E
   
2. 
   . M4 D. `4 Z, u; y# F4 m/ ?
3. plot(xa,ya,'.','markersize',5,'color','red');

复制代码

这部分开始绘制图形，用红色的小点表示湖南省的边界。

1. plot(x,y,'.','markersize',6);

复制代码

这部分继续绘制图形，用蓝色的小点表示湖南省97个县的位置。

1. [xi,yi]=meshgrid(linspace(min(x),max(x),25),linspace(min(y),max(y),25));
   . h( Y: ]0 M/ k. t
2. ( I. k, G! T% v* X8 X- T+ Z
   
3. zi=griddata(x,y,z,xi,yi,'cubic');/ c0 E9 l( |$ n9 ]
   
4. 4 C  t* ~+ p; y
   

复制代码

这里创建了一个网格，然后使用 griddata 进行三维数据插值，得到 zi 表示的湖南省97个县的温度分布。linspace 用于生成指定范围内的等间隔数据。

1. hold on;
   7 y/ J\" {9 b5 T& ~) a2 `3 A! m' v
2. 
   & Z5 v- z$ O, C9 n
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

这部分绘制等值线图，表示湖南省97个县的温度等值线。contour 用于画等值线，xi 和 yi 是插值后的网格，zi 是温度数据，'b-' 表示蓝色实线。

1. clabel(c,h);
   : ^6 P0 V7 R! Z/ O* a) ]
2. $ l4 j. X6 q3 l! o0 H, a6 A\" C
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
) n1 g9 [4 ]0 Z) y7 r  W0 w请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。


, @# H) |  O* V) C/ P1 q& [

最后结果图如下：

6 r! v2 \# P- Y) u3 A0 b8 U% E

具体附件代码如下：

" f4 _9 u% s9 S  p: ]$ j2 O0 f
" Q! }) g" M9 r# _# C: L

5 b1 F+ J+ C3 A1 E7 m