66 lines
2.0 KiB
Matlab
66 lines
2.0 KiB
Matlab
%% This function was automatically generated. When modifying its signature, take care to apply
|
|
%% modifications also to the descriptor files in the repository.
|
|
%% function Geomap_function(Catalog)
|
|
function Figure1 = Geomap_function_4(Catalog, BIN1, BIN2)
|
|
%% GEOMAP_FUNCTION
|
|
%% TODO: Put your application code here
|
|
|
|
% Visualization of events on a Geomap
|
|
% by ML 14 Dec 2021
|
|
%close all; clc
|
|
|
|
|
|
% Set interpreters to correctly display underscore symbols
|
|
set(0, 'DefaultTextInterpreter', 'none')
|
|
set(0, 'DefaultLegendInterpreter', 'none')
|
|
set(0, 'DefaultAxesTickLabelInterpreter', 'none')
|
|
|
|
% Load catalog and strip the file name
|
|
%[baseFileName, folder] = uigetfile('*.mat');
|
|
%fullFileName = fullfile(folder, baseFileName);
|
|
%load(fullFileName);
|
|
%chart_title=baseFileName(1:end-4);
|
|
|
|
% Extraction of longitude and latitude from the Catalog
|
|
lat=Catalog(strcmp('Lat',{Catalog.field})).val;
|
|
long=Catalog(strcmp('Long',{Catalog.field})).val;
|
|
|
|
% Extraction of magnitude from the Catalog which type recognition
|
|
if sum((strcmp('ML',{Catalog.field})))==1
|
|
Magn=Catalog(strcmp('ML',{Catalog.field})).val;
|
|
else
|
|
Magn=Catalog(strcmp('Mw',{Catalog.field})).val;
|
|
end
|
|
|
|
% Assignment of categories to magnitude data
|
|
%a=1;
|
|
%b=2;
|
|
atext=num2str(BIN1);
|
|
btext=num2str(BIN2);
|
|
below='<';
|
|
above='>';
|
|
between='-';
|
|
bin1label=append(below, atext);
|
|
bin2label=append(atext, between, btext);
|
|
bin3label=append(above, btext);
|
|
Category = discretize(Magn,[-100 BIN1 BIN2 100],'categorical', {bin1label, bin2label, bin3label});
|
|
|
|
% Assignment of categories' colors (first category in first row, etc.) in RGB standard
|
|
color_RGB=[0.1 1.0 0.3;...
|
|
0.8 0.8 0.2;...
|
|
1.0 0.0 0.0];
|
|
|
|
% Drawing a Geomap
|
|
gb=geobubble(lat, long, Magn, Category,...
|
|
'Basemap', 'satellite',...
|
|
'BubbleColorList', color_RGB,...
|
|
'ColorLegendTitle', 'Magnitude range',...
|
|
'SizeLegendTitle', 'Event magnitude',...
|
|
'GridVisible', 'on',...
|
|
'ScalebarVisible', 'on');
|
|
gb.BubbleWidthRange = [1 20];
|
|
%title(chart_title);
|
|
%Figure1 = ancestor(gb, 'figure');
|
|
saveas(gcf,'Figure1','png');
|
|
end
|