- 积分
- 543
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-11-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
clear
clc
data=load('C:\Users\HP\Desktop\canshu_data\wind.txt');
lat=data(:,2);
lon=data(:,3);
wind10=data(:,8);
red=find(wind10>18);
red_lon=lon(red);
red_lat=lat(red);
scatter(red_lon,red_lat,'.','r');
hold on
yellow=find(wind10<=18 & wind10>16);
yellow_lon=lon(yellow);
yellow_lat=lat(yellow);
scatter(yellow_lon,yellow_lat,'.','y')
hold on
green=find(wind10<=16 & wind10>14);
green_lon=lon(green);
green_lat=lat(green);
scatter(green_lon,green_lat,'.','g')
hold on
blue=find(wind10<=14);
blue_lon=lon(blue);
blue_lat=lat(blue);
scatter(blue_lon,blue_lat,'.','b')
hold on
title('10年一遇极端风速','Fontsize',15);
h=legend('>18','16-18','14-16','<=14','location','SouthWest')
set(h,'Fontsize',15);
hold on;
shengjie=shaperead('D:\shengjie.shp');
x=[shengjie.X];
y=[shengjie.Y];
m_proj('Equidistant Cylindrical','lon',[114,120],'lat',[29,35]);
set(gca,'fontsize',16,'linewidth',2)
hold on;
m_plot(red_lon,red_lat,'ro','MarkerFaceColor','r')
m_plot(blue_lon,blue_lat,'bo','MarkerFaceColor','b')
m_plot(green_lon,green_lat,'go','MarkerFaceColor','g')
m_plot(yellow_lon,yellow_lat,'yo','MarkerFaceColor','y')
m_plot(x,y,'k','linewidth',2.5)
m_grid
|
-
-
|