IQ(var₁=>param₁, var₂=>param₂, ...)

Image quilting simulation solver as described in Hoffimann et al. 2017.

Parameters

Required

• trainimg - Training image from which to extract tiles
• tilesize - Tuple with tile size for each dimension

Optional

• overlap - Overlap size (default to (1/6, 1/6, ..., 1/6))
• path - Simulation path (:raster (default), :dilation, or :random)
• mapping - Data mapping method (default to NearestMapping())
• inactive - Vector of inactive voxels (i.e. CartesianIndex) in the grid
• soft - A pair (data,dataTI) of geospatial data objects (default to nothing)
• tol - Initial relaxation tolerance in (0,1] (default to 0.1)

Global parameters

Optional

• progress - Whether to show or not the estimated time duration (default to true)
• threads - Number of threads in FFT (default to number of physical CPU cores)
• rng - Random number generator (default to Random.GLOBAL_RNG)

References

• Hoffimann et al 2017. Stochastic simulation by image quilting of process-based geological models.
• Hoffimann et al 2015. Geostatistical modeling of evolving landscapes by means of image quilting.

iqsim(trainimg::AbstractArray{T,N}, tilesize::Dims{N},
      simsize::Dims{N}=size(trainimg);
      overlap::NTuple{N,<:Real}=ntuple(i->1/6,N),
      soft::AbstractVector=[], hard::Dict=Dict(), tol::Real=.1,
      path::Symbol=:raster, nreal::Integer=1,
      threads::Integer=cpucores(), debug::Bool=false,
      progress::Bool=true, rng::AbstractRNG=Random.GLOBAL_RNG)

Performs image quilting simulation as described in Hoffimann et al. 2017.

Parameters

Required

• trainimg is any Julia array
• tilesize is the tile size

Optional

• simsize is the size of the simulation grid (default to training image size)
• overlap is the percentage of overlap (default to 1/6 of tile size)
• soft is a vector of (data,dataTI) pairs (default to none)
• hard is a dictionary mapping coordinates to data values (default to none)
• tol is the initial relaxation tolerance in (0,1] (default to .1)
• path is the simulation path (:raster, :dilation or :random)
• nreal is the number of realizations (default to 1)
• threads is the number of threads for the FFT (default to all CPU cores)
• debug informs whether to export or not the boundary cuts and voxel reuse
• progress informs whether to show or not estimated time duration
• rng is the random number generator (default to Random.GLOBAL_RNG)

The main output reals consists of a list of realizations that can be indexed with reals[1], reals[2], ..., reals[nreal]. If debug=true, additional output is generated:

reals, cuts, voxs = iqsim(..., debug=true)

cuts[i] is the boundary cut for reals[i] and voxs[i] is the associated voxel reuse.