自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(66)
  • 资源 (1)
  • 收藏
  • 关注

原创 object_extract_rate.py

【代码】object_extract_rate.py。

2024-05-29 20:38:14 87

原创 CCD光度校准手册

The CCD Photometric Calibration Cookbook

2024-05-16 11:45:35 319

原创 深度学习相关网站

The latest in Machine Learning | Papers With CodeThe latest in Machine Learning | Papers With Code

2024-01-29 10:57:30 338

原创 t4vvvvv

【代码】t4vvvvv。

2023-06-19 20:43:39 983 1

转载 radon https://scikit-image.org/docs/stable/auto_examples/transform/plot_radon_transform.html#:~:text

https://scikit-image.org/docs/stable/auto_examples/transform/plot_radon_transform.html#:~:text=import%20numpy%20as%20np%20import%20matplotlib.pyplot%20as%20plt,dx%2C%20-dy%2C%20sinogram.shape%20%2B%20dy%29%2C%20aspect%3D%27auto%27%29%20fig.tight_layout%28%

2023-06-19 17:30:23 119

原创 t4vvvvv

【代码】t4vvvvv。

2023-06-15 17:33:57 147

原创 opencv program-t3

【代码】opencv program-t3。

2023-06-14 22:31:20 59

原创 期刊分区查询网站

最新SCI期刊影响因子查询及投稿分析系统(2022-2023年) - LetPubJournal Citation Reports欢迎来到中国科学院文献情报中心期刊分区表https://www.webofscience.com/wos/alldb/basic-search

2023-05-30 11:33:34 143

原创 simulator.py

【代码】simulator.py。

2023-04-03 09:58:39 61

原创 python读取txt文件以空行作为数据的切分处理

python读取txt文件以空行作为数据的切分处理

2022-08-25 15:41:47 1850

原创 清华大学开源软件镜像站

pypi | 镜像站使用帮助 | 清华大学开源软件镜像站 | Tsinghua Open Source Mirror

2022-06-20 17:26:47 144

原创 plot_wcs_center

import numpy as npimport matplotlib.pyplot as pltfrom astropy.io import fitsfrom function import listdirfilepath = './J125035+002000_tmp/'wcslist = listdirfile(path,'.wcs')ls_ra = []ls_dec = []wcslist.sort()for i in wcslist: print(i) wcsf.

2022-05-17 16:34:14 172

原创 azalt_transformto_radec

from astropy.io import fitsimport astropy.units as ufrom astropy.coordinates import SkyCoord,EarthLocationfrom astropy.time import Timehdulistex = fits.open('/media/zyy/Disk_D/望远镜检测/ali20220506/站点-20220505-694低轨-600ms-230407-可见光/站点-20220505-694低轨-600.

2022-05-14 17:34:27 199

原创 check special objects

import osimport respobj = ['4382','40099','40100','41744','41745','37779','43530','48274','49317','49323']path = '/data1/imageprocess/'filetype = '.log'spobj_log = '4691_sp_obj.log'fits4382 = []txt4382 = []fits48274 = []txt48274 = []fits37779 =.

2022-04-26 14:52:20 106

原创 【Python】字符串不足位数补零/补空格

目录一、补零 二、补空格回到顶部一、补零if __name__ == '__main__': str1 = '130000' for a in range(100): print(str1 + str(a).zfill(5))回到顶部二、补空格2.1、rjust:右对齐if __name__ == '__main__': str_1 = ['cscs', 'dad', 'a', 'adasdadadad', 'jdsad'.

2022-04-20 19:47:51 6104

转载 sed 文本修改、替换命令用法

概述sed命令是一个面向字符流的非交互式编辑器,也就是说sed不允许用户与它进行交互操作。sed是按行来处理文本内容的。在shell中,使用sed来批量修改文本内容是非常方便的。sed命令的选项语法:sed [选项] [动作]选项与参数:function:一般function的前面会有一个地址的限制,例如 [地址]function,表示我们的动作要操作的行。删除行//test.txt 内容如下sed '1,2d' test.xx输出:其中1,2d中的d.

2022-04-15 10:42:24 10234

原创 字符串中提取浮点数

# -*- coding: cp936 -*-import re string="A1.45,b5,6.45,8.82"print re.findall(r"\d+\.?\d*",string) # ['1.45', '5', '6.45', '8.82']

2022-04-11 12:07:41 808

原创 read 4681_04_07.log

import ospath = '/data1/gabor_test/'os.listdir(path)# path = '/data1/imageprocess/'filetype = '.log'spobj_log = '4691_sp_obj.log'ls_log = []for file in os.listdir(path): if file[-4:] == filetype and file[0:12]=='4691_2022_03': ls_log.a.

2022-04-09 19:06:02 313

原创 out process log

import osfrom astropy.io import fitsfpath='/data2/data/'date = '2022_03_01/'fitspath = fpath+datecmd = 'ls '+fitspath+'*fits >fit.dat'os.system(cmd)tmppath = './tmp/'solved = open('solved.log','a')failed = open('failed.log','a')nostar = open.

2022-04-08 16:03:05 290

原创 check_result new

import pandas as pdimport os'''modify time:2022-01-19 15:10:28editor:zhaoyuanyuan'''# start_time = time.time()# yesterday = datetime.today()+timedelta(-1)# yesterday_format = yesterday.strftime('%Y_%m_%d')yesterday_format = '2022_03_01'fitspath.

2022-04-08 12:02:43 100

原创 write reg

from astropy.io import fitsimport os# cmd1 = 'ls ./tmp/*589*.2acat >2acat.dat'# cmd2 = 'ls ./tmp/*589*.xy >xy.dat'## os.system(cmd1)# os.system(cmd2)blist = [line.rstrip()for line in open('in_gaborstar.log')]path = './tmp/'for j in blist.

2022-03-31 16:00:16 564

原创 mkdir special_obj

currentHour=$(date + %H)echo $currentHourif [[ $currentHour -le 21 ]];then currentTime=$(date +%s); echo ${currentTime}; tmp=$(($currentTime-3600*21)); today=$(date -d"@$tmp" +%Y_%m_%d); todayrd=$(date +%Y%m%d); today1=$(date +.

2022-03-24 12:22:37 142

原创 gnuplot 安装及官网

gnuplot homepage

2022-03-07 11:52:50 694

原创 启动iraf

cd ~/.iraf xgtermcl

2022-03-07 11:25:44 274

原创 check_result.py

import osimport datetimefrom datetime import timedelta, datetime'''modify time:2022-01-19 15:10:28editor:zhaoyuanyuan'''cmd0 = 'rm process.log'os.system(cmd0)curTime = datetime.now()procDay = str("%04d_%02d_%02d" % (curTime.year, curTime.mont.

2022-02-16 17:51:03 206

原创 shell文件--批量压缩文件(.sh)

for line in $(cat datepath.dat)do echo ${line} p="/data2/data/"${line}"/" cd ${p} pwd find *fits -print | xargs -I {} zip {}.zip {}done

2022-02-16 16:01:07 682

原创 shell程序--批量删除fits文件(.sh)

for line in $(cat datepath.dat)do echo ${line} p="/data2/data/"${line}"/" cp rmfits_zyy.py ${p} cd ${p} pwd python3 rmfits_zyy.pydone

2022-02-16 15:58:59 85

原创 删除fits文件

import oscmd1 = 'ls *fits >fit.dat'os.system(cmd1)alist = [line.rstrip() for line in open('fit.dat')]for i in alist: print(i) f = i+'.zip' print(f) if os.path.exists(f): cmd2 = 'sudo rm '+i os.system(cmd2)...

2022-02-16 10:34:58 224

转载 linux教程

Shell 传递参数 | 菜鸟教程Shell 传递参数 我们可以在执行 Shell 脚本时,向脚本传递参数,脚本内获取参数的格式为:$n。n 代表一个数字,1 为执行脚本的第一个参数,2 为执行脚本的第二个参数,以此类推…… 实例 以下实例我们向脚本传递三个参数,并分别输出,其中 $0 为执行的文件名(包含文件路径): 实例 [mycode4 type='bash'] #!/bin/bash # author:菜鸟教程 # url:www.runo..https://www.runoob.com/linux

2022-02-16 10:02:31 52

原创 批量压缩文件

find *fits -print | xargs -I {} zip {}.zip {}

2022-02-15 21:06:56 113

原创 自动下载gaia dr2星表

import astropy.units as ufrom astropy.coordinates import SkyCoordfrom astroquery.gaia import Gaia'''input:center_ra,center_dec,radius,outfile,out_columns'''center_ra = 280center_dec = -60radius = u.Quantity(0.05, u.deg)outfile = 'cc.csv'out_colu.

2022-02-15 16:22:02 627

原创 numpy 使用手册

NumPy documentation — NumPy v1.22 Manual

2022-02-15 10:25:34 454

转载 查看numpy版本

python下安装numpy后首先进入python shell(如果已经设置python路径,则进入终端后用下面命令) python #进入python shell import numpy numpy.__version__#注意前后都有两横

2022-02-15 10:07:27 8750 1

原创 去掉重复恒星--处理星表

import pandas as pdfrom astropy.table import Tableimport numpy as npfrom multiprocessing import Process,Manager,Pooldef cleargaia (newgaia1,newgaia2): # 读取星表数据(ra,dec,mag) gaiadata = pd.read_csv(newgaia1) print('pre combine', newgaia1, le.

2022-02-08 18:35:07 639

原创 筛选并合并星表

#coding=utf-8import pandas as pdfrom astropy.table import Table, vstack,hstackdef filter_gaia(oldgaia): data = Table.read(oldgaia) ra = data['ra'] dec = data['dec'] # Gmag = data['phot_g_mean_mag'] Gmag = data['Gmag'] id = dat.

2022-02-08 18:33:30 629

原创 RAW to FITS

import numpy as npfrom astropy.io import fitsalist = [line.rstrip()for line in open('raw.dat')]for file in alist: raw = file img=np.fromfile(raw,dtype=np.uint16,offset=1024) b=np.reshape(img,(4096,4096)) fit = file[0:-4]+'.fits' fit...

2021-12-23 17:18:17 688

原创 Python列表中去重的多种方法

怎么快速的对列表进行去重呢,去重之后原来的顺序会不会改变呢?去重之后顺序会改变set去重列表去重改变原列表的顺序了l1 = [1,4,4,2,3,4,5,6,1]l2 = list(set(l1))print(l2) # [1, 2, 3, 4, 5, 6]但是,可以通过列表中索引(index)的方法保证去重后的顺序不变。l1 = [1,4,4,2,3,4,5,6,1]l2 = list(set(l1))l2.sort(key=l1.index)print(..

2021-11-08 18:33:01 470

转载 ubuntu 批量修改文件名

方法一:Rename函数Linux下修改文件名有mv和rename。其中mv命令只能对单个文件重命名,这是mv命令和rename命令的在重命名方面的根本区别。linux的rename命令有两个版本,一个是c语言版本的,一个是perl语言版本的,判断方法:输入man rename 看到第一行是:RENAME(1) Linux Programmer’s Manual RENAME(1)这个就是C语言版本的而如果出现:RENAME(1) Perl Programmers Reference Gu

2021-09-27 15:05:58 4651

原创 python 字符串对齐

场景:字符串对齐python提供非常容易的方法,使得字符串对齐 >>> print("abc".center (30,'-')) -------------abc-------------- >>> print("abc".ljust (30)+'|') abc | >>> print("abc".rjust (30...

2021-09-16 20:22:42 88

原创 按照文件名排序

csv_list = glob.glob('*region.csv')csv_list.sort(key=lambda x:int(x[-12:-10]))print(csv_list)

2021-09-15 21:15:19 84

暗标准星场+50.xlsx

暗标准星场+50.xlsx

2021-09-26

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除