matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');8 T- R7 D( g9 G' J) q; x
   
2. 
   / b6 D1 E' k: y- D' ^
3. hnb=fgetl(fphn);
   1 |! N/ j/ n, e% `; a, Z
4. 
   : y4 r% j- `0 R) [
5. hnmap=fscanf(fphn,'%f %f',[2,59]);$ T3 Y% u4 S; k5 Y7 O
   
6. 
   . w5 d' J# r  Y
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   6 \; E6 y. k, m# }) ~6 f
2. + [9 u$ l\" y- q# u3 C+ O9 ~
   
3. xa=hnmap(:,[1]);
   ! ~. B+ k& Y, I0 u' B) E  ~
4. 
   8 `3 c2 b0 f) I7 _\" y\" T
5. ya=hnmap(:,[2]);4 t6 b: i8 F' {& R
   
6. 
   , r# @; ], w3 Y

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   ( b0 r! g) ^: I) @- C/ T
2.   U% e& @; v' L6 N
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);! H! N( ~: k4 \* J8 ?2 @. O- h2 {
   
4. 
   3 k& F3 K! z+ q4 i
5. fclose(fp);
   & q# Q) @4 D/ w/ J
6. \" N; V$ r! r5 v! Y. g& `
   

复制代码

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

1. LL57=LL57';
   ' ]; E4 r; X% G5 k
2. 
   5 ~. e% ~& ]4 o4 ]( R- Q
3. x=LL57(:,[3])/10;  V: X/ [3 t/ K
   
4. / K4 v( [+ @% I3 r1 U/ p8 {# N
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   ' b9 M* Q  J: x, {6 M
2. 
   6 j6 M0 m7 {. U4 z& O0 s' g
3. ymd57=fscanf(fpy,'%d',[3,1]);
   ( B6 M5 i6 ~( F6 F* v+ ^- d1 q
4. 
   # T3 v: O& j) w
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);' x0 U' K6 B  k* p. j- H0 @2 X
   
6. 7 w9 f4 V( U1 s0 f' _( M
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';  g/ V& W3 f& p
   
2. 2 A/ A6 O\" N7 p) Z: x) J% m
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;6 G  j0 P7 U$ ?! N3 P/ A- e
   
2. 
   : U  Z# d7 c7 C
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));  c$ r3 |! a: j% U7 H# R
   
2. 
   2 L\" ^/ z  ^8 H( T
3. zi=griddata(x,y,z,xi,yi,'cubic');
   / H9 Z, v' V9 q/ y% |% a5 d
4. & A# ^2 i% C; D- ^
   

复制代码

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

1. hold on;! R, r3 |- T/ \\" _( V$ v
   
2. \" `- u! t$ ]5 [8 }  C  t
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   * _- U\" S0 |7 Q
2. 
   5 U* ^3 P9 b- M4 `8 G1 J1 `! w2 X
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
' P4 U* ]. Y9 \# v( M请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
2 y: O/ E; l0 s: q
% l  |+ ~! K, d9 W, E  u
( ]0 }7 B" B: A; ]

最后结果图如下：

! Z0 L* L: Z; l. d6 ?$ g4 @8 Z

3 U5 |9 N, c3 W" r5 }' e% N# d

具体附件代码如下：

0 q3 \& x) }7 A+ ?& k+ c
# N; M0 r+ V8 j5 q7 F% v  G