matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');& c. a! \2 A- z4 l5 L% q( h. h
   
2. 
   0 d! E- a8 v\" E( n
3. hnb=fgetl(fphn);
   . g0 c4 q# j4 b  b/ R
4. & F& ~. \1 f& k  A/ y
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   : X\" S; x7 U: W: R
6. 
   ! [2 G3 q6 Y; n\" M: h& ^% f
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   ) k3 t/ x' m9 n8 ], I
2. 0 `  ^& t. j1 `7 m, |$ d
   
3. xa=hnmap(:,[1]);1 c& ^; Z- o* d, u9 @' q
   
4. 0 @# z7 P5 L; u2 X2 D4 v8 }7 f+ A
   
5. ya=hnmap(:,[2]);
   / H9 V4 _8 L/ R& l# Z8 i
6. % u( s: h5 M4 y2 u4 f2 O; k
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   ; U& ?! \/ I: @* Z% G
2. 8 I7 E7 _2 m' ^4 y9 z4 C- c% R- w# W) D
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   2 @+ t; a. }\" f8 M$ U
4. 
   # {' P+ O5 x- J% B# x8 G, C% @5 _
5. fclose(fp);9 d! _+ @3 E+ }- a; {
   
6. 8 R  {. q4 s* Q1 _* s
   

复制代码

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

1. LL57=LL57';
   & K6 L\" E1 V2 M5 W9 z5 S
2. 3 ^' e8 v- N5 U7 Y* K# S
   
3. x=LL57(:,[3])/10;
   ! S7 N, f0 N1 M0 H\" a$ F
4. / J% E- m8 D* O; t
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');1 ?% p4 ^3 Y6 [6 M
   
2.   Z7 Z1 |/ K0 _1 Y4 f+ W
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   & o; ^. t- s\" Q
4. 
   % P  H4 B0 i- s- x$ I& t' @; u
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);3 L( `3 I) Z( l
   
6. . G/ ^6 ^2 B4 v! h
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';2 A& N& @\" Y& b\" U
   
2. , @* b' `: F' Y* o
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   9 j) V  @3 W& O% S
2. 3 W* ^0 L* ~2 k, f9 L
   
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));0 o4 y+ F& M. i- S0 q) X
   
2. 
   ! j. b( }# M8 H; t\" D5 d
3. zi=griddata(x,y,z,xi,yi,'cubic');
   4 o; P3 F7 T* b8 e+ u7 E
4. , T3 J2 e5 W* U: ^5 ]
   

复制代码

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

1. hold on;4 y% u, x4 u/ ]0 D8 W' X
   
2. - d. M$ Y9 Q) R& S2 N+ J6 T( K, h% v
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   \" G- E0 k7 C! B) v7 O# ^; s
2. ) @& |; \1 q. q3 R
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
! w0 i2 ^7 G( J; K  n* D, w  O
5 Q7 h3 N6 m$ k1 n7 I5 R* o5 @
4 p/ k1 A* b. \& N- [

最后结果图如下：

$ s' f3 T: f% u0 T/ _4 V, r0 w

9 y3 R8 C. G" M& ]) K+ m
7 Y5 v4 ^- _9 V. O( K

具体附件代码如下：

/ q- S  h4 C& z4 S( D% C

6 W* x1 J2 I. ?, t: I

6 k$ b; t2 {8 d7 B2 g" r$ ]* @
/ w$ P6 }( D0 C$ Z& N