matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   5 K* g- ~) p  ]1 L
2. 5 Y: @! R6 i- n4 @
   
3. hnb=fgetl(fphn);1 S6 k; E9 C6 D! S5 F; V
   
4. 7 T. l\" j& }9 k# c
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);; Y( i) r, T; U  Q0 u
   
6. 
   * s7 H  l; F  T2 g3 {- Q  Q
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   & u1 e# z) S1 h
2. . t% T- N) D$ U- u5 Q2 j. C
   
3. xa=hnmap(:,[1]);
   \" z9 T- T1 c$ t* l# F\" _5 U  F. S9 \# M
4. 
   3 ^; W' ]+ {3 |
5. ya=hnmap(:,[2]);
   6 |* @# E$ k- O4 m2 W  w0 k
6. 
   # `1 p0 V. ~  P% j

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   ' \0 j6 A6 U1 c& t/ j
2. 
     j+ m2 i  u5 o' t: g0 q
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   9 w3 }% ?$ {: C8 p8 \
4. . o( F, ]8 z7 _( j8 W0 [) j
   
5. fclose(fp);
   2 R2 Y7 w) b* r: g
6. 
     ^5 ]7 Y+ H! U

复制代码

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

1. LL57=LL57';
   ( K$ p  V/ P5 \) P3 O3 _: n' a
2. 
   6 z* e: }7 G3 S\" V8 v
3. x=LL57(:,[3])/10;1 |\" O$ c! x: [9 ~# Z
   
4. 
   # m1 Q3 d1 ]! p8 x
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');& k# a7 D3 W5 |5 R\" j' B
   
2. - q+ m% W  M4 k\" @4 `
   
3. ymd57=fscanf(fpy,'%d',[3,1]);; H: t# f1 i- J( V  W) x3 U  k2 D
   
4. 9 A9 N( B& [3 [. c* d6 d
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   0 o. O& V7 S4 T\" Y/ a
6. * {, j% O% n7 w) b+ s\" X% B8 @0 [
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   9 d; i. C) b! G; Z+ O
2. 
   4 a& ~% \3 W* U' q* \5 p2 P7 B' I
3. z=yu97(:,[2]);

复制代码

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

1. hold on;% r0 B  `/ x4 y' M
   
2. 5 V- h# j, a7 J% q* g* N- r
   
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));5 x1 p0 l3 ]) K2 n2 x# s4 P1 Z
   
2. 
   , t9 X2 `* C5 h& ^5 e
3. zi=griddata(x,y,z,xi,yi,'cubic');5 ~% q: Q0 ~/ a' w5 \! W
   
4. , P. Z. N0 ^- `1 n1 ^
   

复制代码

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

1. hold on;
   . O7 X# N; E1 z5 C) H0 B
2. ; @- I; B6 ?! @' |
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   ' z9 S+ H4 Z, M1 R. N7 Y
2. 4 {8 x# A% y( Q# I1 r
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
& G0 I# g* x+ A! f" s- K9 y请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。

: N. q6 }5 ]5 _  J

最后结果图如下：

5 H: X( x: w$ U1 \/ O4 T1 s
" A: _0 `9 }* V6 d* H2 W! T  W3 D

具体附件代码如下：

! p- A7 {$ O' E0 m3 }



) g* o7 b9 [3 ^' I( U
; [4 x$ t+ s$ Z' J! U; h- o