R を使用して、イベリア半島から収集されたデータの補間周波数のマップを実行しようとしています。(このようなものhttps://gis.stackexchange.com/questions/147660/strange-spatial-interpolation-results-from-ordinary-kriging )
私の問題は、autokrige 関数の属性 new_data に何らかのエラーがあるため、プロットに補間データが表示されないことです。
https://cran.r-project.org/web/packages/automap/automap.pdf new_data: 予測位置を含む sp オブジェクト。new_data は、ポイント セット、グリッド、またはポリゴンにすることができます。NA を含めてはなりません。このオブジェクトが提供されない場合、デフォルトが計算されます。これは、input_data の凸包を取得し、その凸包に約 5000 個のグリッドセルを配置することによって行われます。
問題は、この new_data 属性を回避すると、クリギング値の適切なプロットが得られるため、R がポリゴンに変換されたマップをうまく読み取れないことだと思います。しかし、イベリア半島の形がよくわかりません。誰か助けてくれませんか？本当にありがたく思います
ここで私のデータを見ることができます: http://pastebin.com/QHjn4qjP
実際のコード: データ座標を UTM 投影に変換するため、エラー メッセージは表示されませんが、最後のプロットは補間されず、マップ全体が 1 つの単色で表示されます :(
setwd("C:/Users/Ruth/Dropbox/R/pruebas")
#Libraries
library(maps)
library(mapdata)
library(automap)
library(sp)
library(maptools)
library(raster)
library(rgdal)
####################MAPA#############
#obtain the grid of the desired country
map<-getData('GADM', country='ESP', level=0)
#convert the grid to projected data
mapa.utm<-spTransform(mapa3,CRSobj =CRS(" +proj=utm +zone=29 +zone=30 +zone=31 +ellps=WGS84"))
###############################Datos#######################
#submit the data
data1<-read.table("FRECUENCIASH.txt",header=T)
head(data1)
attach(data1)
#convert longlat coordinates to UTM
coordinates(data1)<-c("X","Y")
proj4string(data1) = CRS("+proj=longlat +datum=WGS84")
data1.utm=spTransform(data1, CRS("+proj=utm +zone=29 +zone=30 +zone=31 +ellps=WGS84 "))
######################Kriging interpolation #####################
#Performs an automatic interpolation
kriging_result<-autoKrige(Z~1,data1.utm,mapa.utm,data_variogram = data1.utm)
#Plot the kriging
result1<-automapPlot(kriging_result$krige_output,"var1.pred",sp.layout = list("sp.points", data1.utm));result1
result2<-plot(kriging_result);result2