How to add wind direction to image

2 Ansichten (letzte 30 Tage)
Tanmoyee Bhattacharya
Tanmoyee Bhattacharya am 6 Mai 2016
Kommentiert: Chad Greene am 6 Mai 2016
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 Kommentare
Tanmoyee Bhattacharya
Tanmoyee Bhattacharya am 6 Mai 2016

Melden Sie sich an, um zu kommentieren.

Antworten (1)

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

Kategorien

Mehr zu Agriculture finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by