matlabの地図にデータをプロットする

Question

私はMATLAB R2015aを使用しています。基本的に、グリッド点の緯度を第1列に、グリッド点の経度を第2列に格納したグリッドを定義しました。列ベクトルに格納された領域の地震エネルギーのデータがあります各要素は、対応する格子点のエネルギーに対応する。私はこのコードを使用して表面プロットを行っている（ここでe_latとe_longは、それぞれのグリッド行列の1列目と2列目である）： -

これは一例であり

function [b] = cumulative_plot(beam, e_lat,e_long, t_start, t_end) 
%CUMULATIVE_PLOT Plots the cumulative energy of the earthquake 
%% Input Arguments 
% *beam* - Energy for each time increment (columns) for each grid point (rows) 
% 
% *e_lat* - Vector containing Latitudes of the grid points 
% 
% *e_long* - Vector containing Longitudes of the grid points 
% 
% *t_start* - starting time 
% 
% *t_end* - ending time 
% 
% *t_start* and *t_end* define the time window within which the energy is 
% to be considered 

%% Code 

b = []; 
b = sum(beam(:,t_start:t_end)')'; % Adding the energy within the time window 
b = b./max(b); % Normalising 

fn = 'cumulative_energy.txt'; 
f = fopen(fn,'w'); 

for i=1:length(e_lat) 
    fprintf(f,'%f %f %f \n',e_long(i),e_lat(i),b(i)); 
end 

fclose(f); 

energy_surf = fit([e_long,e_lat],b, 'loess'); 
plot(energy_surf,'style','contour'); 
hold on; 
plot3(73.6400 ,34.5239 ,20,'s','MarkerSize',20,'MarkerEdgeColor','k','MarkerFaceColor','k') 
hold on; 
plot3(94.709,23.03,20,'s','MarkerSize',20,'MarkerEdgeColor','b','MarkerFaceColor','b') 
shading interp 
alpha(1) 
view(0,90) 
box off 
colorbar 
title(['Cumu Energy(0.05 - 0.2 Hz) at seconds = ' num2str(t_start)],'FontWeight','bold','FontSize',15,'FontName','Times'); 
xlabel('Long/degree','FontWeight','bold','FontSize',13,'FontName','Times'); 
ylabel('Lat/degree','FontWeight','bold','FontSize',13,'FontName','Times'); 

end 

（Iは処理しています実際のデータ）： -

cumulative_plot(b_corr,e_lat,e_long,1,20); 

私は、指定された地域の地理的地図上にエネルギーデータの等高線図を作りたいです。これは可能ですか？

そして、これは私が達成したいものの一種である - ：

は、より良いアイデアを与えるために、これは私が今持っているものである（紫の円形マーカーや他のものなしだけで基本エネルギー。） ： -

Answer 1

あなたが山の細部のBMP画像をお持ちの場合は、RGB形式でデータを保存し、その強度によって累積エネルギースケールのデータとそれを混ぜます。 Intensityはアルファブレンディング値を提供します。