How to perform clutter removal on this signal?

Hello everyone,

I am trying to perform clutter removal on this signal attached below.

I want to set the signals to zero except a higher peak region in the middle like shown below.

Any suggestions would be appreciated. The higher peak in the middle have values both in positive as well as negative directions. Signals beyond that regions are set to zero. I am attaching the excel data sheet for this signal.

Attaching two more pictures (from literature, not mine). I am trying to doing something like this.

[KALYAN ACHARJYA]

Can you clarify the codition again?

Those y data forced to be zero ....if....

[Adam Danz]

What defines the region of interest?

• A range of x values?
• The location of the tallest peak for each inidividual line?

What should happen when an individual line doesn't have a significant peak in that area?

• Should the entire line be set to 0?
• Should that section be preserved while the rest of the line is set to 0? What defines that section?

Lots of open questions that must be defined before thinking about the code.

[Deepshikha Bhargava]

@KALYAN ACHARJYA Thanks for showing interest.

I want those Y data to be zero everywhere except that higher peak region between 0.005-0.01 xaxis.

Means one signal gets zero values everywhere expcept when it is in that higher peak region.

Am i clear enough? let me know.

[Deepshikha Bhargava]

@Adam Danz Thankyou for the reply. I am trying here to answer your questions. Let me know If I havenot answered them correctly.

What defines the region of interest?

1. No, xaxis data has nothing to do with the higher peak region. I need to work on yaxis data.
2. The line does not necessarily have higher peak in that region. But whatever peak the line has in that region should be preserved and rest set to zero (even the peak is higher outside that region).

What should happen when an individual line doesn't have a significant peak in that area?

1. Whatever peak the line has in that region should be preserved. Overall I want that region of the graph to have higher peak only rest set to zero.
2. yes, that section should be preserved.

[Deepshikha Bhargava]

@KALYAN ACHARJYA @Adam Danz Please see the question again. I added two pictures from literature. This is what I am trying to do.

Many Thanks.

[Adam Danz]

If the x-axis has nothing to do with it then what do you mean by "in that region"? What defines the region?

For example, maybe the region is defined by x=.1 to x=.15 but you're saying the x-axis has nothing to do with it so I'm confused.

What if there are multiple peaks within the same line within the region?

[Deepshikha Bhargava]

@Adam Danz

You are right. xaxis would play a role in preserving that higher peak region. The signal also has some higher peak in the starting (x=0) that we want to set to zero. But want to preserve the peak between x= 0.005 to 0.01 (roughly).

The graphs from literature also has some higher peak in the starting (x=0). They set them to zero and kept the peaks between x= 0.1 to 0.15 (roughly).

[Adam Danz]

Looks like you have some good advice below in the answers section.

Answer 1

[Image Analyst]

Try this:

clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 22;
fprintf('Beginning to run %s.m ...\n', mfilename);

% Read in data.
x=readmatrix('xaxis_data.xlsx')';
y=readmatrix('yaxis_data.xlsx');
[rows, columns] = size(y)
% Define limits of the range we want to view.
lowerLimit = 0.005;
highLimit = 0.010;
% Find the indexes where x is outside the limits because we want to erase there.
badIndexes = x < lowerLimit | x > highLimit;
% Make a copy of the array that will be zeroed out outside of the good range.
cleanedY = y;
cleanedY(badIndexes, :) = 0; % Zero outside the good range.
% Plot the cleaned data.
for k = 1 : size(y, 2)
    plot(x,cleanedY(:,k), '-');
    hold on;
end 
grid on;
xlim([0, 0.015]);
xlabel('x', 'FontSize', 20);
ylabel('y', 'FontSize', 20);
% Put lines on either side of the good range.
darkGreen = [0, 0.7, 0];
xline(lowerLimit, 'Color', darkGreen, 'LineWidth', 3);
xline(highLimit, 'Color', darkGreen, 'LineWidth', 3);
% NOTE: there are no data points exactly at the lower and upper limits.
% so if you want it to be zero EXACTLY up to that point, you'll have to insert those points into the data.
fprintf('Beginning to run %s.m ...\n', mfilename);

NOTE: there are no data points exactly at the lower and upper limits, so if you want it to be zero EXACTLY up to that point, you'll have to insert those points into the data. If you can't figure that out, let me know and I'll do it.

[Image Analyst]

Alright, you just add this:

% Insert points that are zero exactly at the limits.
x = [x; lowerLimit; highLimit];
y(end+1:end+2, :) = 0;
[x, sortOrder] = sort(x, 'ascend');
y = y(sortOrder, :);

The full demo is attached.

Answer 2

[KALYAN ACHARJYA]

@Adam Danz clearly stated, If you know the x range, which is to be set to zero, then you can directly locate the indices of the corresponding x value range. Which I have shown in my code below. If you are looking for both the peaks and the x data range, then you also have to look for the y data. But there should be a clear definitive understanding of which values need to be set to zero.

x=readmatrix('xaxis_data.xlsx');
y=readmatrix('yaxis_data.xlsx');
temp_y=y;
temp_y(x>0.025 & x<0.18,:)=0;
for i=1:size(y,2)
    plot(x,temp_y(:,i));
    hold on;
end 
grid on;

I hope, you can modify it and get your desired result, hence I posted this as an answer.