How to add wind direction to image
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
Tanmoyee Bhattacharya
on 6 May 2016
Tanmoyee Bhattacharya
on 6 May 2016
Tanmoyee Bhattacharya
on 6 May 2016
Tanmoyee Bhattacharya
on 6 May 2016
Tanmoyee Bhattacharya
on 6 May 2016
Answers (1)
Bjorn Gustavsson
on 6 May 2016
0 votes
Perhaps it is as simple as a swap of latitude and longitude coordinates?
HTH
1 Comment
Chad Greene
on 6 May 2016
That's what it looks like to me.
Categories
Find more on Google Earth 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!

