minor: update doc

Author: kongdd

DESCRIPTION

@@ -30,7 +30,7 @@ Suggests:
 License: GPL-3
 Encoding: UTF-8
 LazyData: true
-RoxygenNote: 7.2.3
+RoxygenNote: 7.3.2
 Config/testthat/edition: 3
 Roxygen: list(markdown = TRUE)
 VignetteBuilder: knitr

NAMESPACE

@@ -64,8 +64,10 @@ importFrom(fields,image.plot)
 importFrom(fields,predictSurface)
 importFrom(glue,glue)
 importFrom(hydroGOF,KGE)
+importFrom(magrittr,"%$%")
 importFrom(magrittr,"%<>%")
 importFrom(magrittr,"%>%")
+importFrom(magrittr,set_names)
 importFrom(matrixStats,rowMeans2)
 importFrom(plyr,llply)
 importFrom(plyr,round_any)

R/spInterp_Tps.R

@@ -19,7 +19,7 @@ spInterp_Tps <- function(points, dat, range, res = 1,
   if (nrow(points) != nrow(dat)) stop("Length of points and dat should be equal!")
 
   if (is.null(Z)) ZGrid = NULL
-  
+
   grid <- query_gridInfo(range, res)
   m = Tps(points, dat, Z = Z, ...)  
   pred = predictSurface(m, grid, extrap = TRUE, ZGrid = ZGrid)

R/spInterp_adw.R

@@ -54,7 +54,6 @@ spInterp <- function(points, dat, range, res = 1,
       .[, .(value = weighted.mean(x, w, na.rm = TRUE)), .(lon, lat)] %>% 
       { merge(grid, ., all.x = TRUE)$value }
   }, .parallel = .parallel, .progress = "text") %>% do.call(cbind, .)
-  
   list(weight = weight, coord = grid, predicted = pred) %>% set_class("spInterp")
 }
 

man/spInterp-package.Rd

@@ -32,12 +32,11 @@ dat <- dat_RH[, .(lon, lat, alt, RH)]
 X <- dat[, .(lon, lat, alt)] %>% as.matrix()
 Y <- dat[, .(RH)] %>% as.matrix()
 
-res <- 1
+res   <- 1
 range <- c(70, 140, 15, 55)
 r_dem <- get_chinadem(res = res) %>% round(1)
 f_dem <- "output/dem.asc"
-
-write_dem(r_dem, f_dem, digits = 1)
+# write_dem(r_dem, f_dem, digits = 1)
 ```
 
 # 1. anusplin

scripts/anusplin.Rmd

@@ -1,7 +1,6 @@
 ```{r}
 library(Ipaper)
 library(data.table)
-
 library(tidymet)
 library(lubridate)
 library(terra)
@@ -30,11 +29,11 @@ df_mat2 = df_mat[year(date) >= 1960]
 ## 1. 数据量过大，每年单独处理
 
 ```{r}
-cdd = 500
-# l = cal_weight(X, range, res=cellsize, cdd = cdd, wFUN = "wFUN_adw")
-# weight = do.call(rbind, l) %>% as.data.table()  
-r <- spInterp_adw(X, Y, range, res = cellsize, cdd = cdd)
-print(str(r))
+# cdd = 500
+# # l = cal_weight(X, range, res=cellsize, cdd = cdd, wFUN = "wFUN_adw")
+# # weight = do.call(rbind, l) %>% as.data.table()  
+# r <- spInterp_adw(X, Y, range, res = cellsize, cdd = cdd)
+# print(str(r))
 ```
 
 
@@ -43,14 +42,17 @@ print(str(r))
 > 只做目标区域的插值
 
 ```{r}
-cellsize = 0.01
+# bbox(105.62916666665245, 28.827500000006523, 116.74499999997568, 34.548333333334654)
+range <- c(105, 117, 28, 35)
+# poly <- st_rect(range)
+# write_shp(poly, "湖北_shape.shp")
+cellsize = 0.05
 range = st_range(basin, cellsize)
 
 ra = make_rast(range, cellsize)
 values(ra) = 1:size(ra)
-
-plot(ra)
-plot(vect(basin), add=TRUE)
+# plot(ra)
+# plot(vect(basin), add=TRUE)
 ```
 
 
@@ -75,6 +77,7 @@ print(str(r))
 ## 3. 转为nc文件
 
 ```{r}
+# 之前只做了官山
 res <- r %$% cbind(coord, predicted)
 ra <- rast(res, extent = ext(range))
 
@@ -88,7 +91,6 @@ writeCDF(ra, fout, "prcp", "precipitation", "mm/h",
   overwrite = TRUE)
 ```
 
-
 ### 日降水
 
 ```{r}