目的:
经常需要对信息分析中的结果文件进行分析,有时可能期望提取指定染色体范围内,指定区间内的文件。
代码:
#!/bin/bash echo -e "\033[32m\n地球是圆的,而看似像终点的地方可能也只是起点。\n \033[0m" # 如果不指定输入结果名称,使用默认名称 outfile="result."`date +%Y-%m-%d` function Site() { echo -e "\033[31m输出结果名称为:$outfile\033[0m" less $infile | awk -F "\t" ‘NR < 2{print $0}NR>1{if($2 == ‘‘‘$chr‘‘‘ && $3 > ‘‘‘$newstart‘‘‘ && $3 < ‘‘‘$newend‘‘‘)print $0}‘ > $outfile } function x2t() { echo -e "结果文件进行格式转换:${outfile} ===> ${outfile}.xlsx" python /WORK/Disease/lmt/Code/xls2txt/xls2txt.py ${outfile} ${outfile}.xlsx t2x } while getopts ":i:c:s:e:o:hf" opt;do case $opt in i) infile=$OPTARG echo -e "\033[32m输入文件为:$infile\033[0m" ;; c) chr=$OPTARG echo -e "\033[32m查询的染色体为:$chr\033[0m" ;; s) start=$OPTARG newstart=`expr ${start} - 1` echo -e "\033[32m查询的起始位置为:$start\033[0m" ;; e) end=$OPTARG newend=`expr ${end} + 1` echo -e "\033[32m查询的终止位置为:$end\033[0m" ;; o) outfile=$OPTARG ;; f) ;; h) echo -e ‘‘‘ 使用帮助: 1) site -i infile -c chr -s start -e end 2) site -i infile -c chr -s start -e end -o result # 自己指定输出结果名称 3) site -i infile -c chr -s start -e end [-o result] -f # 是否对生成的结果文件进行格式转换 4) site -o file -f # 对file文件进行格式转换 ‘‘‘ esac done if [ $infile"x" != "x" -a $chr"x" != "x" -a $start"x" != "x" -a $end"x" != "x" ]; then echo "开始提取结果......" Site if [[ $f"x" == "x" ]]; then x2t fi elif [ -s $outfile -a $f"x" == "x" ]; then x2t elif [[ "$@x" == "x" ]]; then # echo "小朋友, 你是不是忘记啥啦......" # echo "查看帮助: site -h" site -h fi
中间有个格式转换的代码,如下:
#!/usr/bin/env python #-*- coding:utf-8 -*- import openpyxl from openpyxl import load_workbook #读取excel文件 from openpyxl import Workbook # 创建xlsx文件 class ChangeType(object): def __init__(self, args): self.source = args["source"] self.target = args["target"] def xls2txt(self): wb = load_workbook(self.source) ws = wb.active out = open(self.target, ‘w‘) for row in ws.rows: for cell in row: value = cell.value out.write("{value}\t".format(**locals())) out.write("\n") def txt2xls(self): wb = Workbook() ws = wb.create_sheet("Result", index=0) with open(self.source, ‘r‘) as fr: for line in fr: linelist = line.strip().split("\t") ws.append(linelist) wb.save(self.target) def main(): demo = ChangeType(args) if args["type"]=="x2t": demo.xls2txt() elif args["type"] == "t2x": demo.txt2xls() else: exit("请输入想要转换的格式") if __name__ == "__main__": import argparse parser = argparse.ArgumentParser() parser.add_argument("source", help="原格式文件") parser.add_argument("target", help="想要转换的格式文件") parser.add_argument("type", help="请输入想要转换的格式", choices=["x2t", "t2x"]) args = vars(parser.parse_args()) main()