Question about analysis of "radial intensity distribution" in single fluorescence dot?

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

Question about analysis of "radial intensity distribution" in single fluorescence dot?

Liu, Dongfang (NIH/NIAID) [F]
        Dear friends,
        I am new user of Image J.
        I am sorry to bother you.
        I am looking for your help about the analysis of "radial intensity distribution" . The radial intensity distribution I(r) of the flourescence dot was fit for each frame with a nonlinear Levenberg-Marquardt routine to the Gaussian:I(r)=Io*exp(-R^2/w^2) +BG,  where R is the distance of each pixel to the center of mass. The fitting parameters are Io, w, and BG, where Io is the peak intensity, BG is the background intensity, and w is the measure of the width (the Gauss width). It is very important for us to get the peak intensity and width. some people suggest me to use Plot Profile(Analyze -> Plot.
        But I think the "plot profile" is not good. Because the concept of "plot profile"---Displays a two-dimensional graph of the intensities of pixels along a line within the image. The x-axis represents distance along the line and the y-axis is the pixel intensity. For rectangular selections, displays a "column average plot", where the x-axis represents the horizontal distance through the selection and the y-axis the vertically averaged pixel intensity
        I was confused that how to calculate the "radial intensity distribution"  in Image Pro Plus or Image J or other software. Could you help me?  Advice would be greatly appreciated.
        I am looking forward to your reply.

Thank you very much.
Best regards

> ----------
> From: Wayne Rasband
> Reply To: List IMAGEJ
> Sent: Thursday, April 13, 2006 2:30 PM
> To: List IMAGEJ
> Subject: Re: XY Coords in a Stack and in Macro Problem with Variable File Location
>
> You have to use string concatenation (the "+" operator) to pass
> variables to commands called using the run() function. There is an
> example at
>
>      http://rsb.info.nih.gov/ij/macros/ArgumentPassingDemo.txt
>
> -wayne
>
> On Apr 13, 2006, at 6:57 AM, Brian Sellar wrote:
>
> > Hello.
> >
> > I'm using the save XY coords on a stack of images in a macro
> > and it works fine but only for a single stack. I'm trying to
> > do a batch job since I have 25000 images per test so
> > splitting the images into groups of 100 or so. Everything
> > seems to be working fine (i don't get any error messages)
> > but the XY coord command doesn't seem to like being given a
> > Variable as its file location so doesn't write any ouput
> > file on each loop of the macro. Any ideas. Am i missing
> > something obvious. Thanks in advance for any help.
> >
> > Heres my macro.........
> >
> >
> > ..............................................................
> >
> > for (n=0; n<1; n++) {
> >
> > run("Image
> > Sequence...",
> > "open=C:\\BGS\\FOOTAGE\\BGSTEST2VIDEO\\IMAGES\\NECESSARY\\Imag
> > e00000.bmp number=10 starting="+(n*100)+2+" increment=1 scale=100
> > file=");
> >
> > makeRectangle(20, 150, 680, 250);
> > run("Crop");
> >
> > for (i=1; i<=nSlices; i++) {
> >               run("Set Slice...", "slice="+i);
> >           run("8-bit");
> > run("Min...", "value=100");
> > setThreshold(240, 255);
> > run("Threshold", "thresholded remaining black");
> > }
> >
> > setThreshold(240, 255);
> > run("Threshold", "stack");
> > run("Skeletonize", "stack");
> >
> > //This works but will just overwrite NECESSARY.txt each time
> > run("Save XY Coordinates...", "process
> > save=C:\\BGS\\FOOTAGE\\BGSTEST2VIDEO\\JPEGRESULTS\\NECESSARY.txt");
> >
> > LABEL0="C:\\BGS\\FOOTAGE\\BGSTEST2VIDEO\\JPEGRESULTS\\";
> > LABEL1="Sequence"+n;
> > LABEL2=LABEL0+LABEL1;
> > LABEL3=LABEL2+".txt";
> > // Prints Something Like
> > C:\\BGS\\FOOTAGE\\BGSTEST2VIDEO\\JPEGRESULTS\\Sequence#(loop
> > number).txt
> > print(LABEL3);
> >
> > //Trying to save XY Coords on each loop of the macro.
> >
> > run("Save XY Coordinates...", "process save=LABEL3");
> >
> > }
> >
>
>