python---太阳黑子

最新推荐文章于 2021-05-10 12:42:59 发布

转载最新推荐文章于 2021-05-10 12:42:59 发布 · 313 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：https://my.oschina.net/zhangdapeng89/blog/55188

文章标签：

#python

使用Python和ReportLab库从NASA网站获取太阳黑子活动数据，并绘制预测图表。该脚本下载了太阳黑子活动的历史记录及未来一年的预测数据，通过折线图展示每月的平均黑子数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

2019独角兽企业重金招聘Python工程师标准>>>

太阳黑子（python2.7.2）

#!/usr/bin/env python
from urllib import urlopen
from reportlab.graphics.shapes import *
from reportlab.graphics.charts.lineplots import LinePlot
from reportlab.graphics.charts.textlabels import Label
from reportlab.graphics import renderPDF

URL = 'http://www.swpc.noaa.gov/ftpdir/weekly/Predict.txt'
COMMENT_CHARS = '#:'

drawing = Drawing(400,200)
data = []
for line in urlopen(URL).readline():
    if not line.isspace() and not line[0] in COMMENT_CHARS:
       data.append([float(n) for n in line.split()])
pred = [row[2] for row in data]
high = [row[3] for row in data]
low = [row[4] for row in data]
times = [row[0] + row[1]/12.0 for row in data]

lp = LinePlot()
lp.x = 50
lp.y = 50
lp.height = 125
lp.width = 300
lp.data = [zip(times.pred),zip(times,high),zip(times,low)]
lp.lines[0].strokeColor = colors.blue
lp.lines[1].strokeColor = colors.red
lp.lines[2].strokeColor = colors.green

drawing.add(lp)

drawing.add(String(250,150,'Sunspots',fontSize=14,fillColor=colors.red))

renderPDF.drawToFile(drawing,'report1,pdf','Sunspots')

转载于:https://my.oschina.net/zhangdapeng89/blog/55188