Binary striplog

Maybe a binary striplog is a special thing? Where the main component property is True or False… then you could allow operations like dilate and erode. (Can’t really see how you could allow this on non-binary logs).

[rgmyr]

Did some more hacky stuff (to get around non-numeric components, etc.) and came up with this for the grey version (with three components mapping to log values of [1,2,3]):

[rgmyr]

This is a little hacky so I won’t PR it, but here’s an example:

from scipy.ndimage import morphology

for iv in s:
    iv.data.update(iv.primary)

def dilate(s, field, height=0.75, step=0.01):
    
    blog, basis, _ = s.to_log(step=step, field='pay', return_meta=True)
    
    blog = morphology.binary_dilation(blog, structure=np.ones(int(height/step)))
    
    return Striplog.from_log(blog, cutoff=0.5, components=comps[::-1], basis=basis)
    

def erode(s, field, height=0.75, step=0.01):
    
    blog, basis, _ = s.to_log(step=step, field='pay', return_meta=True)
    
    blog = morphology.binary_erosion(blog, structure=np.ones(int(height/step)))
    
    return Striplog.from_log(blog, cutoff=0.5, components=comps[::-1], basis=basis)


dilated = dilate(s, 'pay')
for iv in dilated:
    iv.data.update(iv.primary)

eroded = erode(s, 'pay')
for iv in eroded:
    iv.data.update(iv.primary)

opened = dilate(eroded, 'pay')

closed = erode(dilated, 'pay')

I’ve double checked that the blog array comes out right, but for some reason it seems like the components get swapped and I have to use comps[::-1] list in from_log. Don’t know why that’s happening.

In any case, I made the structuring element a bit large for emphasis, but here’s the plot: