Question

我正在尝试使用geoshow命令在mollweide投影中绘制世界地图。但是，我无法根据经过切割的色彩图值修改绘图中的颜色。很可能这是一个如何使用axesm和geoshow命令的问题，请帮我解决这个问题。见参考。代码如下：

G = rand(180,360);
G(1:90,:)=-1*G(1:90,:);
R = georasterref('RasterSize',size(G),...
    'Latlim',[-90 90], 'Lonlim', [-180 180],'ColumnsStartFrom','north');

% ref this link: http://stackoverflow.com/questions/34727526/matlab-plot-raster-map-with-custom-colormap/34740112#34740112
%setting up graphics parameters 
my_colormap = [254  204   92
               253  141   60
               240   59   32
               189    0   38]/255 ;

startval=min(min(G));
endval=max(max(G));
nElements =size(my_colormap,1); 
stepSize=(endval-startval)/(nElements-1);
breaks = startval:stepSize:endval;
labels = textscan(num2str(round((breaks*100))/100),'%s') ;
labels = labels{1};
[~,indices] = histc(G,breaks);

%actual graphics               
figure
hm=axesm ('mollweid', 'Frame', 'on', 'Grid', 'off');
geoshow(G,R);%geoshow(indices,R);

colormap(my_colormap);
set(gca,'color','none');
set(gca,'box','off','xtick',[],'xcolor','none')
set(gca,'box','off','ytick',[],'ycolor','none')
hc=colorbar('location','southoutside');
caxis([breaks(1) breaks(length(breaks))])%caxis([0 length(breaks)-1])
hcP = [.7 .28 .2 .02];
set(hc,'position',hcP);
h.YTickLabel = labels ;

Answer 1

我的一位同事向我指出了解决方案，并将其发布在此处供其他人参考： geoshow（G，R，＆＃39;显示类型＆＃39;＆＃39; texturemap＆＃39）

在matlab

1 个答案: