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 이미지가있는 경우
오, 알겠습니다. 나는 당신의 질문에 어떻게 대답하지 않았는지 알기 때문에 대답을 삭제했습니다. 원하는 것을 위해 맵핑 도구 상자 https://www.mathworks.com/help/map/index.html을 확인할 수 있다고 생각합니다.하지만 가지고 있지 않아 더 이상 당신을 도울 수 없습니다. 죄송합니다. –
예, 고마워요. 이 플롯은 GMT에서 matlab에서 데이터를 처리 한 후에 수행되었습니다. 나는 matlab에서 전체 작업을하고 싶었다. 슬프게도 상사가 도구 상자를 배울 수있는 여분의 시간을주지 않습니다. – Archon
문제에 대해 생각하고 있었고 플롯 할 토폴로지 데이터를 제공하는 경우 매핑 도구 상자없이 이러한 플롯을 만들 수 있습니다. 또는이 정보를 Matlab 읽기 형식으로 쉽게 제공 할 수있는 서비스를 알아 두십시오. –