tfXddlZddlZddlZddlZddlZdZdZdZdZ dZ dZ d dZ y) NcZt|||\}}t||dz|dz\}}||||fS)zbconvert Google-style Mercator tile coordinate to (minlat, maxlat, minlng, maxlng) bounding box)get_lng_lat_from_tile_pos)zoomxyminlngminlatmaxlngmaxlats [/var/www/html/software/conda/envs/higlass/lib/python3.12/site-packages/clodius/tiles/geo.py get_tile_boxrs?/tQ:NFF.tQUAEBNFF FFF ++ctjdtjz|zd|zz z }t|}|dztjz }dd|zd|zz z}||fS)zZconvert Google-style Mercator tile coordinate to (lng, lat) of top-left corner of tilef@gfv@)mathpi gudermannian)rrrlat_radlatlngs r rrsf ggDGG a1955G7#G E/DGG #C 519T * *C :rc|tjzdz }t|}d|z|dzzdz }d|ztj|z zdtjzz }||fS)z[convert lng/lat to Google-style Mercator tile coordinate (x, y) at the given zoom levelrrr)rrinv_gudermannian)rrrrrrs r get_tile_pos_from_lng_latr"sdDGGme#Gw'G T S5[!E)A T TWWw&'1tww;7A q6Mrcdtjtj|ztjdz z SNr)ratanexpr)rs r rr0s+ tyy!% %! 33rctjtj|tjdz zdz Sr)rlogtanr)rs r rr4s+ 88DHHa$''A+o23 44rc tjj|sddiStj|}|j dj }|d|d|d|d|dg|d |d g|dd|dzzd }|S) NerrorzTileset info is not available!zSELECT * FROM tileset_inforrr) zoom_step tile_sizemax_zoommin_posmax_posmax_data_length)ospathisfilesqlite3connectexecutefetchone)filepathdbresos r tileset_infor;8s 77>>( #9::  "B **1 2 ; ; =CVVFFCF#FCF#q6AQK/  A HrcBtj|}|j}t|||\}} } } t|||zdz ||zdz \} } } } d} |j | ||| | | fj }t jt}t|||zD]}t|||zD] }g|||f< |D]} |djd} |djd}t|d|d|\}}t|d|d|\}} tj|d } tj|d }t|||zD]]}t|||zD]I}||dzks ||k\s||dzks||k\s!|||fxx|d|d|d|d|d ||||d gz cc<K_|j!|S#t$r |d}YwxYw#t$r |d}YwxYw#t$rd }YwxYw#t$rd }YwxYw)a Retrieve a contiguous set of tiles from a 2D db tile file. Parameters ---------- db_file: str The filename of the sqlite db file zoom: int The zoom level x: int The x position of the first tile y: int The y position of the first tile width: int The width of the block of tiles to retrieve height: int The height of the block of tiles to retrieve Returns ------- tiles: {pos: tile_value} A set of tiles, indexed by position raO SELECT minLng, maxLng, maxLat, minLat, uid, importance, geometry, properties, intervals.id FROM intervals, position_index WHERE intervals.id=position_index.id AND zoomLevel <= ? AND rMaxLng >= ? AND rMinLng <= ? AND rMinLat <= ? AND rMaxLat >= ? r(zutf-8rrr&r)Nr') xStartxEndyStartyEnd importanceuidgeometry propertiesid)r3r4cursorrr5fetchallcol defaultdictlistrangedecodeAttributeErrorrjsonloads Exceptionclose)db_filerrrwidthheightconnclng_from_lat_fromlng_tolat_toqueryrowsnew_rowsijrrDrGx_starty_startx_endy_endrErFs r get_tilesrhLs0 ??7 #D A+D!Q7Ha1'a%i!mQZ!^LAvq& E 99UT8VXvF G P P RDt$H 1a%i "q!f*% "A!HaV  ""* A$++g&C 1W%B5QqT1Q4F01qtTB u zz!A$'H  AaD)J q!e)$ A1a&j) QU?uzgAo%ST*aV$&'d$%aD&'d$%aD*+A$#&(0*4"$  ) $   3*V JJL OU A$C  1B  H    J  sH G G*!G?:HG'&G'*G<;G<? H  H  HH)rr) rPrr0r3 collectionsrJrrrrrr;rhrr rks8 ,  45 (ir