6

私は、世界地図上の数千の地理的位置をmatlabでマークしようとしています。私はそれらの場所の緯度と経度を持っています。これを行うための良い方法はありますか？ありがとう。

4

3 に答える 3

8

これは、ツールボックスを必要としない例です。

まず、メルカトル図法を使用して経度/緯度の位置を変換する関数を作成します。

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

出力

于 2012-07-25T21:26:18.017 に答える
2

世界をプロットするのに最適な方法です！

以下を変更する必要があります。

imshow(I, 'InitialMag',100, 'Border','tight'), hold on

の中へ

imshow(img, 'InitialMag',100, 'Border','tight'), hold on
于 2012-11-26T22:33:33.223 に答える
2

アムロの答えは私のために働いたが、私はいくつかの変更をしなければならなかった。

私はMatlab7.9を使用しており、imshowは画像処理ツールボックスの一部です。imshow関数を使用せずに地図を表示するために、次の行を置き換えました。

imshow(img, 'InitialMag',100, 'Border','tight')

これで：

image(img)

そしてそれはうまくいった。

于 2013-12-03T20:31:23.170 に答える