gGddlZddlZddlZddlZddlmZmZmZmZm Z m Z m Z m Z ddl mZddlmZddlmZddlmZddlmZddlmZdd lmZmZmZmZmZdd l m!Z!m"Z"dd l#m$Z$dd l%m&Z&m'Z'dd l(m)Z)Gdde*Z+GddeZ,y)N) addressofbyrefc_bufferc_char_pc_doublec_intc_void_p string_at)Path)Driver) GDALException)rasterBandList)GDALRasterBase)GDAL_RESAMPLE_ALGORITHMSVSI_DELETE_BUFFER_ON_READVSI_FILESYSTEM_PREFIXVSI_MEM_FILESYSTEM_BASE_PATHVSI_TAKE_BUFFER_OWNERSHIP)SpatialReference SRSException) json_regex) force_bytes force_str)cached_propertyceZdZddddZfdZedZejdZedZejd ZxZ S) TransformPoint)r)))originscaleskewc|j|j|d}|j|j|d}t| ||g||_||_y)Nrr ) geotransformindicessuper__init___raster_prop)selfrpropxy __class__s }/var/www/python.vincentserveurtest.ovh/public_html/venv/lib/python3.12/site-packages/django/contrib/gis/gdal/raster/source.pyr+zTransformPoint.__init__*s^    T 21 5 6    T 21 5 6 !Q   c |dSNrr.s r3r0zTransformPoint.x1 Awr4c|jj}|||j|jd<||j_yr6r,r(r)r-r.valuegtfs r3r0zTransformPoint.x59ll''+0DLL $Q '($' !r4c |dSNr r7r8s r3r1zTransformPoint.y;r9r4c|jj}|||j|jd<||j_yrAr;r<s r3r1zTransformPoint.y?r?r4) __name__ __module__ __qualname__r)r+propertyr0setterr1 __classcell__r2s@r3rr#spG XX(( XX((r4rceZdZdZej ZddZfdZdZ dZ dZ e dZ edZe d Zed Ze d Ze d Ze d Zej*dZe dZej*dZe dZej*dZe dZe dZe dZe dZe dZddZddZ ddZe dZ xZ!S) GDALRasterz0 Wrap a raster GDAL Data Source object. c |rdnd|_tjt|tr*t j |rtj|}t|ttfr}t |}|jts-tjj|std|z t!j"t%||j|_yt|t*rd|_t-j.|}t1||_tjj5t6t t9j:}t!j<t%|t?|j2|t@ t!j"t%||j|_yt|tDrd|_t|jGdd}|jHdk7r(d|vr$td j)|jHd |vsd |vr td d |vr tdg}|jGdijKD]D\}} dj)|| } |jM| jOjQF|jMdtStU|z|}t!jV|j&t%|jGdd|d |d |jGdtU|jGdg|jGddt?||_tY|jGdgD]\} } |jZ| } d| vr@| d| _.| j\*d| vsd| vsd| vr| j_| j\fd| j_| jGd| jGd| jGd| jGd|jGd |_0d|vr&|d\|jb_2|jb_3d|vr&|d\|jh_2|jh_3d |vr'|d \|jj_2|jj_3yyt|tlr||_ytd!j)to|#t$r }tdj)||d}~wwxYw#t$r*t!jBt%|tdwxYw)"Nr rz(Unable to read raster source input "%s".z+Could not open the datasource at "{}" ({}).z1Failed creating VSI raster from the input buffer.driverMEMnamez5Specify name for creation of raster with driver "{}".widthheightz;Specify width and height attributes for JSON or dict input.sridz$Specify srid for JSON or dict input. papsz_optionsz{}={} nr_of_bandsbandsdatatype nodata_valuedatasizeshape)r r )rZr\offset)rZr[r\r]r$r%r&z%Invalid data source input type: "{}".)8_writer ensure_registered isinstancestrrmatchjsonloadsr startswithrospathexistsr capiopen_dsr_ptrformatbytessys getsizeofr _ds_inputjoinruuiduuid4create_vsi_file_from_mem_bufferrrunlink_vsi_filedictgetrOitemsappendupperencoderlen create_ds enumeraterVrYrZsrsr$r0r1r%r&r type)r.ds_inputwriteerrr[vsi_pathrMrSkeyvaloptioni band_inputbands r3r+zGDALRaster.__init__Ms aa   " h $)9)9()Czz(+H hd ,8}H&&'<=bggnnG$>I  LLX)> L % (DK==*D&h/DNww||$@#djjlBSTH  0 0H%dnn%)   Y LLX)> L $ 'DKHLL59:F{{e#h(>#KRR h&((*B#Q X%#$JKKM$LL"=CCES S1$$V\\^%:%:%<=F   &&M(::]KM HLL45!" ]C Wb0I,JK Z+m$DI"+8<<+D!E :zz!}!Z/(2>(BD%((4j0!Z/"j0 (9(9';6 J #/#/$..1%>>(3 "F* ||F+DH8#/7/A, t{{}("-5g->* djjl!+3F+;( TYY[" ( + DI7>>tH~N a! #AHHSVW 2! Y$$[%:;#$WXX Ys$>.T !.T9 T6T11T693U,c|jr(tjt|jt |yN) is_vsi_basedrirurrOr*__del__)r.r2s r3rzGDALRaster.__del__s-     TYY!7 8 r4c|jSr)rOr8s r3__str__zGDALRaster.__str__s yyr4cDdtt|jzS)zJ Short-hand representation because WKB may be very large. z)hexrrkr8s r3__repr__zGDALRaster.__repr__s'Ytyy-A)BBBr4cp|js tdtj|jy)a  Flush all data from memory into the source file if it exists. The data that needs flushing are geotransforms, coordinate systems, nodata_values and pixel values. This function will be called automatically wherever it is needed. z9Raster needs to be opened in write mode to change values.N)r^r riflush_dsrkr8s r3_flushzGDALRaster._flushs-{{K  dii r4c|jr|jjtsyt }t j t|jt|t}t||jSr) rrOrerrriget_mem_buffer_from_vsi_filerrrr r=)r. out_lengthdats r3 vsi_bufferzGDALRaster.vsi_bufferse   $))"6"67S"TW //  " *  %  j..//r4c\|jxr|jjtSr)rkrOrerr8s r3rzGDALRaster.is_vsi_baseds!yyHTYY112GHHr4cRttj|jS)zi Return the name of this raster. Corresponds to filename for file-based rasters. )rriget_ds_descriptionrkr8s r3rOzGDALRaster.names 00;< Return the GDAL Driver used for this raster. )ri get_ds_driverrkr )r. ds_drivers r3rMzGDALRaster.drivers# &&tyy1 i  r4c@tj|jS)z+ Width (X axis) in pixels. )ri get_ds_xsizerkr8s r3rPzGDALRaster.width   ++r4c@tj|jS)z, Height (Y axis) in pixels. )ri get_ds_ysizerkr8s r3rQzGDALRaster.heightrr4c tj|j}|syt|dS#t$rYywxYw)zF Return the SpatialReference used in this GDALRaster. Nwkt)srs_type)riget_ds_projection_refrkrr)r.rs r3rzGDALRaster.srs&s@  ,,TYY7C#C%8 8  s!1 1 ==ct|tr|}n-t|ttfr t|}n t dt j |j|jj|jy)z Set the spatial reference used in this GDALRaster. The input can be a SpatialReference or any parameter accepted by the SpatialReference constructor. z/Could not create a SpatialReference from input.N) r`rintra ValueErrorriset_ds_projection_refrkrr{r)r.r=rs r3rzGDALRaster.srs3s` e- .C Sz *"5)CNO O ""499cggnn.>? r4c.|jjS)zA Shortcut to access the srid of this GDALRaster. )rrRr8s r3rRzGDALRaster.sridCs xx}}r4c||_y)zE Shortcut to set this GDALRaster's srs from an srid. N)r)r.r=s r3rRzGDALRaster.sridJs r4ctdz}tj|jt |t |S)z Return the geotransform of the data source. Return the default geotransform if it does not exist or has not been set previously. The default is [0.0, 1.0, 0.0, 0.0, 0.0, -1.0]. rX)rriget_ds_geotransformrkrlist)r.r>s r3r(zGDALRaster.geotransformQs2!|   E#J7Cyr4ct|dk7std|Ds tdtdz|}t j |j t||jy)z)Set the geotransform for the data source.rXc3HK|]}t|ttfywr)r`rfloat).0r0s r3 z*GDALRaster.geotransform..`s&SFqz!c5\'BFs "z.Geotransform must consist of 6 numeric values.N) r|allrrriset_ds_geotransformrkrr)r.valuess r3r(zGDALRaster.geotransform]sX v;! 3&SF&S#SMN NQ,(   E&M: r4ct|dS)z3 Coordinates of the raster origin. r$rr8s r3r$zGDALRaster.origings dH--r4ct|dS)z@ Pixel scale in units of the raster projection. r%rr8s r3r%zGDALRaster.scalens dG,,r4ct|dS)z7 Skew of pixels (rotation parameters). r&rr8s r3r&zGDALRaster.skewus dF++r4c|jj|jj|jzz}|jj|jj|j zz}t ||jj}t||jj}t ||jj}t||jj}||||fS)zJ Return the extent as a 4-tuple (xmin, ymin, xmax, ymax). )r$r0r%rPr1rQminmax)r.xvalyvalxminxmaxyminymaxs r3extentzGDALRaster.extent|s {{}}tzz||djj88{{}}tzz||dkk994'4'4'4'T4%%r4ct|Srrr8s r3rVzGDALRaster.bandss ~r4c |jd|j|jd|j|jd|jj|jd|j |jd|j |jd|j|jd|jjd|vr)|jd z|jjz|d<d |vr |jd j|d <|jDcgc]}d |jic}|d <t|d}t|}tj |j"|jj$j'|j"|jj$j'|d|t)t)t) |j+|Scc}w)a  Return a warped GDALRaster with the given input characteristics. The input is expected to be a dictionary containing the parameters of the target raster. Allowed values are width, height, SRID, origin, scale, skew, datatype, driver, and name (filename). By default, the warp functions keeps all parameters equal to the values of the original source raster. For the name of the target raster, the name of the source raster will be used and appended with _copy. + source_driver_name. In addition, the resampling algorithm can be specified with the "resampling" input parameter. The default is NearestNeighbor. For a list of all options consult the GDAL_RESAMPLE_ALGORITHMS constant. rPrQrRr$r%r&rMrO_copy.rWrrYrVTr) setdefaultrPrQrrRr$r%r&rMrOrVrWrYrKrrireproject_imagerkrr{r r)r.r resampling max_errorbndtarget algorithms r3warpzGDALRaster.warps$ GTZZ0Hdkk2FDHHMM2Hdkk2GTZZ0FDII.Hdkk&6&67  !#yy83dkk6F6FFHV  X %#'::a=#9#9#;HZ LP::V:Cnc.>.>?:VHD1-Z8   II HHLL   ! KK JJNN ! ! #    J J J    3Ws8H c |r|}n|jjdk7r'|jdz|jjz}n?tjj t t tj}ttj|jjt||jttt!t!|j"S)z"Return a clone of this GDALRaster.rNrr)rMrOrfrgrqrrarrrsrKricopy_dsrkrrrr r^)r.rO clone_names r3clonezGDALRaster.clones J [[   &X- 0@0@@J&BC DUVJ LL   J'     ++  r4c t|}t|tr|}n-t|ttfr t|}n t d|j |j k(r,|r||jjk(r|j|Stj|j|jjj|jj||t!}t#|}|j |j$|j&|j(j*|j(j,g|j.j*|j.j,g|j0j*|j0j,gd} |r|| d<|r|| d<|j3| ||S)zk Return a copy of this raster reprojected into the given spatial reference system. zETransform only accepts SpatialReference, string, and integer objects.)rRrPrQr$r%r&rMrO)rr)rr`rrra TypeErrorrRrMrOrriauto_create_warped_vrtrkrrr{r rKrPrQr$r0r1r%r&r) r.rrMrOrrr target_srsrrZs r3 transformzGDALRaster.transformse-Z8 c+ ,J c3Z ()#.J  ??dii '6T[[EUEU;U::d# #,, II HHLL   ! NN ! ! #   J  F#OO\\mm}} 8llnnfllnn5[[]]FKKMM2   #DN DLyy* yJJr4c^tj|jdjS)z Return information about this raster in a string format equivalent to the output of the gdalinfo command line utility. N)ri get_ds_infoptrdecoder8s r3infozGDALRaster.infos$ $/6688r4)F)NearestNeighbourrr)NNrr)"rCrDrE__doc__riclose_ds destructorr+rrrrrFrrrrOrMrPrQrrGrRr(r$r%r&rrVrrrrrHrIs@r3rKrKFsJGR C !00 II==!!,, ,,    ZZ    [[   .. -- ,,  & &;z ,UX3Kj99r4rK)-rcrfrnrrctypesrrrrrrr r pathlibr django.contrib.gis.gdal.driverr django.contrib.gis.gdal.errorr "django.contrib.gis.gdal.prototypesrri#django.contrib.gis.gdal.raster.bandr#django.contrib.gis.gdal.raster.baser$django.contrib.gis.gdal.raster.constrrrrrdjango.contrib.gis.gdal.srsrrdjango.contrib.gis.geometryrdjango.utils.encodingrrdjango.utils.functionalrrrrKr7r4r3rsh    17=8>G283 (T (FW9W9r4