Py学习  »  Python

一图胜千言!深度掌握 Python 绘图

数据STUDIO • 5 月前 • 160 次点击  


让数据集改变你的思维

数据可视化是数据科学家传达洞见和讲述数据故事的关键工具。作为 Python 开发者,我们拥有丰富的可视化库和工具,能够创建各种引人入胜的图表。本文将探索一些鲜为人知但实用的可视化类型,如桑基图(Sankey Diagrams)、脊线图(Ridge Plots)、内嵌图(Insets)、雷达图(Radar Chart)和词云图(Word Cloud Plots)。我们将主要使用流行的 Matplotlib、Seaborn 和 Plotly 这些 Python 库来实现这些有趣的可视化效果,让你的数据讲述更加生动有趣。

我们将使用形状、大小、颜色、方向、面积和标记符号等属性,为十种不同的使用案例创建绘图。在每个使用案例中,我们的目标都是创建有效、高效和美观的可视化效果。描述这些词语在图表中的含义如下:

(a) 有效:图表中包含了所有需要传达的信息
(b) 高效:图表中没有多余的数据
(c) 美观:图表以清晰的方式展示数据,吸引注意力

所有图表都是二维图,因为从效率和效果的角度来看,二维图比三维图更清晰易懂,同时更容易描绘距离。文中还将介绍每个用例的代码,并讨论代码和图表中的要点。

用例 1

描述大学之间学生交流流动的桑基图。

这些图展示了资源流动的情况,下面的代码展示了用例的实现。 字符 "A" 代表第一所大学,字符 "B" 代表第二所大学。 数字 3、4、5 分别代表不同的系,即{Statistics, Math, Physics}。第 25 行创建了一个图表,其中 nodelink是字典。 node 使用的 label 对象由唯一的 Depts 院系组成,而 link 使用的两个列表分别由 sending"院系的索引和 acepting 院系的索引组成。

import pandas as pd
import plotly.graph_objects as gr
data = {
    'Sending_Dept':   ['5A''4A''5B''5A''4B''4A''3A''3B''3A''3B''3A''3B'],
    'Accepting_Dept': ['4B''5B''5A''5B''4A''4B''5B''5A''4B''4A''3B''3A'],
    'FlowValue': [134344113253]
}
df = pd.DataFrame(data)
unique_departments = set(df['Sending_Dept']).union(set(df['Accepting_Dept']))
Depts = list(unique_departments)
Dept_indices = {}
for i, dept in enumerate(Depts):
    Dept_indices[dept] = i
sending_indices = []
for dept in df['Sending_Dept']:
    dept_index = Dept_indices[dept]
    sending_indices.append(dept_index)
    print(f"Sending indices are: {sending_indices}")
accepting_indices = []
for dept in df['Accepting_Dept']:
    dept_index = Dept_indices[dept]
    accepting_indices.append(dept_index)
flowvalues = df['FlowValue'].tolist()
# Sankey diagram
fig = gr.Figure(data=[gr.Sankey(
    node=dict( pad=10,thickness=25,line=dict(color="red", width=0.8),label=Depts,),
    link=dict(source=sending_indices,target=accepting_indices,value=flowvalues
    ))])
fig.update_layout(title_text="Sankey Diagram of exchange students flow between University Depts", font_size=12)
fig.show()

生成的"桑基图"图(1)中,节点3A旁的橙色矩形显示了光标放置在节点上时的情况。当光标位于节点"3A"上时,我们可以看到A大学3系接受和派遣交换生的频率。它接受学生1次,派遣学生3次。我们还可以从上面代码片段中的 data 字典推断出这一点,因为"3A"在Sending_Dept列表中出现了3次,在Accepting_Dept列表中出现了1次。节点 "3A" 左边的数字9是它向B大学派出的交换生总数。我们还可以通过在Sending_Dept列表中添加与3A相对应的FlowValues来推断。

我们还注意到,当我们点击节点 "3A" 时,从它发出的箭头会变暗,并显示出与 "3A" 交换学生的其他节点。箭头的粗细与 FlowValues 相对应。总之,桑基图利用箭头的方向和粗细来传递流动信息,并以文字为基础为每个节点形成累积流动。

图 1. 桑基图显示了两所大学各系之间的学生交流流

用例 2

绘制一家房地产中介公司的房屋销售数据。

一位数据科学家在房地产中介公司工作,机构要求绘制上个月售出房屋信息的二维图。每栋售出的房屋需包含房价、距离市中心、方向、代理佣金和销售代理的公司级别(助理、副总裁、合伙人)的信息。二维图形信息量大,可使用复杂对象描述地块上的每栋房屋。具体来说,使用“笑脸表情符号”实现方法的代码片段如下。

import matplotlib.pyplot as plt
import numpy as np
np.random.seed(125)
num_houses = 10
distances = np.random.uniform(030, num_houses)  # distance from city center
prices = np.random.uniform(4002000, num_houses) * 1000  # sale price in thousands
directions = np.random.choice(['N''S''E''W'], num_houses)  # direction from city center
agent_levels = np.random.choice([123], num_houses)  # agent's level
def get_emoji_size(level):
    size_map = {125023803700}
    return size_map.get(level, 120)  # Increased size for better visibility
def get_emoji_color_new(price):
    if price 600000:
        return 'white'  # Light yellow for $400k-$600k
    elif price 800000:
        return 'yellow'        # White for $600k-$800k
    elif price 1000000:
        return 'pink'         # Pink for $800k-$1 million
    else:
        return 'lime'         # Lime for $1 million-$2 million
def rotate_smiley(direction):
    rotation_map = {'N'0'E'270'S'180'W'90}
    return rotation_map.get(direction, 0 )  # default no rotation if direction not found
plt.figure(figsize=(128))
for i in range(num_houses):
    plt.scatter(distances[i], prices[i], s=get_emoji_size(agent_levels[i]),\
    c=get_emoji_color_new(prices[i]), 
                marker='o', edgecolors='black', alpha=0.8)
    plt.text(distances[i], prices[i], "😊", fontsize=agent_levels[i]*10
             rotation=rotate_smiley(directions[i]), ha='center', va='center',\
             fontweight='bold')
plt.xlabel('Distance from City Center (km)')
plt.ylabel('Sale Price ($)')
plt.title('House Sales Data for 10 Houses: Price vs Distance with New Color Scheme')
plt.grid(True)
plt.show()

如上面的代码和下面图 2 中的散点图所示,X 轴和 Y 轴分别对应于与市中心的距离和销售价格。以下是一些要点:

  • 表情符号的大小用于描述销售代理的级别。尺寸越大,代理人的级别越高
  • 表情符号微笑的位置(顶部、底部、左侧、右侧)表示从市中心出发的方向。例如,如果微笑在顶部,房子就在市中心的北面。
  • 表情符号的颜色表示中介的佣金。例如,石灰色的佣金为 6%,粉红色的佣金为 5%(房地产中介的佣金政策是,售价越高,佣金越高)。

作为本使用案例的结论,我们使用散点图和表情符号的形状、颜色和大小,在二维图中表达了已售房屋的五个属性。使用复杂的对象(如笑脸)来表示数据点,确实有助于在图中包含大量信息。

图 2. 房屋销售数据散点图

用例 3.在旭日图中显示某大学各学院和系的组成。还必须表达各学院和部门的规模信息。

这是一种不同组成部分具有层次结构的情况。在这种情况下,旭日图是理想的选择。下面的代码片段展示了实现过程。第一个数组 labels 包含旭日图的名称。第二个数组 parents 包含层次结构,而数组 values 包含段的大小。

import plotly.graph_objects as go
labels = ["PineForest University"
          "College of Engineering""Natural Sciences""Medical School",
          "Chemical Engineering""Mechanical Engineering""Computer Science""Electrical Engineering",\
          "Biomedical Engineering",
          "Astronomy""Biology""Mathematics""Physics","Chemistry","Statistics","EarthSciences",\
          "Emergency Medicine""Neurology""Cardiology""Pediatrics""Pathology","Family Medicine", \
          "Orthopedics""Obstetrics""Anesthesiology""Anatomy""Physiology",\
          "Microbiology","Immunology"]
parents = [""
           "PineForest University""PineForest University""PineForest University",
           "College of Engineering""College of Engineering""College of Engineering",\
           "College of Engineering","College of Engineering",
           "Natural Sciences","Natural Sciences","Natural Sciences","Natural Sciences","Natural Sciences",\
           "Natural Sciences","Natural Sciences",
           "Medical School""Medical School""Medical School""Medical School","Medical School", \
           "Medical School","Medical School""Medical School""Medical School",\
           "Medical School""Medical School""Medical School","Medical School"]
values = [0
          5030200
          85100180163,120,
          90,70 ,104,54,180,100,70,
          70,180,20085,170,89,75,120,76,150,67,56,239]
fig = go.Figure(go.Sunburst(
    labels=labels,
    parents=parents,
    values=values,
))
fig.update_layout(margin=dict(t=0, l=0, r=0, b=0))
fig.show()

图 3 显示了上述算法的输出结果。请注意,每个片段的大小与 "值 "数组中的相应数字成正比。下图显示,当我们点击一个片段时,它的大小就会显示出来(Immunology 239)。

旭日图使用大小、颜色和文本来描述大学不同实体的层次结构。

图 3. 旭日图显示了一所大学不同学院和系的结构

用例 4

我们的房地产客户需要一张二维图,显示上个月售出房屋的信息: (a) 售价,(b) 面积,(c) 海边距离,(d) 火车站距离。必须放大已售出最多房屋的图段。

本用例与用例 2 类似,但数据科学家有新任务:创建一个放大版的最繁忙图段。这就是插图。插图可以放大图表的重要部分,提高效果。下面的代码片段显示了插图的实现。

import seaborn as sns
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
np.random.seed(0)
size_of_house = np.random.uniform(10004000100)  # Size in square feet
price = size_of_house * np.random.uniform(150350) + np.random.normal(050000100)  # Price
distance_from_train = np.random.uniform(0.55100)  # Distance from train in miles
distance_from_ocean = np.random.uniform(0.110100)  # Distance from ocean in miles
df = pd.DataFrame({
    'Size of House': size_of_house,
    'Price': price,
    'Distance from Train': distance_from_train,
    'Distance from Ocean': distance_from_ocean
})
# Sample 10 points for a less cluttered plot
sampled_df = df.sample(10)
# Adding the inset 
# Filtering sampled data for the inset plot
#inset_data = sampled_df[(sampled_df['Size of House'] >= 2000) & (sampled_df['Size of House'] <= 3000) & (sampled_df['Price'] >= 250000) & (sampled_df['Price'] <= 600000)]
inset_data=sampled_df
inset_ax = ax.inset_axes([0.70.050.250.25])  # Inset axes 
# Scatter plot in the inset with filtered sampled data
sns.scatterplot(data=inset_data, x='Size of House', y='Price', ax=inset_ax, size='Distance from Train', sizes=(40200), hue='Distance from Ocean', alpha=0.7, marker='D', legend=False)
# Capturing the part of the original regression line that falls within the bounds
sns.regplot(
    x='Size of House'
    y='Price'
    data=df,  # Using the entire dataset for the trend line
    scatter=False
    color='red'
    ax=inset_ax, 
    lowess=True,
    truncate=True  # Truncating the line within the limits of the inset axes
)
# Adjusting the limits of the inset axes
inset_ax.set_xlim(20003500)
inset_ax.set_ylim(6000001000000)
# Display the plot
plt.show()

插图 4 展示了散点图,其中 X 轴和 Y 轴分别代表房屋面积和售价。

菱形点代表已售出的房屋,颜色表示与海边的距离,大小对应着与火车站的距离。 这个用例展示了如何通过嵌入地块的重要部分来提高地块的有效性。

图 4. 带有插图的房屋销售数据

用例 5

为公众号数据STUDO制作一张吸引眼球的海报

在本例中,我们将使用词云图。这种图表包含不同大小的词语,词语的大小取决于它在原文中出现的频率。图 5 显示了词云图,下面的代码片段显示了生成该图的 Python 代码。

from wordcloud import WordCloud
import matplotlib.pyplot as plt

# Sample from the reviews of a book
reviews = """
🏴‍☠️宝藏级🏴‍☠️ 原创公众号『数据STUDIO』内容超级硬核。公众号以Python为核心语言,垂直于数据科学领域,包括可戳👉 Python|MySQL|数据分析|数据可视化|机器学习与数据挖掘|爬虫 等,从入门到进阶!
"""

font_path="/System/Library/fonts/PingFang.ttc" #解决中文乱码

wordcloud = WordCloud(width=800, height=400, background_color='white', font_path=font_path).generate(reviews)
plt.figure(figsize=(127), dpi=150)
plt.imshow(wordcloud, interpolation='bilinear')
plt.axis('off')  
plt.show()
plt.figure(figsize=(127))
plt.imshow(wordcloud, interpolation='bilinear')
plt.axis('off')  
plt.show()
图 5. 根据公众号简介创建的词云

用例 6

以清晰而有说服力的方式介绍一种新零食产品的优势。

我们的数据科学家的任务是展示市场营销部门进行的一次民意调查的结果。在这次民意调查中,客户被要求对一种高蛋白棒的新旧配方的性价比、营养价值、外观和口味进行评分。

据科学家使用雷达图来显示投票结果。雷达图将不同的属性(口味、外观等)放在不同的坐标轴上,然后连接属于同一实体(本例中为蛋白棒)的属性值,形成一个多边形区域。不同的区域使用不同的颜色,便于查看者掌握产品之间的差异。

雷达图(Radar Chart),又可称为戴布拉图、蜘蛛网图(Spider Chart),每个分类都拥有自己的数值坐标轴,这些坐标轴由中心向外辐射, 并用折线将同一系列的值连接,用以显示独立的数据系列之间,以及某个特定的系列与其他系列的整体之间的关系。

图 6 显示了生成的图形,下面的代码片段显示了生成图形的 Python 代码。下面代码的第 7-9 行有一个有趣的地方。enpoint被设置为False,第一个值被添加到 old_valuesnew_values 数组的末尾,以关闭循环。

import matplotlib.pyplot as plt
import numpy as np
categories = ['Taste''Nutritional value''Appearance''Money Value']
old_values = [3443]
new_values = [5545]
N = len(categories)
angles = np.linspace(02 * np.pi, N, endpoint=False).tolist()
old_values = np.append(old_values, old_values[0])
new_values = np.append(new_values, new_values[0])
angles = np.append(angles, angles[0])
fig, ax = plt.subplots(figsize=(66), subplot_kw=dict(polar=True))
ax.fill(angles, old_values, color='blue', alpha=0.4)
ax.fill(angles, product_values, color='lime', alpha=0.4)
ax.set_yticklabels([])
ax.set_xticks(angles[:-1])
ax.set_xticklabels(categories)
plt.title('Benchmarking: Standard vs New Product', size=20)
plt.legend(['Old Product''New Product'], loc='lower right', bbox_to_anchor=(1.11.1))
plt.show()

如下图所示,人们认为新的蛋白质棒在三个关键方面优于旧产品:营养价值、口味和性价比。

在这个案例中,展示了如何利用颜色和形状在雷达图上区分两种产品。

图 6. 展示两种产品属性差异的雷达图

用例 7

在同一幅图上,显示不同类型图书在这些类型的图书节前后的销售情况。

在这个案例中,我们需用同一图绘制不同类别的时间数据。Python的ridge plot适合这种情况。桥形图以不同类别分布的垂直堆叠图呈现,便于比较它们的异同。图7显示了不同图书类型的销售差异以相应图书节为中心。将分布图放在同一图上,使用不同颜色,可创建一个视觉效果良好的图表,有效比较。

相应代码如下:

(a) ridge plot使用joyplot Python软件包创建
(b) 第6-13行模拟了各种图书流派的分布情况,书籍类型后面的数字是相应分布的中心和标准偏差。
(c) 在第20行,我们将数据放入一个DataFrame中,因为joyplot希望输入一个DataFrame。
(d) 在第30行、第33-36行,我们指定了轴[-1]上x轴的参数,即底图的参数。

import matplotlib.pyplot as plt
import pandas as pd
from joypy import joyplot
import numpy as np
import scipy.stats as stats
book_festivals = {
    'Mystery': [1,1], 
    'Adventure fiction':[6,2],
    'Fantasy': [11,3],
    'Self-help':[10,2],
    'Sci-Fi': [7,2],
    'Historical fiction': [12,1]
}
genre_data = {}
for genre, params in book_festivals.items():
    mu,sigma=params
    data = np.random.normal(mu, sigma, 300)
    genre_data[genre] = data
genre_data_dict = {genre: data.tolist() for genre, data in genre_data.items()}
data=pd.DataFrame(genre_data_dict)
# Now we can plot it using joyplot
fig, axes = joyplot(
    data,
    figsize=(3018),
    colormap=plt.cm.plasma,
    alpha=1.0
)
colors = plt.cm.plasma(np.linspace(01, len(data.columns)))
# Set x-axis labels and limits
axes[-1].set_xlim(112)
month_order = ['January''February''March''April''May''June',
               'July''August''September''October''November''December']
axes[-1].set_xticks(range(113))
axes[-1].set_xticklabels(month_order)
axes[-1].tick_params(axis='x', labelsize='30')
axes[-1].set_xticklabels(month_order)
for ax in axes:
    ax.set_yticks([])
    
legend_handles = [plt.Line2D([0], [0], color=color, lw=6, label=genre) for genre, \
                  color in zip(book_festivals.keys(), colors)]
ax.legend(handles=legend_handles, title='Book Genres', loc='upper right', fontsize='26')
plt.show()
图 7. 图书类型的脊状分布图

用例 8

给定信用卡公司提供的有关年龄和购买金额的数据,在图中突出显示购买金额最高的组别。

二维直方图似乎适合这类数据,与标准直方图不同,标准直方图按一维分类,而二维直方图同时按两个维度分类。以下是代码片段和图8显示的生成的二维直方图。

import numpy as np
import matplotlib.pyplot as plt
np.random.seed(0)
ages = np.random.randint(186120)  # Random ages between 18 and 60 for 20 customers
purchase_amounts = np.random.randint(200100120)  # Random purchase amounts between 200 and 1000 for 20 customers
plt.figure(figsize=(106))
plt.hist2d(ages, purchase_amounts, bins=[55], range=[[1860], [2001000]], cmap='plasma')
plt.colorbar(label='Number of Customers')
plt.xlabel('Age')
plt.ylabel('Purchase Amount ($)')
plt.title('2D Histogram of Customer Age and Purchase Amount')

下面的二维柱状图显示,消费金额最高的是 35-45 岁的人群。总之,二维柱状图有效且直观地显示了购买金额最高的群体。

图 8:显示年龄与购买金额的二维柱状图

用例 9

在一幅图中,显示 1990-2010 年生物、天文、物理和数学系的学生人数。

在这种情况下,我们可以使用堆叠面积图。这种类型的图表用于显示不同类别的时间数据,可以使用堆叠区域图来展示不同类别的时间数据,用户可以看到:

(a) 每个类别随时间的变化;
(b) 它们的相对规模;
(c) 它们随时间的总规模。每个类别用不同的颜色标识。

代码如下,图 9 显示了不同院系学生人数的堆叠面积图。

堆叠面积图与面积图类似,都是在折线图的基础上,将折线与自变量坐标轴之间区域填充起来的统计图表,主要用于表示数值随时间的变化趋势。 而堆叠面积图的特点在于,有多个数据系列,它们一层层的堆叠起来,每个数据系列的起始点是上一个数据系列的结束点。

import matplotlib.pyplot as plt
import numpy as np
years = np.arange(19902011)
biology_students = np.linspace(1001000, len(years))
astronomy_students = np.linspace(50150, len(years))
math_students = np.linspace(80500, len(years))
physics_students = np.linspace(70480, len(years))

plt.figure(figsize=(106))
plt.stackplot(years, biology_students, astronomy_students, math_students,\
              physics_students, labels=['Biology''Astronomy''Math''Physics'],\
              colors=['#76b947''#fcb001''#00a2e8''#ff7f27'])
plt.legend(loc='upper left')
plt.title('University Department Student Numbers (1990-2010)')
plt.xlabel('Year')
plt.ylabel('Number of Students')
plt.tight_layout()

plt.show()
图 9. 不同院系学生人数的堆叠面积图

用例 10

显示一家公司 2010 年至 2020 年的营销成本与收入对比。

与上一个用例一样,我们也可以在这个用例中使用堆叠面积图。然而,堆叠面积图并不是万能的,也不是绘制随时间变化的分类差异的一站式解决方案。对于当前用例,比较营销成本与收入的更有效方法是r收入柱状图,并将营销成本作为趋势线。代码如下所示。将营销成本作为趋势线绘制在收入条形图前,可以更有效地传达营销成本增加的信息。

import matplotlib.pyplot as plt
import numpy as np
years = np.arange(20102021)
revenue = np.linspace(1060, len(years))  
marketing_costs = revenue * 0.4  

plt.figure(figsize=(86))
plt.bar(years, revenue, color='lightblue', label='Revenue')
plt.plot(years, marketing_costs, color='red', marker='o', \
linestyle='-', linewidth=2, markersize=6, label= 'Marketing Costs')
plt.xlabel('Year')
plt.ylabel('Revenue (Millions)')
plt.title('Bar and Line Combination Plot')
plt.legend()
plt.tight_layout()
plt.show()

plt.figure(figsize=(86))
plt.stackplot(years, marketing_costs, revenue, labels=['Marketing Costs''Revenue'], colors=['lightcoral''palegreen'])
plt.xlabel('Year')
plt.ylabel('Total Amount (Millions)')
plt.title('Area Plot')
plt.legend()
plt.tight_layout()
plt.show()
图 10. 以营销成本为趋势线的收入柱状图

图 11. 收入和营销成本的堆叠区域图

结论

在这篇文章中,我们讨论了各种用例和相应的 Python 图,强调了 Python 语言在可视化方面的丰富性。从 Sankey 图的流程到蜘蛛图的轮廓再到山脊图的多峰,展示了 Python 将原始数据转化为引人注目的数据故事的能力。讨论的图只是 Python 广泛工具箱中的一个缩影。


🏴‍☠️宝藏级🏴‍☠️ 原创公众号『数据STUDIO』内容超级硬核。公众号以Python为核心语言,垂直于数据科学领域,包括可戳👉 PythonMySQL数据分析数据可视化机器学习与数据挖掘爬虫 等,从入门到进阶!

长按👇关注- 数据STUDIO -设为星标,干货速递

Python社区是高质量的Python/Django开发社区
本文地址:http://www.python88.com/topic/170622
 
160 次点击