2012-07-30 5 views
1

저는 R에서 초보자이며 R에서 보간을하려고합니다 (불충분 한 질문 인 경우 미안). 나는 그들을 보간 할 필요가있는 몇몇 점 (위, 란 및 각 점에 대한 하나의 정수 값)을 가지고 있습니다. 나는 ARCGIS에서 GIS 사용자이지만 ArcGIS는 내가 아는 한 스플라인에 대해 RMSE를 제공하지 않으므로 R에서 생각하고 R을 배울 생각을했습니다. 음모를 꾸미는 것과 같은 기본 자료를 배웠습니다. ,, 포인트를 보간하고 (RMSE를 계산하기 위해) R 언어로해야 할 일을 단계별로 알려주십시오.어떻게하면 R에서 스플라인 보간을 할 수 있습니까

은 정말 어떤 도움을 주셔서 감사합니다 행크

내가 TPS는 얇은 판 스플라인의 약자 필드 패키지에서 Tps 기능을 살펴 것
+0

'mgcv' 패키지를보세요. – mnel

답변

3

. Tps (the documentation에서) :

불규칙한 간격의 데이터에 얇은 판 스플라인 표면을 맞 춥니 다. 평활화 매개 변수는 일반화 된 교차 유효성 검사에 의해 선택됩니다. 가정 모델은 additive Y = f (X) + e이고, 여기서 f (X)는 d 차원 표면이다. 이것은 공간 프로세스 추정의 특별한 경우입니다.

일부 코드 예제는 Tps 설명서를 참조하십시오.

+1

+1 필드는 훌륭한 패키지입니다. 이 방법은 저에게 효과적이었습니다. 예제를 잠깐 게시했습니다 (방법 3 참조 : http://menugget.blogspot.de/2012/03/xyz-geographic-data-interpolation-part.html). –

+0

@Marc, 링크 끝 부분에 여분의')'이 있습니다. 귀하의 게시물은 매우 멋져 보이는데, 고마워요! –

+0

감사합니다. @Paul. 건배 –

3

스플라인을 사용하는 이유가 있습니까? 소프트웨어에서 가지고있는 것 외에 다른 이유가 있습니까? 아마도 http://cran.r-project.org/web/views/Spatial.html을 읽고 Geostatistics http://spatial-analyst.net/book/에 대한 좋은 책을 고려해야하며 Kriging을 사용하여 보간 작업을 할 수 있습니다.

+2

+1. 나는 지형 통계학을 보는 것이 OP를위한 가치가 있다고 동의한다. 그러나 Kriging은 항상 최선의 선택은 아니지만 많은 변수 (데이터의 양, 데이터의 부드러움 등)에 따라 달라집니다. 보간 방법을 비교하는 보고서를 작성한 것은 OP에 대해 흥미로울 수 있습니다. http://www.numbertheory.nl/files/report_evap.pdf. kriging을 수행하는 쉬운 방법은'automap' 패키지를 사용하는 것입니다 (면책 조항, 필자가 작성했습니다). –

+0

OP에서 언급 한 정수가 약간 걱정 스러웠습니다 ... – Spacedman

관련 문제