rbi!dZddlZddlZddlZddlmZddlmZddlmZddl m Z ddl m Z ddl mZdd lmZdd lmZmZd ed edzfd Zdeee fd edzfdZdeee fdededed ef dZ d%deee fdedzdedzdedzd ef dZded efdZdeee fd efdZdZ d&dZd'deezdzd eee ffd Z d(dedzdedzdedzdedzdeded#edeezdzd e!fd$Z"dS))z Based on the logic in the PROJ projsync CLI program https://github.com/OSGeo/PROJ/blob/9ff543c4ffd86152bc58d0a0164b2ce9ebbb8bec/src/apps/projsync.cpp N)datetime)partial)Path)Any) urlretrieve)get_proj_endpoint)BBox) get_data_dirget_user_data_dircoordsreturnc t|\}}tt|t|t|t|S#t$rYnwxYwd}|D]}t |}|||}t|j|j|_t|j|j|_t|j|j|_t|j |j |_ |S)z/ Get the bounding box from coordinates )westsoutheastnorthN) zipr minmax ValueError_bbox_from_coordsrrrr)r xxxyyy coord_bbox coord_setbboxs [/var/www/html/mdtn/previsions/meteo_cartes/venv/lib/python3.11/site-packages/pyproj/sync.pyrrs <SSSCCQQQQ      J > >  ++ <   JJ!*/49==JO":#3TZ@@J ":#3TZ@@J !*/49==JOO sAA A"!A"geomcVd|vsd|vrdS|d}|ddkrt|Sd}d}g}|D]D}t|}||jdkrd}n |jdkrd}||Ed}|D]}|r$|r"|jdkrd|_|xjd z c_||}-t |j|j|_t |j|j|_t |j|j|_t |j|j|_|S) z4 Get the bounding box from geojson geometry coordinatestypeN MultiPolygonFLTh)rrrappendrrrr)rr found_minus_180found_plus_180bboxescoordinate_setr grid_bboxs r_bbox_from_geomr,-sbD  F$$6$6t}%K F|~%% ---ON F% 00 <  9  "OO Y#  !N dI < <  ~ $)t2C2CDI II II  II ;;IN!)/4:>>IO!)/4:>>IO ;;INN featurer spatial_testinclude_world_coveragec|d}|t|}|dS|j|jz dkrA|j|jz dkr.|sdSt d |_t d|_n6|jdkr+|jdkr |xjd zc_|xjd zc_t|||SdS) zC Filter by the bounding box. Designed to use with 'filter' geometryNFiginfr$r#r%)getr,rrrrfloatgetattr)r.rr/r0r geom_bboxs r _filter_bboxr9Qs ;;z " "D #D))  5 NY^ +c 1 1)/1C77) u#Ell]IN"5\\INN ^c ! !di$&6&6 NNc !NN NNc !NN*wt\**9555 5r- source_id area_of_usefilenamec|d}|sdS|d}|d}|r|sdS|dup||v}|dup||ddv}|dup||v} |r|r| rdSdS) zA Filter by the properties. Designed to use with 'filter' propertiesFnamer:Nr;T)r5) r.r:r;r<r> p_filename p_source_idsource_id__matchedarea_of_use__matchedfilename__matcheds r_filter_propertiesrFks\**J u''J..--K [u"d*Fi;.F&$.+rBB3!D(BH ,B27Ht 5r- grid_namec tt|rdStt jD]'}t||rdS(dS)zX Run through all of the PROJ directories to see if the file already exists. FT)rr existsr splitospathsep)rGdata_dirs r_is_download_neededrNs    ++2244u NN((44 ) $ $ + + - - 55  4r-c|d}|sdS|d}|sdSt|S)zM Filter grids so only those that need to be downloaded are included. r>Fr?)r5rN)r.r>r<s r_filter_download_neededrPsN\**J u~~f%%H u x ( ((r-ctj}t|d5tfddD]}|| dddn #1swxYwY|S)z7 Return sha256 checksum of file given by path. rbc.dS)Ni)read)filesrz_sha256sum..s$))E"2"2r-r-N)hashlibsha256openiterupdate hexdigest) input_filehasherchunkrUs @r _sha256sumr`s ^  F j$  !42222C88 ! !E MM%  !!!!!!!!!!!!!!!!     s,AA"%A"Fc|rtd|t||d} t|||%|t|krt d||t|| t j|dS#t$rYdSwxYw# t j|w#t$rYwwxYwxYw)z. Download resource file from PROJ url z Downloading: z.partNzSHA256 mismatch: ) printrrr` RuntimeErrorreplacerKremoveFileNotFoundError)file_url short_name directoryverboserXtmp_paths r_download_resource_filerls * (h(()))I*33344H Hh'''  &Jx,@,@"@"@?:??@@ @i44555  Ih          DD   Ih         D sBAB+B B('B(+C-CC C CCCtarget_directoryc|td}t|d}|rHtjtj|jz jdkr!ttdd|tj | dS) zR Returns ------- dict[str, Any]: The PROJ grid data list. NTz files.geojsonrz/files.geojson)rgrhrizutf-8)encoding)r rrIrnow fromtimestampstatst_mtimedaysrlrjsonloads read_text)rm local_paths r_load_grid_geojsonrys,T22&88J      (01B1B1KLL LRUVVV)++;;;&&    :j**G*<< = ==r- intersectsTinclude_already_downloadedct|d}||jdkr1|j|jkr!|xjdzc_|xjdzc_n|jdkr1|j|jkr!|xjdz c_|xjdz c_nPt|jdkr8t|jdkr |j|jkr|xjdz c_t t t ||||}t t t||||}|rt|Stt t|S) a= Get a list of transform grids that can be downloaded. Parameters ---------- source_id: str, optional area_of_use: str, optional filename: str, optional bbox: BBox, optional spatial_test: str, default="intersects" Can be "contains" or "intersects". include_world_coverage: bool, default=True If True, it will include grids with a global extent. include_already_downloaded: bool, default=False If True, it will list grids regardless of if they are downloaded. target_directory: str | Path, optional The directory to download the geojson file to. Default is the user writable directory. Returns ------- list[dict[str, Any]]: A list of geojson data of containing information about features that can be downloaded. )rmfeaturesNr$r%r#)rr/r0)r:r;r<) ryrrabsfilterrr9rFtuplerP) r:r;r<rr/r0r{rmr}s rget_transform_grid_listrsqF"3CDDDZPH  9s??ty4944 II II II III Y  $)di"7"7 II II II III ^^c ! !c$)nns&:&:ty49?T?T II II )'=        #    H"X /:: ; ;;r-)NNN)FN)N)NNNNrzTFN)#__doc__rWrurKr functoolsrpathlibrtypingrurllib.requestr pyproj._syncr pyproj.aoir pyproj.datadirr r listrdictstrr,boolr9rFrNrPr`rlryrrr-rrs+  &&&&&&******::::::::dtd{0!$sCx.!TD[!!!!H #s(^#'7:TX 8!"  #s(^TztDj   4 3 4     )T#s(^ ) ) ) ) )   <@*>>td):>d3PS8n>>>>,!"$#'',*.B<B<TzB<tB<DjB< + B<  B< ! B<!%B<Dj4'B< B<B<B<B<B<B