tf<ddlZddlZddlmcmZdZddZddZy)Nc0g}|D]}|jd}t|dkr tdt|d}t|d}t|d}t ||||j d}||t j|fgz }|S)N.zNot enough tile info present)splitlen IndexErrorinttilesreshapehgfoformat_dense_tile) gridtile_ids tile_valuestile_idpartszxy ret_arrays `/var/www/html/software/conda/envs/higlass/lib/python3.12/site-packages/clodius/tiles/npmatrix.py tiles_wrapperrsK F c" u:>;< < aM aM aM$1a(00"6 $"8"8"CDEE  F cd}t|j}tjtj||z tjdz }|dkrdn|}d|z|z}d|z|z}||z }|1|d|dg}|d|dg} | d|dz |z}| d|dz |z}nddg}|j} t |jdkDrt d|j|||| |d|dS) z+ Get the tileset info for the grid rrrrz)Grid's shape is not conducive to plottingfalse) max_width max_width1min_posmax_posmax_zoom mirror_tilesbins_per_dimension)maxshapemathceillogr ValueError) rboundsbin_sizemax_dimr%r!r"scale_upr#r$s r tileset_infor2sH$**oGyy'H"45 CDHqLqhHX (Ih)J7"H !9fQi(!9fQi(QZ'!*,8 aj71:-9 a&** 4::DdjjQQ & rct|j}tjtj||z tjdz }|dkrdn|}d||z z|z}||z} ||z} t |jd| |z} t |jd| |z} d||z z} || | | | f}| tjtd|jd| z z}| tjtd|jd| z z}||jdz }||jdz }t j|d|fd|ffdt jt jf}t j|j|jdd| fd}t j|jj|jdd| dj}t j||dk(<||| | | | f}t j|d|fd|ffdt jt jf}t j|j|jdd| fd}t j|jj|jdd| dj}||z }||jdz }||jdz }t j|d|fd|ffdt jt jfS) aR Return tiles at the given positions. Parameters ----------- grid: np.array An nxn array containing values z: int The zoom level (0 corresponds to most zoomed out) x: int The x tile position y: int The y tile position bin_size: int The number of values per bin rrrconstant)constant_valuesr )axisg) r(r)r*r+r,minnppadnannansumrT)rrrrnan_gridr/r0r% tile_widthx_starty_startx_endy_end num_to_sumdatadivisible_x_widthdivisible_y_widthdivisible_x_paddivisible_y_padabr not_nan_datananb norm_arrayx_pady_pads rrrDs"$**oGyy'H"45 CDHqLqhHx!|$x/J*nG*nG  1 w3 4E  1 w3 4Ex!|$J  wu}, -D #TYYs1djjm/Dz/Q%RR"TYYs1djjm/Dz/Q%RR'$**Q-7O'$**Q-7O  _ ?34(  A !))QWWQZZ89BA !##++aggaj"jAJLLI"$&&Ii3   wu} <= VV  !A#7 8 VVRVV,   YYrzz288A;J"?@q IYYrtt||BHHQKZHqQSS  *  yq) )E yq) )E 66 UaZ (  r)N)Nr) r*numpyr8clodius.tiles.formatrformatrrr2rrrUs" ##(%PSr