超简易!Circos在Python中的实现方法
很高兴再次见到生信技能树的粉丝们,我是技能树VIP小编tsznxx,目前在安德森肿瘤医院工作,记忆力好的小朋友应该对我之前的教程有印象:
用GenePred注释文件进行数据分析
在那里我用了最优雅的方式来解决基因结构的呈现方式,希望对大家胃口。
如何下载注释文件并查看基因结构
对bed格式的基因组区间文件进行基因注释
生物信息学技能面试题(第3题)-探索人类基因组注释文件
生物信息学技能面试题(第5题)-根据GTF画基因的多个转录本结构
这里我将带来完全不一样的教程:
人生苦短
我用Python!
之前用过Perl版的circos,配置文件一堆,安装包一堆,非常痛苦。后来也试过R版的circos,也并不方便。于是自己花时间用Python实现了这个东西。代码一共200多行,大概一半是注释行。还有部分是优化呈现。基本功能的实现极其简单。
1
Packages (git@github.com:tsznxx/PyCircos.git)
- Python版本2.7.
- 用到的只有3个基本包:matplotlib,pandas,numpy。这些都是Pythonb必装包。没有任何其他额外的包!!!
2
实现技巧
2.1 matplotlib 的polar 画图。 说白了就是按照极坐标的形式画常规的matplotlib图。
?
pax = fig.add_axes([0,0,1,1],polar=True)
然后就可以像普通2D的axis对象一样call各种bar,scatter,fill_between, vlines, hlines等图了。参数上有细微的差别,角度(弧度单位)是横坐标,半径方向是纵坐标。 染色体band 用pax.bar() (填充和非填充两种方式) 刻度用pax.vlines() CNV的图用的是pax.fill_between() 2.2 画inter-chrom link 这个用的是SVG里面的路径Path里的CURVE3(这个功能很简单,在Javascript的SVG画图里面很常用,但是由于我很少用JavaScript,我花了好久才找到这种实现方式)。CURVE3的意思就是由3个点确定一条曲线(quadratic Bézier curve),(起点,控制点,终点)。控制点是表示切线方向。
对于一条路径,就是一组点组成的闭合曲线:(起点,控制点1,点2,控制点2,...... ,点N,控制点N,起点)。我们所见的inter-chrom 链接就可以用9个点组成一条闭合曲线(chord):
3
画图
当呈现方式就绪后,下一步就是做图了。
cytoband是用的UCSC的文件。
3.1 画染色体cytoband (barplot,颜色采用的是R里面的染色体配色)
3.2 画染色体空框(barplot, {'edgecolor':'k','linewidth':1,'linestyle':'-','fill':False})
3.3 画 ticks (每50M为间隔,刻度用M表示)
3.4 画染色体名字 (旋转)
3.5 画CNV amplification (fill_between, 低于cutoff用灰色表示)
3.5 画CNV deletion (fill_between, 负值,向内画,低于cutoff用灰色表示)
3.6 因为这里并用不到inter-chromosome interaction,我只是随便画了2个link作为例子。注:染色体内部的控制点我偷懒,直接用了半径,所以靠近染色体的弧线有些平(以后有空了改下)。
最终效果图:
推荐阅读
-
在 ts 中实现类 java hashmap 的简单方法
-
Java 8新特性探究(十三)JavaFX 8新特性以及开发2048游戏-JavaFX历史## 跟java在服务器端和web端成绩相比,桌面一直是java的软肋,于是Sun公司在2008年推出JavaFX,弥补桌面软件的缺陷,请看下图JavaFX一路走过来的改进 从上图看出,一开始推出时候,开发者需使用一种名为JavaFX Script的静态的、声明式的编程语言来开发JavaFX应用程序。因为JavaFX Script将会被编译为Java bytecode,程序员可以使用Java代码代替。 JavaFX 2.0之后的版本摒弃了JavaFX Script语言,而作为一个Java API来使用。因此使用JavaFX平台实现的应用程序将直接通过标准Java代码来实现。 JavaFX 2.0 包含非常丰富的 UI 控件、图形和多媒体特性用于简化可视化应用的开发,WebView可直接在应用中嵌入网页;另外 2.0 版本允许使用 FXML 进行 UI 定义,这是一个脚本化基于 XML 的标识语言。 从JDK 7u6开始,JavaFx就与JDK捆绑在一起了,JavaFX团队称,下一个版本将是8.0,目前所有的工作都已经围绕8.0库进行。这是因为JavaFX将捆绑在Java 8中,因此该团队决定跳过几个版本号,迎头赶上Java 8。 ##JavaFx8的新特性 ## ###全新现代主题:Modena 新的Modena主题来替换原来的Caspian主题。不过在Application的start方法中,可以通过setUserAgentStylesheet(STYLESHEET_CASPIAN)来继续使用Caspian主题。 参考http://fxexperience.com/2013/03/modena-theme-update/ ###JavaFX 3D 在JavaFX8中提供了3D图像处理API,包括Shape3D (Box, Cylinder, MeshView, Sphere子类),SubScene, Material, PickResult, LightBase (AmbientLight 和PointLight子类),SceneAntialiasing等。Camera类也得到了更新。从JavaDoc中可以找到更多信息。 ###富文本 强化了富文本的支持 ###TreeTableView ###日期控件DatePicker 增加日期控件 ###用于 CSS 结构的公共 API
-
挂载在Linux中的作用及实现方法独到详解
-
配置HP FC存储设备多路径驱动的方法在linux和Windows操作系统中即可实现
-
快速傅里叶变换在Python中的实现及其代码解析
-
使用matplotlib在Python中实现正弦信号的时域波形和频谱图示例
-
在JSP中引入Java和Vue的实现方法
-
在Bootstrap中实现水平和垂直方向的div居中方法
-
在HTML中实现div元素垂直居中的三种方法
-
在div中实现文字垂直居中的方法