- 博客(20)
- 收藏
- 关注
原创 重复源的筛选
在天文学中有一些数据中是有很多重复源的,例如LAMOST,大约有20%的重复观测,实际使用中,我们应当去掉重复源,否则在数据拟合或者进行机器学习的时候,会导致部分源的权重偏大,导致模型拟合的不好。对于重复源,一般是两种做法,如果数据中有信噪比,那么我们肯定是选择重复源中信噪比最大的一个数据,如果没有信噪比,那么我们就应该对重复源求一个平均值或者中位数,
2022-11-01 15:56:38 255 1
原创 天文学论文绘图注意事项
1. 需要让刻度线在上下左右都显示,并且刻度向内,如下图所示可以通过如下步骤实现:首先使用如下命令找到matplotlib配置文件存在的位置import matplotlibprint(matplotlib.matplotlib_fname())接着在命令行下使用vim 打开配置文件,找到ticks部分将xtick.top,ytick.right改为True,xtick.direction,ytick.direction改为in2. 绘图时要注意绘图时要导出axes,如果画.
2021-10-19 20:42:41 962
原创 ADS文章查询技巧
当我们写完论文初稿交给导师的时候,导师一般会让我们加一些参考文献,一般导师不会直接告诉我们名字,只会告诉我们第一作者是谁,哪一年发表的,大体意思是什么,这就需要我们自己去ADS查了,而一般输入一个人名,你会搜到浩如烟海的文章,下面我讲一下我自己在查询了若干篇文章总结的一些技巧首先,输入作者名字,然后选择年份,就在下图这里把起止都输入导师告诉你的年份如果你知道这是一篇很著名的文章,那么在排序里按照引用次数排序如果你不知道是不是一篇著名的文章,那就按照第一作者排序,也就是上图的firs.
2021-09-25 19:05:35 483
原创 用0和1记录多重标准筛选时,每个标准的吻合情况
当我们对大量数据进行多重标准的筛选的时候,往往想记录一下每一个点是否吻合每一个标准如果是一个标准,一般在这时引入一个flag量,用0和1来表示是否吻合,但是当有多重标准时怎么操作呢?还是引入一个flag量,初始值设为0,不通过第一个标准的flag+1,不通过第二个标准的flag+2,不通过第三个标准的flag+4,不通过第四个标准的flag+8,按照以此类推,注意要累加,这样才能记录不吻合多个标准的点具体不吻合那些那些标准,最后把他们转化为二进制数字,1就表示不通过,0就表示通过,Python十进制
2021-08-20 16:58:29 117
原创 用AAS 模板写Latex 教程
用AAS 模板写Latex 教程(未完成)首先你要下载他们的模板,地址是https://journals.aas.org/aastexguide/#intro先添加一些必要的部分,说明你用的模板是什么\documentclass[linenumbers]{aastex631}%linenumbers表示显示行号\newcommand{\vdag}{(v)^\dagger}\newcommand\aastex{AAS\TeX}\newcommand\latex{La\TeX}\rece
2021-05-07 13:49:31 2127 4
原创 Python 高斯拟合
Python 高斯拟合通常我们进行高斯拟合的办法是导入scipy的curve_fit 包,不过这需要自己手写一个高斯分布的函数表达式,不是很方便,astropy提供了一个写好的高斯拟合包调包from astropy.modeling import models, fittingimport numpy as npimport matplotlib.pyplot as plt生成一个高斯的数据为了检验拟合结果的好坏,我们先生成一个μ=0.5,σ=0.2的高斯数据,并赋予他一个噪
2020-12-23 20:35:53 15492 6
原创 python计算金星凌日
首先,需要安装几个包jplephem,astropy,numpy直接pip安装即可需要下载NASA的行星历表,jplephem的文档里面讲了怎么下载jplephem的官方文档,下面的例子中使用的是de421代码#计算金星凌日发生的时间(1900-2050)from jplephem.spk import SPKfrom astropy.time import Timeimport numpy as npstart_time = '1900-01-03T00:00:00.0'end.
2020-06-12 10:50:13 1489
原创 pandas 使用
对于操作csv文件,pandas是一个很好的工具,比如分析部分学生的成绩:#读入数据文件data1=pd.read_csv('students_grade1.csv')>>>> Unnamed: 0 学号 语文 数学 英语0 0 20167598 89 89 891 1 20167589 85 85 852 2 20164569 52 52 52
2020-05-25 10:10:06 198
原创 python列表表达式
对于新手来说,生成一个新的列表,他可能会使用for 循环,但是for循环一般来说是要少用的(因为比较慢),那么如何不使用for循环生成一个新的列表呢?那就需要用到列表表达式了,例子如下:#生成一个从0到9的列表a=[x for x in range(10)]print(a)>>>>[0, 1, 2, 3, 4, 5, 6, 7, 8, 9]如果生成的列表每一个元素是二元的呢?a=[(x,y) for x in range(3) for y in range(3)]
2020-05-25 09:15:07 335
原创 获取交叉列表后的索引
交叉列表求并集是一个很简单的操作,使用set函数即可,例子如下a=[1,2,3,4,5,6,7,8,9]b=[4,5,1,2,9,8]print(set(a)&set(b))>>>>{1, 2, 4, 5, 8, 9}但是有些时候是需要求出这个交集在一个列表里面对应的索引的,如何求出对应的索引呢?比如求列表a和b的交集在列表a中的索引,方法是把列表a的索引和数据设置为一个字典,然后使用set函数交叉这个字典与列表b,最后输出对应的字典就行了,例子如下:in
2020-05-25 09:00:54 383
原创 正确的python画图方式
1.通常,我们直接通过plt.plot(),plt.scatter()等命令来画图,但这个方式画图有很多限制,比如你只想看图中的一部分,就只能改变x,y来实现,如果要反转坐标轴,根本无法实现2. 这里我介绍一下正确的画图方式,首先,放置一个fig和ax,也就是画布和坐标系fig1,ax1=plt.subplots()然后再画图,命令变成ax1.plot(x,y)ax1.scatter(...
2020-04-20 18:29:20 1001 1
原创 matplotlib.pyplot.scatter的颜色问题
1.通常在import包时,将matplotlib.pyplot缩写为plt,也就是import matplotlib.pyplot as plt2.plt可用于画散点图,也就是plt.scatter,这个函数常见的可选参数有plt.scatter(x,y,s,c)其中x,y是x轴,y轴的坐标,s是点的大小,如果s是一个数字,那么就表示所有点的大小都是这个,如果s是一个和x,y一样长的...
2020-04-20 18:11:38 5540
原创 pandas dataframe合并时的问题
使用函数:pd.concat([df1,df2])这样就可以将两个dataframe df1,df2合并可选参数:(1)ignore_index,默认是false ,但是会导致合并后的dataframe的索引错误,ser1 = pd.Series(['A', 'B', 'C'])ser2 = pd.Series(['D', 'E', 'F'])pd.concat([ser1, s...
2020-03-23 16:03:30 927
翻译 使用python进行普通日期和儒略历的转换
使用软件包astropy的time包安装pip install astropy使用方法调包import numpy as npfrom astropy.time import Time写入一个时间(为了演示数组的操作,我输入了两个时间)times = ['1999-01-01T00:00:00.123456789', '2010-01-01T00:00:00']t ...
2020-03-18 15:54:03 2105
原创 使用python进行天文学坐标转换
使用软件包astropy的coordinates包使用方法1.调包from astropy import units as u #用于单位转换的包from astropy.coordinates import SkyCoord2.写入数据,比如你现在有两个list,一个是赤经a,一个是赤纬b,先把他们写进来c=SkyCoord(ra=a*u.degree,dec=b*u....
2020-03-12 18:53:30 6769
原创 fits文件的读取与查看
fits文件读写(写给自己,大神勿喷)fits文件的结构一个fits文件包括多个HDU(两个以上),开头一个是Primary HDU,储存着类似于观测日期等的信息,剩下的是Bin HDU,储存着观测数据每个HDU由header 和data组成,data中储存着数据,header表示相应的信息,储存着和Primary HDU一样的信息,并且描述第几列的数据是什么,类型是什么TTYPE1=...
2020-01-19 10:13:11 14106 10
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人