From db7d8cf7b3329de7280c14fcb841df81790b953f Mon Sep 17 00:00:00 2001 From: Michal_Lelonek Date: Thu, 9 Dec 2021 16:47:32 +0100 Subject: [PATCH] Update 'Geomap_function_3.m' --- Geomap_function_3.m | 53 ++++++++++++++++++++++++++++++++++++++++++++++++++--- 1 file changed, 50 insertions(+), 3 deletions(-) diff --git a/Geomap_function_3.m b/Geomap_function_3.m index e0d6715..58fd64d 100644 --- a/Geomap_function_3.m +++ b/Geomap_function_3.m @@ -1,8 +1,55 @@ %% 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_3(Catalog) - -%% GEOMAP_FUNCTION_3 +%% function Geomap_function(Catalog) +function Figure1 = Geomap_function_3(Catalog) +%% GEOMAP_FUNCTION %% TODO: Put your application code here +% Visualization of events on a Geomap +% by ML 25 Nov 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 +Category = discretize(Magn,[-100 1 2 100],'categorical', {'<1', '1-2', '>2'}); + +% 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', 'streets',... + 'BubbleColorList', color_RGB,... + 'ColorLegendTitle', 'Magnitude category',... + 'SizeLegendTitle', 'Event magnitude',... + 'GridVisible', 'on',... + 'ScalebarVisible', 'on'); +gb.BubbleWidthRange = [1 20]; +%title(chart_title); +%Figure1 = ancestor(gb, 'figure'); +saveas(gcf,'Figure1','png'); end