I made some tests and the threshold can not be applied to all scales as a constant.

We should use a relative thresholded approach here. Let's stick to the following version `(Axx + Ayy)*(Axx + Ayy) >= rpc_threshold*(Axx*Ayy - Axy*Axy)` and rename the parameter.