tf)ddlZddlZddlZddlZddlZddlZddl m cm Z ddl mZmZddlmZdZdZdZdZej,eZiZddded<dd Zd Z dd Zidfd Z d Zy)N) abs2genomicget_quadtree_depth)ThreadPoolExecutor significant Significant)namevaluecDtj||}t|d<|S)N range_modes)hgbi tileset_infor )bbpath chromsizestis ^/var/www/html/software/conda/envs/higlass/lib/python3.12/site-packages/clodius/tiles/bigbed.pyrrs#  6: .B#B} Ic |\ }}}}}}}}} |j|} tjtj|| || } |dk(r| | } } n| | } } d}d}i}|D]}|||j|<||z }|dz } g}d}g}| s|S| D]%} |jt|d|dz }'||k\r5||kr0| D])} t|d}|j| ||d+|S||kDrg}||z }ttj|d|z }| D].} t|d}||k\r|j| ||d0t|}||kDr:tjt!||}t#|Dcgc]}|| }}|S#t $r d\} } YUt $r d\} } YewxYw#tt f$r|jtY:wxYw#tt f$r |j| t|dYYwxYw#tt f$r)t|k\r|j| t|dY6wxYwcc}w)Nr)NNrr) chrOffset importancefields)indexftpartialbbifetch_intervals IndexErrorKeyErrorappendint ValueError DEFAULT_SCOREnpquantilelenrandomsamplerangesorted)arbinsizer range_mode min_elements max_elementscidstartendchrom fetch_factory intervals intervals2offset offsetIdx chrOffsetschrSizefinal_intervalsintervals_lengthscoresintervalscorethresholded_intervals desired_percthresholded_scorethresholded_intervals_lengthindicesis r fetch_datarGs     (T+  % 3#6#6ucR  &$1O]_zI$1O]_zIFIJ28 :##I./'Q  O F  ) MM#hqk* + A <',< ,L" H HQK(&&%/%6&+"* ` ; L ( "#&66  FA 4D EF" H HQK(--)00)3E):*/&.  ((++@'A$ ', 6mmE*F$GVGAGQA4Q7QOQ e +!+ : + * :+,J' ) MM- ( )  + &&%/%6&3"* 0 +  $55)00)3E):*7&. RsTAF+ G&G;2+H- I(+G < G  G $G87G8;+H*)H*-4I%$I%c|tj|}|t}|t}tj|}||} t t |||} td5} t | jt| D cgc]} t|| ||||gt | z!c} } ddd Dcgc] }|gk7s | }}|Dcgc] }|D]}| c}}Scc} w#1swY$C#" C(9 C4C4C9#C((C1c d}d}g}|D]}|jddd}|jdd} | jd} ttt| dd} t | dkDr| dnt } | t vr| nd} t|Dcgc]}|jdc}}d |vrt|d }d |vrt|d }||kDr|}|}|}n ||k(r|}|dz}|dkr t}t}|rE|Dcgc]}|d }}|Dcgc]}t|d}}tj|| }nd}d |vr|d }||vr||}nd}| d}| d}|tj|}t|tj}tjd ||z zz}||z}||z}t!|||||| ||}|||fgz }|Scc}wcc}wcc}w)a  Generate tiles from a bigbed file. Parameters ---------- tileset: tilesets.models.Tileset object The tileset that the tile ids should be retrieved from tile_ids: [str,...] A list of tile_ids (e.g. xyx.0.0) identifying the tiles to be retrieved chromsizes_map: {uid: []} A set of chromsizes listings corresponding to the parameters of the tile_ids. To be used if a chromsizes id is passed in with the tile id with the `|cos:id` tag in the tile id chromsizes: [[chrom, size],...] A 2d array containing chromosome names and sizes. Overrides the chromsizes in chromsizes_map Returns ------- tile_list: [(tile_id, tile_data),...] A list of tile_id, tile_data tuples |rNr.:minmax)rcos)r.r/r0)splitrOrPr"r'DEFAULT_RANGE_MODEr dictrLrMpdSeriesrrKr TILE_SIZEr^)rtile_idschromsizes_maprr/r0generated_tilestile_idtile_option_partstile_no_options tile_id_parts tile_position return_valuer.o tile_optionstemp_min_elementsrX chromnames chromlengthschromsizes_to_use chromsizes_idrRtile_pos max_depth tile_sizerSrT tile_values rtilesrs]2LLOC3#MM#.qr2!--,Q/'--c2 SmAa&89: #M 2Q 6M! rs  ##21"     8 $ &3mL M Sv.=b,.$bJ%r