[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