matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   # v. @, }2 s; j# j- F: \$ j
2. . b% y* |& z, L6 Q, m
   
3. hnb=fgetl(fphn);  C& \5 H0 ]% D# a  l, O
   
4. 
   - P& _- g* b9 F! J  P0 |, B* V
5. hnmap=fscanf(fphn,'%f %f',[2,59]);  h1 [9 m; ?: R7 ?* S/ S$ j9 L
   
6. 7 f& @+ m( Q( n/ h
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';: N7 i7 b$ k\" N
   
2. 2 d9 @8 l5 V\" {\" A, R, o& A
   
3. xa=hnmap(:,[1]);
   8 W  t\" d5 T\" u- Q/ k6 z8 K
4. 
   \" ?; e6 ^! [' P2 Z1 h& W
5. ya=hnmap(:,[2]);0 k3 x. z0 N: q\" {
   
6. 6 v! s6 S3 f4 ]5 Z5 d+ z
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   4 |' d! [! ?/ Y7 W6 |% D
2. 3 G7 d) i5 x6 N# z  F& N
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   ( ]' d' c' P7 {/ o  F
4. 
   . e+ O! ^7 E\" \; w- h7 _; p
5. fclose(fp);) B' T. t( Q! O6 Z0 n) }
   
6. 
   7 u/ u: n  O9 {- U9 b

复制代码

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

1. LL57=LL57';
   . J8 d* e\" ]+ w9 T- _& y
2. 
   9 `6 p- u6 i0 R( y( k
3. x=LL57(:,[3])/10;
   . I# O: @3 `6 k- G' J
4. / I! X( `% y' v) g6 W) c) |\" b
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');$ B0 Q# o1 ~' ]) e
   
2. 
   8 X  X5 N( t! {* R! l4 ?4 ?, P
3. ymd57=fscanf(fpy,'%d',[3,1]);1 |  ^% g. g2 e$ b& {8 X* H1 Q
   
4. % D  |6 ?# K5 }) t2 \1 ~! c
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);! X1 X3 s% L6 Y2 i: c0 f( V
   
6. $ |3 h5 v4 W\" x+ o0 }
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   . ], i! {( Z' Q4 b2 U6 I
2. ( ~. n3 B6 R3 [; {
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;. z8 h4 g3 B, ?) g/ `\" a) ~
   
2. 
   7 u& U, Q\" R\" s1 S) Y3 Z( k
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));
   , d) F7 D8 g( |% x$ ^
2. 0 B. y: t( K: m' i! [
   
3. zi=griddata(x,y,z,xi,yi,'cubic');; O' T: }2 c1 L\" q
   
4. 
   1 _4 \* ^\" W) O  f; k. s

复制代码

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

1. hold on;
   % \# M! e7 S# ^8 S
2. 
   $ M\" @. m3 ?0 N% |0 @
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   - `$ Y% [& Y* Q8 L
2. 
   # _; x0 R5 D) d1 e
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
' |: M: J* h! r- o+ b+ D
& X  I. x& M0 {' x8 Q
4 d8 C: @. I2 q, P4 `9 b0 x. e' _, K

最后结果图如下：

, `- \4 u+ w' x3 G$ D1 F6 h# W: j# Q

具体附件代码如下：

% j1 l* D3 v9 Y  E' T  \



- t  S1 S9 r$ d