if (FALSE) {
# Not run:
library(BoneProfileR)
path_Hedgehog <- system.file("extdata", "Erinaceus_europaeus_fem_2-1_small.png",
package = "BoneProfileR")
bone <- BP_OpenImage(file=path_Hedgehog)
bone <- BP_DetectBackground(bone=bone, analysis="logistic")
bone <- BP_DetectForeground(bone=bone, analysis="logistic")
bone <- BP_DetectCenters(bone=bone, analysis="logistic")
bone <- BP_EstimateCompactness(bone, analysis="logistic")
bone <- BP_FitMLCompactness(bone, analysis="logistic")
fittedpar <- BP_GetFittedParameters(bone, analysis="logistic")
bone <- BP_DuplicateAnalysis(bone, from="logistic", to="flexit")
bone <- BP_FitMLCompactness(bone,
fitted.parameters=c(fittedpar, K1=1, K2=1),
fixed.parameters=NULL, analysis="flexit")
compare_AIC(Logistic=BP_GetFittedParameters(bone, analysis="logistic", alloptim=TRUE),
Flexit=BP_GetFittedParameters(bone, analysis="flexit", alloptim=TRUE))
bone <- BP_FitMLRadialCompactness(bone, analysis="logistic")
# Test using the change of orientation using default.angle from BP_EstimateCompactness():
bone <- BP_DuplicateAnalysis(bone, from="logistic", to="logistic_rotation_pi")
# With a pi rotation, the top moves to the bottom and the left moves to the right
bone <- BP_EstimateCompactness(bone, rotation.angle=pi, analysis="logistic_rotation_pi")
bone <- BP_FitMLRadialCompactness(bone, analysis="logistic_rotation_pi")
BP_Report(bone=bone,
analysis=1,
docx=NULL,
pdf=NULL,
xlsx=file.path(getwd(), "report.xlsx"),
author="Marc Girondot",
title=attributes(bone)$name)
BP_Report(bone=bone,
analysis=1,
docx=NULL,
pdf=file.path(getwd(), "report.pdf"),
xlsx=NULL,
author="Marc Girondot",
title=attributes(bone)$name)
BP_Report(bone=bone,
analysis=1,
docx=file.path(getwd(), "report.docx"),
pdf=NULL,
xlsx=NULL,
author="Marc Girondot",
title=attributes(bone)$name)
}
Run the code above in your browser using DataLab