matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ' g, Y. W% v5 Q' x) |5 L4 d
2. 
   2 B4 G# u* b1 c' h& u- L$ f1 t( X
3. hnb=fgetl(fphn);\" l# N- G) j! @) U
   
4. ' e\" y9 V+ C$ j& f: }# R
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);, ], |5 N' k  Y, p4 f( w
   
6. . v5 J% x. O9 q4 G/ f
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';0 _4 g2 B, N- i1 o
   
2. ' e+ b9 @  u! e5 Y. ]/ Z
   
3. xa=hnmap(:,[1]);
   ! t- L- b- O5 o\" k* |+ b/ [& U+ s
4. 
   8 k' R9 C- K& \8 a- Y8 Q\" {( S  W
5. ya=hnmap(:,[2]);
   ! Q* ]8 W( o+ R0 `
6. 8 r( u) _7 {4 V7 O0 L- w# |
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');; C( W. b2 ~! k4 a( ?, W4 ?
   
2. 
   - x, A! Y- p/ t5 p+ c7 d7 H5 }
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   3 l9 o0 K9 k' N/ M/ [$ P
4. 
   ' m5 n- u+ H* F) M
5. fclose(fp);
   2 l: s1 j% X2 P9 O
6. % [6 i, B& U. z  o
   

复制代码

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

1. LL57=LL57';
   % d9 A\" D0 D& _7 |. n6 l, D
2. 
   ) _9 l* w& }. h0 A% a3 O
3. x=LL57(:,[3])/10;
     d1 R+ U8 W$ o1 D3 O. p. u\" `
4. 4 b$ s. M9 T  j: s8 R) }. s. W
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');$ s& [5 G4 ?5 U  x5 x  H) c1 {
   
2. ) e; d+ }3 i5 z: J. o) ]
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   % z) d+ q7 ~$ [3 m( O
4. 
   , M  D* g6 C( C
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);4 |7 y3 h+ q- l6 d+ L# m, Z2 ]
   
6. ) s5 S- j9 b! u& X$ R
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';! [. V- `$ `) S! V0 F
   
2. : L/ R0 |; c) Z) P. i
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;1 v+ T# P5 T2 _2 U
   
2. & C) A# G. L  s) f* [. P0 V$ m3 q
   
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));
   # L3 B; R4 Q1 `  t1 ~6 x! Q& q
2. 
   % L5 w. p+ r\" T0 i) L0 K
3. zi=griddata(x,y,z,xi,yi,'cubic');3 K* K$ q4 h4 X
   
4. ) S$ P& {, u  v! j) m
   

复制代码

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

1. hold on;  p, T( i4 M: h0 X/ b8 A
   
2. . r  H: d' f0 B  Z
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);* |; B$ d) w* h- e\" {( [: v
   
2. 
   6 L: A2 S- _/ h: u# Y$ J  \6 |\" N
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
/ C& R8 I! J/ b& U) {& S2 e请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
! a2 M# }7 w) W' q
" `$ K0 T8 k5 K2 k5 I
' u* N6 f1 s0 a# ?% K

最后结果图如下：

/ G1 f8 h* E; t- e7 {0 i

4 n+ N1 F; \5 r$ w2 b2 {
5 l2 N, q  M: t

具体附件代码如下：

: y* d; j" D  b9 x