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

Waichler, Scott R Scott.Waichler at pnl.gov
Sat Sep 11 17:26:06 CEST 2010


> > 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.
> 
>  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)

You're right, it's simple for a grid where dx and dy are constant.  Thanks.  But if they vary?

Scott Waichler



More information about the R-sig-Geo mailing list