sgspy 1.0.2__cp313-cp313-manylinux_2_28_x86_64.whl

sgspy/utils/vector.py

@@ -0,0 +1,268 @@
+# ******************************************************************************
+#
+#  Project: sgs
+#  Purpose: GDALDataset wrapper for vector operations
+#  Author: Joseph Meyer
+#  Date: June, 2025
+#
+# ******************************************************************************
+
+import sys
+import os
+import site
+import tempfile
+from typing import Optional
+import warnings
+
+import matplotlib.pyplot as plt
+import matplotlib #fpr type checking matplotlib.axes.Axes
+
+from.import plot
+from .plot import plot_vector
+
+#ensure _sgs binary can be found
+site_packages = list(filter(lambda x : 'site-packages' in x, site.getsitepackages()))[0]
+sys.path.append(os.path.join(site_packages, "sgspy"))
+sys.path.append(os.path.join(os.path.dirname(__file__), ".."))
+from _sgs import GDALRasterWrapper
+
+from _sgs import GDALVectorWrapper
+
+try:
+    import geopandas as gpd
+    GEOPANDAS = True
+except ImportError as e:
+    GEOPANDAS = False
+
+PROJDB_PATH = os.path.join(sys.prefix, "sgspy")
+
+##
+# @ingroup user_utils
+# This class represents a spatial vector, and is used as an input to many sgs functions.
+#
+# It has a number of additional uses, including displaying info about the vector, converting 
+# to a GDAL or GeoPandas object. 
+# 
+# Accessing vector info:
+#     
+# vector metadata can be displayed using the info() function. All layers
+# are displayed unless a specific layer is specified. The per-layer info
+# includes: name, number of features, number of fields, geomtype, and bounds.
+# 
+# Public Attributes:
+# --------------------
+# layer_names : list[str] @n
+#     a list of layer names
+# 
+# Public Methods:
+# --------------------
+# info() @n
+#     takes an optional argument specify the band, and prints vector metadata to console
+class SpatialVector:
+
+    def __init__(self, 
+                 image: str | GDALVectorWrapper):
+        """
+        Constructing method for the SpatialVector class.
+
+        Has one required parameter to specify a gdal dataset. The following
+        attributes are populated:
+        self.cpp_vector
+        self.layer_names
+
+        Parameters
+        --------------------
+        image: str | GDALVectorWrapper
+           specifies a path to a vector file or the C++ class object itself
+        """
+        if type(image) is str:
+            self.cpp_vector = GDALVectorWrapper(image, PROJDB_PATH)
+        elif type(image) is GDALVectorWrapper:
+            self.cpp_vector = image
+        else:
+            raise TypeError("'image' parameter to SpatialVector constructor must be of type str or GDALVectorWrapper.")
+
+        self.layers = self.cpp_vector.get_layer_names()
+
+    def print_info(self, 
+                   layer_name: str, 
+                   layer_info: dict):
+        """
+        prints layer information using the layer_info from self.cpp_vector.
+
+        This is an internal function not meant to be used by the end user.
+
+        Parameters
+        --------------------
+        name : str
+            str containing the layer name
+        layer_info : dict
+            dict containing 'feature_count', 'field_count', 'geometry_type', 'xmax', 'xmin', 'ymax', and 'ymin' items
+        """
+        print("{} layer info:".format(layer_name))
+        print("feature count: {}".format(layer_info['feature_count']))
+        print("field count: {}".format(layer_info['field_count']))
+        print("geometry type: {}".format(layer_info['geometry_type']))
+        print("bounds (xmin, xmax, ymin, ymax): ({}, {}, {}, {})".format(
+            layer_info['xmin'], 
+            layer_info['xmax'], 
+            layer_info['ymin'],
+            layer_info['ymax']
+        ))
+        if layer_info['crs']: print("crs: {}".format(layer_info['crs']))
+        print()
+
+    def info(self, 
+             layer: Optional[int | str] = None):
+        """
+        calls self.print_info depending on layer parameter. If no layer is given,
+        print all layers. A layer may be specified by either a str or an int.
+
+        Parameters
+        --------------------
+        layer : str or int
+            specifies the layer to print information on
+        """
+        if layer is not None and type(layer) not in [int, str]:
+            raise TypeError("'layer' parameter, if given, must be of type int or str.")
+
+        if type(layer) == str:
+            self.print_info(layer, self.cpp_vector.get_layer_info(layer))
+        elif type(layer) == int:
+            self.print_info(self.layers[layer], self.cpp_vector.get_layer_info(self.layers[layer]))
+        else:
+            for layer in self.layers:
+                self.print_info(layer, self.cpp_vector.get_layer_info(layer))
+
+    def samples_as_wkt(self):
+        """
+        Calls get_wkt_points on the underlying cpp class, to return
+        the samples as wkt strings. 
+
+        This function requires that there be a layer named 'samples' which
+        is comprised entirely of Points or MultiPoints. These conditions
+        will be satisfied if this SpatialVector is the output of one of the
+        sampling functions in the sgs package.
+        """
+        if "samples" not in self.layers:
+            print("this vector does not have a layer 'samples'")
+        else:
+            return self.cpp_vector.get_wkt_points('samples')
+
+    def plot(self,
+        geomtype: str,
+        ax: Optional[matplotlib.axes.Axes] = None,
+        layer: Optional[int | str] = None, 
+        **kwargs):
+        """
+        Calls plot_vector on self.
+
+        Paramters
+        --------------------
+        ax : matplotlib.axes.Axes
+            axes to plot the raster on
+        geomtype : str
+            the geometry type to try to print
+        layer : None | int | str
+            specification of which layer to print
+        **kwargs
+            any parameter which may be passed ot matplotlib.pyplot.plot
+        """
+
+        if ax is not None: 
+            plot_vector(self, ax, geomtype, layer, **kwargs)
+        else:
+            fig, ax = plt.subplots()
+            plot_vector(self, ax, geomtype, layer, **kwargs)
+            plt.show()
+
+    @classmethod
+    def from_geopandas(cls, obj, layer_name: str=None):
+        """
+        This function is used to convert a geopandas object into an sgspy.SpatialVector. The geopandas object
+        may either by of type GeoDataFrame or GeoSeries.
+
+        If a particular layer name is desired, it can be passed as a parameter.
+
+        Examples:
+
+        gdf = gpd.read_file("access.shp")
+        access = sgspy.SpatialVector.from_geopandas(gdf)
+
+        
+        gs = gpd['geometry'] #geometry column is a geoseries
+        access = sgspy.SpatialVector.from_geopandas(gs)
+
+
+        gdf = gpd.read_file("access.shp")
+        gdf = gdf[gdf == "LineString"]
+        access = sgspy.SpatialVector.from_geopandas(gdf)
+        """
+        if layer_name is not None and type(layer_name) is not str:
+            raise TypeError("layer_name, if given, must be of type 'str'.")
+
+        if not GEOPANDAS:
+            raise RuntimeError("from_geopandas() can only be called if geopandas was successfully imported, but it wasn't.")
+
+        if type(obj) is not gpd.geodataframe.GeoDataFrame and type(obj) is not gpd.geoseries.GeoSeries:
+            raise TypeError("the object passed must be of type geopandas GeoDataFrame or GeoSeries")
+
+        #get a the geopandas object as a geoseries
+        if type(obj) is gpd.geodataframe.GeoDataFrame:
+            if 'geometry' not in obj.columns:
+                raise RuntimeError("'geometry' must be a column in the geodataframe passed")
+            gdf = obj
+            gs = gdf['geometry']
+        else:
+            gs = obj
+            gdf = gpd.GeoDataFrame(gs)
+ 
+        if layer_name is None: layer_name = 'geopandas_geodataframe'
+        projection = gs.crs.to_wkt()
+       
+        # the conversion to geojson may raise a warning about how the projection is unable to be converted to the new format correctly
+        #
+        # however, we have the projection from the geodataframes geoseries crs, so it won't be an issue
+        with warnings.catch_warnings():
+            warnings.simplefilter("ignore")
+            geojson = gdf.to_json().encode('utf-8')
+
+        cpp_vector = GDALVectorWrapper(geojson, projection, layer_name, PROJDB_PATH)
+        return cls(cpp_vector)
+
+    def to_geopandas(self):
+        """
+        This function is used to convert an sgspy.SpatialVector into a geopandas geodataframe.
+
+        Examples:
+
+        access = sgspy.SpatialVector("access.shp")
+        gdf = access.to_geopandas()
+        """
+        if not GEOPANDAS:
+            raise RuntimeError("to_geopandas() can  only be called if geopandas was successfully imported, but it wasn't.")
+
+        tempdir = tempfile.gettempdir()
+        file = os.path.join(tempdir, "temp.geojson")
+
+        #get the projection info
+        projection = self.cpp_vector.get_projection()
+
+        #write the dataset to a tempfile
+        self.cpp_vector.write(file)
+    
+        # This method of writing to a file then reading from that file is definitely clunky,
+        # however it's easy. Theres the possiblity of iterating through every field within
+        # every feature, and needing to then call a different function depending on the data
+        # type of the field (because C++ types are rigid). That may still be done in the future,
+        # but for now this works.
+
+        #have geopandas read the tempfile
+        gdf = gpd.read_file(file)
+        if projection != "": gdf.set_crs(projection, inplace=True, allow_override=True)
+
+        #remove the geojson file
+        os.remove(file)
+
+        return gdf
+        

sgspy-1.0.2.dist-info/METADATA

@@ -0,0 +1,13 @@
+Metadata-Version: 2.1
+Name: sgspy
+Version: 1.0.2
+Requires-Dist: numpy
+Requires-Dist: matplotlib
+Requires-Dist: daal==2025.9.0
+Requires-Dist: tbb==2022.*
+Requires-Dist: mkl==2025.3.0
+Requires-Dist: onemkl-license==2025.3.0
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Requires-Dist: geopandas; extra == "test"
+

sgspy-1.0.2.dist-info/RECORD

@@ -0,0 +1,40 @@
+sgspy/__init__.py,sha256=ONzqVZBLUs0Zm4s-x09HDoVWICfjDeygpzBscn9JHA0,2582
+sgspy/_sgs.cpython-313-x86_64-linux-gnu.so,sha256=TfFOr1UZP8mCbG5KS7nL5opDzWu2WWOdVJeRDPalOdI,86466384
+sgspy/calculate/__init__.py,sha256=CA6GtS3Lym_TtYYlMp6j02uf57DDlq3vGzN8qXxYDuc,331
+sgspy/calculate/pca/__init__.py,sha256=L6jpsKWDSch-Rq4CqJSBO0elALd_yzuvTY2JASAipB8,39
+sgspy/calculate/pca/pca.py,sha256=yaaZitdTQ0p25REkpUQwnIxNi0kiCeOE8mkxWl8aoZg,5422
+sgspy/calculate/representation/__init__.py,sha256=yiDee7Cw6LnAbuATKMWXneVPBwhogV_5rIunK-USh7I,72
+sgspy/calculate/representation/representation.py,sha256=nKNNlWmmF1uFyM9KuD9Y8Qj-YfYGfpSRuPZ_zUpNalE,72
+sgspy/sample/__init__.py,sha256=Z20McdhG5A4J5qtwZ9O7CyH_tbXNKnYA54WINVwU-F0,422
+sgspy/sample/ahels/__init__.py,sha256=2oAo9Fkck1BRqyVNfzr9rbWvQ4Jt2lg72EG7SIVbKzU,45
+sgspy/sample/ahels/ahels.py,sha256=skPsJ3OhQAHBP9wHY3LGVlmnqkRlsidAFkv7CqpgWHM,63
+sgspy/sample/clhs/__init__.py,sha256=vddVyTqOvsBoODikmdValw-pwEfYohlHo4Wlz-MnyQE,42
+sgspy/sample/clhs/clhs.py,sha256=0F3Z4pIUZ0FeOPtK5WRO2ToNt7JF0M0Uh7Yw_yLU6uo,7528
+sgspy/sample/nc/__init__.py,sha256=c5P2cH42Q23-hBguVF81b8KbcaFypg1PBw0TGY7nf6o,36
+sgspy/sample/nc/nc.py,sha256=77sqsBpNYsyqVRoAoEnhAGM_WjQrNBnyHTdSUc67334,60
+sgspy/sample/srs/__init__.py,sha256=0_Qm4f_6iE1Ds4AHQp_8t9UZZ7I0m46xWRvL4o6pKow,39
+sgspy/sample/srs/srs.py,sha256=QqV1xdHo9a-jQv4C4nrtOpZdxVys_xmrDcYhR0LdwKA,7974
+sgspy/sample/strat/__init__.py,sha256=WYiabUGZckk38Qhr1fFXhS8sh9BqyXtGEBPJ9lstpJo,45
+sgspy/sample/strat/strat.py,sha256=KPzwzEj2HACUsOuMDrbUZOTM8GM3q6ldT6EAVDYuw2M,16171
+sgspy/sample/systematic/__init__.py,sha256=frl-x5lN3zY7Q5HzbvwT0zz50Byh7VDoSgYDcG_vKbA,60
+sgspy/sample/systematic/systematic.py,sha256=HPM2VTzFm2Q5Ge2Qhq-3yqL-u-KuiBCr2QzBiRykRVo,8336
+sgspy/stratify/__init__.py,sha256=RULBeQemsCmVjHu57Me2tptDIlhe64M_PtvCsRpXe4Q,399
+sgspy/stratify/breaks/__init__.py,sha256=45VAzOTJVGOn4kfU8w-mA2lq4PP5Z0PZjWbIt6RLj_k,48
+sgspy/stratify/breaks/breaks.py,sha256=ekkanqFzjE3cUhVO9vwnukodJdnaSf3HXzOKBV6zhwo,9099
+sgspy/stratify/kmeans/__init__.py,sha256=PjR6Uh8xhpId2HIl7vbysdduQpXO7l00FDmR9QAWa4M,48
+sgspy/stratify/kmeans/kmeans.py,sha256=unem1jiJAwzdfdGnCRYn-47eeZXJccCaT7nL9d5ViXs,64
+sgspy/stratify/map/__init__.py,sha256=uw8HivY0wg9Tke4XI8vtHCMNMYJOmncHb2Enl3jTg7s,71
+sgspy/stratify/map/map_stratifications.py,sha256=RYdpBsZV61Evo73OzGx2KUXC0IufcgsZ6iChXh2PvWQ,10652
+sgspy/stratify/poly/__init__.py,sha256=CfCIvJBQIqqS8aQG-52Ja9bRxrBN1jbcUBTkbkhD9KY,42
+sgspy/stratify/poly/poly.py,sha256=Dt_wTxshxw02lbKVjLgS1yFWh5qJbVklTnVqPdzrSW0,6233
+sgspy/stratify/quantiles/__init__.py,sha256=_g9MTbV7PIZxdshhEvJtZRLBCS5ehoUqZNCMBHrnOsQ,57
+sgspy/stratify/quantiles/quantiles.py,sha256=ZjMoJNt2WUkpE7Jsbs_nPQcNP9rvf4EMYUyVOmGr9Y4,12195
+sgspy/utils/__init__.py,sha256=XXEYEVdjUa345mCs0d5Bj7lVCyHEwMoDvselH-Deebs,287
+sgspy/utils/plot.py,sha256=9gnW0C8vK5AUDKBi00830V9i5ciXVXROTH7A1TLqN0c,5528
+sgspy/utils/raster.py,sha256=EOSZJOoWZ7IKo_iar4iHiWV74q_IT9HNQ4NRIfVoIQg,22945
+sgspy/utils/vector.py,sha256=vHkIg7v85ljgO6MsRt-cRxLVw_bAZQPbiRto8k4YQko,9576
+sgspy-1.0.2.data/data/sgspy/libonedal.so.3,sha256=sAmzMql47Yy9eO70uraCrAoAkJiYfYLMXs5CrwTr82I,10841849
+sgspy-1.0.2.data/data/sgspy/proj.db,sha256=sGOmBu2wzZwqDJ9V5yu9mez2dvZCrqH4izR1mp9Q0Lk,9359360
+sgspy-1.0.2.dist-info/METADATA,sha256=hiYB5ZO_eU0RMBYvkhNY45vq2MA07RhGMZb3EDhJJwQ,325
+sgspy-1.0.2.dist-info/WHEEL,sha256=22Lftwvjz7BjVNTxxDM8SyAm3GLe_fJvO71HQ9qo0Mw,99
+sgspy-1.0.2.dist-info/RECORD,,

sgspy-1.0.2.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: meson
+Root-Is-Purelib: false
+Tag: cp313-cp313-manylinux_2_28_x86_64
+