matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');& c/ V/ D+ j4 k! v3 u: g2 V
   
2. 
   . E; R2 `* b9 X1 F1 @\" a
3. hnb=fgetl(fphn);# d5 w4 K$ R( V$ L
   
4. 
   ( T( r  L! U2 P0 k
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   1 I+ n2 A4 ?* \: }
6. 
   $ j: k( H! E% B7 Z- b
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   . h5 T* E/ u, ]
2.   s. t/ i, I% ?5 I' ~\" k) k
   
3. xa=hnmap(:,[1]);
   # C0 J# U& R+ j- a( Z% J6 A1 C
4. 
   ( Q& A+ C8 v2 E) [2 H+ o4 k
5. ya=hnmap(:,[2]);
   - M1 v6 P+ W3 V
6. 
   0 o2 @9 q$ w7 F: e% z, G

复制代码

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

1. fp=fopen('LATLON57.txt','r');8 g% o\" h& T9 a$ Z. L+ o  U
   
2. ( M  d2 t5 I5 \* t& y2 v/ D
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   + G9 V\" d/ b) \5 s\" F! q5 d
4. - {$ m& \' D0 O  D: a* D6 k  S/ e
   
5. fclose(fp);
   ; o$ p7 ?- J8 [) |3 k: s
6. , p+ c0 j5 C+ Z6 G9 k7 M  |: P
   

复制代码

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

1. LL57=LL57';2 j5 @7 E  c% b  i3 H& j! }- l
   
2. 
   7 J* ?, B8 y; {8 I
3. x=LL57(:,[3])/10;
   \" {* s! ?/ D$ l  M
4. ; r% N/ Y. P7 k3 y6 J7 A. W( k
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');( M$ s$ Z$ p7 `% P. v6 H! E7 L
   
2. 
   ) ?\" Z7 f( n9 N
3. ymd57=fscanf(fpy,'%d',[3,1]);
   ( V& A2 x9 \9 r4 x2 `
4. , P  R8 c4 n. o8 O\" ?
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);( S' \6 J$ l5 Y5 p- e/ U7 _
   
6. 
   . p8 @; b; V* t4 v5 l\" R
7. fclose(fpy);

复制代码

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

1. yu97=yu97';6 f  N; B) E4 ~
   
2. 
   5 V4 k2 x7 ~2 w
3. z=yu97(:,[2]);

复制代码

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

1. hold on;( h% n1 {: ^6 w
   
2. 
   - N\" r# `7 X$ j) J# R
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));) U: }3 t+ M: L& O2 l
   
2. ( K9 j3 A1 i8 R- V- j
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
     w/ g4 C0 n2 N9 Q
4. 
   9 |4 R' H6 o5 m, A8 H+ |

复制代码

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

1. hold on;
   2 v8 f% _  m$ q7 b4 ~
2. 
   ; h3 `& r. o! v. g' n0 ^0 z
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);3 l- J: T  j2 h+ V$ i
   
2. 
   / G/ j8 a$ _\" Y5 P
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
6 x" P2 `6 |* \2 j$ Q5 R请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。

. R. r- z8 k: K1 V

最后结果图如下：

! ~6 d  F# e3 N5 l2 z

! I! ?2 {% p( v, p8 p

具体附件代码如下：

5 N1 \9 y, U: C0 S
2 J8 [. P/ z- `  L+ [2 a
8 y( y$ M- P$ c3 f. t- b
7 e, G9 k3 V( X* z9 s( M
5 M; e1 N9 n1 [- L: G