jG"ddlZddlZddlZddlZddlmZmZmZmZm 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,dS)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||_||_dS)Nrr ) geotransformindicessuper__init___raster_prop)selfrpropxy __class__s r/var/www/finuniver-perm.ru/html/student/venv/lib/python3.11/site-packages/django/contrib/gis/gdal/raster/source.pyr+zTransformPoint.__init__*sc   T 21 5 6   T 21 5 6 !Q     c|dSNrr.s r3r0zTransformPoint.x1 Awr4cl|jj}|||j|jd<||j_dSr6r,r(r)r-r.valuegtfs r3r0zTransformPoint.x54l'+0DL $Q '($' !!!r4c|dSNr r7r8s r3r1zTransformPoint.y;r9r4cl|jj}|||j|jd<||j_dSrAr;r<s r3r1zTransformPoint.y?r?r4) __name__ __module__ __qualname__r)r+propertyr0setterr1 __classcell__r2s@r3rr#sG XX((X( XX((X(((((r4rceZdZdZejZd dZfdZdZ dZ dZ e dZ ed Ze d Zed Ze d Ze d Ze dZejdZe dZejdZe dZejdZe dZe dZe dZe dZe dZd!dZd"dZ d#dZe dZ xZ!S)$ GDALRasterz0 Wrap a raster GDAL Data Source object. Fc |rdnd|_tjt|tr(t j|rtj|}t|ttfrt |}| ts1tj |std|z t!jt%||j|_dS#t$r(}td||d}~wwxYwt|t*rd|_t-j|}t1||_tj t6t t9j}t!jt%|t?|j|t@ t!jt%||j|_dS#t$r1t!j!t%|tdwxYwt|tDrAd|_t|#dd}|j$dkr+d|vr'td |j$d |vsd |vrtd d |vrtdg}|#di%D]T\}} d|| } |&| '(U|&dtStU|z|}t!j+|jt%|#dd|d |d |#dtU|#dg|#ddt?||_tY|#dgD]\} } |j-| } d| vr=| d| _.| j.)d| vsd| vsd| vr| /| j.fd| /| #d| #d| #d| #d|#d |_0d|vr|d\|j1_2|j1_3d|vr|d\|j4_2|j4_3d |vr!|d \|j5_2|j5_3dSdSt|tlr ||_dStd!to|)"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 'aaa  """ h $ $ ,)9()C)C ,z(++H hd , ,} 8}}H&&'<== bgnnGG $>I  LX)>)> LL     #AHHSVWW % ( (n DK=**D&h//DNw||$@#djllBSBSTTH  0H%%dn%%)     Y LX)>)> LL   Y Y Y$[%:%:;;;#$WXXX Y$ ' 'U DKHLL599::F{e##h(>(>#KRR h&&((*B*B#Q X%%#$JKKKM$LL"==CCEE > >S S11$$V\\^^%:%:%<%<====   & & &&M(:(::]KM HLL4455!" ]C Wb0I0I,J,JKK Z++m$$DI"+8<<+D+D!E!E   :z!}!Z//(2>(BD%(4j00!Z//"j00 (9';6 JJJ #//#//$..11%>>(33  ||F++DH8##/7/A, t{}(""-5g->* djl!!+3F+;( TY[[["! ( + +  DIII7>>tH~~NN s$,D D3 #D..D3:,H((;I#c|jr&tjt|jt dSN) is_vsi_basedrirurrOr*__del__)r.r2s r3rzGDALRaster.__del__sC   9  TY!7!7 8 8 8 r4c|jSr)rOr8s r3__str__zGDALRaster.__str__s yr4cJdtt|jzS)zJ Short-hand representation because WKB may be very large. z)hexrrkr8s r3__repr__zGDALRaster.__repr__s!'Yty-A-A)B)BBBr4cd|jstdtj|jdS)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  di     r4c |jr|jtsdSt }t jt|jt|t}t||j Sr) rrOrerrriget_mem_buffer_from_vsi_filerrrr r=)r. out_lengthdats r3 vsi_bufferzGDALRaster.vsi_buffersz   "&)"6"67S"T"T 4WW /  " " *   %   j.///r4cN|jo|jtSr)rkrOrerr8s r3rzGDALRaster.is_vsi_baseds yHTY112GHHHr4cNttj|jS)zi Return the name of this raster. Corresponds to filename for file-based rasters. )rriget_ds_descriptionrkr8s r3rOzGDALRaster.names 0;;<< Return the GDAL Driver used for this raster. )ri get_ds_driverrkr )r. ds_drivers r3rMzGDALRaster.drivers$ &ty11 i   r4c4tj|jS)z+ Width (X axis) in pixels. )ri get_ds_xsizerkr8s r3rPzGDALRaster.width  +++r4c4tj|jS)z, Height (Y axis) in pixels. )ri get_ds_ysizerkr8s r3rQzGDALRaster.heightrr4c tj|j}|sdSt|dS#t$rYdSwxYw)zF Return the SpatialReference used in this GDALRaster. Nwkt)srs_type)riget_ds_projection_refrkrr)r.rs r3rzGDALRaster.srs&s[  ,TY77C t#C%888 8   44 s00 >>c6t|tr|}n;t|ttfrt|}nt dt j|j|j | dS)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- . . PCC Sz * * P"5))CCNOO O "49cgnn.>.>??? r4c|jjS)zA Shortcut to access the srid of this GDALRaster. )rrRr8s r3rRzGDALRaster.sridCs x}r4c||_dS)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.geotransformQs:!|  E#JJ777Cyyr4c t|dkstd|Dstdtdz|}t j|jt||dS)z)Set the geotransform for the data source.rXc3NK|] }t|ttfV!dSr)r`rfloat).0r0s r3 z*GDALRaster.geotransform..`s0&S&Sqz!c5\'B'B&S&S&S&S&S&Sr4z.Geotransform must consist of 6 numeric values.N) r|allrrriset_ds_geotransformrkrr)r.valuess r3r(zGDALRaster.geotransform]sz v;;!  3&S&SF&S&S&S#S#S MNN NQ,(  E&MM::: r4c"t|dS)z3 Coordinates of the raster origin. r$rr8s r3r$zGDALRaster.origings dH---r4c"t|dS)z@ Pixel scale in units of the raster projection. r%rr8s r3r%zGDALRaster.scalens dG,,,r4c"t|dS)z7 Skew of pixels (rotation parameters). r&rr8s r3r&zGDALRaster.skewus dF+++r4cb|jj|jj|jzz}|jj|jj|jzz}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 {}tz|dj88{}tz|dk994''4''4''4''T4%%r4c t|Srrr8s r3rVzGDALRaster.bandss~~r4NearestNeighbourc |d|j|d|j|d|jj|d|j|d|j|d|j|d|jj d|vr|j d z|jj z|d<d |vr"|j d  |d <d |j D|d <t|d}t|}tj|j|jj|j|jj|d|t't't' ||S)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.rWrc g|] }d|ji SrYr)rbnds r3 z#GDALRaster.warp..s VVVCnc.>?VVVr4rVTrr) setdefaultrPrQrrRr$r%r&rMrOrVrWrKrrireproject_imagerkrr{r r)r.r resampling max_errortarget algorithms r3warpzGDALRaster.warps$ GTZ000Hdk222FDHM222Hdk222GTZ000FDI...Hdk&6777  ! !#y83dk6FFHV  X % %#':a=#9#9#;#;HZ WV4:VVVHD111-Z8   I HL   ! ! K JN ! ! # #    JJ JJ JJ    r4Nc |r|}nk|jjdkr|jdz|jjz}nCtjt t tj}ttj |jj t||j ttt!t!|jS)z"Return a clone of this GDALRaster.rNrr)rMrOrfrgrqrrarrrsrKricopy_dsrkrrrr r^)r.rO clone_names r3clonezGDALRaster.clones  WJJ [  & &X- 0@@JJ&BC DUDUVVJ L  J''      +    r4c t|}t|tr|}n;t|ttfrt|}nt d|j|jkr'|r||jjkr| |Stj |j |j j|j||t!}t#|}|j|j|j|jj|jjg|jj|jjg|jj|jjgd} |r|| d<|r|| d<|| ||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.transforms}-Z8 c+ , , JJ c3Z ( ( )#..JJ  ?di ' ' '6T[EU;U;U::d## #, I HL   ! ! N ! ! # #   JJ   F##O\m} 8lnfln5[]FKM2     $#DN  DLyy* yJJJr4cZtj|jdS)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% $//66888r4)F)rrr)NNrr)"rCrDrE__doc__riclose_ds destructorr+rrrrrFrrrrOrMrPrQrrGrRr(r$r%r&rrVrrrrrHrIs@r3rKrKFsJGGGGR CCC ! ! !00X0 II_I==X=!!_!,,X, ,,X,   X  Z  Z X  [[   X ..X. --X- ,,X,  & &X &X;;;;z    ,UX3K3K3K3Kj99X99999r4rK)-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.functionalrrrrKr7r4r3rs                     111111777777======888888>>>>>>GFFFFFFF22222288888888333333 ( ( ( ( (T ( ( (FW9W9W9W9W9W9W9W9W9W9r4