爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7500|回复: 3

[讨论] matlab画点去掉shp外的点

[复制链接]

新浪微博达人勋

发表于 2017-11-20 19:44:40 | 显示全部楼层 |阅读模式
2金钱
clc;clear all;close all;
data=load('data.txt');
%data(14252,3)=nan;%发现有一个错误的降水数据,值为7.4728e+04
lat=data(:,1);
lon=data(:,2);
rain=data(:,3);
clear data

level0=find(rain==0);
level1=find(rain>0&rain<2);
level2=find(rain>=2&rain<4);
level3=find(rain>=4&rain<6);
level4=find(rain>=6&rain<8);
level5=find(rain>=8&rain<10);
level6=find(rain>=10&rain<20);
level7=find(rain>=20&rain<50);
level8=find(rain>=50);
rain(level0)=nan;
rain(level1)=1;
rain(level2)=2;
rain(level3)=3;
rain(level4)=4;
rain(level5)=5;
rain(level6)=6;
rain(level7)=7;
rain(level8)=8;
cmap=[0,0,0;166,242,142;57,168,3;93,184,255;4,0,249;1,117,72;248,4,252;225,79,0;111,2,0]/255;
set(gcf,'position',[100 100 618 750])

scatter(lon,lat,10,rain)
set(gca,'fontsize',18,'linewidth',1.5,'box','on')
xlim([105 112])
ylim([31 40])
colormap(cmap)
cb=colorbar;
caxis([-0.5 8.5])
set(cb,'ytick',0:1:9,'yticklabel',{'0','0-2mm/h','2-4mm/h','4-6mm/h','6-8mm/h','8-10mm/h','10-20mm/h','20-50mm/h','>=50mm/h'});

map = shaperead('shaanxi.shp');
title('陕西省雷达降水','fontsize', 14);
hold on
bou1_4lx=[map(:).X];
bou1_4ly= [map(:).Y];
plot(bou1_4lx,bou1_4ly); % 画边界线
hold on
用matlab画的图,加载了shp地图文件,怎么讲shp外面的点去掉

rain.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-25 21:11:53 | 显示全部楼层
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12340
参考一下这个帖子~主要是这下面的代码:
  1. isin=inpolygon(lon,lat,yunnan.Lon,yunnan.Lat);
  2. Z(~isin)=NaN;
复制代码


密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2017-12-11 10:00:53 | 显示全部楼层
liuxiaoyue 发表于 2017-11-25 21:11
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12340
参考一下这个帖子~主要是这下面的代码:
...

我试过这个方法,因为data里面的数据特别多,有45W条左右,这个方法画的时候程序会卡死,不知道什么原因造成的
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-12-11 19:24:52 | 显示全部楼层
门中门 发表于 2017-12-11 10:00
我试过这个方法,因为data里面的数据特别多,有45W条左右,这个方法画的时候程序会卡死,不知道什么原因 ...

那可能就没办法了,我当时mask掉中国以外的区域也将近运行了几十秒钟,你最后解决了吗?
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表