各位前輩，
最近遇到一個比較棘手的問題，我現有一個省份的鄉鎮市區的boundary的座標（經緯度），希望通過stata來計算這些鄉鎮市區的面積和中心，請問各位前輩具體的語法應該怎麼寫。我之前在網路查過一個Pro.Nick寫的Program ,試著運行但是數據沒有什麽變化，即是我的data當中并沒有產生一個新的叫area的變項。
以下是語法
*! NJC 1.0.0 15 August 2003
program polyarea, rclass
version 12
syntax varlist(min=2 max=2 numeric) [in]

tokenize `varlist'
args x y

marksample touse, novarlist
qui count if `touse' & (missing(`x') | missing(`y'))
if r(N) {
di as err "missing values in data"
exit 198
}

tempvar order area

gen long `order' = _n
su `order' if `touse', meanonly

qui gen `area' = `touse' * ///
(`x'[_n + 1] - `x') * (`y'[_n + 1] + `y') / 2
qui replace `area' = ///
(`x'[`r(min)'] - `x') * (`y'[`r(min)'] + `y') / 2 in `r(max)'
su `area', meanonly

di as res %9.0g abs(r(sum))
return scalar area = abs(r(sum))
end


謝謝各位前輩，希望斑竹不要再刪除我的帖子。謝謝。