Deleting Specific Points in a Wave

Hi,

Couple of issues with you approach is delete is slow and more problematic is that in each loop that deletes a point and changes the length of the wave.  So when you don't actually test all the points and you since it shifts and you could get an out of range error towards the end of wave as the number of points is reduced.  One way around this potential issue that I have used in the past is to work backwards from the end of the wave to point 0.  That way you you won't get an out of range error.

Another alternative is to use extract.

Extract source, dest, Topwave!=2 && NumType(BotWave)!=2

Then kill the source wave.

Andy

This function extracts the points and places them all together in a separate wave, I need to maintain them in the same index order.  I haven't figured out how to kill the source wave.  Does that delete the extracted points?

Hi,

Give this a try.

Extract TopWave, IsolatedTop, Topwave!=2 && NumType(BotWave)!=2

This will create a new wave IsolatedTop with the points removed and in the same order as in top wave.

Andy

That's what I tried and the new wave doesn't match the index of the source wave.

I hope I'm being clear about what I need.

I essentially want to delete the topwave points when both top and bottom waves have real numbers.  So referring to the image I attached, the isolatedtop wave should have blanks, and then have numbers after the 3rd index.

Deleting a point is not the same as setting the value of the point to blank/NaN.

Which do you want?

Sorry for not being clear.  I want to set the value to a blank/NaN.

When I try the following code, it assigns all of my other values as 0.

How about:

topwave = Numtype(BotWave)==2? Topwave : NaN

If topwave is a NaN then it stays a NaN

if Botwave is a NaN then the topwave remains the same whether it is a real value or a NaN

if Botwave is a real number then it changes topWave at that point to NaN

The result is a topwave with real values and NaNs.

If you want to subsequently remove the NaNs from topwave you can do a wavetransform zapnans.

Andy

j s s wrote:

When I try the following code, it assigns all of my other values as 0.

 

Function IsolateTop(TopWave, BotWave)
    Wave TopWave,BotWave
    Variable i, N = numpnts(Topwave)
    Make/O/N=(N) IsolatedTop
   
    for (i=0; i<N; i+=1)
        if ((NumType(TopWave[i]) == 0) && (NumType(BotWave[i]) == 0))
        IsolatedTop[i] = NaN
        endif
    endfor

end

That is because IsolatedTop starts out as all 0s and is assigned nothing else inside your for loop when the conditions do not match. Why did you use Make instead of Duplicate as before?