我正在尝试使用Matlab在世界地图上标记数千个地理位置。 我有这些位置的纬度和经度。 有没有好的方法可以做到这一点?谢谢。
我正在尝试使用Matlab在世界地图上标记数千个地理位置。 我有这些位置的纬度和经度。 有没有好的方法可以做到这一点?谢谢。
function [x,y] = mercatorProjection(lon, lat, width, height)
x = mod((lon+180)*width/360, width) ;
y = height/2 - log(tan((lat+90)*pi/360))*width/(2*pi);
end
% GPS positions (latitude,longitude) of some markers
data = [
-22.976730, - 43.195080 ;
55.756950, 37.614975 ;
33.605381, - 7.631940 ;
35.670479, 139.740921 ;
51.506325, - 0.127144 ;
40.714550, - 74.007124 ;
-33.869629, 151.206955 ;
-26.204944, 28.040035 ;
37.777125, -122.419644 ;
30.083740, 31.255360 ;
6.439180, 3.423480
];
labels = {
'Rio de Janeiro'
'Moscow'
'Casablanca'
'Tokyo'
'London'
'New York'
'Sydney'
'Johannesburg'
'San Francisco'
'Cairo'
'Lagos'
};
接下来我们从维基百科加载地图,应用投影并叠加标记:
% world map in Mercator projection
fname = 'https://upload.wikimedia.org/wikipedia/commons/thumb/7/74/Mercator-projection.jpg/773px-Mercator-projection.jpg';
img = imread(fname);
[imgH,imgW,~] = size(img);
% Mercator projection
[x,y] = mercatorProjection(data(:,2), data(:,1), imgW, imgH);
% plot markers on map
imshow(img, 'InitialMag',100, 'Border','tight'), hold on
plot(x,y, 'bo', 'MarkerSize',10, 'LineWidth',3)
text(x, y, labels, 'Color','w', 'VerticalAlign','bottom', 'HorizontalAlign','right')
hold off
很好的绘制全球地图的方法!
您只需更改以下内容:
imshow(I, 'InitialMag',100, 'Border','tight'), hold on
进入
imshow(img, 'InitialMag',100, 'Border','tight'), hold on
Amro的答案对我有用,但我需要做一些更改。
我正在使用Matlab 7.9和imshow是图像处理工具箱的一部分。为了显示地图而不使用imshow函数,我替换了这行代码:
imshow(img, 'InitialMag',100, 'Border','tight')
使用这个:
image(img)
它成功了。