Hi David, I need more time to evaluate your latest code. So far I could take apart and understand your previous code almost completely. I have not a final conclusion yet. I can say that the code seems to work as it should and that their algorithm seems ok. However, I'm not too sure about the result, I need to continue investigating. Give me some more days......

This experiment is not confidential, You can ask what you want. This is a situation in which light from a laser is incident on the input fiber with an angle to the surface. The result is that the output light is a ring-shaped projection. The image that you see, is obtained with the help of a beam spliter reflecting the ring projected onto a target at about 1 meter from the output of the fiber. The analysis of the results can give us informations about the polishing angle of the fiber in test.

Thanks a lot.....

Best regards,

Antonio

Hi Antonio,

The code I posted yesterday uses an algorithm like you first described -- it determines a point "centered" with respect to a collection of maxima along radials, and it iterates. It may work.

Regarding the addition of a list of images, that is what it does. You give it a list of file names. It imports those files as images, adds them to make one image, and analyzes it. This can be immediately extended as you describe above. If you partition the file names into a list of lists, where each low level list is a list of files you want added and analyzed, then you will get a list of results. Look at the function named Map. Map[f,{a,b,c}] outputs { f[a], f[b], f[c] }.

So if you have a large list of files which you want to be analyzed in groups, you partition them:

partitionedFiles = { {file1, file2, file3, file4, file5}, {file6, files6, file8, file9, file10} }

And then

Map[ AnalyzeImageGroup, partitionedFiles]  or the shortcut AnalyzeImageGroup/@partitionedFiles

will produce an analysis for each of the sublists as a group. You will get a list of results, one result for each group. And an annotated graph will be saved as a png for each of the groups, using the name of the first file in the group for the name of the png.

Kind regards,

David

Please attach a sample image. It will help in understanding.

Best regards, David

Hi Antonio,

What I understand is that you wish to calculate the centroid of the image, and then perhaps the mean radial distance from the centroid to the annular ring. The attached notebook uses the usual normalized weighted sum to get the centroid, and then reprocesses the data into a purely radial dependency. From that it extracts the radius of the annulus. Note that I deleted the output cells to shrink the file for posting.

Some of the code may be a little difficult at first if you're just starting with Mathematica. But time spent learning Mathematica will be repaid many times. Here are two nice resources: An Elementary Introduction and Virtual Book

If I am not understanding this correctly, please let me know.

Kind regards,

David

Attachments:

David, I'll take some time to fully understand and test your code by myself. I can see clearly que your approach is different of mine.This is something really very liberating.

As I said before, the code needs to repeat the calculation process for 20 pictures or more, always adding each group of five images in sequence. Can you do this repetition?

My big difficulty with Mathematica is to imagine repetitive processes in long codes without use of looping.

Thanks

Cesar

Hello Cesar,

I became fascinated with this problem. Attached you will find a notebook. I deleted the output cells to get it below the 20MB limit for attached files. I also attach a PDF of it, so you can see what it should look like when executed.

The notebook exemplifies several aspects I think are important in using Mathematica.

1. As I discuss above, its goal is to develop a single complicated function which can be applied to an argument, or even mapped onto a list of arguments. (In this case, the argument is the list of files to be added together and analyzed.) But the development process is to first perform the work the function will perform, but in a sequence of individual evaluations. Then this code can be copied to build up a function from proven code.
2. In the code, you will see many instances in which Fortran or C would use a loop to make a calculation on several data items, but Map is used instead. In most of the occasions in which we use a loop in a procedural language, the Nth iteration does not depend on (N-1). For these cases map is the Mathematica way.
3. In many of these cases, what is mapped onto a list is an "anonymous function." This is code like #1^2+#2&. It is worth understanding these. They are very useful.
4. However -- the method of iterating on estimates of the center is in fact a true iteration. This is done using the function FixedPoint, which repeatedly composes a function until the Nth result is the same as the (N-1) result. (Where we can define the meaning of same.) This is an example of another Mathematica truth: There is often a built-in function which will do the heavy lifting.
5. The function finally developed will stand on its own. The first section was only used to develop the code in pieces. It accepts a list of file names , which should be the group of fits files. It loads them, adds them into a single image, and processes it to get a center, a radius to the annulus, and outputs an analysis. In the process, it exports a plot that can be checked for reasonableness. Like any function, you could take your all file names and partition them into a list of lists, each of which contains the names of files which should be analyzed as a set. Then you could Map this function onto the list of lists, to get a list of results, one result for each file group.

This is a lot more that I started out to do, but I really became fascinated with the problem!

Kind regards,

David

Attachments:

Dear Frank, Thanks for your suggestion. I have a tendency to use looping because my only computer training was over 30 years with FORTRAN and looping was a very reasonable option. I will work hard with your suggestion ......... Thanks

Initially I would like to apologize for this question involving a so large code. I am no expert in Mathematica and have been trying to learn from limited resources. I have been involved with this code longer than it should and that makes me really discouraged.From the moment I posted the code here I get some progress to find the center of a single image. This was achieved by using a looping type (Goto). My difficulty now is define another loop to redo the operations for the other images imported on begin of the code. I would really appreciate some help to finish this code to obtain the final value of the center of mass coordinates. Thanks for any help...

Attachments: