matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ! G3 F7 v$ V; u; y5 n1 `
2. % j$ [2 E- X4 B4 p7 B
   
3. hnb=fgetl(fphn);
   7 F7 x2 K\" F: f
4. 5 [! `/ p- y\" i
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   - s( \! v$ K8 c/ x% Z
6. 0 Y: G1 o' v% a2 r; Q
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   % q\" i6 ]7 l; U8 w8 _1 Q- O
2. 0 w5 y. _2 [& n- g' q# y/ f- n( r
   
3. xa=hnmap(:,[1]);2 `/ i; _/ g1 e  j& ~\" N8 b0 D1 k
   
4. 1 I. C, h; j\" I2 q  x7 V# w
   
5. ya=hnmap(:,[2]);( m8 ~\" `6 `& a4 R% ?8 B* A: w
   
6. 
   . ]4 W- {& E; z( [9 T! W

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   / a( ^/ `( b: A% c( g% G# K
2. 
   ! O! c) {% ]; x; c4 K
3. LL57=fscanf(fp,'%d %f %f',[3,97]);\" f8 U6 z% O  k+ s4 l# y) ?: ?4 C
   
4. 
     l) |) R! h8 M& r! {
5. fclose(fp);
   ' {1 @* |1 K8 s. D
6. 2 L! q1 Z0 ^( F) v
   

复制代码

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

1. LL57=LL57';
   \" v) R# f\" b/ r
2. ) s, ~# C+ P& c. {
   
3. x=LL57(:,[3])/10;
   ! F  W9 W9 _2 s, F
4. 1 d$ M. u7 F( l
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   1 h: t2 j' \8 q
2. 
   5 q6 b( V3 p  |2 [
3. ymd57=fscanf(fpy,'%d',[3,1]);/ L2 h) G: k2 ^- P6 @/ s
   
4. 
   0 \/ E% c\" v0 X3 f5 v( ]
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   , a* w- a& |& o2 U
6. . Y3 a( F. k& D1 |; P0 `\" z
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   0 K, j- ]* N5 Z( @$ k
2. 
   9 m2 s$ W\" X, p: u& V8 t5 _
3. z=yu97(:,[2]);

复制代码

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

1. hold on;. t  g7 t0 B8 \\" P% O8 ?
   
2. % D2 Y6 n  \5 b  n! q* C3 o/ G$ _
   
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));
   . a% x/ M8 w# u3 R6 B9 F* E! o
2. + K; L6 p1 j+ O4 P- f
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   0 E* b. Z0 d8 k/ w8 D2 I# W2 T+ `
4.   Y3 f- W, X\" a
   

复制代码

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

1. hold on;% p1 P+ J  u- `1 Z2 t
   
2. ! t9 h) _\" R0 X6 f: s. k3 ?
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
     r( ^3 r  p5 X3 F6 L, r. `6 `
2. 
   9 M- L6 _5 L7 _$ ?+ C\" F
3. hold off;

复制代码

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

9 f' U% Q9 G" T' D4 R! r

最后结果图如下：

& m0 L, c$ O, F; W5 A4 Y: ^

" a2 l0 h% }* j/ w' P

具体附件代码如下：

' }: {1 j; O; D
* x8 n: z. k; s! g% E( k
  r& P$ t* y0 ~, }
/ \5 L' q$ `$ P. d; U! ]0 O3 F1 C  A" }