How to generate a 2D meshgrid based on the coordinates (longitude and latitude) of points

I have several points' coordinates, I want to merge a meshgrid of the points looks like the following figure.

For example, I have five points with coordinates info, the ideal generated meshgrid should be like this (the purple square which is overshadowed).

Thank you very much!

Answer 1

[Voss]

xy = [1 3; 2 9; 6 6; 8 8; 9 1];

% [xx,yy] = meshgrid(sort(xy(:,1)),sort(xy(:,2))) % note: maybe you want to sort x and y first
[xx,yy] = meshgrid(xy(:,1),xy(:,2))

xx = 5×5

1 2 6 8 9 1 2 6 8 9 1 2 6 8 9 1 2 6 8 9 1 2 6 8 9

yy = 5×5

3 3 3 3 3 9 9 9 9 9 6 6 6 6 6 8 8 8 8 8 1 1 1 1 1

plot(xy(:,1),xy(:,2),'ok','MarkerFaceColor','k')
hold on
xline(xx(1,:))
yline(yy(:,1))
xlim([0 10]);
ylim([0 10]);
set(gca(),'Visible','off')

[Wang Zhen]

Hello, impressive! Thanks!

However, the input value I use is longitude and latitude follows:

22.7397 113.4884

22.7670 113.5093

22.7704 113.5073

22.7670 113.5083

22.7670 113.5083

22.7696 113.5113

22.7712 113.5088

22.7712 113.5088

22.7695 113.5128

22.7694 113.5132

22.7680 113.5180

22.7691 113.5203

22.7690 113.5203

22.7690 113.5203

22.7690 113.5203

22.7693 113.5185

22.7692 113.5185

22.7684 113.5182

22.7704 113.5148

22.7708 113.5113

22.7703 113.5152

22.7712 113.5164

22.7692 113.5161

22.7716 113.5113

22.7727 113.5102

22.7727 113.5103

22.7733 113.5084

22.7730 113.5084

22.7681 113.5157

22.7715 113.5128

22.7746 113.5084

The xline command has some error... 'Passing multiple values to ConstantLine is not supported.'

New to MATLAB, Hoping for your further update.

[Voss]

Maybe that error is due to different MATLAB versions, because it runs here:

xy = [
   22.7397  113.4884
   22.7670  113.5093
   22.7704  113.5073
   22.7670  113.5083
   22.7670  113.5083
   22.7696  113.5113
   22.7712  113.5088
   22.7712  113.5088
   22.7695  113.5128
   22.7694  113.5132
   22.7680  113.5180
   22.7691  113.5203
   22.7690  113.5203
   22.7690  113.5203
   22.7690  113.5203
   22.7693  113.5185
   22.7692  113.5185
   22.7684  113.5182
   22.7704  113.5148
   22.7708  113.5113
   22.7703  113.5152
   22.7712  113.5164
   22.7692  113.5161
   22.7716  113.5113
   22.7727  113.5102
   22.7727  113.5103
   22.7733  113.5084
   22.7730  113.5084
   22.7681  113.5157
   22.7715  113.5128
   22.7746  113.5084
   ]; % actually yx

[xx,yy] = meshgrid(xy(:,2),xy(:,1))

xx = 31×31

113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128 113.4884 113.5093 113.5073 113.5083 113.5083 113.5113 113.5088 113.5088 113.5128 113.5132 113.5180 113.5203 113.5203 113.5203 113.5203 113.5185 113.5185 113.5182 113.5148 113.5113 113.5152 113.5164 113.5161 113.5113 113.5102 113.5103 113.5084 113.5084 113.5157 113.5128

yy = 31×31

22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7397 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7704 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7670 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7696 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7712 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7695 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694 22.7694

plot(xy(:,2),xy(:,1),'ok','MarkerFaceColor','k')
hold on
xline(xx(1,:))
yline(yy(:,1))

If you really need to plot the grid lines, there are different ways. Easiest might be to use a couple of loops:

figure
plot(xy(:,2),xy(:,1),'ok','MarkerFaceColor','k')
hold on
for ii = 1:size(xx,2)
    xline(xx(1,ii));
end
for ii = 1:size(yy,1)
    yline(yy(ii,1));
end

[Wang Zhen]

Thank you very much, I'm currently using MATLAB2020b, I think it's time for a update.

[Voss]

You're welcome!

[Wang Zhen]

https://www.mathworks.com/matlabcentral/answers/1743855-how-to-calculate-the-l2-euclidean-and-l1-cityblock-or-we-say-manhattan-distance-distance-based?s_tid=srchtitle

Hi, here is another problem I encountered which is associated with this, please take a look if you are free.

Thank you very much!