Image Moment - Optimizing Code for Speed

[SectorEffector]

Hello 

I'm want to find the moment of inertia of an 2d-array. The array is converted from an image using "IMAQ ImageToArray".

The algorithm I'm using is discribed here:

Wikipedia - Image Moments

I need to calculate this formula with different values for i and j:

i.j = 0.0 - 0.1 - 1.0 - 1.1 - 0.2 - 2.0

I programmed the code shown above/attached VI, but I need the optimize it for speed.

The 2d array can be any size with a maximum of 2048 x 2048 with values varying between 0 and 4095.

My question:

How can I make this code faster?  

Thank you and kudos will be given! 

[TiTou]

I think if you want good performance you should try to avoid the "Image to Array" and see if your processing doesn't exist in the IMAQ functions, open the Vission concepts manual in C:\Program Files\National Instruments\Vision\Documentation and search for moment the fews hits might be of interest to you.

Hope this helps

[SectorEffector]

Hello Titou,

Thank you for your reply.

The Vision palette only has mathematical functions for binary images (0 to 1 instead of 0 to 4095). 

The "IMAQ Particle Analysis" vi is promising and has functions for finding the orientation and moment of inertia,

but only for binary particles.

The only Vision VI that makes use of the weight of the pixels is the "IMAQ Centroid" VI.

[SectorEffector]

Is there someone who can help me optimize this code?

It is to slow for our application at the moment and I would be vey thankful for every optimize!

[Darin.K]

Looking at the code I would try the following:

combine all of the operations that involve the transposed arrays into a single for loop.  same for the original arrays.  (ie. combine M10 and M20 loops as well as M01 and M02 loops)

parallelize the for loops (only the outer one for M11).

[GregSands · Accepted solution]

You can speed it up about 60% just by combining some of your loops -- the calculation of M01 and M02, and M10 and M20, both share some computation.