2012-07-30 36 views
1

我是R中的新手,並試圖在R中進行插值(抱歉,如果它看起來是一個蹩腳的問題)。 我有一些點(Lat,Lan和每個點的一個整數值),我需要插入它們。我是一位GIS用戶,我是在ARCGIS中完成這項工作的,但據我所知,ArcGIS並沒有給我RMSE樣條曲線,所以我想在R中學習並學習一些R.我已經學習了basif的東西,比如plotting ,,, ,, 您能否一步一步告訴我,我應該用R語言來插入點(全局水平)並計算RMSE?我如何在R中做樣條插值

我真的很感謝所有幫助 漢克

+0

看看'mgcv'包。 – mnel 2012-07-30 04:36:44

回答

3

我想看看從田間地頭包,其中Tps內代表薄板樣條Tps功能。 Tps(從the documentation):

將薄板樣條曲面擬合爲不規則間隔的數據。平滑參數通過廣義交叉驗證來選擇。假定模型是加性Y = f(X)+ e,其中f(X)是d維的表面。這是空間過程估計的特例。

有關代碼示例,請參閱Tps文檔。

+1

+1域是一個很棒的包。這種方法對我很有幫助 - 我後來發佈了一個示例(請參閱方法3,在這裏:http://menugget.blogspot.de/2012/03/xyz-geographic-data-interpolation-part.html)。 – 2012-07-30 05:43:22

+0

@Marc,在你的鏈接末尾有一個額外的''',打破它。你的帖子看起來很不錯,謝謝! – 2012-07-30 06:04:04

+0

謝謝@保羅。歡呼聲 – 2012-07-30 16:22:35

3

你有沒有理由使用樣條曲線(除了這是你在軟件中)?您應該閱讀http://cran.r-project.org/web/views/Spatial.html並考慮一本關於地統計學http://spatial-analyst.net/book/的好書,然後您可能需要使用Kriging進行插值。

+2

+1。我同意,尋找地質統計學對OP來說是值得的。但克里格並不總是最好的選擇,但這取決於很多變量(數據量,數據的平滑度等)。我後來寫了一篇報告,比較插值方法,對OP來說可能很有趣:http://www.numbertheory.nl/files/report_evap.pdf。執行克里格的一個簡單方法是使用'automap'包(免責聲明,我寫它)。 – 2012-07-30 06:53:32

+0

OP提到的整數確實讓我有些擔心...... – Spacedman 2012-07-30 07:32:14