登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
xb=sdbj(:,1);
yb=14-sdbj(:,2);
plot(xb,yb,'k')
hold on;
xz=sdzd03(:,1);
yz=14-sdzd03(:,2);
plot(xz,yz,'k.')
text(xz(1),yz(1),['德州'],'FontSize',7,'FontWeight','normal')
%text(xz(2),yz(2),['惠民'],'FontSize',7,'FontWeight','normal')
%text(xz(3),yz(3),['乐陵'],'FontSize',7,'FontWeight','normal')
%text(xz(4),yz(4),['沾化'],'FontSize',7,'FontWeight','normal')
text(xz(2),yz(2),['东营'],'FontSize',7,'FontWeight','normal')
%text(xz(6),yz(6),['博兴'],'FontSize',7,'FontWeight','normal')
%text(xz(7),yz(7),['莱州'],'FontSize',7,'FontWeight','normal')
%text(xz(8),yz(8),['龙口'],'FontSize',7,'FontWeight','normal')
%text(xz(9),yz(9),['福山'],'FontSize',7,'FontWeight','normal')
%text(xz(10),yz(10),['文登'],'FontSize',7,'FontWeight','normal')
%text(xz(11),yz(11),['临清'],'FontSize',7,'FontWeight','normal')
%text(xz(12),yz(12),['冠县'],'FontSize',7,'FontWeight','normal')
%text(xz(13),yz(13),['禹城'],'FontSize',7,'FontWeight','normal')
%text(xz(14),yz(14),['平阴'],'FontSize',7,'FontWeight','normal')
text(xz(3),yz(3),['济南'],'FontSize',7,'FontWeight','normal')
text(xz(4),yz(4),['泰安'],'FontSize',7,'FontWeight','normal')
text(xz(5),yz(5),['淄博'],'FontSize',7,'FontWeight','normal')
%text(xz(18),yz(18),['寿光'],'FontSize',7,'FontWeight','normal')
%text(xz(19),yz(19),['沂源'],'FontSize',7,'FontWeight','normal')
%text(xz(20),yz(20),['平度'],'FontSize',7,'FontWeight','normal')
text(xz(6),yz(6),['潍坊'],'FontSize',7,'FontWeight','normal')
%text(xz(22),yz(22),['安丘'],'FontSize',7,'FontWeight','normal')
%text(xz(23),yz(23),['高密'],'FontSize',7,'FontWeight','normal')
%text(xz(24),yz(24),['诸城'],'FontSize',7,'FontWeight','normal')
%text(xz(25),yz(25),['莱阳'],'FontSize',7,'FontWeight','normal')
%text(xz(26),yz(26),['即墨'],'FontSize',7,'FontWeight','normal')
%text(xz(27),yz(27),['乳山'],'FontSize',7,'FontWeight','normal')
text(xz(7),yz(7),['荷泽'],'FontSize',7,'FontWeight','normal')
%text(xz(29),yz(29),['梁山'],'FontSize',7,'FontWeight','normal')
text(xz(8),yz(8),['济宁'],'FontSize',7,'FontWeight','normal')
%text(xz(31),yz(31),['泗水'],'FontSize',7,'FontWeight','normal')
%text(xz(32),yz(32),['新太'],'FontSize',7,'FontWeight','normal')
%text(xz(33),yz(33),['平邑'],'FontSize',7,'FontWeight','normal')
%text(xz(34),yz(34),['沂水'],'FontSize',7,'FontWeight','normal')
text(xz(9),yz(9),['临沂'],'FontSize',7,'FontWeight','normal')
%text(xz(36),yz(36),['胶南'],'FontSize',7,'FontWeight','normal')
text(xz(10),yz(10),['日照'],'FontSize',7,'FontWeight','normal')
%text(xz(38),yz(38),['单县'],'FontSize',7,'FontWeight','normal')
text(xz(11),yz(11),['枣庄'],'FontSize',7,'FontWeight','normal')
%text(xz(40),yz(40),['郯城'],'FontSize',7,'FontWeight','normal')
X1=0:0.5:20
Y1=0:0.5:14
[X2,Y2]=meshgrid(X1,Y1)
Z=griddata(xz,yz,upm10(:,3),X2,Y2)
[c,h]=contour(X2,Y2,Z,'k')
clabel(c,h)
hold on
%[c,h]=contour(X2,Y2,Z,[0 0.2 0.33 0.43 0.52 0.6 0.7 0.8],'k')
%clabel(c,h)
%hold on
%[c,h]=contour(X2,Y2,Z,[-0.1 -0.2 -0.33 -0.43 -0.52 -0.6],':k')
%clabel(c,h)
|