f!.dZddlZddlZddlZddlmZddlmZddlmZddl m Z m Z ddl m Z ddlmZdd lmZdd lmZmZd ed e efd Zdeee fd e efdZdeee fdededed ef dZ d%deee fde ede ede ed ef dZded efdZdeee fd efdZdZ d&dZ d'd eee ffdZ! d(de ede ede ede ededed"ed#e ed 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)AnyOptional) 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 E/var/www/html/Alfredo/env/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/rr0r1r geom_bboxs r _filter_bboxr:Ns ;;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)r6) r/r;r<r=r? p_filename p_source_idsource_id__matchedarea_of_use__matchedfilename__matcheds r_filter_propertiesrGhs\**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)rHdata_dirs r_is_download_neededrOs    ++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@)r6rO)r/r?r=s r_filter_download_neededrQsN\**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_filehasherchunkrVs @r _sha256sumras ^  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: ) printrrra RuntimeErrorreplacerLremoveFileNotFoundError)file_url short_name directoryverboserYtmp_paths r_download_resource_filerms * (h(()))I*33344H Hh'''  &Jx,@,@"@"@?:??@@ @i44555  Ih          DD   Ih         D sBAB+B B('B(+C-CC C CCCc|td}t|d}|rHtjtj|jz jdkr!ttdd|tj | dS) zR Returns ------- dict[str, Any]: The PROJ grid data list. NTz files.geojsonrz/files.geojson)rhrirjzutf-8)encoding)r rrJrutcnow fromtimestampstatst_mtimedaysrmr jsonloads read_text)target_directory local_paths r_load_grid_geojsonrzs,T22&88J        X3JOO4E4E4NOO OU    )++;;;&&    :j**G*<< = ==r. intersectsTinclude_already_downloadedrxct|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) aI 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: Union[str, Path, None], 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. )rxfeaturesNr%r&r$)rr0r1)r;r<r=) rzrrabsfilterrr:rGtuplerQ) r;r<r=rr0r1r|rxr~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)NNNNr{TFN)$__doc__rXrurLr functoolsrpathlibrtypingrrurllib.requestr pyproj._syncr pyproj.aoir pyproj.datadirr r listrdictstrr-boolr:rGrOrQrarmrzrrr.rrs   &&&&&&******::::::::dx~,!$sCx.!Xd^!!!!H #s(^#'7:TX 8 $!%"  #s(^}#sm   4 3 4     )T#s(^ ) ) ) ) )   <@*>>c3h>>>>. $!%"$#'',&*B<B<}B<#B<smB< 4. B<  B< ! B<!%B<smB< B<B<B<B<B<B