matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');1 n* \: h  y5 c* i3 S
   
2. 
   4 {* Z# c; v6 d3 C& s
3. hnb=fgetl(fphn);( z2 y* `$ G: C, q
   
4. 
   ( e\" X/ Q! D5 w/ x# S* [/ _6 `
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   - R9 ^, _* {- d. O
6. 
   ! v* O+ O; ?8 ^) U& t+ x
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   8 l! ^$ X* S3 n\" W3 c
2. 3 I4 e- h! Q$ y; I
   
3. xa=hnmap(:,[1]);1 l3 n' y) g1 H; C
   
4. 
   2 ?9 d7 b3 h\" }7 R) X4 A
5. ya=hnmap(:,[2]);
   . L% t, [4 o; i
6. 2 O% X# R- i9 H; j4 p
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   , a$ @( x* a3 D: w
2. ( V5 n/ N$ R( T
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);6 y6 C  h6 }1 K$ B\" P/ H5 X4 t
   
4. 
   2 c! s) [2 B; t. F: s
5. fclose(fp);
   5 J4 E/ ]& \2 n4 X
6. 
   & F, Y1 ?- ~) l& G. U; m7 j

复制代码

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

1. LL57=LL57';
   9 m+ t6 C% D2 o) p  w0 K
2. % q' @1 g2 J8 P# c
   
3. x=LL57(:,[3])/10;' |# d: q( c0 p, x
   
4. : z! t4 f; d- R$ {
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   ' |' Z6 Y3 Q% E' d4 d- A' Y: J
2. 3 S\" b) H- }* `1 |1 T
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   # s; M% ?\" L1 C9 S8 O2 B
4. 9 E, M8 C7 Z4 t7 p
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);; v. A0 h. V6 q  A& ^* A! G* C: t
   
6. 
   9 h$ W/ R  u/ C  f0 o$ ^( ~
7. fclose(fpy);

复制代码

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

1. yu97=yu97';  P  F0 G. D$ f8 m: L/ d& t
   
2. 
   3 A, I, M9 A* j( x
3. z=yu97(:,[2]);

复制代码

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

1. hold on;$ U5 A  h\" U  x
   
2. 
   4 r3 }6 s5 c* J, [* V$ Z7 E% h
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));
   : ?1 w7 y  `: C4 s  ~7 T
2. 
   4 O5 ^6 H2 l6 f
3. zi=griddata(x,y,z,xi,yi,'cubic');
   \" \1 Z) O8 z+ k- E
4. 
     T& s- D3 _% E+ m

复制代码

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

1. hold on;) K* s- J- P* E3 G
   
2. 5 |- d- B9 e0 i
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);* n* J: ?. c# f( y( {$ [
   
2. ; B5 H8 [& `7 L. T# S0 p! G  ~
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
, H+ E6 s  _- G9 |- J
6 T, {9 z, U3 x' X) R
% F9 k6 b5 e+ e- `; \

最后结果图如下：

# R% T" a9 J0 {, ]+ y9 d6 n2 y$ w

具体附件代码如下：

& v& S. L- H3 j  o! q

: Z6 o5 l' P  U7 k7 o- V
" m1 F# e, i7 J7 V