Description Usage Arguments Details Value Examples
Convert observations stored in the data.frame format into the array format.
1 |
contrast |
the dataset containing the observations in rows and the contrast parameters in columns. vector or data.frame. REQUIRED. |
coords |
the spatial coordinates of the observations. matrix with a number of rows equal to the number of rows of |
format |
the format of the output. Can be |
default_value |
the element used to fill the missing observations. numeric. |
range.coords |
the maximum coordinate in each dimension to be considered. numeric vector with length equal to the number of columns of |
FUNCTION:
If contrast
contains several parameters, they are treated one at a time and the result is returned in the form of a list.
If range.coords
is NULL
then the maxima coordinates are those of the coords
argument.
If only one parameter is specified and the format is set to "any"
then a vector is returned.
a list containing :
[[contrast]]
: list containing the new contrast in the new format.
[[coords]]
: a data.frame containing the coordinates of each observation.
[[unique_coords.group]]
: two list containing the possibles coordinates in each dimension.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 | #### 1- with simulated data ####
## simulate
set.seed(10)
n <- 4
Y <- rnorm(n^2)
## conversion
res1 <- dt2array(contrast = Y, coords = expand.grid(1:n + 0.5, 1:n + 0.5))
res2 <- dt2array(contrast = Y, coords = expand.grid(1:n, 1:n), format = "matrix")
res3 <- dt2array(contrast = Y, coords = expand.grid(2 * (1:n), 2 * (1:n)))
res4 <- dt2array(contrast=cbind(Y ,Y, Y), coords = expand.grid(2 * (1:n), 2 * (1:n)),
range.coords = c(10,10))
## display
par(mfrow = c(2,2), mar = rep(2,4), mgp = c(1.5,0.5,0))
fields::image.plot(unique(res1$coords[,1]), unique(res1$coords[,2]), res1$contrast[[1]],
xlab = "", ylab = "")
fields::image.plot(unique(res2$coords[,1]), unique(res2$coords[,2]), res2$contrast,
xlab = "", ylab = "")
fields::image.plot(res3$contrast[[1]])
fields::image.plot(res4$contrast[[2]])
#### 2- with MRIaggr data ####
## load a MRIaggr object
data("MRIaggr.Pat1_red", package = "MRIaggr")
carto <- selectContrast(MRIaggr.Pat1_red, param = "DWI_t0", format = "vector")
coords <- selectCoords(MRIaggr.Pat1_red)
coords[,1] <- coords[,1] + 30
coords[,2] <- coords[,2] + 15
## converion 1
array.DWI_t0 <- dt2array(carto, coords = coords, default_value = 1000)$contrast[[1]]
# display
fields::image.plot(min(coords[,1]):max(coords[,1]), min(coords[,2]):max(coords[,2]),
array.DWI_t0[,,1], xlab = "i", ylab = "j")
## conversion 2
array.DWI_t0 <- dt2array(contrast=carto, coords = coords, default_value = 1000,
range.coords = c(128,128,3))$contrast[[1]]
# display
fields::image.plot(1:128, 1:128, array.DWI_t0[,,1], xlab = "i", ylab = "k")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.