tilupy 0.1.4__py3-none-any.whl → 1.0.0__py3-none-any.whl

tilupy/utils.py

@@ -11,45 +11,47 @@ import matplotlib.pyplot as plt
 import shapely.geometry as geom
 import shapely.ops
 
+
 def CSI(pred, obs):
-    
-    ipred = pred>0
-    iobs = obs>0
-    
-    TP = np.sum(ipred*iobs)
-    FP = np.sum(ipred*~iobs)
-    FN = np.sum(~ipred*iobs)
-    
-    return TP/(TP+FP+FN)
-
-def diff_runout(x_contour, y_contour, point_ref, section=None,
-                orientation='W-E'):
-    
+    ipred = pred > 0
+    iobs = obs > 0
+
+    TP = np.sum(ipred * iobs)
+    FP = np.sum(ipred * ~iobs)
+    FN = np.sum(~ipred * iobs)
+
+    return TP / (TP + FP + FN)
+
+
+def diff_runout(
+    x_contour, y_contour, point_ref, section=None, orientation="W-E"
+):
     npts = len(x_contour)
-    contour = geom.LineString([(x_contour[i],y_contour[i])
-                               for i in range(npts)])
+    contour = geom.LineString(
+        [(x_contour[i], y_contour[i]) for i in range(npts)]
+    )
     point = geom.Point(point_ref)
     if section is None:
         return point.distance(contour)
     elif isinstance(section, np.ndarray):
         section = geom.LineString(section)
-        
-    assert(isinstance(section, geom.LineString))
+
+    assert isinstance(section, geom.LineString)
     section = revert_line(section, orientation)
     intersections = section.intersection(contour)
     if isinstance(intersections, geom.MultiPoint):
         intersections = geom.LineString(section.intersection(contour))
     intersections = np.array(intersections.coords)
-    if orientation == 'W-E':
-        i = np.argmax(intersections[:,0])
-    if orientation == 'E-W':
-        i = np.argmin(intersections[:,0])
-    if orientation == 'S-N':
-        i = np.argmax(intersections[:,1])
-    if orientation == 'N-S':
-        i = np.argmin(intersections[:,1])
-    intersection = geom.Point(intersections[i,:])
-    
+    if orientation == "W-E":
+        i = np.argmax(intersections[:, 0])
+    if orientation == "E-W":
+        i = np.argmin(intersections[:, 0])
+    if orientation == "S-N":
+        i = np.argmax(intersections[:, 1])
+    if orientation == "N-S":
+        i = np.argmin(intersections[:, 1])
+    intersection = geom.Point(intersections[i, :])
+
     #######
     # plt.figure()
     # cont = np.array(contour.coords)
@@ -61,71 +63,97 @@ def diff_runout(x_contour, y_contour, point_ref, section=None,
     #           pt[:,0], pt[:,1], 'o',
     #           inter[:,0], inter[:,1],'x')
     #######
-    
+
     return section.project(intersection) - section.project(point)
-    
-def revert_line(line, orientation='W-E'):
-    
+
+
+def revert_line(line, orientation="W-E"):
     pt_init = line.coords[0]
     pt_end = line.coords[-1]
-    if orientation == 'W-E':
-        if pt_init[0]>pt_end[0]:
+    if orientation == "W-E":
+        if pt_init[0] > pt_end[0]:
             line = shapely.ops.substring(line, 1, 0, normalized=True)
-    elif orientation == 'E-W':
-        if pt_init[0]<pt_end[0]:
+    elif orientation == "E-W":
+        if pt_init[0] < pt_end[0]:
             line = shapely.ops.substring(line, 1, 0, normalized=True)
-    if orientation == 'S-N':
-        if pt_init[1]>pt_end[1]:
+    if orientation == "S-N":
+        if pt_init[1] > pt_end[1]:
             line = shapely.ops.substring(line, 1, 0, normalized=True)
-    elif orientation == 'N-S':
-        if pt_init[1]<pt_end[1]:
+    elif orientation == "N-S":
+        if pt_init[1] < pt_end[1]:
             line = shapely.ops.substring(line, 1, 0, normalized=True)
-    
+
     return line
 
-def get_contour(x,y,z,zlevels,indstep=1,maxdist=30, closed_contour=True):
-    
-    #Add sup_value at the border of the array, to make sure contour
-    #lines are closed
+
+def get_contour(x, y, z, zlevels, indstep=1, maxdist=30, closed_contour=True):
+    # Add sup_value at the border of the array, to make sure contour
+    # lines are closed
     if closed_contour:
         z2 = z.copy()
         ni = z2.shape[0]
         nj = z2.shape[1]
         z2 = np.vstack([np.zeros((1, nj)), z2, np.zeros((1, nj))])
-        z2 = np.hstack([np.zeros((ni+2, 1)), z2, np.zeros((ni+2, 1))])
-        dxi = x[1]-x[0]
-        dxf = x[-1]-x[-2]
+        z2 = np.hstack([np.zeros((ni + 2, 1)), z2, np.zeros((ni + 2, 1))])
+        dxi = x[1] - x[0]
+        dxf = x[-1] - x[-2]
         dyi = y[1] - y[0]
         dyf = y[-1] - y[-2]
-        x2 = np.insert(np.append(x, x[-1]+dxf), 0, x[0]-dxi)
-        y2 = np.insert(np.append(y, y[-1]+dyf), 0, y[0]-dyi)
+        x2 = np.insert(np.append(x, x[-1] + dxf), 0, x[0] - dxi)
+        y2 = np.insert(np.append(y, y[-1] + dyf), 0, y[0] - dyi)
     else:
         x2, y2, z2 = x, y, z
-    
+
     backend = plt.get_backend()
-    plt.switch_backend('Agg')
+    plt.switch_backend("Agg")
     plt.figure()
-    ax=plt.gca()    
-    cs=ax.contour(x2, y2, np.flip(z2, 0), zlevels)
-    nn1=1
-    v1=np.zeros((1,2))
-    xcontour={}
-    ycontour={}
+    ax = plt.gca()
+    cs = ax.contour(x2, y2, np.flip(z2, 0), zlevels)
+    nn1 = 1
+    v1 = np.zeros((1, 2))
+    xcontour = {}
+    ycontour = {}
     for indlevel in range(len(zlevels)):
         for p in cs.collections[indlevel].get_paths():
-            if (p.vertices.shape[0]>nn1):
-                v1=p.vertices
-                nn1=p.vertices.shape[0]
-        xc = [v1[::indstep,0]]
-        yc = [v1[::indstep,1]]
+            if p.vertices.shape[0] > nn1:
+                v1 = p.vertices
+                nn1 = p.vertices.shape[0]
+        xc = [v1[::indstep, 0]]
+        yc = [v1[::indstep, 1]]
         if maxdist is not None and not closed_contour:
-            ddx=np.abs(v1[0,0]-v1[-1,0])
-            ddy=np.abs(v1[0,1]-v1[-1,1])
-            dd=np.sqrt(ddx**2+ddy**2)
-            if (dd>maxdist) :
-                xc[0]=None
-                yc[0]=None
-        xcontour[zlevels[indlevel]]=xc[0]
-        ycontour[zlevels[indlevel]]=yc[0]
+            ddx = np.abs(v1[0, 0] - v1[-1, 0])
+            ddy = np.abs(v1[0, 1] - v1[-1, 1])
+            dd = np.sqrt(ddx**2 + ddy**2)
+            if dd > maxdist:
+                xc[0] = None
+                yc[0] = None
+        xcontour[zlevels[indlevel]] = xc[0]
+        ycontour[zlevels[indlevel]] = yc[0]
     plt.switch_backend(backend)
-    return xcontour,ycontour
+    return xcontour, ycontour
+
+
+def format_path_linux(path):
+    """
+    Change a Windows-type path to a path formatted for Linux. \\ are changed
+    to /, and partitions like "C:" are changed to "/mnt/c/"
+
+    Parameters
+    ----------
+    path : string
+        String with the path to be modified.
+
+    Returns
+    -------
+    path2 : string
+        Formatted path.
+
+    """
+    if path[1] == ":":
+        path2 = "/mnt/{:s}/".format(path[0].lower()) + path[2:]
+    else:
+        path2 = path
+    path2 = path2.replace("\\", "/")
+    if " " in path2:
+        path2 = '"' + path2 + '"'
+    return path2

{tilupy-0.1.4.dist-info → tilupy-1.0.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: tilupy
-Version: 0.1.4
+Version: 1.0.0
 Summary: Thin-layer models unified processing tool
 Author-email: Marc Peruzzetto <m.peruzzetto@brgm.fr>
 License: CeCILL-C FREE SOFTWARE LICENSE AGREEMENT
@@ -529,6 +529,7 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: seaborn
 Requires-Dist: requests
+Requires-Dist: scipy
 Provides-Extra: dev
 Requires-Dist: pipreqs ; extra == 'dev'
 Provides-Extra: gis
@@ -559,7 +560,7 @@ It contains one submodule per thin-layer model for writing and reading raw input
 Outputs are then easily compared between different simulations / models. The models themselves are not part of this package and must
 be installed separately.
 
-Note that `tilupy` is still under development, thus only minimal documentation is available at the moment and testing is underway.
+Note that `tilupy` is still under development, thus only minimal documentation is available at the moment, and testing is underway.
 Contributions are feedback are most welcome. Reading and writing is available for the `SHALTOP` model (most commonly used by the author) and `r.avaflow`
 (only partly maintained).
 

tilupy-1.0.0.dist-info/RECORD

@@ -0,0 +1,26 @@
+tilupy/__init__.py,sha256=4r-aUNht-u_rUoYIyz3BOFEew11xCURWpSRYEfc8Dis,515
+tilupy/calibration.py,sha256=AGZq1YMI60d0RYzto1xdXglZ8DfSRCfo5gCAHzvwfwg,3318
+tilupy/cmd.py,sha256=Ly5HArNtgpZJ1IBBwjgpgJjcvMc_Pf9iQz2azEZZtk8,4046
+tilupy/compare.py,sha256=5PkDmMen_-yJImowvFJYB3IV0x0xs9WSFAYhX81wFLg,6294
+tilupy/download_data.py,sha256=D_R9l5XtgsmxEoKOXXPlWFV6q3e9SMa7asRRy-blQvw,1057
+tilupy/initdata.py,sha256=3vk-pw5IzCOKimBrGiT1iBxtqp2aJ9uqAY3GN7LFPYw,2579
+tilupy/make_mass.py,sha256=XM4HzwVLho_QDZzoXxVQCVJwNeMWawg-nRaDxiN9vW0,3442
+tilupy/make_topo.py,sha256=6fby5Zkl657pNgPFNWFF37-sORQvyBKsYn9atKESB_I,14302
+tilupy/notations.py,sha256=nSv49NbhlhD-zMBWFE6vk1UojCMH1R05oqrV8_6bCqk,10390
+tilupy/plot.py,sha256=O4eF71mYfX0YmzqsBvwIQTFVLhg4HZCFq4NedeqFyxs,21252
+tilupy/raster.py,sha256=4d_-RWGHHO0ZXkjwaHPEOBxPjb36HW8ISq_7z32sgQg,3605
+tilupy/read.py,sha256=Z93_zbaJ1bhnSibyQHA_aCWWge62zZy9WPTGKM5TPNA,27178
+tilupy/utils.py,sha256=N3qGhjOvVoQl9ySqxNsUzpE0W0gin1v6z94v-2Mya0Q,4705
+tilupy/models/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tilupy/models/ravaflow/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tilupy/models/ravaflow/initsimus.py,sha256=VCr_k2HuYecbX7KYS54FtY109OCzlslJ9pY1J1_iWNQ,5104
+tilupy/models/ravaflow/read.py,sha256=xKuMFRHyC0SV6QubMOcQ9EpJ3c843g2RKViDq7g-yHY,7900
+tilupy/models/shaltop/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tilupy/models/shaltop/initsimus.py,sha256=lmAr_rKozsv1CW0hfHsIi7me_X0ZuwG36F1cBdN0emE,8597
+tilupy/models/shaltop/read.py,sha256=w6c-1eh-SEQ73_W9Gr2jUOGf0L2ElV6nRxncVsI4O3w,13161
+tilupy-1.0.0.dist-info/LICENSE,sha256=bNrodROEaIAk5t7X7Db_7d2CjDfHBp_FjweStZO2gZk,21863
+tilupy-1.0.0.dist-info/METADATA,sha256=_7us2zgoWLJ0LI9gSApA4AEHZnqNUqSEIEu6-oS8R18,41875
+tilupy-1.0.0.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
+tilupy-1.0.0.dist-info/entry_points.txt,sha256=jf71F8J1P5s-zVt5BIxKZDmCIaoUqryMCtT7tc3BkkA,104
+tilupy-1.0.0.dist-info/top_level.txt,sha256=jhJY_A8Vzmqub_dVy0wPSfelUqFPjHMZOaepGVcrpxA,7
+tilupy-1.0.0.dist-info/RECORD,,

{tilupy-0.1.4.dist-info → tilupy-1.0.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.40.0)
+Generator: bdist_wheel (0.42.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

tilupy-0.1.4.dist-info/RECORD

@@ -1,24 +0,0 @@
-tilupy/__init__.py,sha256=4r-aUNht-u_rUoYIyz3BOFEew11xCURWpSRYEfc8Dis,515
-tilupy/calibration.py,sha256=O3UT-H9-eebaDttAN7GCQ4QDXynCb4IvwXSs8dHsL0E,3424
-tilupy/cmd.py,sha256=A_Xv9rq0e6vx4QIBQx_EunjDBKITXrZpQjQpjWu6tKA,4260
-tilupy/compare.py,sha256=5PkDmMen_-yJImowvFJYB3IV0x0xs9WSFAYhX81wFLg,6294
-tilupy/download_data.py,sha256=D_R9l5XtgsmxEoKOXXPlWFV6q3e9SMa7asRRy-blQvw,1057
-tilupy/initdata.py,sha256=UmUMff9zIu8df0s5prV6DVBD-KzBqdTUfMDJiDDzmQI,1334
-tilupy/notations.py,sha256=s4p9jsz0n2nXNxxVG5SU14j37b2kqYoIcTKvSehNOOI,3509
-tilupy/plot.py,sha256=bD1MeZyTH3hOYwbqto2mn1v49aUI2dJvagGy2qKwUwY,16552
-tilupy/raster.py,sha256=4d_-RWGHHO0ZXkjwaHPEOBxPjb36HW8ISq_7z32sgQg,3605
-tilupy/read.py,sha256=vyhbqBeh3ade8sXcCZK8ZZIjjVX7Dr9E7h1JhwY8mZw,12602
-tilupy/utils.py,sha256=P4OyIGggrUVC0aKqTlZuwGhUN-nyORI870S9QFWWJAA,4118
-tilupy/models/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-tilupy/models/ravaflow/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-tilupy/models/ravaflow/initsimus.py,sha256=VCr_k2HuYecbX7KYS54FtY109OCzlslJ9pY1J1_iWNQ,5104
-tilupy/models/ravaflow/read.py,sha256=xKuMFRHyC0SV6QubMOcQ9EpJ3c843g2RKViDq7g-yHY,7900
-tilupy/models/shaltop/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-tilupy/models/shaltop/initsimus.py,sha256=X0FSwwDSHWh-ODEu_zTZMCpB4cR4yFQpiDAF5ufNXxM,4600
-tilupy/models/shaltop/read.py,sha256=LetbNT19KvOJx0gdMOnlHxr0KlAg8TN3FmfBB9NIct0,9209
-tilupy-0.1.4.dist-info/LICENSE,sha256=bNrodROEaIAk5t7X7Db_7d2CjDfHBp_FjweStZO2gZk,21863
-tilupy-0.1.4.dist-info/METADATA,sha256=uP3gHp-Qq6dV7ydTSALDWiLh5r9AEJZ4XsTCblVw1QM,41853
-tilupy-0.1.4.dist-info/WHEEL,sha256=pkctZYzUS4AYVn6dJ-7367OJZivF2e8RA9b_ZBjif18,92
-tilupy-0.1.4.dist-info/entry_points.txt,sha256=jf71F8J1P5s-zVt5BIxKZDmCIaoUqryMCtT7tc3BkkA,104
-tilupy-0.1.4.dist-info/top_level.txt,sha256=jhJY_A8Vzmqub_dVy0wPSfelUqFPjHMZOaepGVcrpxA,7
-tilupy-0.1.4.dist-info/RECORD,,