Calculating center of mass of a array IDL

最新推荐文章于 2023-09-12 22:56:59 发布

juliosun

最新推荐文章于 2023-09-12 22:56:59 发布

阅读量356

点赞数

分类专栏： IDL

IDL 专栏收录该内容

51 篇文章 0 订阅

订阅专栏

   FUNCTION Centroid, array
   s = Size(array, /Dimensions)
   totalMass = Total(array)
   xcm = Total( Total(array, 2) * Indgen(s[0]) ) / totalMass
   ycm = Total( Total(array, 1) * Indgen(s[1]) ) / totalMass                                                                        
   RETURN, [xcm, ycm]
   END

3D

Reasoning by analogy to the 2D case, this should work, I think:

   xcm = Total( Total(Total(array,3),2) * Indgen(s[0])) / totalMass
   ycm = Total( Total(Total(array,3),1) * Indgen(s[1])) / totalMass
   zcm = Total( Total(Total(array,2),1) * Indgen(s[2])) / totalMass