How to add wind direction to image
    3 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 Mapping Toolbox 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!

