Trying to use AB Join to compare two well log curves called Gamma Ray (GR) and pink major trends with similar magnitudes about stumpy HOT 2 CLOSED

@Philliec459
Thank you for your question, and welcome to the STUMPY community!

I love this library.

Great! Glad to hear it!!

However, it picks up similar trends but does not match the magnitudes. Is there a way to do this?

If magnitude is important in your data, then you need to set the parameter normalize to False.

w10_mp = stumpy.stump(
    T_A = df['GR'],
    m = m,
    T_B = df2['GR'],
    ignore_trivial = False,
    normalize=False,
)

The parameter normalize has the default value of True, which z-normalizes subsequences prior to computing distances. Read more in the documentation

Note that your current code uses the default value True for the normalize parameter. Therefore, after finding the closest match, you need to z-normalize each of those before plotting them against each other. In case that matters, each subsequence S can be z-normalized as follows:

import stumpy
z_normalized_S = stumpy.core.z_norm(S)

from stumpy.

Thank you Nima, I did not realize the normalize was default True, and the documentation does point that out. I am on the learning curve with this library. I see huge potential here once I figure it out. I also find that using Panel (or ipywidgets in notebook) is very useful.

This is now working much better with a much better Matrix Profile. I am getting a pretty good match everywhere, but the end of the data (bottom of the well) the GR does increase on both wells, but the program matches my hight GR interval at the bottom of Well 10 with a high Matrix Profile number to a low GR interval with a low Matrix Profile number in Well 31.

I would think it would find the similar, hight GR pattern at the bottom of Well 31 ?

Also, as a new user I have no idea of how to implement your second comment:

"Therefore, after finding the closest match, you need to z-normalize each of those before plotting them against each other. In case that matters, each subsequence S can be z-normalized as follows:"
z_normalized_S = stumpy.core.z_norm(S)

from stumpy.