Mh ddlZddlZejejdejeZdZdZdZ dZ dZ d Z dS) Nz)%(asctime)s - %(levelname)s - %(message)s)levelformatc||d|dz |z }|d|dz |z }|d|d|dd| gS)zI This function computes the resolution based on data dimensions. r)extentnlinesncolsresxresys ./home/jovyan/projeto_inpe/dados_merge/utils.pygetGeoTrsO 1Iq !U *D 1Iq !V +D 1ItQq Au 55c||dz |dz }||dz |dz }t|t|fS)a; This function returns the grid point associated to a spatial location. - param x The spatial x-coordiante. - param y The spatial y-coordiante. - param geoTransform A list of 6 coefficients describing an affine transformation to georeference a grid. - return The grid point line and column. rrr)int)xygeoTlincols rgeo2gridrsG tAw;$q' !C tAw;$q' !C s88SXX rcZ||dz|dz}||dz|dz}||fS)a! This function returns the spatial location of a grid point. - param geoTransform A list of 6 coefficients describing an affine transformation to georeference a grid. - param col The grid point column. - param row The grid point row. - return The spatial location. rrrrr )rrrrrs rgrid2geor#s9 d1g QA d1g QA a4Krc|d}|d|d|dzz}|d|d|dzz}|d}||||fS)zo This function returns the extent based on the given GDAL geo-transform parameters and dimensions. rrrrr )gtshapellxllyurxurys r getExtentr$/s[ Q%C Q%E!H1% &C Q%E!HA& 'C Q%C c3 rct|||}t|d5}dD]0}|d||1 ddddS#1swxYwYdS)Nw)rrrrrz{} )ropenwriter)pathr r r rfileis rcreateWorldFiler-:s 665 ) )D dC/D# / /A JJv}}T!W-- . . . . ///////////////////s4A$$A(+A(c tj|stj|dSdS#t $r+}td|d|Yd}~dSd}~wwxYw)z. It Create the directory if it does not exist.TFz%Directory created or already exists: zError creating directory z: N)osr*existsmakedirsloggerinfoOSErrorerror) directoryr5s rcreate_directory_if_not_existsr7Bs w~~i((  K " " "45  EEEeEEFFFuuuuus39 A. A))A.) loggingr/ basicConfigINFO getLogger__name__r2rrrr$r-r7r rrr=s ',/Z[[[[  8 $ $666          ///     r