Additional comments
add cr reject settings to m2fs_merge so I don't keep editing merge.py
add documentation
Use sexconvert
scattered light work. Works but doesn't compute estimation regions on its own.
Tidy __init__ still not sure I'm doing this right....
Create order location routine.

Smooth image with guassian kernel and then pick out peaks