Как решить ошибку кригинга с внешним дрейфом в R? - PullRequest
Новая
       53

Как решить ошибку кригинга с внешним дрейфом в R?

0 голосов
/ 03 апреля 2020

Я пытаюсь использовать кригинг с внешним дрейфом, используя hydrostm и hydrokrige. Пока, пытаясь интерполировать, не могу хорошо интерполировать. Я использую дем 30 м данных. Код ниже: 

library(zoo)
library(hydroTSM)
library(xts)
library(sp)
library(automap)
library(gstat)
library(sp)
library(rgdal)
library(raster)
bena.ts<- read.csv("EbroPPtsDaily.csv")
bena.gis = read.csv("EbroPPgis.csv")
bena.rs <- raster("DEM-UTM.tif")
bena.sa = readOGR("subBASIN-UTM.shp")
# Covert Raster to spatial grid data frame 
y <- as(bena.rs, "SpatialGridDataFrame")
class(y)
# Make Projection according to study area
p4s <- CRS("+proj=utm +zone=43 +datum=WGS84 +units=m +no_defs +ellps=WGS84
+towgs84=0,0,0")
bena.ts1 <- as.numeric(bena.ts[180, 2:ncol(bena.ts)])
names(bena.ts1) = colnames(bena.ts[180,2:ncol(bena.ts)])
# Computing OK, over the spatial grid defined by the DEM
## Not run: 
x.ok <- hydrokrige(x.ts= bena.ts1, x.gis=bena.gis, 
                   X="EAST_ED50", Y="NORTH_ED50", sname="ID", 
                   formula=value~1,
                   subcatchments= bena.sa,
                   p4s= p4s, 
                   predictors=y,
                   ColorRamp= "Precipitation", 
                   main= "OK Precipitation on the Ebro",
                   arrow.plot= TRUE, 
                   arrow.offset= c(900000,4750000), arrow.scale= 20000,
                   scalebar.plot= TRUE, 
                   sb.offset= c(400000,4480000), sb.scale= 100000)

1 Ответ

0 голосов
/ 07 апреля 2020

bena.ts показывает суточные значения осадков на разных станциях bena.gis показывает координаты в UTM, высоту точек и название суббассейна bena.sa рассказывает о шейп-файле исследования область в UTM bena.rs показывает растровый файл DEM 30 м

Моя проблема заключается в том, что я могу интерполировать его с использованием отметки. Я также конвертирую свой растровый файл DEM 30 м в пространственный массив данных сетки, но все равно он не показывает никаких результатов. Пожалуйста, помогите мне.

...