Question

我创建了以下地图，该地图有一个均匀的灰色网格，间隔为1°，为经线和纬线：

我还希望每隔5°间隔（同时保持1°网格）使经线和平行线变粗和变黑，以便网格线与纬度和经度标签相匹配，如下所示：

我知道MATLAB对于标准2D图有major and minor grids，我过去曾使用它们。但是，据我所知，地图没有此功能。

我认为我想要做的就是通过访问地图对象属性（使用gcm或getm）并为经络和纬线的特定子集指定黑色属性来实现（使用setm）。也许函数gridm或axesm可以处理这个问题，但我不确定。

在实践中，我不知道该怎么做，因为我没有任何地图经验。我真的很感激你的帮助。

代码：

注意：此代码需要Mapping Toolbox。

% Read vector features and attributes from shapefile.
landareas = shaperead('landareas.shp', 'UseGeoCoords', true);

% Define map axes and set map properties.
axesm ('lambert',...
    'MapLonLimit', [-70 10],...
    'MapLatLimit', [30 70],...
    'MapParallels', [38.00555556 71.01111111],...
    'Frame', 'on',...
    'FLineWidth', 1,...
    'Grid', 'on',...
    'GLineStyle', '-',...
    'GLineWidth', 0.1,...
    'GColor', [.7 .7 .7]);

% Display map latitude and longitude data.
geoshow(landareas, 'FaceColor', [1 1 .5], 'EdgeColor', [.3 .3 .3]);

% Toggle and control display of graticule lines.
gridm('MLineLocation', 1,...
    'MLabelLocation', 5,...
    'PLineLocation', 1,...
    'PLabelLocation', 5);

% Toggle and control display of meridian labels.
mlabel on;

% Toggle and control display of parallel labels.
plabel on;

axis off;

Answer 1

这是一个非常简单的技巧：在同一个axes中制作第二个'Position'并在那里制作所需的主要网格。

从您的代码开始进行一些修改（我合并了axesm和gridm）：

landareas = shaperead('landareas.shp', 'UseGeoCoords', true);

% make the first axes and get the handle of it
ha = axes;

axesm ('lambert',...
    'MapLonLimit', [-70 10],...
    'MapLatLimit', [30 70],...
    'MapParallels', [38.00555556 71.01111111],...
    'Grid', 'on',...
    'GLineStyle', '-',...
    'GLineWidth', 0.1,...
    'GColor', [.7 .7 .7], ...
    'MLineLocation', 1,...
    'PLineLocation', 1);

geoshow(landareas, 'FaceColor', [1 1 .5], 'EdgeColor', [.3 .3 .3]);

axis off;

然后制作第二个轴：

% get the position of the first axes
L = get(ha, 'Position');

% make a new axes in the same position
axes('Position', L)

axesm ('lambert',...
    'MapLonLimit', [-70 10],...
    'MapLatLimit', [30 70],...
    'MapParallels', [38.00555556 71.01111111],...
    'Frame', 'on',...
    'FLineWidth', 2,...
    'Grid', 'on',...
    'GLineStyle', '-',...
    'GLineWidth', 2,...
    'GColor', [.2 .2 .2], ...
    'MLineLocation', 5,...
    'MLabelLocation', 5,...
    'PLineLocation', 5,...
    'PLabelLocation', 5);

mlabel on;
plabel on;
axis off;

你会得到这个结果：

甚至不需要获取和设置位置，因为它们都将在默认位置创建。

Answer 2

<强> TL; DR

使用您稍微修改过的代码版本：

% Read vector features and attributes from shapefile.
landareas = shaperead('landareas.shp', 'UseGeoCoords', true);

% Define map axes and set map properties.
hAx = axesm ('lambert',...  <<<<<<<<<<<<<<<<<<<<<<<<< Change #1
    'MapLonLimit', [-70 10],...
    'MapLatLimit', [30 70],...
    'MapParallels', [38.00555556 71.01111111],...
    'Frame', 'on',...
    'FLineWidth', 1,...
    'Grid', 'on',...
    'GLineStyle', '-',...
    'GLineWidth', 0.1,...
    'GColor', [.7 .7 .7]);

% Display map latitude and longitude data.
geoshow(landareas, 'FaceColor', [1 1 .5], 'EdgeColor', [.3 .3 .3]);

% Toggle and control display of graticule lines.
hGrid = gridm(... <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Change #2
  'MLineLocation', 1,... 
  'PLineLocation', 1);

% Change #3:
hGridCopy = copyobj(hGrid,hAx); 
set(hGridCopy,'Tag','oldGrid');  % We use this to hide the grid from handlem,
gridm(...                          so that it doesn't get found and deleted.
  'MLineLocation', 5,...
  'MLabelLocation', 5,...
  'PLineLocation', 5,...
  'PLabelLocation', 5,...
  'GColor',[0,0,0],...
  'GLineWidth',3);

% Toggle and control display of meridian labels.
mlabel on;

% Toggle and control display of parallel labels.
plabel on;

axis off;

要了解这个hacky解决方案的来源，请参阅以下部分。

当查看axesm生成的句柄时，我们可以了解一下地图的结构：

% Change from this:
axesm ('lambert',...
% To this:
hAx = axesm ('lambert',...

然后，hAx.Children给出了以下内容：

ans = 

  30×1 graphics array:

  Text     (PLabel)
  ... 
  Text     (MLabel)
  Line     (Parallel)
  Line     (Meridian)
  Group
  Patch    (Frame)

由此我们了解到，每个行集合，即Meridian和Parallel，实际上都是单个对象，这意味着某个类型的所有行都必须具有完全相同的样式 - 不像＆＃34;传统＆＃34; 次要和主要 grid行（可能不同）。如果我们希望某些行不同，我们必须创建一个具有所需属性的新对象（行的子集）。

现在，创建这些极地网格线的最简单方法是再次调用gridm，但是，这会删除已经绘制的所有网格。它是如何做到的？ gridm在内部调用handlem('grid')，查找Tag Parallel或Meridian的图形对象，找到的对象稍后会被删除并替换为具有新面貌。那么我们该怎么办？我们＆＃34;隐藏＆＃34;当然，来自handlem的网格（请参阅Change #3;通过更改这些对象的Tag属性）！然后我们只创建另一个网格，我们就完成了。

地图的主要和次要刻度？

代码：

2 个答案: