“My IDL program”的版本间差异

来自Shiyin's note
跳到导航 跳到搜索
 
(未显示同一用户的10个中间版本)
第1行: 第1行:
==常用==
==天文==
*spherematch_kth
===binxy===
:两个样本之间找到天球上第k个靠近的 [[file:Spherematch_kth.pro]]
*将x,y分bin(x上),然后显示每个bin里面的平均(或者median)的x和y的数值


==函数==
===Sch_LF===
===Sch_LF===
*Schechter function,微分和数值积分
*Schechter function,微分和数值积分


===Gauss_fun===
*高斯分布函数
*有内禀的Gauss_pdf,Gauss_cvf,Gaussint等程序



==SDSS==
*Conv_specobjID: 将长的SpecObjID转换为plate,mjd,fiber[[file:Con_specObjID.pro.txt]]

==统计相关==
===小程序===
*conf_lev: 用于计算Chi^2分布中的置信区间。

===binxy===
*将x,y分bin(x上),然后显示每个bin里面的平均(或者median)的x和y的数值


===subsample===
*找到两个样本之间共同的子样本。[[文件:Subsample.pro.txt]]
a=randomn(seed,1000)
b=randomn(seed,1000)-0.5
subsample,a,b,suba,subb,0.02,count=NS,KS=ks
sdssidl_plothist,a,bin=0.2
sdssidl_plothist,b,bin=0.2,linestyle=1,/overplot
sdssidl_plothist,a[suba],bin=0.2,linestyle=1,color=cgcolor('red'),/overplot
sdssidl_plothist,b[subb],bin=0.2,linestyle=1,color=cgcolor('blue'),/overplot


==文件读写==
==文件读写==
第11行: 第40行:


===read_col ===
===read_col ===
我自己写的一个read_col 程序,可以只读一个文件的指定列。
*读一个ascii文件的指定列。
例:读入filename文件的第2列,格式是双精度
例:读入filename文件的第2列,格式是双精度
arr=read_col('filename',2,type='dbl')
arr=read_col('filename',2,type='dbl')
第17行: 第46行:
*列的计数从1开始
*列的计数从1开始
*关键词type的默认是字符格式
*关键词type的默认是字符格式

===my_read_csv===
对read_csv程序的修改,用来读取[[CSV文件]],可以实现header命名Tag的功能


==画图==
==画图==

2015年9月25日 (五) 06:56的最新版本

天文

  • spherematch_kth
两个样本之间找到天球上第k个靠近的 文件:Spherematch kth.pro

函数

Sch_LF

  • Schechter function,微分和数值积分


Gauss_fun

  • 高斯分布函数
  • 有内禀的Gauss_pdf,Gauss_cvf,Gaussint等程序


SDSS

统计相关

小程序

  • conf_lev: 用于计算Chi^2分布中的置信区间。

binxy

  • 将x,y分bin(x上),然后显示每个bin里面的平均(或者median)的x和y的数值


subsample

a=randomn(seed,1000)
b=randomn(seed,1000)-0.5
subsample,a,b,suba,subb,0.02,count=NS,KS=ks
sdssidl_plothist,a,bin=0.2
sdssidl_plothist,b,bin=0.2,linestyle=1,/overplot
sdssidl_plothist,a[suba],bin=0.2,linestyle=1,color=cgcolor('red'),/overplot
sdssidl_plothist,b[subb],bin=0.2,linestyle=1,color=cgcolor('blue'),/overplot

文件读写

read_txt

我自己写了一个read_txt文件:Read txt.pro.txt程序,可以智能的将一个N行M列的数据文件(空格分隔)读为一个data[N,M]的flota array。

read_col

  • 读一个ascii文件的指定列。
例:读入filename文件的第2列,格式是双精度
arr=read_col('filename',2,type='dbl')
  • 文件要求空格做分隔符
  • 列的计数从1开始
  • 关键词type的默认是字符格式

my_read_csv

对read_csv程序的修改,用来读取CSV文件,可以实现header命名Tag的功能

画图

aitoff2和aitoff_grid2

  • astro_lib中的atioff,atioff_grid可以实现天球投影,中心在0。通过修改使得X坐标从0,360投影。