
How to add wind direction to image
2 views (last 30 days)
Show older comments
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 Comments
Answers (1)
Bjorn Gustavsson
on 6 May 2016
Perhaps it is as simple as a swap of latitude and longitude coordinates?
HTH
See Also
Categories
Find more on Weather and Atmospheric Science in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!