Find Peaks, Max and Min

Share icon

Hi there, I try to automatize the calculation of the weight difference of a silo for certain time periods. The signal of the silo weight is shown in the attached image. There are always time periods of charging (weight increase: almost a step function) and discharging the silo (slow weight decrease). I'm interested in how much of the charged matter was used between two charges of the silo. Until now I made this manually. Then I tried the PeakFind. Min and Max etc. functions in IGOR, but didn't get a good result.
Does somebody have an idea?
Thank you...K.
An approach I suggest is to create a wave of "charging" times or points by detecting the large step-like increases. You may want to set a minimum increase threshold for this. The 'Differentiate' operation might be an easy way to assist this. Then to find the used matter, simply integrate ('Integrate' operation) the original data wave between adjacent charging points, taking into account the time scaling.

July 29, 2010 at 03:18 am - Permalink
Finding the charging points by the differentiate operation is a good approach.
Maybe the following function does what you want:
function detUsage(silo)
	wave silo
	make/n=0/o/free usage_tmp,usageX_tmp
	variable i=0,mmin=pnt2x(silo,0),mmax=pnt2x(silo,numpnts(silo)-1),V_max2
	Differentiate silo/D=silo_diff
	findpeak/q/r=(mmin,mmax)/m=(0.1*silo_diff(mmin)) silo_diff
	V_max2=wavemax(silo,mmin,V_LeadingEdgeLoc)
	do
		if(V_flag!=0)
			killwaves/z silo_diff
			extract/o usageX_tmp,usageX,usage_tmp>0
			extract/o usage_tmp,usage,usage_tmp>0
			return 0
		endif
		wavestats/q/r=(mmin,V_LeadingEdgeLoc) silo
		usage_tmp[i]={V_max2-V_min}
		usageX_tmp[i]={V_minloc}
		wavestats/q/r=(V_LeadingEdgeLoc,V_TrailingEdgeLoc) silo
		V_max2=V_max
		mmin=V_TrailingEdgeLoc
		findpeak/q/r=(mmin,mmax)/m=(0.1*silo_diff(V_maxloc)) silo_diff
		i+=1
	while(1)
end

July 29, 2010 at 07:08 am - Permalink