In [1]:

    

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
import welly
welly.__version__


    

    
Out[1]:





'0.3.3'

In [2]:

    

import os
env = %env


    

In [3]:

    

import csv
fname = env['HOME'] + '/Dropbox/dev/recipes/data/Nova_Scotia_wells_2015.csv'
reader = csv.DictReader(open(fname))

wells = []

for line in reader:
    params = {'header': {},
              'location': {},
             }
    
    for key, value in line.items():
        sect, item = key.split('.')
        params[sect][item] = value
        
    # Deal with surface location nonsense.
    
    # First the numerics.
    try:
        lat = float(params['location'].get('latitude'))
        lon = float(params['location'].get('longitude'))
    except:
        lat, lon = 0, 0
        
    try:
        x = float(params['location'].get('x'))
        y = float(params['location'].get('y'))
    except:
        x, y = 0, 0

    # Then the strings.
    datum = params['location'].get('datum')  # Empty string if missing
    utm = params['location'].get('datum')    # Empty string if missing
    
    # STRICT RULES
    # Either a well location has what it needs, or it doesn't. 
    
#     # If we have (lat, lon) then chuck everything else out.
#     if lat and lon:
#         del params['location']['x']
#         del params['location']['y']
#         del params['location']['utm']
#         if not datum:
#             del params['location']['latitude']
#             del params['location']['longitude']
            
#     # Otherwise, keep fully qualified (x, y)s
#     else:
#         del params['location']['latitude']
#         del params['location']['longitude']
#         if (not x and y) or (not datum) or (not utm):
#             del params['location']['x']
#             del params['location']['y']
#             del params['location']['datum']
#             del params['location']['utm']

    # Build the well and add it to the list.
    h = welly.Header(params['header'])
    l = welly.Location(params['location'])
    w = welly.Well({'header': h, 'location': l})
    wells.append(w)


    

In [4]:

    

[w.header.license for w in wells].index('P-129')


    

    
Out[4]:





10

In [5]:

    

wells[10].add_curves_from_las('../../recipes/data/las/P-129_out.LAS')


    

In [6]:

    

wells[10].location.crs_from_epsg(4269)  # NAD83


    

In [7]:

    

wells[10].location.crs


    

    
Out[7]:





CRS({'init': 'epsg:4269', 'no_defs': True})

In [8]:

    

wells[10]


    

    
Out[8]:




Kennetcook #2
td
1935.0
crs
CRS({'init': 'epsg:4269', 'no_defs': True})
basin
Windsor Basin
gl
90.3
kb
94.8
latitude
45.20951028
longitude
63.75679444
data
CALI, DEPT, DPHI_DOL, DPHI_LIM, DPHI_SAN, DRHO, DT, DTS, GR, HCAL, NPHI_DOL, NPHI_LIM, NPHI_SAN, PEF, RHOB, RLA1, RLA2, RLA3, RLA4, RLA5, RM_HRLT, RT_HRLT, RXOZ, RXO_HRLT, SP

In [9]:

    

wells_to_export = wells[:20]  # Most of these will fail.


    

In [10]:

    

my_crs = wells[10].location.crs.data  # Most of these will fail.


    

In [11]:

    

import fiona
from fiona import crs
from shapely.geometry import Point, LineString, mapping

point_schema = {'geometry': 'Point',
                'properties': {'name': 'str',
                               'uwi': 'str',
                               'td': 'int',
                               'gr': 'float'
                               }
                }

fname = 'well.shp'

with fiona.open(fname, "w",
                driver="ESRI Shapefile",
                crs=my_crs,
                schema=point_schema) as f:

    count = 0
    for w in wells_to_export:
        try:
            p = Point(w.location.latitude, w.location.longitude)
            f.write({'geometry': mapping(p),
                     'properties': {'name': w.header.name,
                                    'uwi': w.uwi,
                                    'td': w.location.td,
                                    'gr': float(w.data['GR'].mean())
                                   }
                     })
            count += 1
        except:
            pass
print('Exported {} well(s)'.format(count))


    

    




---------------------------------------------------------------------------
ModuleNotFoundError                       Traceback (most recent call last)
<ipython-input-11-a7e51a270f45> in <module>()
----> 1 import fiona
      2 from fiona import crs
      3 from shapely.geometry import Point, LineString, mapping
      4 
      5 point_schema = {'geometry': 'Point',

ModuleNotFoundError: No module named 'fiona'

In [13]:

    

import pprint
with fiona.open(fname, "r", driver="ESRI Shapefile") as c:
    for i in c:
        pprint.pprint(i)


    

    




{'geometry': {'coordinates': (45.20951028, 63.75679444), 'type': 'Point'},
 'id': '0',
 'properties': {'gr': 78.9863535887685,
                'td': 1935,
                'uwi': None,
                'name': 'Kennetcook #2'},
 'type': 'Feature'}

In [ ]: