Probably you would want to be consistent with what we do in the m/z dimension, which is to use the weighted mean of the signal in the extraction range. This often ends up being just a single value in centroided data, but it also works well in messier profile mode data, like profile-mode TOF m/z. I think it likely works better when resolution decreases and an apex value is not clearly the center of the distribution. This method has the ability to give a value that is between two buckets with relatively similar signal.