[R-sig-Geo] center of mass for variable values in a grid

Barry Rowlingson b.rowlingson at lancaster.ac.uk
Sat Sep 11 10:14:08 CEST 2010


On Sat, Sep 11, 2010 at 12:49 AM, Waichler, Scott R
<Scott.Waichler at pnl.gov> wrote:
> Is there an existing function to compute the coordinates for the center of mass of a variable in a 2D or 3D grid, where dx, dy, and dz are known for all cells?  The mass in each cell varies by location.
>
> Thanks,
> Scott Waichler

 If your data is in a data frame of coordinates x, y, (and maybe z),
and variable W for all cells, and all your cells are the same size,
then its just the weighted average of the coordinate (assuming your
dataframe is the complete grid). So:

 xm = sum(d$x*d$W)/sum(d$W)
 ym = sum(d$y*d$W)/sum(d$W)

Or have I missed something?

Barry



More information about the R-sig-Geo mailing list