sciviews / aurelhy

Hydrometeorological interpolation using the AURELHY method

'SciViews::R' Dialect for Data Processing and Visualization

{SciViews} mainly provides the SciViews::R dialect through the function of the same name. It loads a series of tidyverse and SciViews packages in order to supplement base R with functions to implement that dialect. See, for instance the books Science des Données Biologiques I and Science des Données Biologiques II (in French) for extensive examples of the use of SciViews::R.

Installation

{SciViews} is available from CRAN, but it is an old version. You should install it from the SciViews R-Universe. To install this package and its dependencies, run the following command in R:

install.packages('SciViews', repos = c('https://sciviews.r-universe.dev',
  'https://cloud.r-project.org'))

An older version of {SciViews} can be installed from CRAN:

install.packages("SciViews")

You can also install the latest development version. Make sure you have the {remotes} R package installed:

install.packages("remotes")

Use install_github() to install the {SciViews} package from GitHub (source from main branch will be recompiled on your machine):

remotes::install_github("SciViews/SciViews")

R should install all required dependencies automatically, and then it should compile and install {SciViews}.

Further explore {SciViews}

You can get further help about this package this way: Make the {SciViews} package and all the other packages required by the SciViews::R dialect available in your R session:

SciViews::R()

Get help about this package:

library(help = "SciViews")
help("SciViews-package")
vignette("SciViews") # None is installed with install_github()

For further instructions, please, refer to these help pages at https://www.sciviews.org/SciViews/.

Code of Conduct

Please note that the {SciViews} package is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.

Note to developers

This package used to be developed on R-Forge in the past. However, the latest R-Forge version was moved to this Github repository on 2018-01-05 (SVN version 569). Please, do not use R-Forge anymore for SciViews development, use this Github repository instead.

aurelhy's Issues

Ununderstood problem trying to apply the AURELHY method

Hello,

I am trying to use this interpolation technique to achieve the same result as Morocco, in my case it is the French Pyrénées. You will find my code below :
library('aurelhy')

#read ascii grid and convert it into a geotm object

dem_path <- 'D:/These_PLVG-IMFT-ENIT/Prog/SpatialisationPluies_BVGP/PERFECT_RASTER/intdem30.asc'
GP_dem <- read.geotm(file=dem_path, type = "ascii")

#Read ASCII grid into and convert it as a geomask object

mask_path <- 'D:/These_PLVG-IMFT-ENIT/Prog/SpatialisationPluies_BVGP/PERFECT_RASTER/mask30.asc'
GP_mask <- read.geomask(file=mask_path, type = "ascii")

#Read station data into a dataframe

rain_path <- 'D:/These_PLVG-IMFT-ENIT/Prog/SpatialisationPluies_BVGP/POSTES_BVGP_EV2012.csv'
rain_df <- read.csv(rain_path)

#Rain interpolation
maurelhy <- aurelhy(GP_dem, GP_mask)

And the error that I get is the following one :
Error in par["x2"] <- par["x1"] + nrow(x) * par["size"] :
replacement has length zero

I don't really understand this error message, could you please tell me what does it mean ?
You will find attached the files I used to apply the aurelhy method.

Thanks in advance for your time and help,
Rabab
GP_AUREHYINTERPOLATION.zip

Mask for aurelhy

Hi there,

Just trying to use this interpolation technique to achieve the same result as Morocco, in my case it is New Zealand.

`library('aurelhy')

Read ASCII grid into a geotm object

dem_path <- './Work/Interpolation/Data/otago005_1.asc'
nz_dem <- read.geotm(file=dem_path, type = "ascii")

Read ASCII grid into a geomask object

mask_path <- './Work/Interpolation/Data/otago_mask_3.asc'
nz_mask <- read.geomask(file=mask_path, type = "ascii")

Read station data into a dataframe

rain_path <- './Work/Interpolation/Data/NZ_query_selection.csv'
rain_df <- read.csv(rain_path)

Interpolate

maurelhy <- aurelhy(nz_dem, nz_mask)

Error in apply(mask, 1, any) : dim(X) must have a positive length `

Note:

1. I have no hint why I get the error at aurelhy function call. What does this mean?
2. Also assuming the TRUE values in the mask as targeted area and FALSE values as not interested?
3. The way I created the mask is not correct as you can see from the terrain image. Just tried to extract a small portion with a certain cell value.

Files attached:
NZ_query_selection.zip

Thanks for the time and effort,

Vijay Paul