爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3942|回复: 3

[程序设计] matlab用分辨率较小的数据画大范围风场怎么设置风矢量的密度

[复制链接]

新浪微博达人勋

发表于 2018-4-5 17:01:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
在用分辨率较小的数据画大范围风场的时候,风场几乎铺满了整个屏幕。。。完全看不到箭头,如下图,请问有没有什么办法可以把风矢量密度设置的大一点?以下是我的程序:
clear all;clc;close all;
map=shaperead('D:\matlab\biyesheji\landareas.shp');
X=[map.X];Y=[map.Y];
plot(X,Y,'color',[0.7 0.7 0.7]')
hold on;
filename='D:\biyesheji\ECMWF\hgt7.nc';
lon=ncread(filename,'longitude')-180;
lat=ncread(filename,'latitude');
lev=ncread(filename,'level');
u1=ncread(filename,'u');
v1=ncread(filename,'v');
[x1,y1]=meshgrid(lon,lat);
u=u1(:,:,2,77);%u,v都是四维变量,分别为水平方向,垂直方向和时间
v=v1(:,:,2,77);
quiver(x1,y1,u',v','marker','.','ShowArrowHead','on','LineWidth',1.5,'color','b');
axis equal;
axis([90 180 0 80]);
hold on;
titlename=('7月19日06时风场图');
title(titlename,'fontname','宋体','fontsize',15,'fontweight','b');
str=get(gca,'xticklabel');
strtxt1=strcat(str,'°E');
set(gca,'xticklabel',strtxt1);
str=get(gca,'yticklabel');
strtxt2=strcat(str,'°N');
set(gca,'yticklabel',strtxt2);

捕获.PNG
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2018-4-5 17:22:21 | 显示全部楼层
本帖最后由 二爷名声在外 于 2018-4-5 17:25 编辑

插值成2.5*2.5的就会好很多的
PS这个是分辨率较大,不是较小……
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-4-5 18:31:07 | 显示全部楼层
嗷嗷,想要问下该怎么插值呢?不会插值。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-4-9 17:45:18 | 显示全部楼层
试试隔几个点读一个数呢,比如u=u1(1:count:size(lat),1:count:size(lon),2,77);你可以自己试试看count取多大比较合适~当然lat,lon,和v也要做相应的修改。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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