matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');/ k; g\" G! k5 s9 r9 t2 o
   
2. ! B* z% B! Y# w6 o3 x
   
3. hnb=fgetl(fphn);
   9 w$ ~9 i9 @+ ^. t: }% k6 D2 z
4. 
   7 k9 ?) m7 d- {; t( T5 [
5. hnmap=fscanf(fphn,'%f %f',[2,59]);4 X# D* t- K% Z( }3 H. S. s5 m0 T
   
6. ) i! K/ M. Y% Q; r7 f
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   6 I- u2 a0 j; r: V+ K- `# z
2. 
   7 s' N: F& C* o1 g. W
3. xa=hnmap(:,[1]);
   / M* r9 V; {0 O
4. # Q( m5 t% A. j; F9 k; w: f5 k
   
5. ya=hnmap(:,[2]);$ }+ O  P) k3 G. p0 @
   
6. 
   . O( p: _. L3 E: k

复制代码

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

1. fp=fopen('LATLON57.txt','r');9 w+ F4 R$ }: ^0 D3 n0 d( }! R+ f
   
2. ! i9 h3 f# N% A, v* r
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);- Y' T* r4 J/ h: ^( Y
   
4. 
   % O: h- |3 ?; J1 a
5. fclose(fp);
   ; |. k/ G6 i! H& B
6. . |& ]0 R$ @1 O) x) e
   

复制代码

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

1. LL57=LL57';
   % B$ y+ A\" z* x+ ?3 h8 e! O3 p
2. 
   7 R$ S, @2 E% D
3. x=LL57(:,[3])/10;6 e2 j# u5 m5 A3 [* @6 ~. ?
   
4. 4 Y\" H% k0 H. q. z
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   3 H5 G, T/ Z$ T/ o: H
2. 
   * R5 q. [0 R\" X/ _5 O9 q
3. ymd57=fscanf(fpy,'%d',[3,1]);
   , ~8 j3 ?, p! I
4. 
   ; f, j2 @! ~0 r! P, k$ h, b8 K
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   & ?( Z0 V9 }3 Y5 j4 f
6. 
   . ?- \7 s: x4 n% |0 Q
7. fclose(fpy);

复制代码

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

1. yu97=yu97';0 B: a) b8 G4 e$ Y
   
2. 
   $ L# ^% J- Q' L* |
3. z=yu97(:,[2]);

复制代码

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

1. hold on;3 f5 p, F( R( ?! U1 s8 D/ O, Y4 K
   
2. ( ^% K- F  I4 A4 f# M
   
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));2 q\" G: ?) p$ e2 {4 z
   
2. 4 q2 G3 o( ^: M- ]0 h' K
   
3. zi=griddata(x,y,z,xi,yi,'cubic');5 F- N7 Y/ u/ t; X\" H' X7 C: s
   
4. 
   # r/ \. C$ }3 x* q6 _3 p( h

复制代码

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

1. hold on;
   - i+ M5 o6 j0 [8 h$ _0 ?
2. . s5 ~+ I' {) c2 Z4 `
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   * Z/ G) h2 ^. i6 n; V% D- r$ n
2. 
   ' q\" Z/ S' x0 {9 K0 B2 ?. D
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
* o, H& C- ?- d% O! A& g  o8 W7 j  }

& |2 @- B# V$ A# a4 q7 A( F; V, I

最后结果图如下：

' c' `+ u. Z( P! u! J, e

3 x8 [! K) `* U

具体附件代码如下：

" `. V: c& s( m



2 o% s: T* U2 W% m& z
" I6 C* |! Z3 q7 v! H; A2 q: e; A