DUf4ddlZddlZddlZddlmZddlmZddlmZe dZ e dZ d Z d Z d Zejejd d ejdeejddedejdddddejdddejde d ejd!d"ejgd#d$d%&ejd'd(ddd)ejd*ed+ejd,ed-ejd.d/d01ejd2d3dd45d6ZdS)7N)util)Cooler)cligחAgcAc||d||dz }||d||dz }||zS)Nrr)extent)c row_region col_regionnrowsncolss L/var/www/html/software/conda/lib/python3.11/site-packages/cooler/cli/show.pyget_matrix_sizersd HHZ  #ahhz&:&:1&= =E HHZ  #ahhz&:&:1&= =E 5=c|||||}|dkrtj|}n|dkrtj|}|S)N)balancefieldlog2log10)matrixfetchnprr)r r r rbalancedscalemats r load_matrixrs\ ((85( 1 1 7 7 J O OC gcll '  hsmm Jrc ddlm d d   fd}tjdjj g d jd| dS)NrcD|jj\}}}}||z}||z}||||fS)N)viewLimbounds)axxstartystartxdeltaydeltaxendyends r get_extentzinteractive..get_extent's6)+):&tVT))rc ttj|d|z |z}ttj|d|z |z}ttj|d|z |z}ttj|d|z |z}t d|}t d|}t |ttj||z |z}t |ttj||z |z}||||fS)Nrrr)intrfloorceilmaxmin)r binsize row_chrom_len col_chrom_lenr#r'r$r(s rround_trim_extentz&interactive..round_trim_extent-sRXfQi'122W<==276!9w.//'9::RXfQi'122W<==276!9w.//'9::QQ4RXmg&=>>HIIJJ4RXmg&=>>HIIJJtT6))rc d}|  }|dkrdS|d< t|dt|df}t|dt|df}jd}t ||}|tkr|t jddddft jzdsN d dd gdddgd d \} d d d d d j }||gd<ncdr4d d4|t || ||jjdS)NF prev_extentrrr+r placeholderskg?)r lwz@The requested region is too large to display at this resolution.center)horizontalalignmentverticalalignment transform)set_autoscale_onr,imagesrMAX_MATRIX_SIZE_INTERACTIVEset_dataronesnanplottext transAxespopremover set_extentfigurecanvas draw_idle)eventr new_col_regionnew_row_regionimnelemboxtxtr"rr1r col_chromr3rr) plotstatepltr4 row_chromr2rs rupdate_heatmapz#interactive..update_heatmap:s E"""B""67M=QQ Y}- - - F#) - "Cq NNCq NNB"Cq NNCq NNB Yr]>>BB / / / KK 111d7+bf4 5 5 5^, 7 }aMJ q!]MJ " hh5(0&. l.1#J .)N+ 9.)--//66888N+ 9 KKA~~uhPUVV    f ""$$$$$rzbin-size)r8r6button_release_event) matplotlib.pyplotpyplotr,info chromsizesgcagcfrL mpl_connectshow)r"r rXrUrrrrYr1r3r)rVrWr4r2s``````` @@@@@@@r interactivercs$##### *** * * *-%-%-%-%-%-%-%-%-%-%-%-%-%-%-%-%-%-%^!&$%%GL+ML+M!#JJswwyy4I4IJJIGGII  !7HHHHHJJJJJrcool_uri COOL_PATH)metavarrange)typez--range2z-r2zThe coordinates of a genomic region shown along the column dimension. If omitted, the column range is the same as the row range. Use to display asymmetric matrices or trans interactions.)rhhelpz --balancedz-bTFzQShow the balanced contact matrix. If not provided, display the unbalanced counts.)is_flagdefaultriz--outz-ozSave the image of the contact matrix to a file. If not specified, the matrix is displayed in an interactive window. The figure format is deduced from the extension of the file, the supported formats are png, jpg, svg, pdf, ps and eps.)riz--dpiz*The DPI of the figure, if saving to a filez--scalez-slinearrrzNScale transformation of the colormap: linear, log2 or log10. Default is log10.r)rhrirkz--forcez-fzkForce display very large matrices (>=10^8 pixels). Use at your own risk as it may cause performance issues.z--zminzThe minimal value of the color scale. Units must match those of the colormap scale. To provide a negative value use a equal sign and quotes, e.g. -zmin='-0.5'z--zmaxzThe maximal value of the color scale. Units must match those of the colormap scale. To provide a negative value use a equal sign and quotes, e.g. -zmax='-0.5'z--cmapYlOrRdzThe colormap used to display the contact matrix. See the full list at http://matplotlib.org/examples/color/colormaps_reference.html)rkriz--fieldcountzPixel values to display.)rk show_defaultric P ddl} || dddlm} n?#t$r2t dt jt jdYnwxYwt|}|j }|}||n|}tj ||\}}}tj ||\}}}t|||tkr1|s/t dt jt jd| d | d | d | t)|||| ||d ||||g|| | | |d| |d| }|dddd||r| ||dSt5| |||| ||dS)z Display and browse a cooler in matplotlib. COOL_PATH : Path to a COOL file or Cooler URI. RANGE : The coordinates of the genomic region to display, in UCSC notation. Example: chr1:10,000,000-11,000,000 rNAggz%Install matplotlib to use cooler show)filerzThe matrix of the selected region is too large. Try using lower resolution, selecting a smaller region, or use the '--force' flag to override this safety limit.) )figsizezContact matrixnone) interpolationr vminvmaxcmapz coordinatezrelative contact frequencyz$log 2 ( relative contact frequency )z%log 10 ( relative contact frequency )rl)dpi) matplotlibuser[r\ ImportErrorprintsysstderrexitrr^r parse_regionrMAX_MATRIX_SIZE_FILErKget_current_fig_managerset_window_titletitleimshowrylabelxlabelcolorbar set_labelsavefigrcr_)rdrgrange2routr}rforcezminzmaxr|rmplrWr r^r r rXrow_lorow_hirUcol_locol_hicbs rrbrbqs^    ? GGENNN'''''''  5CJGGGG   xAJJ%~6J $ 1*j I IIvv $ 1*j I IIvv :z226JJJ K  @      JJxJ   !!223CDDDIIbMMMJJAz:uhFF/    JJ)((()))JJ)((())) BLL2:<      P CS !!!!!CGGIIq)YxOOOOOs!$9A A )rclicknumpyrrwrapirrr,rrArrrccommandargumentstroptionChoicefloatrbrrrs s3xx!c#hh OOOd  K000c"""   @    6   @gC&RSSS   111 2 2       ?   Q  Q  Y  wT8RDPDP TS#"10FDPDPDPr