quiverm creates wrong vectors?

Question

Sandra 2012-6-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41860-quiverm-creates-wrong-vectors

评论： Chad Greene 2014-7-28

Hello,

I am trying to plot wind vectors with quiverm. For some reason, and I would really like to know why, the vectors I get point in the wrong direction, distorted in the x and y direction, compared to when I just plot them with quiver.

This is what I have and use:

ywindlat  meshgrid (17x69) (latitude increasing by rows)
xwindlon meshgrid (17x69) (longitude increasing by columns)
u wind (17x69)
v wind (17x69)

then I use the following commands:

   load coast
     axesm('mercator', 'MapLatLimit' ,[-20 20], 'MapLonLimit' ,[120 -70],...
     'ParallelLabel', 'on', 'MeridianLabel', 'on', 'MLabelParallel', 'south', 'Grid', 'on', ...
     'PLabelLocation', 10, 'MLabelLocation', 20, 'MLineLocation', 20,'PLineLocation', 10)
    plotm(lat,long, 'Color',grey,'linewidth',1)
    quiverm(ywindlat,xwindlon, wind_u, wind_v,'black',0.2 )

that gives me different vectors than if I just plotted:

quiver(flipud(wind_u), flipud(wind_v))

What am I missing?

Any help is appreciated a lot

Sandra

5 个评论
显示 3更早的评论隐藏 3更早的评论

Sandra 2012-6-25

Hi Walter,

I think I have a slight idea why it is distorting them. And it's somewhat embarrassing :-) My U and V input vectors are wind speeds in m/s and not in in lat/long degree units as requested by the quiverm function. Obviously this is not going to work out...(ahhh). Do you know of a way that I can use the quiverm but with actual wind values?

Have a great day

Sandra

Walter Roberson 2012-6-26

I have not worked with the Mapping toolkit. Looking around the documentation, I'm thinking that perhaps you need to convert your velocity vectors using reckon() as discussed near the bottom of http://www.mathworks.com/help/toolbox/map/bragdwh.html . Divide your meters/s by 1000 to get the values in kilometers. But to get the arclen parameter for the distance I'm thinking you probably have to find the magnitude of your wind-speed vector... I guess that would be sqrt(u.^2 + v.^2). And I guess the az would be atan2(u,v). Might be easier to do a cart2pol() to calculate the magnitudes and angles.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Chad Greene 2014-5-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41860-quiverm-creates-wrong-vectors#answer_137894

I recently ran into the same problem. The ncquiverref function scales the vectors correctly, and correctly identifies u as a zonal component of wind and v as the meridional component.