I want to turn this code into led light instead of red ones

Hi, I want to turn this code into led light instead of red ones. How may I do that; which part should I change?

Here is my code;

data =imread('IMG_4236.JPG');


% Now to track red objects in real time
% we have to subtract the red component 
% from the grayscale image to extract the red components in the image.
diff_im = imsubtract(data(:,:,1), rgb2gray(data));
%Use a median filter to filter out noise
diff_im = medfilt2(diff_im, [3 3]);
% Convert the resulting grayscale image into a binary image.
diff_im = im2bw(diff_im,0.18);

% Remove all those pixels less than 300px
diff_im = bwareaopen(diff_im,300);

% Label all the connected components in the image.
bw = bwlabel(diff_im, 8);

% Here we do the image blob analysis.
% We get a set of properties for each labeled region.
stats = regionprops(bw, 'BoundingBox', 'Centroid');

% Display the image
imshow(data)

hold on

%This is a loop to bound the red objects in a rectangular box.
for object = 1:length(stats)
    bb = stats(object).BoundingBox;
    bc = stats(object).Centroid;
    rectangle('Position',bb,'EdgeColor','r','LineWidth',2)
    plot(bc(1),bc(2), '-m+')
    a=text(bc(1)+15,bc(2), strcat('X: ', num2str(round(bc(1))), '    Y: ', num2str(round(bc(2)))));
    set(a, 'FontName', 'Arial', 'FontWeight', 'bold', 'FontSize', 12, 'Color', 'yellow');
end

and here is the image

Answer 1

[Image Analyst]

I don't know what this means: "turn this code into led light instead of red ones".

You can't turn code into a light. The light you have in your scene is white. If you want a red LED instead of a white one, then put a red LED into your circuit. It has nothing to do with code. It has to do with electrical components and perhaps a soldering iron or a breadboard or other components to build circuits with.

The bright light in the image is not red. Plus your algorithm is just overly simplistic (even if the light were red) and basically no good. If you want to track a colored object in an image, see my attached demo where I track a green Sharpie marker in a video.

[Batuhan Istanbullu]

What I mean by turning the code into white led light, is my ignorance sorry for that!, All I wanted to say that; I want to capture only that led light color and disable all other parts.

[Image Analyst]

Is it always in the same place in the image? Or is it always the brightest thing in the image? Or is it always brighter than some certain intensity level? Is it always of some certain size? Or color? How can you describe the LED to make its attributes as unique as possible such that they could not describe any other object in the image? I suggested some possibilities but I'm wondering what you would use to describe it.

[Batuhan Istanbullu]

It is not always in the same place, its location may vary cause after I found its HSV value or any other type to find that led I want to add several images in this code so that I want to see its position while in changing with time I suppose. Size is same. Not brighter same level. and same color yes.

[Batuhan Istanbullu]

channel1Min = 0.188;

channel1Max = 0.145;

% Define thresholds for channel 2 based on histogram settings

channel2Min = 0.000;

channel2Max = 0.142;

% Define thresholds for channel 3 based on histogram settings

channel3Min = 0.961;

channel3Max = 1.000;

These are the HSV values of that led.

[Image Analyst]

Can we just find any blob brighter than, say, 250 gray levels.

rgbImage = imread(fileName);
grayImage = rgb2rgay(rgbImage);
mask = grayImage > 250;
% Take the largest one only.
mask = bwareafilt(mask, 1);
% Measure centroid and bounding box.
props = regionprops(mask, grayImage, 'Centroid', 'BoundingBox');
xCentroid = props.Centroid(1)
yCentroid = props.Centroid(2)

It could be made more sophisticated to handle different images, like images with colored lights, different size lights (some that you want to keep and some that you want to exclude). You'd have to try it with some of those images and improve the algorithm if it's a problem, but it should work for the one image you posted.

[Batuhan Istanbullu]

but when I want to find another location of the led, for instance like this; your suggestion gave me the same result of the previous one?

[Batuhan Istanbullu]

I hope that you understand what I meant, I want to track that led just like your video about the greensharpie. But instead of video which I do have the video version. I wanted to try out with the several images.

[Batuhan Istanbullu]

I mean I only want to do this but with another video which I have HSV values like this;

% % Define thresholds for channel 1 based on histogram settings

% channel1Min = 0.188;

% channel1Max = 0.145;

%

% % Define thresholds for channel 2 based on histogram settings

% channel2Min = 0.000;

% channel2Max = 0.142;

%

% % Define thresholds for channel 3 based on histogram settings

% channel3Min = 0.961;

% channel3Max = 1.000;

[Image Analyst]

No, it will not give the same location for a different image when the LED is in a different location. Let me do a full demo for the 2 or 3 images you've posted and I'll show you. Give me an hour or so.

[Batuhan Istanbullu]

Okay sir, I will wait. In the meantime I will show you my Hsv values and the video of that object changing its location. By the way, thank you for effort!. It means a lot to me to my school project.

[Image Analyst]

Here is the code. It works for both your images and does not give the same location because the LED is in different locations in each image.

clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 16;
%===============================================================================
% Get the name of the image the user wants to use.
baseFileName = 'IMG_4236.JPG';
baseFileName = 'distance.JPEG';
folder = pwd;
fullFileName = fullfile(folder, baseFileName);
% Check if file exists.
if ~exist(fullFileName, 'file')
    % The file doesn't exist -- didn't find it there in that folder.
    % Check the entire search path (other folders) for the file by stripping off the folder.
    fullFileNameOnSearchPath = baseFileName; % No path this time.
    if ~exist(fullFileNameOnSearchPath, 'file')
        % Still didn't find it.  Alert user.
        errorMessage = sprintf('Error: %s does not exist in the search path folders.', fullFileName);
        uiwait(warndlg(errorMessage));
        return;
    end
end
%=======================================================================================
% Read in image.
rgbImage = imread(fullFileName);
% Get the dimensions of the image.
[rows, columns, numberOfColorChannels] = size(rgbImage)
% Display image.
subplot(2, 2, 1);
imshow(rgbImage, []);
impixelinfo;
axis on;
caption = sprintf('Original RGB Image\n%s', baseFileName);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.
% Set up figure properties:
% Enlarge figure to full screen.
g = gcf;
g.WindowState = "maximized";
% Get rid of tool bar and pulldown menus that are along top of figure.
% set(gcf, 'Toolbar', 'none', 'Menu', 'none');
% Give a name to the title bar.
set(gcf, 'Name', 'Demo by ImageAnalyst', 'NumberTitle', 'Off')
drawnow;
%=========================================================================================================================
% Threshold the image to create a mask of the objects
% [mask,maskedRGBImage] = createMask(rgbImage);
% Get bright regions.  Gray scale is fine because the LED in the image is not colored -- it's white.
mask = rgb2gray(rgbImage) > 250;
% Get rid of particles less than 1000 pixels.
% props = regionprops(mask, 'Area'); % Measure areas of our initial mask to determine areas of LED and possible noise.
% allAreas = [props.Area]
mask = bwareaopen(mask, 1000);
% Fill any holes in the blobs.
mask = imfill(mask, 'holes');
% Display the mask image.
subplot(2, 2, 2);
imshow(mask, []);
impixelinfo;
axis('on', 'image');
title('Mask Image', 'FontSize', fontSize);
drawnow;
%=========================================================================================================================
% Get location of the LED.
props = regionprops(mask, 'Centroid');
allCentroids = vertcat(props.Centroid) % works for multiple LEDs
% Display the centroid over the image as cross hairs
subplot(2, 2, 3);
imshow(rgbImage, []);
impixelinfo;
hold on;
markerSize = 50;
for k = 1 : length(props)
    plot(allCentroids(k, 1), allCentroids(k, 2), 'r+', 'MarkerSize', markerSize, 'LineWidth', 2);
    caption = sprintf('%d LEDs found.  This one at (x,y) = (%.1f, %.1f)', ...
        numel(allAreas), allCentroids(k, 1), allCentroids(k, 2))
    title(caption, 'FontSize', fontSize);
end
axis('on', 'image');
drawnow;

See how the crosshair in the lower left image is in a different location for each image?

I did not use the Color Thresholder code because there was no need. Your LED light is bright white, not a pure color so we didn't need to do color segmentation.

[Batuhan Istanbullu]

Sir, just thank you for all your effort! I really appreciated it. I do have one last question; I will send it by the winrar file.

[Batuhan Istanbullu]

why can't I read the HSV values that I determined in here like you?, I mean I determined the HSV values of the led with color threshholder app but I dont know why

[Batuhan Istanbullu]

and when I do want to run it I encounter with this problem

[Batuhan Istanbullu]

I do not know where I'm doing wrong when I'm trying to implement your code to my values and my video

Answer 2

[Matt J]

Is this what you want?

IMG=load('Image').A;

c=regionprops(all(IMG>250,3),'Centroid').Centroid; %LED centroid/position

imshow(IMG); hold on
plot(c(1),c(2),'rx'); hold off

[Batuhan Istanbullu]

Yes sir, but I want to know its position like the rest of them and only want to know leds position and show it like others.

[Batuhan Istanbullu]

Hopefully, I could of explained

[Matt J]

As indicated in the comments, the code does give you the LED position in the variable 'c'.

[Batuhan Istanbullu]

Sir, you have no idea how much you helped me about my project! I reall appreciated what you've done for me. Thank you for effort and kindness!

[Matt J]

You're welcome, but please Accept-click the answer as well as answers to other questions you've posted which are now resolved.

[Batuhan Istanbullu]

Yes, I will do that immediately.