matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
     s) x: i4 j) R( @
2. 6 l8 K9 J6 l  B* G
   
3. hnb=fgetl(fphn);
   8 r  t) E1 G& o3 f
4. 7 w* m$ v- B& N
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);- i& |! H3 e1 @: B# J+ g* z/ K0 ]! c
   
6. + X% u0 ]/ ?6 O) j+ V
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';- E$ X6 k7 k4 M; Q0 o
   
2. , ^; T. u6 |1 v3 S: f
   
3. xa=hnmap(:,[1]);! a- }3 m, |) h3 n# a0 O  s
   
4. 9 O% ]+ `\" ^4 G( Q# v+ w
   
5. ya=hnmap(:,[2]);
   4 O/ M  b/ X1 Y# V% Z5 k
6. 
   1 f# ^- e: L1 a* y

复制代码

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

1. fp=fopen('LATLON57.txt','r');% H( K' W% ^+ b7 c+ a- `( Y
   
2. + m4 ]# D7 L) M; ?- A0 T+ V
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);4 r% y8 `% A2 @8 u- o% \0 R\" c
   
4. ! j# H\" t/ P7 F0 g6 M9 k# P( x
   
5. fclose(fp);
   , H  j% |5 t# p- t' C2 R  G
6. % \( c1 ]' p2 k8 j
   

复制代码

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

1. LL57=LL57';- @& a! h+ y- h6 ?
   
2. 
   8 X2 v3 h$ ^! C9 O
3. x=LL57(:,[3])/10;  o\" c  P+ ~\" P/ R4 H) J
   
4. 
   8 K( V' M6 O, S# j3 j- m0 `
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');3 R6 q1 `\" \5 n\" \
   
2. 
   & B; I+ l& q, b2 g
3. ymd57=fscanf(fpy,'%d',[3,1]);8 `( B( O/ J9 y
   
4. # g7 E: f' r0 Z7 i( m4 j
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);# V; \9 s7 i7 d, x0 C
   
6. ; c: E1 F: S0 l
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   ) \- E8 u8 Z. G1 G
2. 
   \" M% `% ?8 a6 d; P& m
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   $ T2 e\" c- u( [7 M4 q. r
2. - s& I- N% D  @) H: t
   
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));* i9 x. y+ e4 Z% ^
   
2. 
   ' S- J7 N: b, U1 p
3. zi=griddata(x,y,z,xi,yi,'cubic');$ F/ f' l: S1 y! m$ P3 e4 P
   
4. . ?1 P% o4 c: u: Z6 G
   

复制代码

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

1. hold on;
   & [: ]3 _9 D6 P5 x8 `! X7 W5 d
2. 
   : u2 [  s, a6 M4 k\" v9 y% p
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   1 U1 o# g3 X% A& y9 a
2. 8 A: D( R9 D* L' X7 ?/ c' k: ~
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
0 J; g2 m9 x4 u0 x
3 c+ @  j' `' O
- h: ~- _. N* M1 A. I9 c+ L1 u

最后结果图如下：

5 ?; t6 Z2 w9 W* ?/ @: Z
" d6 D0 G+ Q7 Q$ G  _( _0 ~

具体附件代码如下：

) l1 m  i. f5 d0 y% P. q. p


& W4 s3 c0 g# S5 _8 T
8 V* N% t& G; R  z+ T, M9 b3 H7 u