画等温线

cdq872

clear;close all;
fphn=fopen('hunan_zh.txt','r');
hnb=fgetl(fphn);
hnmap=fscanf(fphn,'%f %f',[2,1449]); % It has 59 rows now.湖南省界经纬度
fclose(fphn);
hnmap=hnmap';
xa=hnmap(:,[1]);
ya=hnmap(:,[2]);
fp=fopen('LATLON57.txt','r');
LL57=fscanf(fp,'%d %f %f',[3,97]); % It has 97 rows now.湖南省97县名称号码，经纬度
fclose(fp);
LL57=LL57';
x=LL57(:,[3])/10;
y=LL57(:,[2])/10;

fpy=fopen('etw00100.txt','r');
%ymd57=fscanf(fpy,'%d',[3,1]);%实在不懂这句是什么意思，并且后面也没有用到这个变量
yu97=fscanf(fpy,'%d %f %f',[3,97]); % It has 97 rows now.湖南省97县温度
fclose(fpy);
yu97=yu97';
z=yu97(:,[2]);%湖南省97县温度
hold on;
%plot(xa,ya,'.','markersize',5,'color','red');%湖南省界
plot(xa,ya);
hold on;
plot(x,y,'.','markersize',6);%湖南省97县位置
[xi,yi]=meshgrid(linspace(min(x),max(x),25),linspace(min(y),max(y),25)); %meshigrid（x，y）的作用是产生一个以向量x为行，向量y为列的矩阵，从数组a和b产生网格
zi=griddata(x,y,z,xi,yi,'cubic');%湖南省97县温度立方插值，以三角形为基础的三次方程内插
hold on;
[c,h]=contour(xi,yi,zi,'b-');%温度等值线，注意是平面的等值线
%clabel(c,h);  %C——标签矩阵、h——等高线句柄  %标签到等高线图像中
clabel(c,h,'FontSize',15,'Color','r','Rotation',0); %指定字体和颜色
hold off;


地图数据自己弄啊。

mofangbao

谢谢分享源程序，多多交流哈

无径之林

很不错！

tintin

不错

zxm456

如何沿省边界裁掉等值线。谢谢啦

njzqxt

楼主真是厉害。

哈斯巴人

好东西，给力

七夕路

实话说，没看懂，先收藏。谢谢楼主分享。楼主辛苦

cui111yy

谢谢分享源程序，多多交流哈

莫莫泡沫

解决问题了，赞