tests/testthat/test-hhcf.R

In thomasgredig/nanoscopeAFM: Load and Analyze Atomic Force Microscopy Data Files

test_that("Height height correlation function: check roughness", {
  filename = AFM.getSampleImages(type='tiff')
  a = AFM.import(filename)
  a = AFM.flatten(a)
  h1 = AFM.hhcf(a, numIterations = 1e5, degRes=1000, dataOnly = TRUE)
  # AFM.hhcf(a, numIterations = 4e4, degRes=1000, addFit=TRUE)

  expect_equal(mean(h1$g) / (2*AFM.math.params(a)$Rq^2), 1, tolerance = 2e-2)
})