matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');; I6 p! G! c9 s6 y( W/ J
   
2. 
   ; |2 t, I' r  ^5 R2 W7 x: r/ r
3. hnb=fgetl(fphn);
   1 R$ c7 X2 i. M3 U4 {' ^4 \: t
4. 7 N\" U2 R7 F6 D5 m7 j9 E3 h$ [
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);, W' ^+ p  Y2 C5 A4 P
   
6. 
   * k6 d! y4 G  \: n! c
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';* _- ]' W$ U/ Y- }: H2 U
   
2. & R& I( b% o\" e% h  L6 v' `3 a
   
3. xa=hnmap(:,[1]);) G' \: }0 o\" T' u0 X  P4 Y
   
4. 
   2 w4 @! \# `( y& S' d; k
5. ya=hnmap(:,[2]);
   2 Z: Q; C) @5 t5 V+ K
6. 
   2 m# D& Z* j* |; N9 s2 s\" O+ o

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   : U% J2 R+ Z* T2 o/ b& x4 S! H7 b
2. ' U0 ~/ ^$ ^, C! X. \# ?; x
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   ' _0 a/ ^# A\" |' `! E
4. 
   9 b+ O0 M( e: c7 C, y
5. fclose(fp);! a% S; ~' Z2 u
   
6. 
   0 I! q+ J- ~$ e2 s6 d& Q: _! N! u

复制代码

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

1. LL57=LL57';$ q% Y; S4 f9 H) V) V
   
2. 3 o\" ]. h9 N( @4 x
   
3. x=LL57(:,[3])/10;& i) b+ u; J; K* Z6 T, H: N
   
4. : c. z8 c- E) e2 R7 q2 o
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   * D9 L& {5 L# V6 \! x- F\" ?
2. 7 w2 _' C3 @  ?; j* E9 Y
   
3. ymd57=fscanf(fpy,'%d',[3,1]);, u$ `4 n6 q) I0 ]
   
4. 
   - ^: l! m1 X+ C3 ~, I5 j
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   9 z6 D- |4 ?: O! c
6. 
   3 c% h* a\" q+ n2 w& T8 \5 R- R
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   , G2 K6 _\" v& Q' x1 x3 d0 w. n
2. ! k, O( K. A/ K. c& r# j
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   + J# O5 O3 N) P; t2 c) N
2. 6 K$ [, L4 N+ C/ y* X# @- k
   
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));- N0 N4 d! F- B' w/ f# |
   
2. ' C2 `: o) v\" c' f0 w
   
3. zi=griddata(x,y,z,xi,yi,'cubic');  `* S/ x! Y1 [$ L5 M
   
4. 
   , v# \# A- V& T5 H7 C0 P# i

复制代码

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

1. hold on;# Z; P* Q% ]. B  }' I
   
2. # X' A8 P2 Q% G5 L
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);: q1 q3 z; c9 s! }8 L3 j4 w
   
2. 
   8 M4 M* x+ Z+ b+ m+ Z
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
6 P& y$ f) r7 O, _0 F( `& w请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。

/ R5 Q% t+ z6 W+ j' Q

最后结果图如下：

2 k% G5 R$ t6 t- w) J

具体附件代码如下：

  c0 k  o1 Z0 B0 Z2 ^) J$ _7 |; ?