Visualisation de la logistique avec Python

Qu'est-ce que c'est

Un exemple pour un aperçu rapide des tendances des données de débit.

données brutes

Tiré de Recensement de la distribution du ministère des Affaires étrangères [Volume d'écoulement inter-préfectoral (par type de produit) -Weight-](http: // www. mlit.go.jp/sogoseisaku/transport/sosei_transport_fr_000074.html) est utilisé. Chaque feuille est les données de poids du flux physique par préfecture en 2010. Il y a 10 feuilles au total, eau agricole, foresterie, exploitation minière, machines métalliques, chimie, travaux légers, travaux divers, décharges et spéciaux.

Environnement d'exécution

Nous avons préparé une image Docker (tsutomu7 / graphviz) afin que vous puissiez utiliser graphviz rapidement, afin que vous puissiez exécuter ce qui suit.

bash

firefox http://localhost:8888 &
docker run -it --rm -p 8888:8888 tsutomu7/graphviz

Si vous voulez construire le vôtre, faites "conda install graphviz" et "pip install graphviz" après avoir installé Anaconda. Installez le corps principal de graphgviz avec conda et le wrapper avec pip.

Essayez-le avec Python

Lire les données

Lisez toutes les feuilles dans la variable a à la fois avec read_excel. Pour a, 0-9 est la clé et DataFrame est la valeur. a [0] est le DataFrame «total» pour toutes les industries.

python3

import numpy as np, pandas as pd
cat = 'Total Agriculture, eau et sylviculture Exploitation minière Métallurgie Machines chimiques Travaux légers Travaux divers Émissions spéciales'.split()
rng = list(range(len(cat)))
a = pd.read_excel('http://www.mlit.go.jp/sogoseisaku/transport/butsuryucensus/T9-010301.xls',
    rng, skip_footer=1, skiprows=8, header=None, index_col=0, parse_cols=np.arange(1,49))
a[0].ix[:3, :7]

Il représente une préfecture dont la ligne est From et une préfecture dont la colonne est To, et est une matrice 47 x 47.

Normaliser

Trier par industrie par ordre décroissant de volume.

python3

prefs = a[0].index.map(lambda x: x.replace('\u3000', ''))
b = [pd.DataFrame([(prefs[i], prefs[j], a[h].iloc[i, j]) for i in range(47) for j in range(47)
    if i != j], columns=['From', 'To', 'Val']).sort_values('Val', ascending=False) for h in rng]
b[0][:3]

Dessinez un diagramme avec les 5 principaux débits de chaque industrie

La figure est sortie sous la forme "fig_industry.png ". Le nombre de débit est de 1000 tonnes / an.

python3

from graphviz import Digraph
from IPython.display import display
for h, c in zip(rng, cat):
    g = Digraph(format='png')
    g.attr('graph', label=c, labelloc='t')
    g.node_attr['fontsize'] = '10'
    for _, r in b[h][:5].iterrows():
        g.edge(r.From, r.To, label='%d'%(r.Val//1000))
    g.render('fig_%s'%c)
    display(g)

c'est tout