{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Analysing the source of stressors (flow matrix)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "To calculate the source (in terms of regions and sectors) of a certain stressor or impact driven by consumption, one needs to diagonalize this stressor/impact. This section shows how to do this based on the \n", "small test mrio included in pymrio. The same procedure can be use for any other MRIO, but keep in mind that\n", "diagonalizing a stressor dramatically increases the memory need for the calculations." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Basic example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First we load the test mrio:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pymrio\n", "\n", "io = pymrio.load_test()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The test mrio includes several extensions:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['factor_inputs', 'emissions']" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "list(io.get_extensions())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For the example here, we use 'emissions' - 'emission_type1':" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
regionreg1reg2...reg5reg6
sectorfoodminingmanufactoringelectricityconstructiontradetransportotherfoodmining...transportotherfoodminingmanufactoringelectricityconstructiontradetransportother
stressorcompartment
emission_type1air1848064.80986448.09023613787.0028139100.002584141.804132656.321766987.07842090.61697937.30347378.150...4229931910773826.015777996.06420955.5113172450.056022534.04861838.51819562147046542.021632868
emission_type2water139250.4722343.295763569.18273981.55317396.511254477.81012999.12449178.0204835.4429463.944...41998417191006.34826108.11865625.112700193.0753213.72699288.3138923138765784.316782553
\n", "

2 rows × 48 columns

\n", "
" ], "text/plain": [ "region reg1 \\\n", "sector food mining manufactoring electricity \n", "stressor compartment \n", "emission_type1 air 1848064.80 986448.090 23613787.00 28139100.00 \n", "emission_type2 water 139250.47 22343.295 763569.18 273981.55 \n", "\n", "region \\\n", "sector construction trade transport other \n", "stressor compartment \n", "emission_type1 air 2584141.80 4132656.3 21766987.0 7842090.6 \n", "emission_type2 water 317396.51 1254477.8 1012999.1 2449178.0 \n", "\n", "region reg2 ... reg5 \\\n", "sector food mining ... transport other \n", "stressor compartment ... \n", "emission_type1 air 1697937.30 347378.150 ... 42299319 10773826.0 \n", "emission_type2 water 204835.44 29463.944 ... 4199841 7191006.3 \n", "\n", "region reg6 \\\n", "sector food mining manufactoring electricity \n", "stressor compartment \n", "emission_type1 air 15777996.0 6420955.5 113172450.0 56022534.0 \n", "emission_type2 water 4826108.1 1865625.1 12700193.0 753213.7 \n", "\n", "region \n", "sector construction trade transport other \n", "stressor compartment \n", "emission_type1 air 4861838.5 18195621 47046542.0 21632868 \n", "emission_type2 water 2699288.3 13892313 8765784.3 16782553 \n", "\n", "[2 rows x 48 columns]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.emissions.F" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "et1_diag = io.emissions.diag_stressor((\"emission_type1\", \"air\"), name=\"emtype1_diag\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The parameter name is optional, if not given the name is set to the stressor name + '_diag'" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The new emission matrix now looks like this:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
regionreg1reg2...reg5reg6
sectorfoodminingmanufactoringelectricityconstructiontradetransportotherfoodmining...transportotherfoodminingmanufactoringelectricityconstructiontradetransportother
regionsector
reg1food1848064.80.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
mining0.0986448.090.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
manufactoring0.00.0023613787.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
electricity0.00.000.028139100.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
construction0.00.000.00.02584141.80.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
trade0.00.000.00.00.04132656.30.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
transport0.00.000.00.00.00.021766987.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
other0.00.000.00.00.00.00.07842090.60.00.00...0.00.00.00.00.00.00.00.00.00.0
reg2food0.00.000.00.00.00.00.00.01697937.30.00...0.00.00.00.00.00.00.00.00.00.0
mining0.00.000.00.00.00.00.00.00.0347378.15...0.00.00.00.00.00.00.00.00.00.0
manufactoring0.00.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
electricity0.00.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
construction0.00.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
trade0.00.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
transport0.00.000.00.00.00.00.00.00.00.00...0.00.00.00.00.00.00.00.00.00.0
\n", "

15 rows × 48 columns

\n", "
" ], "text/plain": [ "region reg1 \\\n", "sector food mining manufactoring electricity \n", "region sector \n", "reg1 food 1848064.8 0.00 0.0 0.0 \n", " mining 0.0 986448.09 0.0 0.0 \n", " manufactoring 0.0 0.00 23613787.0 0.0 \n", " electricity 0.0 0.00 0.0 28139100.0 \n", " construction 0.0 0.00 0.0 0.0 \n", " trade 0.0 0.00 0.0 0.0 \n", " transport 0.0 0.00 0.0 0.0 \n", " other 0.0 0.00 0.0 0.0 \n", "reg2 food 0.0 0.00 0.0 0.0 \n", " mining 0.0 0.00 0.0 0.0 \n", " manufactoring 0.0 0.00 0.0 0.0 \n", " electricity 0.0 0.00 0.0 0.0 \n", " construction 0.0 0.00 0.0 0.0 \n", " trade 0.0 0.00 0.0 0.0 \n", " transport 0.0 0.00 0.0 0.0 \n", "\n", "region \\\n", "sector construction trade transport other \n", "region sector \n", "reg1 food 0.0 0.0 0.0 0.0 \n", " mining 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.0 0.0 0.0 \n", " construction 2584141.8 0.0 0.0 0.0 \n", " trade 0.0 4132656.3 0.0 0.0 \n", " transport 0.0 0.0 21766987.0 0.0 \n", " other 0.0 0.0 0.0 7842090.6 \n", "reg2 food 0.0 0.0 0.0 0.0 \n", " mining 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.0 0.0 0.0 \n", " construction 0.0 0.0 0.0 0.0 \n", " trade 0.0 0.0 0.0 0.0 \n", " transport 0.0 0.0 0.0 0.0 \n", "\n", "region reg2 ... reg5 reg6 \\\n", "sector food mining ... transport other food mining \n", "region sector ... \n", "reg1 food 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " mining 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " construction 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " trade 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " transport 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " other 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", "reg2 food 1697937.3 0.00 ... 0.0 0.0 0.0 0.0 \n", " mining 0.0 347378.15 ... 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " construction 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " trade 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", " transport 0.0 0.00 ... 0.0 0.0 0.0 0.0 \n", "\n", "region \\\n", "sector manufactoring electricity construction trade transport \n", "region sector \n", "reg1 food 0.0 0.0 0.0 0.0 0.0 \n", " mining 0.0 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.0 0.0 0.0 0.0 \n", " construction 0.0 0.0 0.0 0.0 0.0 \n", " trade 0.0 0.0 0.0 0.0 0.0 \n", " transport 0.0 0.0 0.0 0.0 0.0 \n", " other 0.0 0.0 0.0 0.0 0.0 \n", "reg2 food 0.0 0.0 0.0 0.0 0.0 \n", " mining 0.0 0.0 0.0 0.0 0.0 \n", " manufactoring 0.0 0.0 0.0 0.0 0.0 \n", " electricity 0.0 0.0 0.0 0.0 0.0 \n", " construction 0.0 0.0 0.0 0.0 0.0 \n", " trade 0.0 0.0 0.0 0.0 0.0 \n", " transport 0.0 0.0 0.0 0.0 0.0 \n", "\n", "region \n", "sector other \n", "region sector \n", "reg1 food 0.0 \n", " mining 0.0 \n", " manufactoring 0.0 \n", " electricity 0.0 \n", " construction 0.0 \n", " trade 0.0 \n", " transport 0.0 \n", " other 0.0 \n", "reg2 food 0.0 \n", " mining 0.0 \n", " manufactoring 0.0 \n", " electricity 0.0 \n", " construction 0.0 \n", " trade 0.0 \n", " transport 0.0 \n", "\n", "[15 rows x 48 columns]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "et1_diag.F.head(15)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And can be connected back to the system with:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "io.et1_diag = et1_diag" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finally we can calulate the all stressor accounts with:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.calc_all()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This results in a square footprint matrix. In this matrix, every column respresents the amount of stressor occuring in each region - sector driven by the consumption stated in the column header. Conversly, each row states where the stressor impacts occuring in the row are distributed due (from where they are driven)." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
regionreg1reg2...reg5reg6
sectorfoodminingmanufactoringelectricityconstructiontradetransportotherfoodmining...transportotherfoodminingmanufactoringelectricityconstructiontradetransportother
regionsector
reg1food609347.99874734.9632231.987631e+057.678755e+022.873371e+035.603158e+041.448778e+035.225312e+044.826852e+040.453097...74.449012222.95328925525.5163926.5576001.382548e+0538.345126408.9378914.425472e+0463.9402154.876200e+02
mining2527.44944161271.2496391.232716e+055.406781e+041.632967e+051.459774e+044.916876e+034.975184e+044.311608e+02288.903690...261.326848829.871069268.4067931990.8400547.400739e+046327.871035299.9166991.732632e+04242.8744444.038056e+03
manufactoring1199.04153038.2366794.686837e+068.108636e+021.816229e+047.713610e+032.825229e+031.634290e+043.205155e+022.886944...424.008556517.344586145.90043972.1950393.424703e+0694.046696129.1270597.826193e+03241.1851869.015651e+02
electricity148505.09190212764.7842971.519466e+061.167804e+075.265922e+051.307806e+064.851957e+054.308489e+061.908031e+04240.730963...9294.38155128937.5690167947.002552476.1068971.072302e+0638647.049939495.1593809.712017e+053850.1338269.804940e+03
construction49.0184796.0534594.019302e+023.081457e+022.568908e+067.291250e+025.853107e+029.718253e+034.408230e+000.037972...1.6797914.4610352.2115300.2139512.898921e+021.68297712.6673705.303592e+020.9889793.859558e+00
trade138.0413553.1395961.880042e+038.852940e+011.420349e+032.390871e+064.371531e+022.383210e+032.637857e+010.095853...36.43364150.39205142.0989451.1416821.240785e+034.98978814.1926231.726377e+0645.6334644.276151e+01
transport521.216924122.5049681.585636e+041.428149e+036.467383e+032.741408e+041.018383e+074.388313e+049.410163e+012.099087...2854.767892222.900909568.39399655.9169121.130829e+04285.596244127.3880072.198931e+048734.7269731.309008e+03
other537.06267924.6880207.752597e+031.170368e+038.587431e+031.322798e+044.448616e+037.516913e+066.837072e+010.767113...56.726997435.70335936.6409533.8566485.346664e+0320.75094515.1524251.029493e+0458.2780761.356728e+03
reg2food234.8701080.0412821.213308e+033.241119e+009.719899e+001.004924e+014.822982e-011.281313e+011.694248e+060.041642...3.0264720.47502480.1392060.0251412.183449e+020.4725227.0111101.605455e+020.4551753.092008e+00
mining215.5627621690.1866705.892279e+041.862733e+037.746514e+027.570607e+021.582065e+022.302825e+031.787512e+0310287.905967...200.67883058.08336766.1856181351.8093331.430140e+042650.00818442.2300728.909866e+0379.9819463.551366e+03
manufactoring75.6213464.2294637.185296e+066.763336e+017.171199e+024.429766e+021.914508e+021.076869e+039.941057e+025.361056...352.701574125.16558844.65983819.1575941.140699e+0619.22046150.9827843.009360e+0374.8371902.460570e+02
electricity4.9121834.3922706.448813e+031.947867e+021.552407e+012.864877e+011.010699e+019.682693e+011.143809e+0323.763711...225.0203232.6141120.7405253.1624661.079151e+038.6011510.2546392.214440e+031.1863021.102287e+01
construction0.1792740.2014904.711030e+028.221412e-011.045043e+021.884894e+001.624084e+001.924495e+011.001309e+021.098256...130.3900274.8040880.2302080.1526997.708346e+010.54543347.7502323.449526e+021.0235205.208169e+00
trade9.4878020.4428512.547850e+035.931388e+003.093165e+012.414790e+028.294245e+005.515915e+013.073687e+021.313536...159.34029669.42112025.2782970.7832725.438996e+022.7954868.8522731.174931e+0620.0994092.543311e+01
transport30.16791711.6917031.095783e+047.704484e+012.369344e+021.171461e+034.962901e+031.299130e+039.083313e+0215.109529...789093.398671117.822465101.0208957.8168962.491430e+0351.07975628.6095828.615321e+032234.5858454.125342e+02
other4.7101520.9996563.487999e+031.504331e+016.607788e+018.958452e+013.266330e+011.104173e+032.621005e+023.670327...376.984093213.7006967.1789141.5231516.818453e+027.5761987.3806013.072026e+0337.6951427.442119e+02
reg3food79.4879950.0124202.707179e+033.212251e-011.118965e+008.545640e+003.848642e-012.342262e+016.321926e+010.001120...0.6603781.09557996.4270890.0604301.350698e+030.1365240.2537325.444486e+020.3302141.007920e+02
mining1.6608269.2831442.805174e+034.683637e+003.721404e+001.903501e+011.984584e+001.130052e+021.043411e+001.482323...1.5792974.7600431.0285779.0185131.370585e+0319.2088360.6092437.625285e+021.2349526.075141e+02
manufactoring256.95078712.4969661.951101e+071.576456e+021.369020e+039.774252e+024.651207e+029.975075e+032.274549e+022.093210...856.5214231696.838781304.748914109.5397169.493385e+0671.802017260.5840836.010179e+04587.4916713.817698e+04
electricity346.61894475.1661703.907669e+068.377082e+031.955704e+033.618932e+031.930657e+035.954611e+059.697584e+0233.309324...1833.1413856869.6691391750.8057011639.9271361.908654e+0644575.821242628.0610596.364065e+064888.4724243.661473e+06
\n", "

20 rows × 48 columns

\n", "
" ], "text/plain": [ "region reg1 \\\n", "sector food mining manufactoring electricity \n", "region sector \n", "reg1 food 609347.998747 34.963223 1.987631e+05 7.678755e+02 \n", " mining 2527.449441 61271.249639 1.232716e+05 5.406781e+04 \n", " manufactoring 1199.041530 38.236679 4.686837e+06 8.108636e+02 \n", " electricity 148505.091902 12764.784297 1.519466e+06 1.167804e+07 \n", " construction 49.018479 6.053459 4.019302e+02 3.081457e+02 \n", " trade 138.041355 3.139596 1.880042e+03 8.852940e+01 \n", " transport 521.216924 122.504968 1.585636e+04 1.428149e+03 \n", " other 537.062679 24.688020 7.752597e+03 1.170368e+03 \n", "reg2 food 234.870108 0.041282 1.213308e+03 3.241119e+00 \n", " mining 215.562762 1690.186670 5.892279e+04 1.862733e+03 \n", " manufactoring 75.621346 4.229463 7.185296e+06 6.763336e+01 \n", " electricity 4.912183 4.392270 6.448813e+03 1.947867e+02 \n", " construction 0.179274 0.201490 4.711030e+02 8.221412e-01 \n", " trade 9.487802 0.442851 2.547850e+03 5.931388e+00 \n", " transport 30.167917 11.691703 1.095783e+04 7.704484e+01 \n", " other 4.710152 0.999656 3.487999e+03 1.504331e+01 \n", "reg3 food 79.487995 0.012420 2.707179e+03 3.212251e-01 \n", " mining 1.660826 9.283144 2.805174e+03 4.683637e+00 \n", " manufactoring 256.950787 12.496966 1.951101e+07 1.576456e+02 \n", " electricity 346.618944 75.166170 3.907669e+06 8.377082e+03 \n", "\n", "region \\\n", "sector construction trade transport other \n", "region sector \n", "reg1 food 2.873371e+03 5.603158e+04 1.448778e+03 5.225312e+04 \n", " mining 1.632967e+05 1.459774e+04 4.916876e+03 4.975184e+04 \n", " manufactoring 1.816229e+04 7.713610e+03 2.825229e+03 1.634290e+04 \n", " electricity 5.265922e+05 1.307806e+06 4.851957e+05 4.308489e+06 \n", " construction 2.568908e+06 7.291250e+02 5.853107e+02 9.718253e+03 \n", " trade 1.420349e+03 2.390871e+06 4.371531e+02 2.383210e+03 \n", " transport 6.467383e+03 2.741408e+04 1.018383e+07 4.388313e+04 \n", " other 8.587431e+03 1.322798e+04 4.448616e+03 7.516913e+06 \n", "reg2 food 9.719899e+00 1.004924e+01 4.822982e-01 1.281313e+01 \n", " mining 7.746514e+02 7.570607e+02 1.582065e+02 2.302825e+03 \n", " manufactoring 7.171199e+02 4.429766e+02 1.914508e+02 1.076869e+03 \n", " electricity 1.552407e+01 2.864877e+01 1.010699e+01 9.682693e+01 \n", " construction 1.045043e+02 1.884894e+00 1.624084e+00 1.924495e+01 \n", " trade 3.093165e+01 2.414790e+02 8.294245e+00 5.515915e+01 \n", " transport 2.369344e+02 1.171461e+03 4.962901e+03 1.299130e+03 \n", " other 6.607788e+01 8.958452e+01 3.266330e+01 1.104173e+03 \n", "reg3 food 1.118965e+00 8.545640e+00 3.848642e-01 2.342262e+01 \n", " mining 3.721404e+00 1.903501e+01 1.984584e+00 1.130052e+02 \n", " manufactoring 1.369020e+03 9.774252e+02 4.651207e+02 9.975075e+03 \n", " electricity 1.955704e+03 3.618932e+03 1.930657e+03 5.954611e+05 \n", "\n", "region reg2 ... reg5 \\\n", "sector food mining ... transport \n", "region sector ... \n", "reg1 food 4.826852e+04 0.453097 ... 74.449012 \n", " mining 4.311608e+02 288.903690 ... 261.326848 \n", " manufactoring 3.205155e+02 2.886944 ... 424.008556 \n", " electricity 1.908031e+04 240.730963 ... 9294.381551 \n", " construction 4.408230e+00 0.037972 ... 1.679791 \n", " trade 2.637857e+01 0.095853 ... 36.433641 \n", " transport 9.410163e+01 2.099087 ... 2854.767892 \n", " other 6.837072e+01 0.767113 ... 56.726997 \n", "reg2 food 1.694248e+06 0.041642 ... 3.026472 \n", " mining 1.787512e+03 10287.905967 ... 200.678830 \n", " manufactoring 9.941057e+02 5.361056 ... 352.701574 \n", " electricity 1.143809e+03 23.763711 ... 225.020323 \n", " construction 1.001309e+02 1.098256 ... 130.390027 \n", " trade 3.073687e+02 1.313536 ... 159.340296 \n", " transport 9.083313e+02 15.109529 ... 789093.398671 \n", " other 2.621005e+02 3.670327 ... 376.984093 \n", "reg3 food 6.321926e+01 0.001120 ... 0.660378 \n", " mining 1.043411e+00 1.482323 ... 1.579297 \n", " manufactoring 2.274549e+02 2.093210 ... 856.521423 \n", " electricity 9.697584e+02 33.309324 ... 1833.141385 \n", "\n", "region reg6 \\\n", "sector other food mining manufactoring \n", "region sector \n", "reg1 food 222.953289 25525.516392 6.557600 1.382548e+05 \n", " mining 829.871069 268.406793 1990.840054 7.400739e+04 \n", " manufactoring 517.344586 145.900439 72.195039 3.424703e+06 \n", " electricity 28937.569016 7947.002552 476.106897 1.072302e+06 \n", " construction 4.461035 2.211530 0.213951 2.898921e+02 \n", " trade 50.392051 42.098945 1.141682 1.240785e+03 \n", " transport 222.900909 568.393996 55.916912 1.130829e+04 \n", " other 435.703359 36.640953 3.856648 5.346664e+03 \n", "reg2 food 0.475024 80.139206 0.025141 2.183449e+02 \n", " mining 58.083367 66.185618 1351.809333 1.430140e+04 \n", " manufactoring 125.165588 44.659838 19.157594 1.140699e+06 \n", " electricity 2.614112 0.740525 3.162466 1.079151e+03 \n", " construction 4.804088 0.230208 0.152699 7.708346e+01 \n", " trade 69.421120 25.278297 0.783272 5.438996e+02 \n", " transport 117.822465 101.020895 7.816896 2.491430e+03 \n", " other 213.700696 7.178914 1.523151 6.818453e+02 \n", "reg3 food 1.095579 96.427089 0.060430 1.350698e+03 \n", " mining 4.760043 1.028577 9.018513 1.370585e+03 \n", " manufactoring 1696.838781 304.748914 109.539716 9.493385e+06 \n", " electricity 6869.669139 1750.805701 1639.927136 1.908654e+06 \n", "\n", "region \\\n", "sector electricity construction trade transport \n", "region sector \n", "reg1 food 38.345126 408.937891 4.425472e+04 63.940215 \n", " mining 6327.871035 299.916699 1.732632e+04 242.874444 \n", " manufactoring 94.046696 129.127059 7.826193e+03 241.185186 \n", " electricity 38647.049939 495.159380 9.712017e+05 3850.133826 \n", " construction 1.682977 12.667370 5.303592e+02 0.988979 \n", " trade 4.989788 14.192623 1.726377e+06 45.633464 \n", " transport 285.596244 127.388007 2.198931e+04 8734.726973 \n", " other 20.750945 15.152425 1.029493e+04 58.278076 \n", "reg2 food 0.472522 7.011110 1.605455e+02 0.455175 \n", " mining 2650.008184 42.230072 8.909866e+03 79.981946 \n", " manufactoring 19.220461 50.982784 3.009360e+03 74.837190 \n", " electricity 8.601151 0.254639 2.214440e+03 1.186302 \n", " construction 0.545433 47.750232 3.449526e+02 1.023520 \n", " trade 2.795486 8.852273 1.174931e+06 20.099409 \n", " transport 51.079756 28.609582 8.615321e+03 2234.585845 \n", " other 7.576198 7.380601 3.072026e+03 37.695142 \n", "reg3 food 0.136524 0.253732 5.444486e+02 0.330214 \n", " mining 19.208836 0.609243 7.625285e+02 1.234952 \n", " manufactoring 71.802017 260.584083 6.010179e+04 587.491671 \n", " electricity 44575.821242 628.061059 6.364065e+06 4888.472424 \n", "\n", "region \n", "sector other \n", "region sector \n", "reg1 food 4.876200e+02 \n", " mining 4.038056e+03 \n", " manufactoring 9.015651e+02 \n", " electricity 9.804940e+03 \n", " construction 3.859558e+00 \n", " trade 4.276151e+01 \n", " transport 1.309008e+03 \n", " other 1.356728e+03 \n", "reg2 food 3.092008e+00 \n", " mining 3.551366e+03 \n", " manufactoring 2.460570e+02 \n", " electricity 1.102287e+01 \n", " construction 5.208169e+00 \n", " trade 2.543311e+01 \n", " transport 4.125342e+02 \n", " other 7.442119e+02 \n", "reg3 food 1.007920e+02 \n", " mining 6.075141e+02 \n", " manufactoring 3.817698e+04 \n", " electricity 3.661473e+06 \n", "\n", "[20 rows x 48 columns]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.et1_diag.D_cba.head(20)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The total footprints of a region - sector are given by summing the footprints along rows:" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "sector\n", "food 2.056183e+06\n", "mining 1.794235e+05\n", "manufactoring 9.749300e+07\n", "electricity 1.188759e+07\n", "construction 3.342906e+06\n", "trade 3.885884e+06\n", "transport 1.075027e+07\n", "other 1.582152e+07\n", "dtype: float64" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.et1_diag.D_cba.sum(axis=0).reg1" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
sectorfoodminingmanufactoringelectricityconstructiontradetransportother
stressorcompartment
emission_type1air2.056183e+06179423.5358939.749300e+071.188759e+073.342906e+063.885884e+061.075027e+071.582152e+07
emission_type2water2.423103e+0525278.1920861.671240e+071.371303e+053.468292e+057.766205e+054.999628e+058.480505e+06
\n", "
" ], "text/plain": [ "sector food mining manufactoring \\\n", "stressor compartment \n", "emission_type1 air 2.056183e+06 179423.535893 9.749300e+07 \n", "emission_type2 water 2.423103e+05 25278.192086 1.671240e+07 \n", "\n", "sector electricity construction trade \\\n", "stressor compartment \n", "emission_type1 air 1.188759e+07 3.342906e+06 3.885884e+06 \n", "emission_type2 water 1.371303e+05 3.468292e+05 7.766205e+05 \n", "\n", "sector transport other \n", "stressor compartment \n", "emission_type1 air 1.075027e+07 1.582152e+07 \n", "emission_type2 water 4.999628e+05 8.480505e+06 " ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.emissions.D_cba.reg1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The total stressor in a sector corresponds to the sum of the columns:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "sector\n", "food 1848064.80\n", "mining 986448.09\n", "manufactoring 23613787.00\n", "electricity 28139100.00\n", "construction 2584141.80\n", "trade 4132656.30\n", "transport 21766987.00\n", "other 7842090.60\n", "dtype: float64" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.et1_diag.D_cba.sum(axis=1).reg1" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
sectorfoodminingmanufactoringelectricityconstructiontradetransportother
stressorcompartment
emission_type1air1848064.80986448.09023613787.0028139100.002584141.804132656.321766987.07842090.6
emission_type2water139250.4722343.295763569.18273981.55317396.511254477.81012999.12449178.0
\n", "
" ], "text/plain": [ "sector food mining manufactoring \\\n", "stressor compartment \n", "emission_type1 air 1848064.80 986448.090 23613787.00 \n", "emission_type2 water 139250.47 22343.295 763569.18 \n", "\n", "sector electricity construction trade transport \\\n", "stressor compartment \n", "emission_type1 air 28139100.00 2584141.80 4132656.3 21766987.0 \n", "emission_type2 water 273981.55 317396.51 1254477.8 1012999.1 \n", "\n", "sector other \n", "stressor compartment \n", "emission_type1 air 7842090.6 \n", "emission_type2 water 2449178.0 " ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.emissions.F.reg1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Aggregation of source footprints" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If only one specific aspect of the source is of interest for the analysis, the footprint matrix can easily be aggregated with the standard pandas groupby function. \n", "\n", "For example, to aggregate to the source region of stressor, do:" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
regionreg1reg2...reg5reg6
sectorfoodminingmanufactoringelectricityconstructiontradetransportotherfoodmining...transportotherfoodminingmanufactoringelectricityconstructiontradetransportother
region
reg17.628249e+0574265.6198826.554229e+061.173668e+073.296308e+063.818391e+061.068369e+071.199973e+076.829377e+04535.974719...1.300377e+043.122120e+043.453617e+042.606829e+034.727453e+064.542033e+041.502541e+032.799800e+061.323776e+041.794454e+04
reg25.755115e+021712.1853847.269345e+062.227236e+031.955463e+032.743145e+035.365729e+035.967041e+031.699751e+0610338.264024...7.905415e+055.920865e+023.254335e+021.384431e+031.160092e+062.740299e+031.930713e+021.201258e+062.449865e+034.998925e+03
reg31.054578e+03215.6545912.382082e+079.500254e+037.044809e+031.285847e+042.265449e+042.844767e+061.740953e+0354.246936...2.026441e+043.336623e+045.910854e+032.129253e+031.160211e+074.640647e+042.506927e+032.295789e+074.824482e+041.756144e+07
reg41.147382e+032323.7920841.219510e+073.931814e+036.717898e+032.272241e+032.088695e+031.054927e+048.986664e+025812.809417...1.818417e+031.930517e+032.402317e+031.257375e+068.856003e+067.878447e+032.121123e+036.469320e+065.546359e+031.756339e+04
reg51.283812e+066596.7139071.588478e+071.642030e+041.234768e+041.598377e+042.387931e+042.906402e+041.828518e+041121.255607...4.126106e+071.129572e+074.599517e+047.759982e+031.519859e+071.879346e+048.076635e+037.701012e+064.094084e+042.125918e+04
reg66.769053e+0394309.5700453.176873e+071.188346e+051.853198e+043.363451e+041.259087e+049.314423e+054.368062e+031283.054207...8.360619e+032.378125e+041.508319e+077.406276e+042.990651e+073.671043e+071.822296e+061.286404e+064.794367e+071.839977e+07
\n", "

6 rows × 48 columns

\n", "
" ], "text/plain": [ "region reg1 \\\n", "sector food mining manufactoring electricity construction \n", "region \n", "reg1 7.628249e+05 74265.619882 6.554229e+06 1.173668e+07 3.296308e+06 \n", "reg2 5.755115e+02 1712.185384 7.269345e+06 2.227236e+03 1.955463e+03 \n", "reg3 1.054578e+03 215.654591 2.382082e+07 9.500254e+03 7.044809e+03 \n", "reg4 1.147382e+03 2323.792084 1.219510e+07 3.931814e+03 6.717898e+03 \n", "reg5 1.283812e+06 6596.713907 1.588478e+07 1.642030e+04 1.234768e+04 \n", "reg6 6.769053e+03 94309.570045 3.176873e+07 1.188346e+05 1.853198e+04 \n", "\n", "region reg2 \\\n", "sector trade transport other food mining \n", "region \n", "reg1 3.818391e+06 1.068369e+07 1.199973e+07 6.829377e+04 535.974719 \n", "reg2 2.743145e+03 5.365729e+03 5.967041e+03 1.699751e+06 10338.264024 \n", "reg3 1.285847e+04 2.265449e+04 2.844767e+06 1.740953e+03 54.246936 \n", "reg4 2.272241e+03 2.088695e+03 1.054927e+04 8.986664e+02 5812.809417 \n", "reg5 1.598377e+04 2.387931e+04 2.906402e+04 1.828518e+04 1121.255607 \n", "reg6 3.363451e+04 1.259087e+04 9.314423e+05 4.368062e+03 1283.054207 \n", "\n", "region ... reg5 reg6 \\\n", "sector ... transport other food mining \n", "region ... \n", "reg1 ... 1.300377e+04 3.122120e+04 3.453617e+04 2.606829e+03 \n", "reg2 ... 7.905415e+05 5.920865e+02 3.254335e+02 1.384431e+03 \n", "reg3 ... 2.026441e+04 3.336623e+04 5.910854e+03 2.129253e+03 \n", "reg4 ... 1.818417e+03 1.930517e+03 2.402317e+03 1.257375e+06 \n", "reg5 ... 4.126106e+07 1.129572e+07 4.599517e+04 7.759982e+03 \n", "reg6 ... 8.360619e+03 2.378125e+04 1.508319e+07 7.406276e+04 \n", "\n", "region \\\n", "sector manufactoring electricity construction trade transport \n", "region \n", "reg1 4.727453e+06 4.542033e+04 1.502541e+03 2.799800e+06 1.323776e+04 \n", "reg2 1.160092e+06 2.740299e+03 1.930713e+02 1.201258e+06 2.449865e+03 \n", "reg3 1.160211e+07 4.640647e+04 2.506927e+03 2.295789e+07 4.824482e+04 \n", "reg4 8.856003e+06 7.878447e+03 2.121123e+03 6.469320e+06 5.546359e+03 \n", "reg5 1.519859e+07 1.879346e+04 8.076635e+03 7.701012e+06 4.094084e+04 \n", "reg6 2.990651e+07 3.671043e+07 1.822296e+06 1.286404e+06 4.794367e+07 \n", "\n", "region \n", "sector other \n", "region \n", "reg1 1.794454e+04 \n", "reg2 4.998925e+03 \n", "reg3 1.756144e+07 \n", "reg4 1.756339e+04 \n", "reg5 2.125918e+04 \n", "reg6 1.839977e+07 \n", "\n", "[6 rows x 48 columns]" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.et1_diag.D_cba.groupby(level=\"region\", axis=0).sum()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In addition, the [aggregation function](../notebooks/aggregation_examples.ipynb) of pymrio also work on the diagonalized footprints. Here as example together with the [country converter coco](https://github.com/IndEcol/country_converter):" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import country_converter as coco\n", "\n", "io.aggregate(\n", " region_agg=coco.agg_conc(\n", " original_countries=io.get_regions(),\n", " aggregates={\n", " \"reg1\": \"World Region A\",\n", " \"reg2\": \"World Region A\",\n", " \"reg3\": \"World Region A\",\n", " },\n", " missing_countries=\"World Region B\",\n", " )\n", ")" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
regionWorld Region AWorld Region B
sectorfoodminingmanufactoringelectricityconstructiontradetransportotherfoodminingmanufactoringelectricityconstructiontradetransportother
regionsector
World Region Afood6.413682e+065.952471e+016.070321e+059.326086e+023.306995e+035.987980e+043.514385e+035.762977e+043.437446e+047.789050e+003.648174e+056.373523e+011.339737e+034.505792e+041.557188e+021.163124e+03
mining6.832129e+032.421509e+064.487266e+051.936672e+051.984123e+052.126768e+041.676197e+047.595762e+044.373682e+023.805693e+032.419141e+051.079248e+041.073022e+042.708989e+048.090257e+021.908491e+04
manufactoring1.575255e+044.337974e+035.857218e+075.217160e+031.108166e+051.787260e+044.264953e+048.572782e+041.196817e+032.389005e+025.612155e+071.442696e+036.427781e+037.159629e+043.128223e+039.631846e+04
electricity1.148908e+068.329886e+051.095357e+075.960881e+071.529502e+061.771262e+063.062426e+069.812779e+061.717518e+042.384122e+031.167721e+073.179723e+051.151092e+047.346021e+062.360035e+048.544805e+06
construction1.287094e+033.530581e+037.855368e+034.364648e+031.082799e+071.811093e+031.193404e+044.796826e+048.117404e+001.278994e+009.064738e+032.977577e+011.183350e+031.057865e+041.704971e+025.278625e+04
trade1.302177e+043.358650e+031.320568e+055.581957e+037.080932e+044.827615e+063.359931e+046.189364e+041.772741e+033.095513e+011.747324e+055.736669e+022.618051e+031.843441e+074.414173e+039.440231e+04
transport1.661673e+041.013917e+042.026243e+051.186026e+049.463786e+046.605838e+046.900978e+072.981319e+053.570924e+032.342901e+022.223755e+052.064082e+034.978304e+033.583977e+058.792139e+053.666801e+05
other4.973365e+045.119490e+044.157135e+053.141636e+042.451423e+057.347058e+043.463245e+053.192544e+071.468825e+032.988605e+025.284564e+053.388346e+037.111777e+036.871237e+058.921976e+033.158392e+07
World Region Bfood1.331074e+065.840707e+011.158034e+067.157090e+022.277605e+032.977910e+041.009061e+036.741954e+032.366136e+071.599984e+027.393735e+053.069128e+032.404678e+041.717301e+053.012820e+047.478749e+04
mining1.120813e+041.223669e+056.244537e+063.478494e+053.383509e+041.050007e+053.107988e+041.118387e+053.852121e+041.612039e+066.627586e+061.172268e+065.475534e+051.556577e+058.352660e+042.849191e+05
manufactoring1.165415e+041.112145e+031.403366e+084.253396e+031.644812e+044.518924e+041.350120e+043.504348e+047.406372e+043.204186e+031.254308e+083.173150e+042.217880e+058.865177e+041.225773e+051.739262e+05
electricity6.624049e+032.433130e+045.205461e+061.790907e+051.709844e+042.343087e+061.540263e+048.184191e+057.886875e+052.316855e+046.039318e+061.225545e+084.700901e+052.298746e+051.006105e+062.680870e+06
construction2.441345e+033.338858e+026.498032e+041.067936e+032.843476e+041.094280e+042.562950e+031.372537e+041.229142e+042.193007e+034.813644e+043.446390e+041.097151e+073.431794e+046.586141e+043.153222e+05
trade7.186649e+021.172395e+028.340493e+043.915868e+027.298241e+021.845264e+071.527523e+031.037099e+041.802717e+044.401980e+029.207769e+044.792922e+033.300157e+041.960537e+074.424292e+044.933951e+04
transport1.118710e+041.051599e+033.382896e+053.512873e+034.336422e+036.536434e+043.522339e+062.451637e+043.616496e+042.977747e+032.337329e+053.716876e+049.823204e+042.505787e+051.102146e+082.615115e+05
other4.082899e+021.690002e+024.761138e+045.086643e+027.504197e+022.293492e+042.199852e+034.906621e+065.528324e+038.635489e+024.987069e+048.685131e+032.601768e+043.556939e+043.218539e+044.074021e+07
\n", "
" ], "text/plain": [ "region World Region A \\\n", "sector food mining manufactoring \n", "region sector \n", "World Region A food 6.413682e+06 5.952471e+01 6.070321e+05 \n", " mining 6.832129e+03 2.421509e+06 4.487266e+05 \n", " manufactoring 1.575255e+04 4.337974e+03 5.857218e+07 \n", " electricity 1.148908e+06 8.329886e+05 1.095357e+07 \n", " construction 1.287094e+03 3.530581e+03 7.855368e+03 \n", " trade 1.302177e+04 3.358650e+03 1.320568e+05 \n", " transport 1.661673e+04 1.013917e+04 2.026243e+05 \n", " other 4.973365e+04 5.119490e+04 4.157135e+05 \n", "World Region B food 1.331074e+06 5.840707e+01 1.158034e+06 \n", " mining 1.120813e+04 1.223669e+05 6.244537e+06 \n", " manufactoring 1.165415e+04 1.112145e+03 1.403366e+08 \n", " electricity 6.624049e+03 2.433130e+04 5.205461e+06 \n", " construction 2.441345e+03 3.338858e+02 6.498032e+04 \n", " trade 7.186649e+02 1.172395e+02 8.340493e+04 \n", " transport 1.118710e+04 1.051599e+03 3.382896e+05 \n", " other 4.082899e+02 1.690002e+02 4.761138e+04 \n", "\n", "region \\\n", "sector electricity construction trade \n", "region sector \n", "World Region A food 9.326086e+02 3.306995e+03 5.987980e+04 \n", " mining 1.936672e+05 1.984123e+05 2.126768e+04 \n", " manufactoring 5.217160e+03 1.108166e+05 1.787260e+04 \n", " electricity 5.960881e+07 1.529502e+06 1.771262e+06 \n", " construction 4.364648e+03 1.082799e+07 1.811093e+03 \n", " trade 5.581957e+03 7.080932e+04 4.827615e+06 \n", " transport 1.186026e+04 9.463786e+04 6.605838e+04 \n", " other 3.141636e+04 2.451423e+05 7.347058e+04 \n", "World Region B food 7.157090e+02 2.277605e+03 2.977910e+04 \n", " mining 3.478494e+05 3.383509e+04 1.050007e+05 \n", " manufactoring 4.253396e+03 1.644812e+04 4.518924e+04 \n", " electricity 1.790907e+05 1.709844e+04 2.343087e+06 \n", " construction 1.067936e+03 2.843476e+04 1.094280e+04 \n", " trade 3.915868e+02 7.298241e+02 1.845264e+07 \n", " transport 3.512873e+03 4.336422e+03 6.536434e+04 \n", " other 5.086643e+02 7.504197e+02 2.293492e+04 \n", "\n", "region World Region B \\\n", "sector transport other food \n", "region sector \n", "World Region A food 3.514385e+03 5.762977e+04 3.437446e+04 \n", " mining 1.676197e+04 7.595762e+04 4.373682e+02 \n", " manufactoring 4.264953e+04 8.572782e+04 1.196817e+03 \n", " electricity 3.062426e+06 9.812779e+06 1.717518e+04 \n", " construction 1.193404e+04 4.796826e+04 8.117404e+00 \n", " trade 3.359931e+04 6.189364e+04 1.772741e+03 \n", " transport 6.900978e+07 2.981319e+05 3.570924e+03 \n", " other 3.463245e+05 3.192544e+07 1.468825e+03 \n", "World Region B food 1.009061e+03 6.741954e+03 2.366136e+07 \n", " mining 3.107988e+04 1.118387e+05 3.852121e+04 \n", " manufactoring 1.350120e+04 3.504348e+04 7.406372e+04 \n", " electricity 1.540263e+04 8.184191e+05 7.886875e+05 \n", " construction 2.562950e+03 1.372537e+04 1.229142e+04 \n", " trade 1.527523e+03 1.037099e+04 1.802717e+04 \n", " transport 3.522339e+06 2.451637e+04 3.616496e+04 \n", " other 2.199852e+03 4.906621e+06 5.528324e+03 \n", "\n", "region \\\n", "sector mining manufactoring electricity \n", "region sector \n", "World Region A food 7.789050e+00 3.648174e+05 6.373523e+01 \n", " mining 3.805693e+03 2.419141e+05 1.079248e+04 \n", " manufactoring 2.389005e+02 5.612155e+07 1.442696e+03 \n", " electricity 2.384122e+03 1.167721e+07 3.179723e+05 \n", " construction 1.278994e+00 9.064738e+03 2.977577e+01 \n", " trade 3.095513e+01 1.747324e+05 5.736669e+02 \n", " transport 2.342901e+02 2.223755e+05 2.064082e+03 \n", " other 2.988605e+02 5.284564e+05 3.388346e+03 \n", "World Region B food 1.599984e+02 7.393735e+05 3.069128e+03 \n", " mining 1.612039e+06 6.627586e+06 1.172268e+06 \n", " manufactoring 3.204186e+03 1.254308e+08 3.173150e+04 \n", " electricity 2.316855e+04 6.039318e+06 1.225545e+08 \n", " construction 2.193007e+03 4.813644e+04 3.446390e+04 \n", " trade 4.401980e+02 9.207769e+04 4.792922e+03 \n", " transport 2.977747e+03 2.337329e+05 3.716876e+04 \n", " other 8.635489e+02 4.987069e+04 8.685131e+03 \n", "\n", "region \\\n", "sector construction trade transport \n", "region sector \n", "World Region A food 1.339737e+03 4.505792e+04 1.557188e+02 \n", " mining 1.073022e+04 2.708989e+04 8.090257e+02 \n", " manufactoring 6.427781e+03 7.159629e+04 3.128223e+03 \n", " electricity 1.151092e+04 7.346021e+06 2.360035e+04 \n", " construction 1.183350e+03 1.057865e+04 1.704971e+02 \n", " trade 2.618051e+03 1.843441e+07 4.414173e+03 \n", " transport 4.978304e+03 3.583977e+05 8.792139e+05 \n", " other 7.111777e+03 6.871237e+05 8.921976e+03 \n", "World Region B food 2.404678e+04 1.717301e+05 3.012820e+04 \n", " mining 5.475534e+05 1.556577e+05 8.352660e+04 \n", " manufactoring 2.217880e+05 8.865177e+04 1.225773e+05 \n", " electricity 4.700901e+05 2.298746e+05 1.006105e+06 \n", " construction 1.097151e+07 3.431794e+04 6.586141e+04 \n", " trade 3.300157e+04 1.960537e+07 4.424292e+04 \n", " transport 9.823204e+04 2.505787e+05 1.102146e+08 \n", " other 2.601768e+04 3.556939e+04 3.218539e+04 \n", "\n", "region \n", "sector other \n", "region sector \n", "World Region A food 1.163124e+03 \n", " mining 1.908491e+04 \n", " manufactoring 9.631846e+04 \n", " electricity 8.544805e+06 \n", " construction 5.278625e+04 \n", " trade 9.440231e+04 \n", " transport 3.666801e+05 \n", " other 3.158392e+07 \n", "World Region B food 7.478749e+04 \n", " mining 2.849191e+05 \n", " manufactoring 1.739262e+05 \n", " electricity 2.680870e+06 \n", " construction 3.153222e+05 \n", " trade 4.933951e+04 \n", " transport 2.615115e+05 \n", " other 4.074021e+07 " ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "io.et1_diag.D_cba" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }