How to add wind direction to image
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (1)
Bjorn Gustavsson
le 6 Mai 2016
Perhaps it is as simple as a swap of latitude and longitude coordinates?
HTH
1 commentaire
Voir également
Catégories
En savoir plus sur Mapping Toolbox dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!


