matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   / N* g6 _& A1 K* Y5 x
2. 7 d7 j( ?$ i# `' A5 b! m5 y
   
3. hnb=fgetl(fphn);
   , Z' E& R% ]: o+ i
4. ' F7 f* K2 g/ E  [6 d- @
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);- K* D/ f* S3 t
   
6. 4 ^7 n4 T- g/ e2 }/ q
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';9 k, l- @6 M0 W: ], c
   
2. 0 C6 q+ p5 u2 i6 }: l/ `/ a4 O: r
   
3. xa=hnmap(:,[1]);3 |) D5 O8 m/ Q% {7 N* h7 K
   
4. : `* L8 X' _$ D7 `$ O7 P
   
5. ya=hnmap(:,[2]);; [* h; ^: }\" ?! O, Q
   
6. 
   2 U$ b; {6 P( v- r/ J

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   ( T/ X- o( C& j7 O2 u
2. 6 l/ n( J* b2 n
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   9 \& q9 `\" [' {1 H
4. 
   ) U/ l: q' V7 Q! x
5. fclose(fp);
   / C0 V' Z, C7 f* H: w* W. |2 u
6. 
   8 `; @3 Z4 l' G* \) G

复制代码

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

1. LL57=LL57';
   ! S' ~% l* g) D2 a1 f! [\" i
2. 
   - P3 R+ J3 @) K- V
3. x=LL57(:,[3])/10;
   / ^* f) c  m- F0 b/ {( L* a# K
4. % K& v7 `7 t: F& H
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');. [( m! t+ K1 O: x
   
2. 7 ^0 O' r9 ], z6 Q( B
   
3. ymd57=fscanf(fpy,'%d',[3,1]);! F8 s6 l/ z% g) {4 Y. K
   
4. 
   6 j+ u+ _\" @0 N! ^
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   ' S! g8 u9 y* Z6 B6 x+ E
6. 5 l* Q; o2 j1 J8 ]8 C9 h
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';& l* W# H' n0 O7 |3 ~
   
2. ( z' v* u+ w) f
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   2 \, E8 W; U$ g8 A( S
2. 4 u& ]6 s! g* d* _1 I) U
   
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));. X. Z\" D+ s' T( O
   
2. 
   / `; P8 M0 q$ b' e6 n9 d. n! v; R
3. zi=griddata(x,y,z,xi,yi,'cubic');+ ~# Q  t\" F/ i7 Y$ V
   
4. # x; n* f8 j/ L% v, D
   

复制代码

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

1. hold on;
   ) {1 M6 f. c/ p; ^* p
2. 
   - \% z7 B, S9 }, ?9 O8 X
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);$ ^\" r+ Z% h% e
   
2. 
   3 C/ I. t0 i9 `: U' ^
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
3 j7 C! D6 L, I" p& G. t$ r
0 |* T  o1 z) v9 h! S; W5 m
) y1 V+ q, V% z  {( u6 \

最后结果图如下：

9 _- B3 R& g( |  R7 |
  J0 _! |" }8 o! ?

具体附件代码如下：