dash_fft/app.py at main · geonda/dash_fft · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
import dash
import dash_core_components as dcc
import dash_html_components as html
import dash_bootstrap_components as dbc

from dash.dependencies import Input, Output, State
from plotly.subplots import make_subplots

from scipy.signal import find_peaks
from pyorbital.orbital import Orbital

import plotly.graph_objects as go
import plotly.io as pio
import numpy as np
import datetime
import plotly


plotly_template = pio.templates["plotly_dark"]
pio.templates.default = "plotly_dark"


satellite = Orbital('TERRA')

app = dash.Dash(__name__, external_stylesheets=[
    "https://codepen.io/chriddyp/pen/bWLwgP.css",
    "/assets/style.css",dbc.themes.DARKLY
])


app.layout = html.Div(
    html.Div([
        html.H1("watch Terra satellite",
        style={"margin-left": "10px",
                "margin-top": "40px",
                "margin-bottom": "40px"}),

        dcc.Graph(id='live-update-graph'),

        dcc.Interval(
            id='interval-component',
            interval=10*1000, # in milliseconds
            n_intervals=0
        ),
        html.Div(id='live-update-text',
                style={"margin-left": "10px",
                        "margin-top": "40px"}),
        html.Div(id='live-update-text-2',
                style={"margin-left": "10px",
                        "margin-top": "20px"}),
    ])
)

server = app.server

@app.callback(Output('live-update-text', 'children'),
              [Input('interval-component', 'n_intervals')])
def update_metrics(n):
    lon, lat, alt = satellite.get_lonlatalt(datetime.datetime.now())
    style = {'padding': '5px', 'fontSize': '25px'}
    return [
        html.Span('Longitude: {0:.2f}'.format(lon), style=style),
        html.Span('Latitude: {0:.2f}'.format(lat), style=style),
        html.Span('Altitude: {0:0.2f}'.format(alt), style=style)
    ]

@app.callback([Output('live-update-graph', 'figure'),Output('live-update-text-2', 'children')],
              [Input('interval-component', 'n_intervals')])
def update_graph_live(n):
    satellite = Orbital('TERRA')
    data = {
        'time': [],
        'Latitude': [],
        'Longitude': [],
        'Altitude': []
    }

    # Collect some data
    for i in range(360):
        time = datetime.datetime.now() - datetime.timedelta(seconds=i*100)
        lon, lat, alt = satellite.get_lonlatalt(
            time
        )
        data['Longitude'].append(lon)
        data['Latitude'].append(lat)
        data['Altitude'].append(alt)
        data['time'].append(time)

    # Create the graph with subplots
    fig = make_subplots(rows=1, cols=2, column_widths=[0.7, 0.3])
    fig['layout']['margin'] = {
        'l': 10, 'r': 10, 'b': 10, 't': 10
    }
    fig['layout']['legend'] = {'x': 1, 'y': 1, 'xanchor': 'right'}

    fig.append_trace({
        'x': data['time'],
        'y': data['Altitude'],
        'name': 'Altitude',
        'mode': 'lines+markers',
        'type': 'scatter'

    }, 1, 1)
    fig.update_xaxes(title_text="time, hh:mm:ss", row=1, col=1)
    fig.update_yaxes(title_text="Altitude, m? ",  row=1, col=1)


    fftfreq = np.fft.fftfreq(np.linspace(0,360*100,360).shape[-1])
    fft = abs(np.imag(np.fft.fft(np.array(data['Altitude']))))
    min =  0
    max = int(len(fftfreq)/2)

    fig.append_trace({
        'x': fftfreq[min:max] ,
        'y': fft[min:max],
        'text': data['time'],
        'name': 'FFT(Altitude)',
        'mode': 'lines+markers',
        'type': 'scatter'
    }, 1, 2)
    fig.update_xaxes(title_text="frequency, ?", range=[0, 0.1], row=1, col=2)
    fig.update_yaxes(title_text="FFT(Altitude)", row=1, col=2)

    peaks, _ = find_peaks(fft[min:max], height=200)
    style = {'padding': '5px', 'fontSize': '25px'}

    f1,f2=fftfreq[peaks]

    return fig,[html.Span('First Frequency: {0:.2f}'.format(f1), style=style),
        html.Span('Second Frequency: {0:.2f}'.format(f2), style=style)]


if __name__ == "__main__":
    app.run_server(debug=True)