From 34fe1888251a7c262fab49492c98b43c8e9a8a29 Mon Sep 17 00:00:00 2001 From: Elias Carvalho Date: Thu, 26 Oct 2023 14:47:37 -0300 Subject: [PATCH] Apply suggestions --- ext/GeoStatsProcessesImageQuiltingExt.jl | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/ext/GeoStatsProcessesImageQuiltingExt.jl b/ext/GeoStatsProcessesImageQuiltingExt.jl index aa68512..50fb1ee 100644 --- a/ext/GeoStatsProcessesImageQuiltingExt.jl +++ b/ext/GeoStatsProcessesImageQuiltingExt.jl @@ -10,7 +10,8 @@ using Meshes using GeoTables using GeoStatsBase: initbuff -using ImageQuilting: iqsim, lin2cart +using ImageQuilting: iqsim + using GeoStatsProcesses: IQ, RandSetup import GeoStatsProcesses: randprep, randsingle @@ -19,6 +20,7 @@ function randprep(::AbstractRNG, process::IQ, setup::RandSetup) # retrieve domain info sdomain = setup.domain simsize = size(sdomain) + lin2cart = CartesianIndices(simsize) Dim = embeddim(sdomain) # retrieve process paramaters @@ -52,7 +54,7 @@ function randprep(::AbstractRNG, process::IQ, setup::RandSetup) # create hard data object linds = findall(mask[var]) - cinds = [lin2cart(simsize, ind) for ind in linds] + cinds = [lin2cart[ind] for ind in linds] hvals = view(buff[var], linds) hard = Dict{CartesianIndex{Dim},Real}() for (ind, val) in zip(cinds, hvals)