請問在Stata中如何根據座標計算任意多邊形的面積和中心

****************************************************************************
clear
import excel "D:\temp\data_sample.xlsx", sheet("Sheet1") firstrow clear

sum
drop if x==.
polyarea  x y
return list
di r(area)

******************************************************************************








运行结果如下




. clear

. import excel "D:\temp\data_sample.xlsx", sheet("Sheet1") firstrow clear

.
. sum

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
         _ID |      9487    12.67208    5.652406          1         20
           x |      9467     120.368    .1514226   120.1335    120.656
           y |      9467    23.72252    .0543959    23.5866   23.84686
          id |      9487        4744    2738.805          1       9487

. drop if x==.
(20 observations deleted)

. polyarea  x y
.1077965

. return list

scalars:
               r(area) =  .107796460899408

. di r(area)
.10779646