matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   / X/ u2 [1 X- D& @+ m
2. 
   / X- O! \3 q, ^& [
3. hnb=fgetl(fphn);
   6 Q5 D( }9 f( ]8 f
4. 
   ! C1 _% u) J( J1 D\" U
5. hnmap=fscanf(fphn,'%f %f',[2,59]);6 ~0 b% ]9 ?, z+ }6 Q% y! }
   
6. - w4 G! K5 R3 `
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';: [6 v' R% {' d- S
   
2. : b! v; ]7 m7 X% \0 b8 M
   
3. xa=hnmap(:,[1]);
     ], ~& N' h: z- }+ S& t( U7 O2 o
4. 9 ^: d( R5 ?4 s' Q+ [: g
   
5. ya=hnmap(:,[2]);
   / t. ^! V8 C# Q7 L8 L
6. 
   : F1 e. ~2 E) k. C

复制代码

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

1. fp=fopen('LATLON57.txt','r');3 `# y& k' }' x
   
2. ) h9 N: y/ C; z; D( x* E! x
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);+ `9 G& q; `* `# C0 j% j2 C9 O
   
4. 
   1 j+ m; j/ y! k7 F* `5 v
5. fclose(fp);/ P3 ?% V. @- N! s\" J2 v! `
   
6. ) E\" Y+ Z( t# F\" |
   

复制代码

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

1. LL57=LL57';. T4 Z2 Z6 q3 l- [) {5 [# c' |) P
   
2. 
   7 `# A$ I3 j( X. M6 X
3. x=LL57(:,[3])/10;
   * a* m3 }. q- G  X  V3 c
4. 
   9 U0 H9 Z0 G; J$ l
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');\" M7 j4 P\" D/ S1 o$ W6 z! F
   
2. 
   # G5 U+ _& k\" r5 l& r3 ]
3. ymd57=fscanf(fpy,'%d',[3,1]);4 K) P1 E  i' s% b1 E6 c- O
   
4. & W2 p+ {  ~7 g* ^$ Z* S( M
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);0 H( Z\" y2 G7 N6 O+ m
   
6. 
   8 r; _( m+ `0 l! V, E+ U
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   # U7 P2 o, O- w5 c; T( s
2. \" a6 a* I' D, \' V: J
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;0 i8 {# p  P5 E% x# M
   
2. 3 A4 n2 I& i* ]\" z  y
   
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));
   4 f& t# S: V* ~, V
2. 
   . j- m; j( j  `7 A6 B7 U
3. zi=griddata(x,y,z,xi,yi,'cubic');
   9 K4 `\" g, f! u\" w) b$ z& a
4. % u! p( h# W4 p, \1 j& _
   

复制代码

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

1. hold on;
   # E0 X; c5 M' ^' Q0 I5 A' |\" X# F
2. 5 L1 h6 ?' X# l+ W! H  l9 t1 z
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   + Z: C$ h+ B$ e. G5 \9 @( K: l
2. 
   ( T' ~* L& F' I8 A7 d6 k
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
! ^* K  E, v4 L请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。

- g3 W/ p# h' G. {4 X0 @

最后结果图如下：

9 O, i" S8 f4 y. r5 w! o8 S* K

% i) e, x- ~. S, }& z

具体附件代码如下：

0 H6 X' H) e) h. S; |2 \4 L3 m9 Q


' r  N; X, ~% J5 l+ M, \1 Z7 r
) x" r! f( a$ ?0 U1 |* Z/ y" O. C