Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (16)

Search Parameters:
Keywords = horizontal visibility graph

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
21 pages, 11315 KiB  
Article
Distinction of Chaos from Randomness Is Not Possible from the Degree Distribution of the Visibility and Phase Space Reconstruction Graphs
by Alexandros K. Angelidis, Konstantinos Goulas, Charalampos Bratsas, Georgios C. Makris, Michael P. Hanias, Stavros G. Stavrinides and Ioannis E. Antoniou
Entropy 2024, 26(4), 341; https://doi.org/10.3390/e26040341 - 17 Apr 2024
Cited by 2 | Viewed by 2230
Abstract
We investigate whether it is possible to distinguish chaotic time series from random time series using network theory. In this perspective, we selected four methods to generate graphs from time series: the natural, the horizontal, the limited penetrable horizontal visibility graph, and the [...] Read more.
We investigate whether it is possible to distinguish chaotic time series from random time series using network theory. In this perspective, we selected four methods to generate graphs from time series: the natural, the horizontal, the limited penetrable horizontal visibility graph, and the phase space reconstruction method. These methods claim that the distinction of chaos from randomness is possible by studying the degree distribution of the generated graphs. We evaluated these methods by computing the results for chaotic time series from the 2D Torus Automorphisms, the chaotic Lorenz system, and a random sequence derived from the normal distribution. Although the results confirm previous studies, we found that the distinction of chaos from randomness is not generally possible in the context of the above methodologies. Full article
Show Figures

Figure 1

Figure 1
<p>The time series used in the proof of Proposition 1 is presented with the corresponding Natural Visibility Graph. The resulting graph is the graph <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>K</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msub> </mrow> </semantics></math> which is non-planar.</p>
Full article ">Figure 2
<p>The first row illustrates how a time series is transformed into a network using the natural visibility algorithm, while the second row depicts the same time series transformed using the horizontal visibility algorithm. Notice that the HVG is a subgraph of NVG [<a href="#B23-entropy-26-00341" class="html-bibr">23</a>].</p>
Full article ">Figure 2 Cont.
<p>The first row illustrates how a time series is transformed into a network using the natural visibility algorithm, while the second row depicts the same time series transformed using the horizontal visibility algorithm. Notice that the HVG is a subgraph of NVG [<a href="#B23-entropy-26-00341" class="html-bibr">23</a>].</p>
Full article ">Figure 3
<p>Example of (<b>a</b>) a time series with 10 data values and (<b>b</b>) its corresponding LPHVG with ρ = 1, where every node corresponds to a time series data. The limited penetrable horizontal visibility lines between data points define the links connecting nodes in the graph. Black lines generate the HVG, and red lines are those added to generate the LPHVG for ρ = 1 [<a href="#B44-entropy-26-00341" class="html-bibr">44</a>].</p>
Full article ">Figure 4
<p>The time series used in the proof of Proposition 2 is presented with the corresponding LPHVG. The resulting graph is the graph <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>K</mi> </mrow> <mrow> <mn>5</mn> </mrow> </msub> </mrow> </semantics></math>, which is non-planar.</p>
Full article ">Figure 5
<p>Degree distribution of the Natural Visibility Graphs of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 6
<p>Degree distribution of the Horizontal Visibility Graphs of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 7
<p>Degree distribution of the Limited Penetrable Horizontal Visibility Graphs with limited penetrable distance ρ = 1 of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 7 Cont.
<p>Degree distribution of the Limited Penetrable Horizontal Visibility Graphs with limited penetrable distance ρ = 1 of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 8
<p>Degree distribution of the Limited Penetrable Horizontal Visibility Graphs with limited penetrable distance ρ = 2 of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 9
<p>Degree distribution of the Phase Space Reconstruction Graphs of two Torus Automorphisms corresponding to the chaotic time series with <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> (first row) and <math display="inline"><semantics> <mrow> <mi>h</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> (second row), in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 10
<p>Degree distribution of the Natural Visibility Graph of the Lorenz System in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 11
<p>Degree distribution of the Horizontal Visibility Graph of the Lorenz System in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 12
<p>Degree distribution of the Limited Penetrable Horizontal Visibility Graphs with limited penetrable distance ρ = 1 (first raw) and ρ = 2 (second raw) of the Lorenz System in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 13
<p>Degree distribution of the Phase Space Reconstruction Graph of the Lorenz System in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 14
<p>Degree distribution of the Natural Visibility Graph of the Random Sequence in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 15
<p>Degree distribution of the Horizontal Visibility Graph of the Random Sequence in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 16
<p>Degree distribution of the Limited Penetrable Horizontal Visibility Graphs with limited penetrable distance ρ = 1 (first raw) and ρ = 2 (second raw) of the Random Sequence, in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">Figure 17
<p>Degree distribution of the Phase Space Reconstruction Graph of the Random Sequence in lin–lin scale (left column), in log–log scale (middle column), and in lin–log scale (right column).</p>
Full article ">
23 pages, 12259 KiB  
Article
Exploring the Entropy-Based Classification of Time Series Using Visibility Graphs from Chaotic Maps
by J. Alberto Conejero, Andrei Velichko, Òscar Garibo-i-Orts, Yuriy Izotov and Viet-Thanh Pham
Mathematics 2024, 12(7), 938; https://doi.org/10.3390/math12070938 - 22 Mar 2024
Cited by 2 | Viewed by 1447
Abstract
The classification of time series using machine learning (ML) analysis and entropy-based features is an urgent task for the study of nonlinear signals in the fields of finance, biology and medicine, including EEG analysis and Brain–Computer Interfacing. As several entropy measures exist, the [...] Read more.
The classification of time series using machine learning (ML) analysis and entropy-based features is an urgent task for the study of nonlinear signals in the fields of finance, biology and medicine, including EEG analysis and Brain–Computer Interfacing. As several entropy measures exist, the problem is assessing the effectiveness of entropies used as features for the ML classification of nonlinear dynamics of time series. We propose a method, called global efficiency (GEFMCC), for assessing the effectiveness of entropy features using several chaotic mappings. GEFMCC is a fitness function for optimizing the type and parameters of entropies for time series classification problems. We analyze fuzzy entropy (FuzzyEn) and neural network entropy (NNetEn) for four discrete mappings, the logistic map, the sine map, the Planck map, and the two-memristor-based map, with a base length time series of 300 elements. FuzzyEn has greater GEFMCC in the classification task compared to NNetEn. However, NNetEn classification efficiency is higher than FuzzyEn for some local areas of the time series dynamics. The results of using horizontal visibility graphs (HVG) instead of the raw time series demonstrate the GEFMCC decrease after HVG time series transformation. However, the GEFMCC increases after applying the HVG for some local areas of time series dynamics. The scientific community can use the results to explore the efficiency of the entropy-based classification of time series in “The Entropy Universe”. An implementation of the algorithms in Python is presented. Full article
(This article belongs to the Special Issue Mathematical Modeling, Optimization and Machine Learning, 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>The workflow diagram of the proposed method of assessing the global efficiency of entropy.</p>
Full article ">Figure 2
<p>Illustrative example of the natural visibility graph representation for a time series (<b>left</b>) and the horizontal visibility graph representation for the same time series (<b>right</b>). The arrows in the images explains the projection of the visibility horizon when constructing the graph.</p>
Full article ">Figure 3
<p>Main steps of NNetEn calculation [<a href="#B1-mathematics-12-00938" class="html-bibr">1</a>]. The figure shows the main stages of calculation NNetEn based on the reservoir neural network LogNNet, where the reservoir is filled with the time series under study, and the entropy value is proportional to the classification metric of the reference database.</p>
Full article ">Figure 4
<p>Section of the buffering diagram of the logistic map, on which two adjacent sets of series are highlighted corresponding to <span class="html-italic">r<sub>j</sub></span><sub>−1</sub> = 3.634 and <span class="html-italic">r<sub>j</sub></span> = 3.636 (<b>a</b>), series (<span class="html-italic">x</span><sub>1</sub>, …, <span class="html-italic">x</span><sub>300</sub>) for <span class="html-italic">r<sub>j</sub></span><sub>−1</sub> = 3.634 (<b>b</b>), series (<span class="html-italic">x</span><sub>1</sub>, …, <span class="html-italic">x</span><sub>300</sub>) for <span class="html-italic">r<sub>j</sub></span> = 3.636 (<b>c</b>), and FuzzyEn values for 100 time series for two classes (MCC = 1) (<b>d</b>). The figure explains the method for calculating the classification metric for the time series of a discrete map for neighboring sets corresponding to two neighboring partitions by <span class="html-italic">r</span>.</p>
Full article ">Figure 5
<p>Distribution of FuzzyEn in Classes 1 and 2 with <span class="html-italic">r<sub>j</sub></span><sub>−1</sub> = 3.688 and <span class="html-italic">r<sub>j</sub></span> = 3.69 (MCC~0.45). The figure shows an example of entropy distribution for poorly separable classes and MCC~0.45.</p>
Full article ">Figure 6
<p>Bifurcation diagrams for the logistic map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>), and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure 7
<p>MCC(<span class="html-italic">r</span>) dependences for FuzzyEn before and after HVG transformation, as well as their difference in ΔMCC (<b>a</b>); MCC(<span class="html-italic">r</span>) dependences for NNetEn before and after HVG transformation, as well as their difference in ΔMCC (<b>b</b>). Calculations were made for the logistic map.</p>
Full article ">Figure 8
<p>ΔMCC(<span class="html-italic">r</span>) dependences for FuzzyEn and NeNetEn. Calculations were made for the logistic map.</p>
Full article ">Figure 9
<p>Bifurcation diagrams for the TMBM map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>); and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure 10
<p>MCC(<span class="html-italic">r</span>) dependences for FuzzyEn before and after HVG transformation, as well as their difference in ΔMCC (<b>a</b>); MCC(<span class="html-italic">r</span>) dependences for NNetEn before and after HVG transformation, as well as their difference in ΔMCC (<b>b</b>). Calculations were made for the TMBM map.</p>
Full article ">Figure 11
<p>ΔMCC(<span class="html-italic">r</span>) dependences for FuzzyEn and NeNetEn. Calculations were made for the TMBM map.</p>
Full article ">Figure A1
<p>Bifurcation diagrams for sine map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>); and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure A1 Cont.
<p>Bifurcation diagrams for sine map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>); and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure A2
<p>MCC(<span class="html-italic">r</span>) dependences for FuzzyEn before and after HVG transformation, as well as their difference in ΔMCC (<b>a</b>); MCC(<span class="html-italic">r</span>) dependences for NNetEn before and after HVG transformation, as well as their difference in ΔMCC (<b>b</b>). Calculations were made for sine map.</p>
Full article ">Figure A3
<p>Bifurcation diagrams for Planck map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>); and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure A3 Cont.
<p>Bifurcation diagrams for Planck map (<b>a</b>); the dependence of entropy on the parameter <span class="html-italic">r</span> for NNetEn_AV (<b>b</b>); and FuzzyEn_AV before and after HVG transformation (<b>c</b>). The figures show changes in the dynamics and irregularity of time series depending on the parameter.</p>
Full article ">Figure A4
<p>MCC(<span class="html-italic">r</span>) dependences for FuzzyEn before and after HVG transformation, as well as their difference in ΔMCC (<b>a</b>); MCC(<span class="html-italic">r</span>) dependences for NNetEn before and after HVG transformation, as well as their difference in ΔMCC (<b>b</b>). Calculations were made for Planck map.</p>
Full article ">
18 pages, 7562 KiB  
Article
Graph- and Machine-Learning-Based Texture Classification
by Musrrat Ali, Sanoj Kumar, Rahul Pal, Manoj K. Singh and Deepika Saini
Electronics 2023, 12(22), 4626; https://doi.org/10.3390/electronics12224626 - 12 Nov 2023
Cited by 2 | Viewed by 1971
Abstract
The analysis of textures is an important task in image processing and computer vision because it provides significant data for image retrieval, synthesis, segmentation, and classification. Automatic texture recognition is difficult, however, and necessitates advanced computational techniques due to the complexity and diversity [...] Read more.
The analysis of textures is an important task in image processing and computer vision because it provides significant data for image retrieval, synthesis, segmentation, and classification. Automatic texture recognition is difficult, however, and necessitates advanced computational techniques due to the complexity and diversity of natural textures. This paper presents a method for classifying textures using graphs; specifically, natural and horizontal visibility graphs. The related image natural visibility graph (INVG) and image horizontal visibility graph (IHVG) are used to obtain features for classifying textures. These features are the clustering coefficient and the degree distribution. The suggested outcomes show that the aforementioned technique outperforms traditional ones and even comes close to matching the performance of convolutional neural networks (CNNs). Classifiers such as the support vector machine (SVM), K-nearest neighbor (KNN), decision tree (DT), and random forest (RF) are utilized for the categorization. The suggested method is tested on well-known image datasets like the Brodatz texture and the Salzburg texture image (STex) datasets. The results are positive, showing the potential of graph methods for texture classification. Full article
Show Figures

Figure 1

Figure 1
<p>Natural visibility graph from time series data.</p>
Full article ">Figure 2
<p>Horizontal visibility graph from time series data.</p>
Full article ">Figure 3
<p>Sample texture images from the Brodatz texture image dataset.</p>
Full article ">Figure 4
<p>Sample texture images from the Salzburg texture image dataset.</p>
Full article ">Figure 5
<p>Flowchart of the proposed method.</p>
Full article ">Figure 6
<p>Accuracy of (<b>a</b>) INVG, (<b>b</b>) IHVG, and (<b>c</b>) combined graphs on Brodatz dataset.</p>
Full article ">Figure 7
<p>Accuracy of (<b>a</b>) degree distribution, (<b>b</b>) clustering coefficient, and (<b>c</b>) combined features on Brodatz dataset.</p>
Full article ">Figure 8
<p>Accuracy of (<b>a</b>) degree distribution, (<b>b</b>) clustering coefficient, and (<b>c</b>) combined features on STex dataset.</p>
Full article ">Figure 9
<p>Accuracy of (<b>a</b>) INVG, (<b>b</b>) IHVG, and (<b>c</b>) combined graphs on STex dataset.</p>
Full article ">
12 pages, 4173 KiB  
Article
Entropy Generation for MHD Peristaltic Transport of Non-Newtonian Fluid in a Horizontal Symmetric Divergent Channel
by Kinda Abuasbeh, Bilal Ahmed, Azmat Ullah Khan Niazi and Muath Awadalla
Symmetry 2023, 15(2), 359; https://doi.org/10.3390/sym15020359 - 29 Jan 2023
Cited by 3 | Viewed by 1571
Abstract
The analysis in view is proposed to investigate the impacts of entropy in the peristaltically flown Ree–Eyring fluid under the stress of a normally imposed uniform magnetic field in a non-uniform symmetric channel of varying thickness. The administering equations of the present flow [...] Read more.
The analysis in view is proposed to investigate the impacts of entropy in the peristaltically flown Ree–Eyring fluid under the stress of a normally imposed uniform magnetic field in a non-uniform symmetric channel of varying thickness. The administering equations of the present flow problem are switched into the non-dimensional form and then reduced by the availing of long wavelengths and creeping flow regime restrictions. The analytical treatment for the developed problem is performed to attain closed-form solutions which are further displayed as graphs of velocity, pressure, temperature, and entropy distribution. The trapping phenomenon has also been an area of our current examination. The role of relevant pronounced parameters such as the Brinkmann number, Hartmann number, and Ree–Eyring parameter for throwing vivid impacts are also concerned. It has been inferred that both the Brinkmann number and Ree–Eyring parameter with rising values inflate temperature and entropy profiles. The velocity profile shows the symmetric nature due to the horizontally assumed symmetric channel of varying thickness. The circulation of streamlines and bolus formations is visibly reduced in response to the increasing Hartmann number. Full article
(This article belongs to the Special Issue Asymmetric and Symmetric Study with PDE)
Show Figures

Figure 1

Figure 1
<p>Schematic flow diagram for peristaltic transport of Ree–Eyring fluid.</p>
Full article ">Figure 2
<p>Comparison of pressure profile based on obtained results with that of Bhatti et al. (2017) [<a href="#B35-symmetry-15-00359" class="html-bibr">35</a>].</p>
Full article ">Figure 3
<p>Variation in velocity for different values of (<b>a</b>) Hartmann number <math display="inline"><semantics> <mi>M</mi> </semantics></math> and (<b>b</b>) Casson parameter <math display="inline"><semantics> <mrow> <mi>ζ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Pressure rises per wavelength for different values of (<b>a</b>) Hartmann number <math display="inline"><semantics> <mi>M</mi> </semantics></math> and (<b>b</b>) Casson parameter <math display="inline"><semantics> <mrow> <mi>ζ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>(<b>a</b>) Temperature distribution for different values of Hartmann number <math display="inline"><semantics> <mi>M</mi> </semantics></math> (<b>b</b>) Different values of Casson parameter <math display="inline"><semantics> <mrow> <mi>ζ</mi> </mrow> </semantics></math> and (<b>c</b>) Different values of Brinkmann number <math display="inline"><semantics> <mrow> <msub> <mi>B</mi> <mi>r</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Entropy generation for different values of (<b>a</b>) Hartmann number <math display="inline"><semantics> <mi>M</mi> </semantics></math> and (<b>b</b>) Casson parameter <math display="inline"><semantics> <mrow> <mi>ζ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Entropy generation for different values of (<b>a</b>) Brinkmann umber <math display="inline"><semantics> <mrow> <msub> <mi>B</mi> <mi>r</mi> </msub> </mrow> </semantics></math> and (<b>b</b>) different values of temperature difference paramter <math display="inline"><semantics> <mrow> <mi>Λ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>Streamlines for different values of <math display="inline"><semantics> <mi>M</mi> </semantics></math> with fixed values of <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mo> </mo> <mi>ζ</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mover accent="true"> <mi>Q</mi> <mo>¯</mo> </mover> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>2</mn> <mo>,</mo> </mrow> </semantics></math> (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>3</mn> <mo>,</mo> <mo> </mo> <mi>and</mi> </mrow> </semantics></math> (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Streamlines for different values of <math display="inline"><semantics> <mi>ζ</mi> </semantics></math> with fixed values of <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mo> </mo> <mi>ϕ</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mover accent="true"> <mi>Q</mi> <mo>¯</mo> </mover> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>ζ</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> </mrow> </semantics></math> (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>ζ</mi> <mo>=</mo> <mn>2.0</mn> <mo>,</mo> <mrow> <mo> </mo> <mi>and</mi> </mrow> </mrow> </semantics></math> (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>ζ</mi> <mo>=</mo> <mn>5.0</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Streamlines for different values of <math display="inline"><semantics> <mi>ϕ</mi> </semantics></math> with fixed values of <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mo> </mo> <mi>ζ</mi> <mo>=</mo> <mn>0.5</mn> <mo>,</mo> <mover accent="true"> <mi>Q</mi> <mo>¯</mo> </mover> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math> (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.4</mn> <mo>,</mo> </mrow> </semantics></math> (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.6</mn> <mo>,</mo> <mo> </mo> <mi>and</mi> <mo> </mo> </mrow> </semantics></math> (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>ϕ</mi> <mo>=</mo> <mn>0.8</mn> </mrow> </semantics></math>.</p>
Full article ">
12 pages, 5618 KiB  
Article
Assigning Degrees of Stochasticity to Blazar Light Curves in the Radio Band Using Complex Networks
by Belén Acosta-Tripailao, Walter Max-Moerbeck, Denisse Pastén and Pablo S. Moya
Entropy 2022, 24(8), 1063; https://doi.org/10.3390/e24081063 - 2 Aug 2022
Cited by 4 | Viewed by 1958
Abstract
We focus on characterizing the high-energy emission mechanisms of blazars by analyzing the variability in the radio band of the light curves of more than a thousand sources. We are interested in assigning complexity parameters to these sources, modeling the time series of [...] Read more.
We focus on characterizing the high-energy emission mechanisms of blazars by analyzing the variability in the radio band of the light curves of more than a thousand sources. We are interested in assigning complexity parameters to these sources, modeling the time series of the light curves with the method of the Horizontal Visibility Graph (HVG), which allows us to obtain properties from degree distributions, such as a characteristic exponent to describe its stochasticity and the Kullback–Leibler Divergence (KLD), presenting a new perspective to the methods commonly used to study Active Galactic Nuclei (AGN). We contrast these parameters with the excess variance, which is an astronomical measurement of variability in light curves; at the same time, we use the spectral classification of the sources. While it is not possible to find significant correlations with the excess variance, the degree distributions extracted from the network are detecting differences related to the spectral classification of blazars. These differences suggest a chaotic behavior in the time series for the BL Lac sources and a correlated stochastic behavior in the time series for the FSRQ sources. Our results show that complex networks may be a valuable alternative tool to study AGNs according to the variability of their energy output. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) A graphical description for modeling a time series with a horizontal visibility graph. Each data correspond to a longitudinal node. (<b>b</b>) According to the visibility of each node, we can calculate the degrees <math display="inline"><semantics> <msub> <mi>k</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>out</mi> </msub> </semantics></math> for directed HVG and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>ud</mi> </msub> </semantics></math> for undirected HVG: that is, how many times the node establishes a connection as a function of time direction and independently of this one.</p>
Full article ">Figure 2
<p>(<b>a</b>) Example radio light curve for the blazar PKS 1502+106 [<a href="#B34-entropy-24-01063" class="html-bibr">34</a>], a FSRQ source. (<b>b</b>) Semilog plot of degree distributions from DHVG, <math display="inline"><semantics> <msub> <mi>P</mi> <mi>in</mi> </msub> </semantics></math> in magenta and <math display="inline"><semantics> <msub> <mi>P</mi> <mi>out</mi> </msub> </semantics></math> in red, and from UHVG, <span class="html-italic">P</span> in black dots with their fit on the dashed line. The value of excess variance <math display="inline"><semantics> <msubsup> <mi>σ</mi> <mi>rms</mi> <mn>2</mn> </msubsup> </semantics></math> of the light curve, <span class="html-italic">D</span>, from the distance between degree distributions <b>in</b> and <b>out</b>, and <math display="inline"><semantics> <mi>γ</mi> </semantics></math>, from the exponential behavior <math display="inline"><semantics> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>∼</mo> <msup> <mi mathvariant="normal">e</mi> <mrow> <mo>−</mo> <mi>γ</mi> <mi>k</mi> </mrow> </msup> </mrow> </semantics></math>, are shown.</p>
Full article ">Figure 3
<p>Probability density function (PDF) of (<b>a</b>) <math display="inline"><semantics> <mi>γ</mi> </semantics></math> and (<b>b</b>) <math display="inline"><semantics> <mrow> <mo form="prefix">log</mo> <mo>(</mo> <mi>D</mi> <mo>)</mo> </mrow> </semantics></math> values for the different subclasses of blazars. Yellow for FSRQ and blue for BL Lac. There are 843 FSRQ and 380 BL Lac sources. With the dotted line in (<b>a</b>), we mark the limit <math display="inline"><semantics> <mrow> <msub> <mi>γ</mi> <mi>un</mi> </msub> <mo>≈</mo> <mn>0.405</mn> </mrow> </semantics></math> between correlated stochastic and chaotic time series [<a href="#B25-entropy-24-01063" class="html-bibr">25</a>].</p>
Full article ">Figure 4
<p>Scatter plots of square root of the excess variance vs. (<b>a</b>) <math display="inline"><semantics> <mi>γ</mi> </semantics></math> and (<b>b</b>) <span class="html-italic">D</span>. Yellow for FSRQ and blue for BL Lac. There are 843 FSRQ and 380 BL Lac sources.</p>
Full article ">Figure A1
<p>Examples of selected probability distribution fits to FSRQ light curves. Degree distributions <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </semantics></math> in blue dots with their fit on the dashed black line.</p>
Full article ">Figure A2
<p>Examples of selected probability distribution fits to BL Lac light curves. Degree distributions <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </semantics></math> in blue dots with their fit on the dashed black line.</p>
Full article ">Figure A3
<p>Examples of discarded probability distribution fits for FSRQ light curves. Degree distributions <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </semantics></math> in blue dots with their fit on the dashed black line.</p>
Full article ">Figure A4
<p>Examples of discarded probability distribution fits for BL Lac light curves. Degree distributions <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </semantics></math> in blue dots with their fit on the dashed black line.</p>
Full article ">
8 pages, 1486 KiB  
Proceeding Paper
Main Climatic Characteristics of the International Airport “Antonio Maceo” during the Period 2017–2021
by Beatriz Valdés Díaz, Amanda Maria De Armas Echevarria and Patricia Coll-Hidalgo
Environ. Sci. Proc. 2022, 19(1), 45; https://doi.org/10.3390/ecas2022-12855 - 25 Jul 2022
Cited by 1 | Viewed by 866
Abstract
The climatological characterization of variables allows us to understand the average behavior of atmospheric conditions, detect extremes and fluctuations, and the relationships of variables with geographic physical factors; it is presented as another aide for weather forecasting. The aim of this research is [...] Read more.
The climatological characterization of variables allows us to understand the average behavior of atmospheric conditions, detect extremes and fluctuations, and the relationships of variables with geographic physical factors; it is presented as another aide for weather forecasting. The aim of this research is characterizing the behavior of the meteorological variables at the Antonio Maceo International Airport in Santiago de Cuba for the period 2017–2021. Antonio Maceo International Airport has the particularity of being located in a complex relief, exposed to marked breeze influences, and a significant number of wind shear pilot reports. The characterization was based on the concepts and graphs of descriptive statistics. The mean monthly distribution of the variables: temperature, relative humidity, atmospheric pressure, and wind speed and direction was obtained and is represented. The distributions of the maximum monthly and annual accumulations, the cloud cover, and the ranges of the reduction of the horizontal and vertical visibility were analyzed. Full article
(This article belongs to the Proceedings of The 5th International Electronic Conference on Atmospheric Sciences)
Show Figures

Figure 1

Figure 1
<p>Location of the “Antonio Maceo” International Airport in Santiago de Cuba.</p>
Full article ">Figure 2
<p>Mean distribution of mean, maximum, and minimum air temperatures.</p>
Full article ">Figure 3
<p>The behavior of (<bold>a</bold>) the mean air temperature and (<bold>b</bold>) relative humidity in MUCU during the period 2017–2021.</p>
Full article ">Figure 4
<p>Annual behavior of the mean, minimum, and maximum atmospheric pressure in MUCU during the period 2017–2021.</p>
Full article ">Figure 5
<p>Behavior (<bold>a</bold>) monthly and (<bold>b</bold>) annual of the accumulated precipitation in MUCU during the period 2017–2021.</p>
Full article ">Figure 6
<p>Rose of the wind direction and velocity for the period 2017–2021.</p>
Full article ">
12 pages, 538 KiB  
Article
The North–South Asymmetry of Sunspot Relative Numbers Based on Complex Network Technique
by Hengyu Xu, Yu Fei, Chun Li, Jiajuan Liang, Xinan Tian and Zhongjie Wan
Symmetry 2021, 13(11), 2228; https://doi.org/10.3390/sym13112228 - 22 Nov 2021
Cited by 6 | Viewed by 1792
Abstract
Solar magnetic activity exhibits a complex nonlinear behavior, but its dynamic process has not been fully understood. As the complex network technique can better capture the dynamics of nonlinear system, the visibility graphs (VG), the horizontal visibility graphs (HVG), and the limited penetrable [...] Read more.
Solar magnetic activity exhibits a complex nonlinear behavior, but its dynamic process has not been fully understood. As the complex network technique can better capture the dynamics of nonlinear system, the visibility graphs (VG), the horizontal visibility graphs (HVG), and the limited penetrable visibility graphs (LPVG) are applied to implement the mapping of sunspot relative numbers in the northern and southern hemispheres. The results show that these three methods can capture important information of nonlinear dynamics existing in the long-term hemispheric sunspot activity. In the presentation of the results, the network degree sequence of the HVG method changes preferentially to the original data series as well as the VG and the LPVG, while both the VG and the LPVG slightly lag behind the original time series, which provides some new ideas for the nonlinear dynamics of the hemispheric asymmetry in the two hemispheres. Meanwhile, the use of statistical feature-skewness values and complex network visibility graphs can yield some complementary information for mutual verification. Full article
Show Figures

Figure 1

Figure 1
<p>The daily time series of sunspot relative numbers from the NAOJ/Mitaka observatory (<b>upper panel</b>) along with the International Sunspot Number (ISSN; <b>lower panel</b>) from 1 March 1939 to 30 September 2019.</p>
Full article ">Figure 2
<p>Monthly distribution of sunspot relative numbers in the northern and southern hemispheres during the time interval from March 1939 to September 2019.</p>
Full article ">Figure 3
<p>Monthly hemispheric sunspot relative numbers (<b>top panel</b>), the absolute difference (<b>middle panel</b>), and the normalized difference (<b>bottom panel</b>) with a running width of <math display="inline"><semantics> <mrow> <mi>ω</mi> <mo>=</mo> <mn>270</mn> </mrow> </semantics></math> months.</p>
Full article ">Figure 4
<p>Probability distribution functions of sunspot relative numbers in the northern (<b>upper panel</b>) and southern (<b>lower panel</b>) hemispheres, respectively.</p>
Full article ">Figure 5
<p>Upper panel: the yearly skewness values of hemispheric sunspot relative numbers. Lower panel: the difference of the yearly skewness of hemispheric sunspot relative numbers.</p>
Full article ">Figure 6
<p>The degree sequences of the sunspot relative numbers for the northern and southern hemispheres (<b>top panel</b>), and the absolute (<b>middle panel</b>) and normalized (<b>bottom panel</b>) hemispheric differences of the degree sequences.</p>
Full article ">Figure 7
<p>Similar to <a href="#symmetry-13-02228-f006" class="html-fig">Figure 6</a>, but for the HVG.</p>
Full article ">Figure 8
<p>Similar to <a href="#symmetry-13-02228-f006" class="html-fig">Figure 6</a>, but for the LPVG.</p>
Full article ">
31 pages, 1660 KiB  
Article
A Comparative Eye Tracking Study of Usability—Towards Sustainable Web Design
by Mihai Țichindelean, Monica Teodora Țichindelean, Iuliana Cetină and Gheorghe Orzan
Sustainability 2021, 13(18), 10415; https://doi.org/10.3390/su131810415 - 18 Sep 2021
Cited by 20 | Viewed by 7482
Abstract
Websites are one of the most frequently used communication environments, and creating sustainable web designs should be an objective for all companies. Ensuring high usability is proving to be one of the main contributors to sustainable web design, reducing usage time, eliminating frustration [...] Read more.
Websites are one of the most frequently used communication environments, and creating sustainable web designs should be an objective for all companies. Ensuring high usability is proving to be one of the main contributors to sustainable web design, reducing usage time, eliminating frustration and increasing satisfaction and retention. The present paper studies the usability of different website landing pages, seeking to identify the elements, structures and designs that increase usability. The study analyzed the behavior of 22 participants during their interaction with five different landing pages while they performed three tasks on the webpage and freely viewed each page for one minute. The stimuli were represented by five different banking websites, each of them presenting the task content in a different mode (text, image, symbol, graph, etc.).; the data obtained from the eye tracker (fixations location, order and duration, saccades, revisits of the same element, etc.), together with the data from the applied survey lead to interesting conclusions: the top, center and right sides of the webpage attract the most attention; the use of pictures depicting persons increase visibility; the scanpaths follow a vertical and horizontal direction; numerical data should be presented through graphs or tables. Even if a user's past experience influences their experience on a website, we show that the design of the webpage itself has a greater influence on webpage usability. Full article
Show Figures

Figure 1

Figure 1
<p>Heat Map BRD bank website in free-viewing conditions.</p>
Full article ">Figure 2
<p>Scan path map showing the first 10 viewed elements on the website of BRD bank.</p>
Full article ">Figure 3
<p>Recommendation for landing page layout and structure.</p>
Full article ">
20 pages, 1475 KiB  
Article
Visibility Graph Analysis of IOTA and IoTeX Price Series: An Intentional Risk-Based Strategy to Use 5G for IoT
by Alberto Partida, Regino Criado and Miguel Romance
Electronics 2021, 10(18), 2282; https://doi.org/10.3390/electronics10182282 - 17 Sep 2021
Cited by 11 | Viewed by 3274
Abstract
The transformation of time series into complex networks through visibility graphs is an innovative way to study time-based events. In this work, we use visibility graphs to transform IOTA and IoTeX price volatility time series into complex networks. Our aim is twofold: first, [...] Read more.
The transformation of time series into complex networks through visibility graphs is an innovative way to study time-based events. In this work, we use visibility graphs to transform IOTA and IoTeX price volatility time series into complex networks. Our aim is twofold: first, to better understand the markets of the two most capitalised Internet of Things (IoT) platforms at the time of writing. IOTA runs on a public directed acyclic graph (DAG) and IoTeX on a blockchain. Second, to suggest how 5G can improve information security in these two key IoT platforms. The analysis of the networks created by the natural and horizontal visibility graphs shows, first, that both IOTA and IoTeX are still at their infancy in their development, with IoTex seemingly developing faster. Second, both IoT tokens form communities in a hierarchical structure, and third, 5G can accelerate their development. We use intentional risk management as a lever to understand the impact of 5G on IOTA and IoTeX. Our results lead us to provide a set of design recommendations that contribute to improving information security in future 5G-based IoT implementations. Full article
(This article belongs to the Section Networks)
Show Figures

Figure 1

Figure 1
<p>Daily price volatility data for IOTA and IoTeX. Subplots (<b>a</b>,<b>b</b>) display the daily volatility time series for IOTA and IoTeX, respectively. Subplots (<b>c</b>,<b>d</b>) display their components: highest and lowest prices in USD per day for IOTA and IoTeX.</p>
Full article ">Figure 2
<p>Visibility and horizontal visibility graphs for IOTA and IoTeX price volatility. Subplots (<b>a</b>,<b>b</b>) display the visibility graph derived from the daily volatility time series for IOTA and IoTeX, respectively (last 20 days of the dataset). Subplots (<b>c</b>,<b>d</b>) display the horizontal visibility graph derived from the daily volatility time series for IOTA and IoTeX, respectively (last 20 days of the dataset). The depiction of these graphs is the outcome of our own Python code.</p>
Full article ">Figure 3
<p>Number of nodes for each degree in the networks stemming from the visibility and horizontal visibility graphs for IOTA and IoTeX price volatility. Subplots (<b>a</b>,<b>b</b>) display the degree of the network stemming from the visibility graph and the best power law fit that the function <span class="html-italic">curve_fit</span> provides together with the corresponding best <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>σ</mi> </semantics></math>. Subplots (<b>c</b>,<b>d</b>) display the degree of the network stemming from the horizontal visibility graph and the best power law fit that the function <span class="html-italic">curve_fit</span> provides and the corresponding best <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>σ</mi> </semantics></math>.</p>
Full article ">Figure 4
<p>Both empirical and fit probability density functions (PDF) and complementary cumulative distribution functions (CCDF) using the <span class="html-italic">powerlaw</span> module by Alstott et al. [<a href="#B55-electronics-10-02282" class="html-bibr">55</a>].</p>
Full article ">Figure 5
<p>Best power law fit using <span class="html-italic">curve_fit</span> with the average of clustering coefficients per degree for (<b>a</b>) IOTA VG, (<b>b</b>) IoTeX VG, (<b>c</b>) IOTA HVG and (<b>d</b>) IoTeX HVG.</p>
Full article ">Figure 6
<p>Communities identified by the <span class="html-italic">networkx</span> module <span class="html-italic">community API</span> in (<b>a</b>) IOTA VG, (<b>b</b>) IoTeX VG, (<b>c</b>) IOTA HVG and (<b>d</b>) IoTeX HVG throughout the timeline in network graph format.</p>
Full article ">Figure 7
<p>Communities identified by the <span class="html-italic">networkx</span> module <span class="html-italic">community API</span> in IOTA, (<b>a</b>,<b>c</b>), and IoTeX, (<b>b</b>,<b>d</b>). Graphs (<b>a</b>,<b>b</b>) use a coloured-coded continuous line and graphs (<b>c</b>,<b>d</b>) a coloured-coded scattered plot.</p>
Full article ">Figure 8
<p>Infographic summary of “Visibility graph analysis of IOTA and IoTeX price series: An intentional risk-based strategy to use 5G for IoT”.</p>
Full article ">
15 pages, 1722 KiB  
Article
Assessment of Arrow-of-Time Metrics for the Characterization of Underwater Explosions
by Ramón Miralles, Guillermo Lara, Alicia Carrión and Manuel Bou-Cabo
Sensors 2021, 21(17), 5952; https://doi.org/10.3390/s21175952 - 4 Sep 2021
Cited by 2 | Viewed by 2237
Abstract
Anthropogenic impulsive sound sources with high intensity are a threat to marine life and it is crucial to keep them under control to preserve the biodiversity of marine ecosystems. Underwater explosions are one of the representatives of these impulsive sound sources, and existing [...] Read more.
Anthropogenic impulsive sound sources with high intensity are a threat to marine life and it is crucial to keep them under control to preserve the biodiversity of marine ecosystems. Underwater explosions are one of the representatives of these impulsive sound sources, and existing detection techniques are generally based on monitoring the pressure level as well as some frequency-related features. In this paper, we propose a complementary approach to the underwater explosion detection problem through assessing the arrow of time. The arrow of time of the pressure waves coming from underwater explosions conveys information about the complex characteristics of the nonlinear physical processes taking place as a consequence of the explosion to some extent. We present a thorough review of the characterization of arrows of time in time-series, and then provide specific details regarding their applications in passive acoustic monitoring. Visibility graph-based metrics, specifically the direct horizontal visibility graph of the instantaneous phase, have the best performance when assessing the arrow of time in real explosions compared to similar acoustic events of different kinds. The proposed technique has been validated in both simulations and real underwater explosions. Full article
(This article belongs to the Special Issue Underwater Acoustic Sensors and Applications)
Show Figures

Figure 1

Figure 1
<p>Main wave reflections produced by an UNDEX event received by a passive acoustic recorder moored at the bottom of the sea under constant speed of sound assumption.</p>
Full article ">Figure 2
<p>Example of how ARIMA(10,1,10) achieves better modeling (smaller residuals) than ARMA(10,10) in signals having a trend similar to that of an UNDEX event.</p>
Full article ">Figure 3
<p>Illustration of the way the horizontal visibility graph is obtained for a time series <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </semantics></math> according to the criteria in (<a href="#FD5-sensors-21-05952" class="html-disp-formula">5</a>).</p>
Full article ">Figure 4
<p>Example of simulated non-UNDEX and UNDEX events (although they are very different, distinguishing between them may be challenging when they are severely clipped).</p>
Full article ">Figure 5
<p>Evolution of all the presented metrics in simulated events when the SNR decreases (Y axis in arbitrary units). The results were obtained for 500 Monte Carlo runs and for unclipped acoustic events.</p>
Full article ">Figure 6
<p>Evolution of all the presented metrics in simulated events when the SNR decreases (Y axis in arbitrary units). The results were obtained for 500 Monte Carlo runs for events clipped to 50% of their maximum value.</p>
Full article ">Figure 7
<p>Cartagena region, showing the tuna farm and the UNDEX sites (red squares) as well as the hydrophone location (green star). The distance from the UNDEX to the hydrophone was approximately 1.2 km.</p>
Full article ">Figure 8
<p>Example of a non-UNDEX event and a real UNDEX event from the database over the duration of the first 0.25 s.</p>
Full article ">Figure 9
<p>Assessment of the arrow-of-time metrics in the database of real events (Y axis in arbitrary units). The colors indicate the human expert classification: red corresponds to UNDEX, blue corresponds to non-UNDEX, and yellow corresponds to UNSURE categories. The threshold (dashed green line) was set to obtain a fixed detection of UNDEX events equal to 90%.</p>
Full article ">
16 pages, 1821 KiB  
Article
Automatic Reconstruction of Multi-Level Indoor Spaces from Point Cloud and Trajectory
by Gahyeon Lim and Nakju Doh
Sensors 2021, 21(10), 3493; https://doi.org/10.3390/s21103493 - 17 May 2021
Cited by 20 | Viewed by 3208
Abstract
Remarkable progress in the development of modeling methods for indoor spaces has been made in recent years with a focus on the reconstruction of complex environments, such as multi-room and multi-level buildings. Existing methods represent indoor structure models as a combination of several [...] Read more.
Remarkable progress in the development of modeling methods for indoor spaces has been made in recent years with a focus on the reconstruction of complex environments, such as multi-room and multi-level buildings. Existing methods represent indoor structure models as a combination of several sub-spaces, which are constructed by room segmentation or horizontal slicing approach that divide the multi-room or multi-level building environments into several segments. In this study, we propose an automatic reconstruction method of multi-level indoor spaces with unique models, including inter-room and inter-floor connections from point cloud and trajectory. We construct structural points from registered point cloud and extract piece-wise planar segments from the structural points. Then, a three-dimensional space decomposition is conducted and water-tight meshes are generated with energy minimization using graph cut algorithm. The data term of the energy function is expressed as a difference in visibility between each decomposed space and trajectory. The proposed method allows modeling of indoor spaces in complex environments, such as multi-room, room-less, and multi-level buildings. The performance of the proposed approach is evaluated for seven indoor space datasets. Full article
(This article belongs to the Section Remote Sensors)
Show Figures

Figure 1

Figure 1
<p>Registered point cloud of an indoor space with complex multiple floors environment (C2 in our dataset). (<b>a</b>) The bird’s-eye view, (<b>b</b>) top view, and (<b>c</b>) side view of the indoor space. The combination of first and 0.5 floors is denoted by red, first floor is denoted by green, second floor is denoted by blue, and a combination of first and second floors is denoted by yellow. The conventional horizontal slicing approach is not suitable for complex multi-level building environments.</p>
Full article ">Figure 2
<p>Mesh model built using our proposed method of indoor structure modeling includes inter-room (doors or openings) and inter-floor connections (stairs), denoted by the red and yellow boxes, respectively.</p>
Full article ">Figure 3
<p>Flowchart of the proposed method.</p>
Full article ">Figure 4
<p>Structural points (grey) and the object points (red) that are segmented from the registered point cloud of our dataset A1.</p>
Full article ">Figure 5
<p>Top view of (<b>a</b>) piece-wise planar segments and (<b>b</b>) 3D cell complex of our dataset A1.</p>
Full article ">Figure 6
<p>Illustration of an indoor space consisting of three rooms and two inter-room connections expressed in 2D for visualization. The (<b>a</b>) structural points (black lines) and trajectory points (black dots) and (<b>b</b>) decomposed 3D spaces. (<b>b</b>) illustrates centers of the cells (colored crosses), visible points in the cells (colored dots), and trajectory points with the smallest difference in visibility (colored circles) for each cell (colored boxes). The cells marked in red, green, and blue select the trajectory point with the smallest difference in visibility but a cell marked in purple, which is outside the indoor space, has large differences in visibility for all trajectory points.</p>
Full article ">Figure 7
<p>Illustration of two pairs of adjacent cells of an indoor space depicted in <a href="#sensors-21-03493-f006" class="html-fig">Figure 6</a>a expressed in 2D for visualization. The area occupied by structural points located on the adjacent face of red and green cells is close to zero, while that of red and blue cells is close to the area of face shared by the two cells.</p>
Full article ">Figure 8
<p>(<b>a</b>) Robot system and (<b>b</b>) backpack system employed for data acquisition.</p>
Full article ">Figure 9
<p>Structural points (grey) and the trajectory points (red) of (<b>a</b>) dataset A1 and (<b>b</b>) dataset C3. The trajectory is continuous through multiple rooms or multiple levels.</p>
Full article ">Figure 10
<p>Results of the generated water-tight meshes. The bird’s-eye view of the entire spaces (<b>left</b>) and the fist floor without ceilings (<b>right</b>) are displayed.</p>
Full article ">Figure 10 Cont.
<p>Results of the generated water-tight meshes. The bird’s-eye view of the entire spaces (<b>left</b>) and the fist floor without ceilings (<b>right</b>) are displayed.</p>
Full article ">Figure 11
<p>Detailed images of dataset C2. (<b>a</b>) inter-floor and three inter-room connections, (<b>b</b>) inter-floor and two inter-room connections, and (<b>c</b>) inter-floor connection.</p>
Full article ">Figure 12
<p>Detailed images of dataset C3. (<b>a</b>) inter-floor connection between underground and first floor, (<b>b</b>) inter-floor connection between first and second floor, and (<b>c</b>) top view of inter-floor connection between third and fourth floor.</p>
Full article ">Figure 13
<p>Point cloud of large-scale and room-less indoor space with curved surfaces and cylindrical pillars.</p>
Full article ">
12 pages, 1478 KiB  
Article
Applying the Horizontal Visibility Graph Method to Study Irreversibility of Electromagnetic Turbulence in Non-Thermal Plasmas
by Belén Acosta-Tripailao, Denisse Pastén and Pablo S. Moya
Entropy 2021, 23(4), 470; https://doi.org/10.3390/e23040470 - 16 Apr 2021
Cited by 17 | Viewed by 3120
Abstract
One of the fundamental open questions in plasma physics is the role of non-thermal particles distributions in poorly collisional plasma environments, a system that is commonly found throughout the Universe, e.g., the solar wind and the Earth’s magnetosphere correspond to natural plasma physics [...] Read more.
One of the fundamental open questions in plasma physics is the role of non-thermal particles distributions in poorly collisional plasma environments, a system that is commonly found throughout the Universe, e.g., the solar wind and the Earth’s magnetosphere correspond to natural plasma physics laboratories in which turbulent phenomena can be studied. Our study perspective is born from the method of Horizontal Visibility Graph (HVG) that has been developed in the last years to analyze time series avoiding the tedium and the high computational cost that other methods offer. Here, we build a complex network based on directed HVG technique applied to magnetic field fluctuations time series obtained from Particle In Cell (PIC) simulations of a magnetized collisionless plasma to distinguish the degree distributions and calculate the Kullback–Leibler Divergence (KLD) as a measure of relative entropy of data sets produced by processes that are not in equilibrium. First, we analyze the connectivity probability distribution for the undirected version of HVG finding how the Kappa distribution for low values of κ tends to be an uncorrelated time series, while the Maxwell–Boltzmann distribution shows a correlated stochastic processes behavior. Subsequently, we investigate the degree of temporary irreversibility of magnetic fluctuations that are self-generated by the plasma, comparing the case of a thermal plasma (described by a Maxwell–Botzmann velocity distribution function) with non-thermal Kappa distributions. We have shown that the KLD associated to the HVG is able to distinguish the level of reversibility that is associated to the thermal equilibrium in the plasma, because the dissipative degree of the system increases as the value of κ parameter decreases and the distribution function departs from the Maxwell–Boltzmann equilibrium. Full article
Show Figures

Figure 1

Figure 1
<p>Construction of Horizontal Visibility Graph. Top, a time series where the degree <math display="inline"><semantics> <msub> <mi>k</mi> <mi>in</mi> </msub> </semantics></math> for in-going links and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>out</mi> </msub> </semantics></math> for out-going links of each of the <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> nodes are detailed. Bottom, probability distribution <span class="html-italic">P</span> in relation to degree <span class="html-italic">k</span>, where <math display="inline"><semantics> <msub> <mi>n</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>n</mi> <mi>out</mi> </msub> </semantics></math> correspond to the frequency of appearance of the degrees <math display="inline"><semantics> <msub> <mi>k</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>k</mi> <mi>out</mi> </msub> </semantics></math>, respectively, defining <math display="inline"><semantics> <msub> <mi>P</mi> <mi>in</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>P</mi> <mi>out</mi> </msub> </semantics></math>.</p>
Full article ">Figure 2
<p>(<b>Left</b>) Average magnetic field energy density fluctuations <math display="inline"><semantics> <msup> <mrow> <mo>(</mo> <mi>δ</mi> <mi>B</mi> <mo>/</mo> <msub> <mi>B</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </semantics></math> as a function of time obtained from Particle In Cell (PIC) simulations for Maxwell–Boltzmann (where MB represents <math display="inline"><semantics> <mrow> <msub> <mi>κ</mi> <mi>e</mi> </msub> <mo>→</mo> <mo>∞</mo> </mrow> </semantics></math>) and Kappa distributions considering different values of the <math display="inline"><semantics> <msub> <mi>κ</mi> <mi>e</mi> </msub> </semantics></math> parameter. (<b>Right</b>) Detrended average magnetic field energy density magnitude.</p>
Full article ">Figure 3
<p>Semi-log plot of the degree distributions of HVG associated to Kappa and Maxwell–Boltzmann distribution. There is an exponential behavior <math display="inline"><semantics> <mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>∼</mo> <mo form="prefix">exp</mo> <mfenced separators="" open="(" close=")"> <mo>−</mo> <mi>γ</mi> <mi>k</mi> </mfenced> </mrow> </semantics></math> and the <math display="inline"><semantics> <mi>γ</mi> </semantics></math> value is shown for each distribution. The left panel corresponds to the results for the magnetic field of the trend data from <a href="#entropy-23-00470-f002" class="html-fig">Figure 2</a> (left), while the right panel for the detrended data from <a href="#entropy-23-00470-f002" class="html-fig">Figure 2</a> (right).</p>
Full article ">Figure 4
<p>KL-Divergence (<span class="html-italic">D</span>) of magnetic field for different Kappa distributions. (<b>Left</b>) Horizontal Visibility Graph (HVG) method applied on the original data. (<b>Right</b>) HVG on the detrended data. The technique used to determine whether the data represent a reversible process consists of applying the HVG algorithm to randomly disordered copies of the data, obtaining the standard deviation <math display="inline"><semantics> <mi>σ</mi> </semantics></math> around the average divergence computed using the disordered data (black dot and vertical lines).</p>
Full article ">Figure 5
<p>Temporal evolution of the KL-divergence considering a moving window that covers 8000 data overlapping every 1000 data on the magnetic time series. (<b>Left</b>) HVG method applied on the original data and (<b>Right</b>) on the detrended data.</p>
Full article ">
21 pages, 2283 KiB  
Article
EEG-Based Emotion Recognition Using an Improved Weighted Horizontal Visibility Graph
by Tianjiao Kong, Jie Shao, Jiuyuan Hu, Xin Yang, Shiyiling Yang and Reza Malekian
Sensors 2021, 21(5), 1870; https://doi.org/10.3390/s21051870 - 7 Mar 2021
Cited by 19 | Viewed by 4076
Abstract
Emotion recognition, as a challenging and active research area, has received considerable awareness in recent years. In this study, an attempt was made to extract complex network features from electroencephalogram (EEG) signals for emotion recognition. We proposed a novel method of constructing forward [...] Read more.
Emotion recognition, as a challenging and active research area, has received considerable awareness in recent years. In this study, an attempt was made to extract complex network features from electroencephalogram (EEG) signals for emotion recognition. We proposed a novel method of constructing forward weighted horizontal visibility graphs (FWHVG) and backward weighted horizontal visibility graphs (BWHVG) based on angle measurement. The two types of complex networks were used to extract network features. Then, the two feature matrices were fused into a single feature matrix to classify EEG signals. The average emotion recognition accuracies based on complex network features of proposed method in the valence and arousal dimension were 97.53% and 97.75%. The proposed method achieved classification accuracies of 98.12% and 98.06% for valence and arousal when combined with time-domain features. Full article
(This article belongs to the Special Issue Emotion Monitoring System Based on Sensors and Data Analysis)
Show Figures

Figure 1

Figure 1
<p>The block diagram of the proposed method for EEG emotion recognition.</p>
Full article ">Figure 2
<p>Horizontal visibility graph of a time series. (<b>a</b>) The histogram of time series; (<b>b</b>) Its corresponding HVG.</p>
Full article ">Figure 3
<p>Graphical illustration of FWHVG of the time series. (<b>a</b>) Angle measurement of FWHVG; (<b>b</b>) Corresponding FWHVG of the time series.</p>
Full article ">Figure 4
<p>Graphical illustration of BWHVG of the time series. (<b>a</b>) Angle measurement of BWHVG; (<b>b</b>) Corresponding BWHVG of the time series.</p>
Full article ">Figure 5
<p>The adjacency matrix of networks based on VG. (<b>a</b>) EEG series with low valence; (<b>b</b>) EEG series with high valence.</p>
Full article ">Figure 6
<p>The adjacency matrix of networks based on HVG. (<b>a</b>) EEG series with low valence; (<b>b</b>) EEG series with high valence.</p>
Full article ">Figure 7
<p>Visualization normalized weight matrices of forward and backward weighted complex networks. (<b>a</b>) FWVG; (<b>b</b>) BWVG; (<b>c</b>) FWHVG; (<b>d</b>) BWHVG.</p>
Full article ">Figure 8
<p>Box plots of the average weighted degree feature of 32 channels. (<b>a</b>) HVG; (<b>b</b>) DWHVG.</p>
Full article ">Figure 9
<p>Classification accuracies of valence.</p>
Full article ">Figure 10
<p>Classification accuracies of arousal.</p>
Full article ">
12 pages, 3113 KiB  
Article
Analysis of Air Mean Temperature Anomalies by Using Horizontal Visibility Graphs
by Javier Gómez-Gómez, Rafael Carmona-Cabezas, Elena Sánchez-López, Eduardo Gutiérrez de Ravé and Francisco José Jiménez-Hornero
Entropy 2021, 23(2), 207; https://doi.org/10.3390/e23020207 - 8 Feb 2021
Cited by 6 | Viewed by 2094
Abstract
The last decades have been successively warmer at the Earth’s surface. An increasing interest in climate variability is appearing, and many research works have investigated the main effects on different climate variables. Some of them apply complex networks approaches to explore the spatial [...] Read more.
The last decades have been successively warmer at the Earth’s surface. An increasing interest in climate variability is appearing, and many research works have investigated the main effects on different climate variables. Some of them apply complex networks approaches to explore the spatial relation between distinct grid points or stations. In this work, the authors investigate whether topological properties change over several years. To this aim, we explore the application of the horizontal visibility graph (HVG) approach which maps a time series into a complex network. Data used in this study include a 60-year period of daily mean temperature anomalies in several stations over the Iberian Peninsula (Spain). Average degree, degree distribution exponent, and global clustering coefficient were analyzed. Interestingly, results show that they agree on a lack of significant trends, unlike annual mean values of anomalies, which present a characteristic upward trend. The main conclusions obtained are that complex networks structures and nonlinear features, such as weak correlations, appear not to be affected by rising temperatures derived from global climate conditions. Furthermore, different locations present a similar behavior and the intrinsic nature of these signals seems to be well described by network parameters. Full article
(This article belongs to the Section Complexity)
Show Figures

Figure 1

Figure 1
<p>Meteorological stations located over the Iberian Peninsula (Spain) selected for this study.</p>
Full article ">Figure 2
<p>The left margin: (<b>a</b>,<b>b</b>) Plots that illustrate mean temperature time series of Valencia and Sevilla in the year 2019, respectively. (<b>c</b>,<b>d</b>) The corresponding anomalies for the same stations and period. The right margin: (<b>e</b>) Example of application of the horizontal visibility graph (HVG) algorithm to the first ten values of Sevilla anomalies in 1960. (<b>f</b>) Network obtained from the previous plot.</p>
Full article ">Figure 3
<p>Annual average of daily mean temperature anomalies for every location.</p>
Full article ">Figure 4
<p>The left margin: (<b>a</b>,<b>b</b>) Degree distributions of Valencia and Sevilla stations in year 2019, respectively. Red lines are the least-square fits of values. (<b>c</b>,<b>d</b>) Annual evolution of slopes obtained from the previous linear fits (<math display="inline"><semantics> <mi>γ</mi> </semantics></math> exponent) for Valencia and Sevilla, respectively. (<b>e</b>,<b>f</b>) Annual evolution of mean values of degree for the same stations. Red lines in every case represent the least-squares fits of curves. The right margin: (<b>g</b>) Normalized histogram of <math display="inline"><semantics> <mi>γ</mi> </semantics></math> exponent obtained for all locations and years. Dashed line represents the theoretical value for an uncorrelated random series (<math display="inline"><semantics> <mrow> <msub> <mi>γ</mi> <mrow> <mi>u</mi> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mi>ln</mi> <mfenced> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </mfenced> </mrow> </semantics></math>). (<b>h</b>) Normalized histogram of mean degree for all locations and years.</p>
Full article ">Figure 5
<p>(<b>a</b>) Annual evolution of global clustering coefficient (<math display="inline"><semantics> <mi>C</mi> </semantics></math> along the years for Valencia station. Red line represents the least-squares fit of the curve. (<b>b</b>) The same plot for Sevilla station. (<b>c</b>) Normalized histogram of <math display="inline"><semantics> <mi>C</mi> </semantics></math> for all years and locations. (<b>d</b>) Scatter plot (grey dots) of <math display="inline"><semantics> <mi>C</mi> </semantics></math> vs. the characteristic degree exponent (<math display="inline"><semantics> <mi>γ</mi> </semantics></math>). Correlation coefficient and <span class="html-italic">p</span>-value for testing the null hypothesis that <math display="inline"><semantics> <mi>C</mi> </semantics></math> and <math display="inline"><semantics> <mi>γ</mi> </semantics></math> are not correlated. This <span class="html-italic">p</span>-value is smaller than the 95% significance level (less than 0.05), thus the correlation is statistically significant. Blue stars are window average values obtained from seven bins of equal size in <math display="inline"><semantics> <mi>γ</mi> </semantics></math> axis and error bars are standard deviations. Red line is the least-squares fit of these average values.</p>
Full article ">
18 pages, 2680 KiB  
Article
Photophysiology and Spectroscopy of Sun and Shade Leaves of Phragmites australis and the Effect on Patches of Different Densities
by Dimitris Stratoulias and Viktor R. Tóth
Remote Sens. 2020, 12(1), 200; https://doi.org/10.3390/rs12010200 - 6 Jan 2020
Cited by 10 | Viewed by 4834
Abstract
Remote sensing of vegetation has largely been revolving around the measurement of passive or active electromagnetic radiation of the top of the canopy. Nevertheless, plants hold a vertical structure and different processes and intensities take place within a plant organism depending on the [...] Read more.
Remote sensing of vegetation has largely been revolving around the measurement of passive or active electromagnetic radiation of the top of the canopy. Nevertheless, plants hold a vertical structure and different processes and intensities take place within a plant organism depending on the environmental conditions. One of the main inputs for photosynthesis is photosynthetic active radiation (PAR) and a few studies have taken into account the effect of the qualitative and quantitative changes of the available PAR within the plants canopies. Mostly large plants (trees, shrubs) are affected by this phenomena, while signs of it could be observed in dense monocultures, too. Lake Balaton is a large lake with 12 km2 dense reed stands, some of which have been suffering from reed die-back; consequently, the reed density and stress condition exhibit a vertical PAR variability within the canopy due to the structure and condition of the plants but also a horizontal variability attributed to the reedbed’s heterogeneous density. In this study we investigate the expression of photosynthetic and spectroscopic parameters in different PAR conditions. We concentrate on chlorophyll fluorescence as this is an early-stage indicator of stress manifestation in plants. We first investigate how these parameters differ across leaf samples which are exposed to a higher degree of PAR variability due to their vertical position in the reed culm (sun and shade leaves). In the second part, we concentrate on how the same parameters exhibit in reed patches of different densities. We then look into hyperspectral regions through graphs of coefficient of determination and associate the former with the physiological parameters. We report on the large variability found from measurements taken at different parts of the canopy and the association with spectral regions in the visible and near-infrared domain. We find that at low irradiance plants increase their acclimation to low light conditions. Plant density at Phragmites stands affects the vertical light attenuation and consequently the photophysiological response of basal leaves. Moreover, the hyperspectral response from the sun and shade leaves has been found to differ; charts of the coefficient of determination indicate that the spectral region around the red-edge inflection point for each case of sun and shade leaves correlate strongly with ETRmax and α. When analysing the data cumulatively, independent of their vertical position within the stand, we found correlations of R2 = 0.65 (band combination 696 and 651) and R2 = 0.61 (band combination 636 and 642) for the ETRmax and α, respectively. Full article
(This article belongs to the Special Issue Remote Sensing of Estuarine, Lagoon and Delta Environments)
Show Figures

Figure 1

Figure 1
<p>Map of the study area in Lake Balaton (inset) and field sampling points’ location in the Kerekedi Bay (main).</p>
Full article ">Figure 2
<p>Light curve of photophysiological parameters of the sun (black circles) and shade (open circles) leaves of <span class="html-italic">Phragmites australis</span>. The ± Standard Deviation (SD) (<span class="html-italic">n</span> = 15) is also indicated for each measurement.</p>
Full article ">Figure 3
<p>The maximum quantum yield for whole chain electron transport (α), the maximum electron transport capacity (ETR<sub>max</sub>) and the theoretical light saturation intensity (I<sub>k</sub>) for the sun and shade leaves of <span class="html-italic">Phragmites australis</span> in dense patches. Boxes encompass the 25% and 75% quartiles of all the data, the central solid and dashed lines represents the median and the average, bars extend to the 95% confidence limits, and dots represent outliers. For all bars <span class="html-italic">n</span> = 15. Differences were tested with the Mann-Whitney Rank Sum Test; *** indicates <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 4
<p>The quantum yield of PSII (Y(II)), photochemical (qP) and non-photochemical (qN) quenching for the sun and shade leaves of <span class="html-italic">Phragmites australis</span> in dense patches. Boxes encompass the 25% and 75% quartiles of all the data, the central solid and dashed lines represents the median and the average, bars extend to the 95% confidence limits, and dots represent outliers. For all bars <span class="html-italic">n</span> = 15. Differences were tested with the Mann–Whitney Rank Sum Test; ** indicates <span class="html-italic">p</span> &lt; 0.01, while *** indicates <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 5
<p>Mean (curve) and standard deviation (shaded area) of relative reflectance (ρ) of <span class="html-italic">Phragmites australis</span> sun (solid line) and shade (dashed line) leaves (upper inset) and the first derivative (lower inset). Each line is an average of 12 measurements.</p>
Full article ">Figure 6
<p>PCA of the spectral response features of sun (open circles) and shade (black circles) leaves of <span class="html-italic">Phragmites australis</span> in Lake Balaton. Axis 1 explains 81.3% of the variation, axis 2—14.5%.</p>
Full article ">Figure 7
<p>Coefficient of determination (R<sup>2</sup>) between the two photophysiological parameters (ETRmax (left) and α (right)) and the NDSI (Ri, Rj) measured on <span class="html-italic">Phragmites australis</span> plants at Kerekedi Bay, Lake Balaton on 14–15 August. The NDSI was calculated using the exhaustive combinations of the hyperspectrum of reflectance of two wavebands at i and j. R<sup>2</sup> corresponding to spectral pairs yielding <span class="html-italic">p</span> &gt; 0.05 (statistical insignificant) have been masked out (blank pixels).</p>
Full article ">Figure 8
<p>Boxplots of major vegetation indexes of sun and shade leaves of <span class="html-italic">Phragmites australis</span>. Data from different leaves were compared with Mann-Whitney Rank Sum Test: ns—Not significant, ** <span class="html-italic">p</span> &lt; 0.01, *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">
Back to TopTop