- 积分
 - 544
 
	- 贡献
 -  
 
	- 精华
 
	- 在线时间
 -  小时
 
	- 注册时间
 - 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 
 
 
 |   
- 
加图层前 
 
 
 
- 
加图层后 
 
 
 
 
 
 
 
 |