i>/ddlZddlmZmZmZddl m Z ddZ dZ ddZdZGd d Zdd d d d dZdddd dZdZdZdZdZddZdZdZdZdZdS)N) as_bin_strtake endian_token)wkbTcL|t||dS)am Dump GeoJSON-like `dict` to GeoPackage binary and write it to the `dest_file`. :param dict obj: A GeoJSON-like dictionary. It must at least the keys 'type' and 'coordinates'. :param dest_file: Open and writable file-like object. :param bool big_endian: specify endianess of the dumped object. :return: N)writedumps)obj dest_file big_endians c/Users/user/workspace/sujinbaek/cqa-test-app/venv/lib/python3.11/site-packages/geomet/geopackage.pydumprs&OOE#z**+++++cDt|S)a Load a GeoJSON `dict` object from a ``source_file`` containing GeoPackage (as a byte string). :param source_file: Open and readable file-like object. :return: A GeoJSON `dict` representing the geometry read from the file. )loadsread) source_files r loadr+s !!## $ $$rc\t|| }tj||d}||zS)a Dump a GeoJSON-like dict to a GeoPackage bytestring. If the dict contains a top-level 'meta' key like so: ``` 'meta': {'srid': 4326} ``` then the srid will be added to the geopackage header, but *not* to the WKB geometry header. If the dict contains a top-level 'bbox' key like so: ``` 'bbox': [0, 0, 3, 3] ``` Then an envelope will be added to the geopackage header with this information. If the geometry's coordinates are empty (an empty list) then the geopackage header's "empty" flag will be set, denoting that this geometry has no coordinates. Please note that while this library can parse geopackages with a mixed byte-order in the header, it will only produce blobs with consistent byte order (albeit properly marked as such). That means you cannot product a geopackage with e.g. little-endian header and big-endian WKB geometry. :param dict obj: The geojson geometry to dump :param bool big_endian: if True, the geopackage binary will use big-endian byte order, little-endian otherwise. :return bytes: bytestring representing the geometry in geopackage format. F) include_meta)_build_geopackage_header_wkb_dumps)r r headerresults r r r 9s6X&cz> : :F [ju = = =F F?rct|}tttj|}t |t |\}}}}}}}t|} | tjz } tt| |} | rt|| |} tj |} |r!dt|i| d<ddd|zid| d<| r| | d<| S)a Construct a GeoJSON `dict` from geopackage (string). This function strips the geopackage header from the string and passes the remaining WKB geometry to the `geomet.wkb.loads` function. The envelope, if present, is added to the GeoJSON as a key called 'bbox' as per the GeoJSON spec, [1]. If an SRID is specified in the geopackage header AND the wkb header, the SRID in the geopackage header will take precedence and will replace that SRID in the returned dict. [1] https://tools.ietf.org/html/rfc7946#section-5 :param bytes string: geopackage byte string. :return dict: GeoJSON represented the parsed geopackage binary. sridmetanamezEPSG%s)type propertiescrsbbox) iter _as_bin_str_take _GeoPackage HEADER_LEN_check_is_valid _parse_header_get_wkb_offset_parse_enveloperrint)stringrgpversionemptyenvelope_indicatoris_little_endianr wkb_offset left_to_take envelope_dataenvelopers r rrjs.&\\F {5v>> ? ?FFfEAq'5,.>!!344J!77LlF ; ;<eZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zd S)r'zy Much more information on geopackage structure can be found here: http://www.geopackage.org/spec/#gpb_format GPrBBBBI 0@N)__name__ __module__ __qualname____doc__MAGIC1MAGIC2VERSION1r(HEADER_PACK_FMTENVELOPE_2D_LENENVELOPE_3D_LENENVELOPE_4D_LEN ENVELOPE_MASKEMPTY_GEOM_MASKENDIANNESS_MASKrr r'r'sQ F FHJOOOOM O OOOrr'r<)rrCrSrCrU)rrSrTr<ct|dd\}}}}}}}|tjks|tjkrdS|tjkrdS|dks|dkrdSdS) a( Check if the data represents a valid geopackage geometry. Input can be either the full geometry or just the header. :param bytes data: bytes representing the geopackage binary. :return (bool, str): Is the geopackage valid, if not, string describing why Nr<)Fz%Missing Geopackage header magic bytes)FzGeopackage version must be 0rrS)Fz&Envelope indicator must be between 0-4)T)r*r'rHrIrJ)datar/r0r1_r3s r is_validr[s2?tBQBx1H1H.Aq'1(!Q [ Q+*<%<%<==+&&&44Q$6$:$:>> 8rc^tjd|dd\}|tjzS)a Check to see if the header is encoded as little endian or big endian. Either the entire binary blob or just the header can be passed in. :param bytes header: geopackage header or binary blob :return bool: is the header little endian BrVrS)_structunpackr'rQ)rflagss r _header_is_little_endianras-~c6!A#;//HU ;. ..rct|}t|tjz}t j||dtj\}}}}}t|\}} } ||||| | |fS)ac Unpack all information from the geopackage header, including "magic" GP bytes. Returns all of them so we can confirm that this geopackage is validly formed. Can also accept the full binary blob. :param header: the header or the full geometry. :return 7-tuple: all attributes stored in the binary header. N)ra _endian_tokenr'rKr^r_r( _parse_flags) rr4fmtr/r0r1r`rr2r3 endiannesss r r*r*s077 ( ) )K,G GC!( V+[++ ,""Aq'5$-9,?,?)E z a%!3Z EErcr|tjz}|tjzdz }|tjzdz }|||fS)a Parse the bits in the "flags" byte of the geopackage header to retrieve useful information. We specifically parse the endianness, the envelope indicator, and the "empty" flag. Much more info can be found in the documentation [1]. [1] http://www.geopackage.org/spec/#gpb_format :param byte flags: The "flags" byte of a geopackage header. :return tuple: rCrS)r'rQrOrP)r`rfr3r2s r rdrd sD"44J+";;A [0 0Q 6E $j 00rc8d}|r|dzdz}|r||z}|dz|zS)aU Create the "flags" byte which goes into the geopackage header. Much more info can be found in the documentation [1]. [1] http://www.geopackage.org/spec/#gpb_format :param int empty: 0 or 1 indicating whether the geometry is empty. True and False also work as expected. :param int envelope_indicator: indicates the dimensionality of the envelope. :param int is_little_endian: 0 or 1 (or False / True) indicating whether the header should be little-endian encoded. :return byte: geopackage header flags rrCrVrR)r2r3r4r`s r _build_flagsri$sB* E !q +** QJ* **rcBt|ddkrdnd}|dg}|didd} tt|}n#t$rt dwxYwt jt jt jt||||g}t|t j z}|dt|zz }| |tj|g|RS) ax Create the geopackage header for the input object. Looks for a 'bbox' key on the geometry to use for an envelope, and a 'meta' key with an SRID to encode into the header. :param dict obj: a geojson object :param bool is_little_endian: which endianness to use when encoding the data. :return bytes: geopackage header. coordinatesrrCr#rrzoBounding box must be of length 2*n where n is the number of dimensions represented in the contained geometries.d)lenget_dim_to_indicatorKeyError ValueErrorr'rHrIrJrircrK_indicator_to_dimextendr^pack)r r4r2r8rr3 pack_argspack_fmts r rrBs2 S'((A--AA1Ewwvr""H 7762   " "61 - -D9.s8}}= 999899 99  U.0@AA  I-..1LLH ();<<=H X < -9 - - --s A::BcRt|\}}|std|zdS)z| Raise if the header is not valid geopackage. :param bytes data: Geopackage data or header. :return None: z6Could not read Geopackage geometry because of errors: N)r[rq)rYvalidreasons r r)r)rsDTNNME6 9/17899 999rc@tj}|t|z|zS)a Get the full byte offset at which the WKB geometry lies in the geopackage geometry. :param int envelope_indicator: indicates the dimensionality of the envelope. :return int: number of bytes until the beginning of the WKB geometry. )r'r(rr)r3base_lens r r+r+s#%H ();< < HHrcpt|}|dt|zz }tj||S)an Parse a geopackage envelope bytestring into an n-tuple of floats. :param int envelope_indicator: indicates the dimensionality of the envelope. :param bytes envelope: Bytestring of the envelope values. :param bool is_little_endian: how to pack the bytes in the envelope. :return tuple[float]: Geometry envelope. rl)rcrrr^r_)r3r8r4rvs r r,r,s:-..H ();<<=H >(H - --r)T)rC)structr^ geomet.utilrr%rr&rrcgeometrrrrr rr'rrror[rar*rdrirr)r+r,rRrr rs ,,,,$ % % %....b555p!!!!!!!!2    ,///"FFF0111.++++<-.-.-.` 9 9 9III".....r