def runtest(test): p = subprocess.Popen(test, shell=True) retval = p.wait() return retval
def checkmsh(meshname): p = subprocess.Popen('[ -f grids/' + str(meshname) + '.msh ]', shell=True) if p.wait() == 0: p = subprocess.Popen('[ grids/' + str(meshname) + '.geo -ot grids/' + str(meshname) + '.msh ]', shell=True) return p.wait() else: return 1
def domesh(meshname): p = subprocess.Popen('gmsh -2 $(pwd)/grids/' + str(meshname) + '.geo ]', shell=True) retval = p.wait() return retval
from skimage.color import rgb2gray import numpy as np import cv2 import matplotlib.pyplot as plt from scipy import ndimage image = plt.imread('../images/subabul.jpg') image.shape plt.imshow(image) plt.show(image) plt.wait(0) gray = rgb2gray(image) plt.imshow(gray, cmap='gray') gray = rgb2gray(image) gray_r = gray.reshape(gray.shape[0] * gray.shape[1]) for i in range(gray_r.shape[0]): if gray_r[i] > gray_r.mean(): gray_r[i] = 3 elif gray_r[i] > 0.5: gray_r[i] = 2 elif gray_r[i] > 0.25: gray_r[i] = 1 else: gray_r[i] = 0 gray = gray_r.reshape(gray.shape[0], gray.shape[1]) plt.imshow(gray, cmap='gray')