Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #по каждой каротажке получаем интервалы где есть измерение (начало исследования, и длина непрерывной записи)
- def get_measured_intervals(well_logs, mnimonic, depth_key):
- ints = []
- pair = []
- nan = 'nan'#'-9999.0'#'nan'
- #добавил 1
- if(len(well_logs[mnimonic].dropna()) == len(well_logs[mnimonic])):
- a = np.array(well_logs[depth_key])[0]
- pair.append(a)
- b = np.array(well_logs[depth_key][-1:], dtype = np.float64)[0]
- pair.append(b-a)
- ints.append(tuple(pair))
- pair = []
- if (str(np.array(well_logs[mnimonic])[0]) == nan):
- for i in range(len(well_logs)-1):
- if(str(np.array(well_logs[mnimonic])[i]) == nan and str(np.array(well_logs[mnimonic])[i+1]) != nan):#i+1 начало интервала
- a = np.array(well_logs[depth_key])[i+1]
- pair.append(a)
- if(str(np.array(well_logs[mnimonic])[i]) != nan and str(np.array(well_logs[mnimonic])[i+1]) == nan):#i конец интервала
- b = np.array(well_logs[depth_key])[i]
- pair.append(b-a)
- ints.append(tuple(pair))
- pair = []
- #добавил
- if(str(np.array(well_logs[mnimonic])[i]) != nan and str(np.array(well_logs[mnimonic])[i+1]) != nan and i == len(well_logs)-2):#i конец интервала
- b = np.array(well_logs[depth_key])[i]
- pair.append(b-a)
- ints.append(tuple(pair))
- pair = []
- if(str(np.array(well_logs[mnimonic])[0]) != nan):
- a = np.array(well_logs[depth_key])[0]
- pair.append(a)
- for i in range(1,len(well_logs)-1):
- if(str(np.array(well_logs[mnimonic])[i]) != nan and str(np.array(well_logs[mnimonic])[i+1]) == nan):#i конец интервала
- b = np.array(well_logs[depth_key])[i]
- pair.append(b-a)
- ints.append(tuple(pair))
- pair = []
- if(str(np.array(well_logs[mnimonic])[i]) != nan and str(np.array(well_logs[mnimonic])[i+1]) != nan and i == len(well_logs)-2):#i конец интервала
- b = np.array(well_logs[depth_key])[i]
- pair.append(b-a)
- ints.append(tuple(pair))
- pair = []
- if(str(np.array(well_logs[mnimonic])[i]) == nan and str(np.array(well_logs[mnimonic])[i+1]) != nan):
- a = np.array(well_logs[depth_key])[i+1]
- pair.append(a)
- return(ints)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement