How to add wind direction to image

2 次查看(过去 30 天)
Sir I have done a plotting of soil moisture map for a region and plot its shapefile upon it.My code is
if true
% code
end
figure;
axesm('MapProjection','eqdcylin','maplatlimit',[26.0 38.0],'maplonlimit',[70.0 85.0],...'ParallelLabel','on','PlabelMeridian','west','MeridianLabel','on','MLabelParallel','south',...'FontSize',6,'FontWeight','bold','PLineLocation',20,'MLineLocation',20);
latitudes = 26:0.25:38;
longitudes = 70:0.25:85;
[latGrid,lonGrid]=meshgrat(latitudes,longitudes);
geoshow(latGrid,lonGrid,double(a),'DisplayType','texturemap');
a=ncread('C:\Users\Tanmoyee\Desktop\GFS\datafiles_RH_20Apr_03May16\gdas1.fnl0p25.2016042000.f00.grib2.thakur159756.nc','RH_P0_L103_GLL0');
colormap(flipud(hot(21)));
chandle = colorbar('Location','EastOutside','FontSize',6,'FontWeight','bold');
set(get(chandle,'ylabel'),'String','Rain Rate(mm/hr)','FontSize',10,'FontWeight','Bold');
set(chandle,'YTick',0:5);
% s = geoshape(shaperead('landareas','UseGeoCoords',true));
% geoshow(s,'DefaultFaceColor','none','DefaultEdgeColor','k');
% tightmap
hold on;
k=shaperead('C:\Users\Tanmoyee\Desktop\GFS\wind2\HP_UK.shp');
plot([k.X],[k.Y], 'k:','linewidth','8');
The image coming is attached.
Noe me adding wind vector upon it by quiver(U,V).Then the image is not falling on previous image.How can I solve it.
  5 个评论

请先登录,再进行评论。

回答(1 个)

Bjorn Gustavsson
Bjorn Gustavsson 2016-5-6
Perhaps it is as simple as a swap of latitude and longitude coordinates?
HTH

类别

Help CenterFile Exchange 中查找有关 Agriculture 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by