Intensity along a line with a "pixel-relative" coordiantes

[itayfay]

Hi,

I am looking for a way to find the intensity along a line (yellow line circled in blue), but instead of the (X,Y) coordinates that correspond to the distance from the line beginning (blue arrows), I wish to obtain the coordinates that correspond to the pixels of the entire image with a given origin. For example, if my image is a 4X4 ( raws X columns ) matrix, and I draw a line from (3,1) to (3,3), I wish to get the following Table:

X; Y: intensity
1; 3; somthing
2; 3; somthing
3; 3; somthing

I noticed that there is a similar post:

Function to get pixel x,y coordinates for a line ROI

Sorry, I have looked around a lot for this I am unable to find the solution. I am stuck on this part: Could you please help me write a function which will output x,y coordinates in raw pixels of a line drawn using line tool (including for oblique lines). Many Thanks Harry

However, if I use the given script, I need to insert the coordinates manually every single time (sorry if I did not understand the code correctly, I am not familiar with the syntax)

[LeChat (Adrien)]

You can trace the line, then to the measurements (ctrl+m) and then in the Results window, you go to the menu Results>Set Measurements and check the boxes “mean gray value” (from what I understand, this is the value of the intensity you are interested in), and as well as “Bounding rectangle” (this will give you BX, BY, Width and Height) from which you basically have the coordinates of the line you plotted.

Hope this helps.

[Mathew]

Hi
@itayfay
To complete Adrien’s [ @LeChat ] info, I did a little macro-test.
It gives you the results in two tables. Personally I will go through a spreadsheet, copying / pasting. It’s easier to deal with.
Can you tell me about your results, if you don’t mind.
Greetings

// Step 1
newImage("test", "8-bit ramp", 512, 512, 1);
makeLine(34, 200, 434, 200);
run("Interpolate", "interval=1");
//run("Interpolate", "interval=1 adjust");
getSelectionCoordinates(X,Y);
print("Number of points=", X.length);
//print( lengthOf(xCoordinates) );
Array.show("Selection Coordinates",X,Y);
print("Table title=",Table.title);

// Step2
p = getProfile();
Plot.create("gray values", "xInterval", "grayvalue", p);
Plot.show;
Plot.getValues( xInterval, grayvalue );
Array.show("Value",xInterval,grayvalue);
print("Table title=",Table.title);
print("Table size=",Table.size);

[eljonco]

Just add this code, to the @Mathew’s code so everything is in 1 table.

for(pt=0;pt<grayvalue.length;pt++){
	setResult("X",pt,X[pt]);
	setResult("Y",pt,Y[pt]);
	setResult("Greyvalue",pt,grayvalue[pt]);
}

and maybe make the second line a bit more exciting:
makeLine(34, 300, 434, 200);

[Mathew]

@eljonco
First of all thank you, because although I knew that I had to go through setResult I did not see how.
A small step forward for me, thanks to you.
Greetings