create grid of squares over map using mapping toolbox

Question

Amir 2016-3-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/270783-create-grid-of-squares-over-map-using-mapping-toolbox

评论： Jos (10584) 2016-3-2

map_grid.jpg

Hi,

I am working on a project where i have specific locations (latitude/longitude) and need to create a type of heat-map which represents how populated each region.

I would like to create a grid of squares over the map (which i have created using the mapping toolbox) and color each square based on the number of points registered inside it. I have done a similar thing using circles, where i verify if the point is inside each circle using the "inpoly" function. However, i am finding it challenging to do the same using squares. I guess as long as i know the corner points of each square, i would still be able to follow the same approach.

I have attached a sample image to show you the type of grid i am looking for. Ideally, i would like to be able to specify the size of each of the squares in for example km. The map was produced using the following code (the grids were overlaid in paint):

figure(1); 
f=worldmap([33 68],[-15 37]); 
geoshow('landareas.shp', 'FaceColor', [1 1 1],'DefaultEdgeColor', 'b')

Thank you very much for you help in advance.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jos (10584) 2016-3-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/270783-create-grid-of-squares-over-map-using-mapping-toolbox#answer_211813

在 MATLAB Online 中打开

You can add a grid like this:

f=worldmap([33 68],[-15 37]); 
geoshow('landareas.shp', 'FaceColor', [1 1 1],'DefaultEdgeColor', 'b') 
AxisHandle = gca ;
GridHandle = gridxy(linspace(AxisHandle.XLim(1),AxisHandle.XLim(2),25), ...
      linspace(AxisHandle.YLim(1),AxisHandle.YLim(2),25)) ;
set(GridHandle,'color',[.6 .6 .6],'LineWidth',.5) ;
uistack(GridHandle,'top')

My function GRIDXY is available at the Matlab File Exchange and can be downloaded here: http://www.mathworks.com/matlabcentral/fileexchange/9973

4 个评论
显示 2更早的评论隐藏 2更早的评论

Amir 2016-3-1

Jos,

Yes it does provide the points, but i would need to have the coordinates in latitude/longitude. Also, the "25" in that lines specifies the number of squares generated over the plot. However, i need to specify the size of each of the squares instead. so for example, the side length of each square to be 25km.

Jos (10584) 2016-3-2

在 MATLAB Online 中打开

I am not sure about the units of the axes, but you can create these points using colon notation (start:step:end) So if the units are kilometers, this would do, I think:

Xpoints = AxisHandle.XLim(1):25:AxisHandle.XLim(2)

请先登录，再进行评论。