/*accept the arguments
&s oddem05 = [response 'enter your dem' oddem05]
&s odlu05 = [response 'enter your lu' odlu05]
&s fc = [response '田間持水量狀態下的土壤含水量' 0.25]
&s sat = [response '飽和狀態下的土壤含水量' 0.25]
&s sw = [response '降雨前土壤含水量' 0.2]
&s rday = [response '日降雨量(mm)' 200]
/*define the internal var
&s slp = slp
&s cn2 = cn2
&s cn3 = cn3
&s cn2s = cn2s
&s cn3s = cn3s
&s cn1s = cn1s
&s **ax = **ax
&s **in = **in
&s w2 = w2
&s w1 = w1
&s s = s
/*define the result val
&s qsurf = qsurf
/* the process of the computing
%slp% = slope(oddem05, percentrise)
items %odlu05%.vat
list %odlu05%.vat
if(%odlu05% == 1)
%cn2% = 83
else if(%odlu05% == 2)
%cn2% = 61
else if(%odlu05% == 3)
%cn2% = 48
else if(%odlu05% == 4)
%cn2% = 55
else
%cn2% = -9999
endif
%cn3% = %cn2% * exp(0.00673 * (100 - %cn2%))
%cn2s% = (%cn3% - %cn2%) / 3 * (1 - 2 * exp(- 13.86 * %slp%))
%cn3s% = %cn2s% * exp(0.00673 * (100 - %cn2s%))
%cn1s% = %cn2s% - 20 * (100 - %cn2s%) / (100 - %cn2s% + exp(2.533 - 0.0636 * (100 - %
cn2s%)))
%**ax% = 25.4 * (1000 / %cn1s% - 10)
%**in% = 25.4 * (1000 / %cn3s% - 10)
%w2% = (ln(%fc% / (1 - %**ax% / %**in%)) - ln(%sat% / (1 - 2.54 / %**ax%) - %sat%)) / (%sat%
- %fc%)
%w1% = ln(%fc% / (1 - %**ax% / %**in%) - %fc%) + %w2% * %fc%
%s% = %**ax% * (1 - %sw%) / (%sw% + exp(%w1% - %w2% * %sw%))
%qsurf% = (%rday% - 0.2 * %s%) / (%rday% + 0.8 * %s%)
/* delete all internal layers
kill %slp% all
kill %cn2% all
kill %cn3% all
kill %cn2s% all
kill %cn3s% all
kill %cn1s% all
kill %**ax% all
kill %**in% all
kill %w2% all
kill %w1% all
kill %s% all
緩和曲線 迴旋曲線的計算
下午打算計算驗證一下迴旋曲線公式，想偷個懶從網上down個 對比驗證，但未檢索到現成 關於迴旋曲線離散點 x,y 的計算分析眾多，但介紹並不詳細，甚至公式和計算方法也都不同。為不枉費半天的查詢 驗證工作，也為方便同仁參考，本文介紹基礎的迴旋線計算過程及 本次計算未涉及座標變換。本部分可忽略不看，直接...
Dubins曲線的計算與實現
dubins曲線計算與實現 1.dubins曲線簡介 dubins曲線是在滿足曲率約束和規定的始端和末端的切線 進入方向 的條件下，連線兩個二維平面的最短路徑，而且限制目標只能向前行進。下圖給出乙個dubins曲線的例子，其中a表示起始位置，該出的箭頭表示當前的速度方向，b表示目標位置，箭頭表示結束...
cocos2d x做圓滑曲線的動畫
參考 1.cccardinalsplineby 2.cccatmullromto cccatmullromby 3.使用cccardinalsplineby自定義路徑 3.貝塞爾曲線生成演算法 這裡先介紹另乙個經典的曲線逼近方法，稱作bezier曲線。想必學過圖形影象的都應該知道啦，所以概念性問題就...