python绘制散点图

文章目录

1.实验目的

掌握Python绘图库matplotlib库

2.需求

3. 代码以及资源文件

python 复制代码
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd

# 读取NPZ文件
data = np.load('../Files/国民经济核算季度数据.npz', allow_pickle=True)

# 获取时间和第一到第三产业生成总值当季值数据
plt.rcParams['font.sans-serif'] = ['SimHei']  # 设置为黑体字体


# 获取数组中所有第二列数据 即
# print(data['values'][:, 1])

def getAssignedData(arr):

    newArr = []

    for i in range(len(arr)):
        if 4 * i <= len(arr):
            newArr.append(arr[4 * i])

    return newArr


time = data['values'][:, 1]
gdp1 = data['values'][:, 3].astype(float)
gdp2 = data['values'][:, 4].astype(float)
gdp3 = data['values'][:, 5].astype(float)

# 创建散点图
plt.scatter(time, gdp1, label='第一产业', marker='o', c='r')
plt.scatter(time, gdp2, label='第二产业', marker='D', c='b')
plt.scatter(time, gdp3, label='第三产业', marker='o', c='#FFFF01')

plt.ylabel('生产总值(亿元)')
plt.title('2000-2017年各产业季度生产总值散点图')
plt.xticks(rotation=45, fontsize=4)
# 选择性显示横轴坐标标签
labels = getAssignedData(time)
arr = [' ']
for label in labels:
    arr.append(label)

plt.xticks(range(4, 73, 4), labels=getAssignedData(time))
start = plt.xticks()[0]
print(start)
arr1 = [-1]
start = start.tolist()
for ele in start:
    arr1.append(ele)

print(arr1)
plt.xticks(np.array(arr1),arr)
# 添加图例 左上角小图
plt.legend(loc='upper left')

# 保存散点图图片
plt.savefig('散点图.png')

# 创建DataFrame并保存为Excel文件
df = pd.DataFrame({'时间': time, '第一产业': gdp1, '第二产业': gdp2, '第三产业': gdp3})

# 显示散点图
plt.show()

资源文件下载地址):

链接:https://pan.baidu.com/s/1_h-q5sBbMqtRohL7qW0jzw?pwd=r02f

提取码:r02f

--来自百度网盘超级会员V5的分享

4.实验结果

相关推荐
AC赳赳老秦4 分钟前
OpenClaw SEO写手Agent实操:生成结构化文章,适配CSDN搜索规则
大数据·人工智能·python·搜索引擎·去中心化·deepseek·openclaw
廖圣平12 分钟前
从零开始,福袋直播间脚本研究【八】《策略模式》
开发语言·python·bash·策略模式
2301_8042154115 分钟前
Python类型提示(Type Hints)详解
jvm·数据库·python
第一程序员24 分钟前
非专业转码心路历程与Rust学习规划
python·github
一帧一画25 分钟前
基于 CNN+LSTM 的全国年度降水预测实践
python·cnn·lstm
骄阳似火_201825 分钟前
Anaconda的详细安装步骤
python
Thomas.Sir27 分钟前
第八章:Python3 之 异常与文件处理【从基础入门到底层原理+项目实战】
python·ai·文件处理·异常
2301_8166512228 分钟前
用户认证与授权:使用JWT保护你的API
jvm·数据库·python
Sunshine for you29 分钟前
Python单元测试(unittest)实战指南
jvm·数据库·python
wefly201735 分钟前
m3u8live.cn:免安装 HLS 在线播放器,流媒体调试效率神器
开发语言·javascript·python·django·ecmascript·hls.js 原理·m3u8 解析