tilupy 0.1.5__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.5.dist-info → tilupy-1.0.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: tilupy
-Version: 0.1.5
+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

{tilupy-0.1.5.dist-info → tilupy-1.0.0.dist-info}/RECORD

@@ -1,26 +1,26 @@
 tilupy/__init__.py,sha256=4r-aUNht-u_rUoYIyz3BOFEew11xCURWpSRYEfc8Dis,515
 tilupy/calibration.py,sha256=AGZq1YMI60d0RYzto1xdXglZ8DfSRCfo5gCAHzvwfwg,3318
-tilupy/cmd.py,sha256=A_Xv9rq0e6vx4QIBQx_EunjDBKITXrZpQjQpjWu6tKA,4260
+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=s4p9jsz0n2nXNxxVG5SU14j37b2kqYoIcTKvSehNOOI,3509
-tilupy/plot.py,sha256=awffgtd_X2ChltgnucIngVSeTNrwZ2YALLEhxYqxHdw,18044
+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=vyhbqBeh3ade8sXcCZK8ZZIjjVX7Dr9E7h1JhwY8mZw,12602
-tilupy/utils.py,sha256=P4OyIGggrUVC0aKqTlZuwGhUN-nyORI870S9QFWWJAA,4118
+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=bezToIzFE25pk6qHRyGcnMQ5Du-s_MVP_KjVZI8KVwI,4601
-tilupy/models/shaltop/read.py,sha256=LetbNT19KvOJx0gdMOnlHxr0KlAg8TN3FmfBB9NIct0,9209
-tilupy-0.1.5.dist-info/LICENSE,sha256=bNrodROEaIAk5t7X7Db_7d2CjDfHBp_FjweStZO2gZk,21863
-tilupy-0.1.5.dist-info/METADATA,sha256=LGpV1Uqymnq-BNbSGsZt9VHxkBP4heel40BObjEtzxA,41875
-tilupy-0.1.5.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-tilupy-0.1.5.dist-info/entry_points.txt,sha256=jf71F8J1P5s-zVt5BIxKZDmCIaoUqryMCtT7tc3BkkA,104
-tilupy-0.1.5.dist-info/top_level.txt,sha256=jhJY_A8Vzmqub_dVy0wPSfelUqFPjHMZOaepGVcrpxA,7
-tilupy-0.1.5.dist-info/RECORD,,
+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,,