2017-02-28 1 views
0

MATLAB R2015a를 사용하고 있습니다. 기본적으로 첫 번째 열의 격자 점의 위도와 두 번째 열의 격자 점의 경도를 저장하는 행렬을 정의했습니다. 열 벡터에 저장된 영역의 지진 에너지에 대한 데이터가 있습니다. 여기서 각 요소는 해당 격자 점의 에너지에 해당합니다. 이 코드를 사용하여 표면 플롯했을 (여기 e_lat 및 e_long 각각 격자 행렬의 제 1 및 제 2 열이다) -matlab에서지도에 데이터 플로팅

이것은 일례이다
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); 

이 에너지 데이터의 윤곽 플롯을 지정된 영역의 지리적지도에 만들고 싶습니다. 이것이 가능한가?

그리고 이것은 내가 무엇을 달성하고자하는의 종류 - :

이 더 나은 아이디어를 제공하기 위해이 내가 지금 무엇을 가지고 (보라색 원형 마커와 다른 것없이 그냥 기본 에너지.) - 당신이 산의 세부 사항의 BMP 이미지가있는 경우

IMG_20170228_173241_HDR (1).jpg

+0

오, 알겠습니다. 나는 당신의 질문에 어떻게 대답하지 않았는지 알기 때문에 대답을 삭제했습니다. 원하는 것을 위해 맵핑 도구 상자 https://www.mathworks.com/help/map/index.html을 확인할 수 있다고 생각합니다.하지만 가지고 있지 않아 더 이상 당신을 도울 수 없습니다. 죄송합니다. –

+0

예, 고마워요. 이 플롯은 GMT에서 matlab에서 데이터를 처리 한 후에 수행되었습니다. 나는 matlab에서 전체 작업을하고 싶었다. 슬프게도 상사가 도구 상자를 배울 수있는 여분의 시간을주지 않습니다. – Archon

+1

문제에 대해 생각하고 있었고 플롯 할 토폴로지 데이터를 제공하는 경우 매핑 도구 상자없이 이러한 플롯을 만들 수 있습니다. 또는이 정보를 Matlab 읽기 형식으로 쉽게 제공 할 수있는 서비스를 알아 두십시오. –

답변

0

는 RGB 형식으로 데이터를 저장, 그 강도에 의해 누적 에너지 규모의 데이터로 섞는다. 강도는 알파 블렌딩 값을 제공합니다.

+0

사실 나는 이것을 가능한 한 일반적인 것으로 만들고 싶었습니다. 지역 지형도의 BMP를 얻으려면 다른 매핑 도구를 사용해야합니다. 제공된 툴박스를 사용하여 MATLAB에서 모든 작업을 수행하고 싶었습니다. – Archon

+0

기본적으로 이미지 오버레이를 만들고 싶습니다. 배경은 Google이나 기타 소스 등의 이미지 여야하는 지리적지도입니다. 오버레이는 데이터 컨투어 플롯으로 알파 벤딩을 사용하여 트랜스 펜션에 표시합니다. – bhamadicharef