Stack Overflow
Stack Overflow

Questions tagged [procrustes]

12 questions
2
votes
1 answer

Method for comparing 2 3-dimensional networks (igraphs) in R

I'm looking for a way to compare 2 (or more) igraph objects in R. These are trajectories in 3-dim which are a network of nodes and edges but are not necessarily the same number of either, just that they have a corresponding coordinate in 3-dim. I…
No Future
  • 41
  • 3
2
votes
0 answers

How can I rotate a 2d image using a target image, landmark coordinates, the least squares approach, and a rotation matrix?

I have two 2d images, one is the source image and the other is a target image; I need to rotate the source image to match the target image using python (scikit & numpy). I have 3 landmark coordinates for each image, as follows: image1_points =…
user5779866
  • 43
  • 4
2
votes
0 answers

vegan protest() finds a high correlation between significantly different data sets?

Assume that we have the following data set.seed(123) Xmat <- matrix(rnorm(2000), ncol = 200, nrow = 10) Ymat <- Xmat[,1:6] Where Ymat is subsampled from Xmat. Running protest() we get a high correlation between the two…
Antsushi
  • 167
  • 10
2
votes
1 answer

Linear projection matrix by solving the procrustes problem

I have two matrices: target = np.array([[1, 1, 1, 1, 1], [2, 2, 2, 2, 2], [3, 3, 3, 3, 3]]) source = np.array([[11, 11, 11, 11, 11], [22, 22, 22, 22, 22], [33, 33, 33, 33, 33]]) and I…
Minions
  • 5,104
  • 5
  • 50
  • 91
1
vote
0 answers

Issue opening geomorph package

Currently having issue opening geomorph in R. I recently updated R to the newest version. I also updated the packages "RRPP" and "rgl" after finding a thread from 2020 (Geomorph not functioning after update) where someone also had issues opening…
Kailie Batsche
  • 11
  • 2
1
vote
0 answers

Weird behavior with orthogonal procrustes

Given two matrices mat1 and mat2, each one representing a set of 3 points in a 3 dimensional space, I want to rotate mat2 so that it matches mat1 using orthogonal procrustes : mat1 <- t(matrix(c(1.16, 0.21, 0.11, -0.78, -0.02, -0.73, -0.37, -0.18,…
mhaddad
  • 119
  • 1
  • 10
1
vote
1 answer

R (SensoMineR) is not producing all plots when using pmfa function

I am trying to perform a Procrustean multiple factor analysis on my dataset. The dataset was created using a technique called 'napping' that is used in sensory science. The problem that I have is that one of the plots pmfa is supposed to produce is…
KatP
  • 45
  • 7
0
votes
1 answer

How do I align 3D landmark data to a Procrustes fit in R?

I'm using the r package called Morpho to create predicted soft tissue nose shapes for populations using the underlying hard tissue (bone). https://cran.r-project.org/web/packages/Morpho/Morpho.pdf Using guidelines provided by the PhD of author of…
Catherine Sulzmann
  • 3
  • 2
0
votes
0 answers

Solving the Orthogonal Procrustes problem in Python

I have two sets of points in 3D space, and I want to find the optimal rotation matrix that aligns one set of points to the other, in a least-squares sense (Orthogonal Procrustes Problem). I've read that this can be done using the Orthogonal…
Simon
  • 1
  • 2
0
votes
1 answer

Get the Transformation Matrix From the SciPy Procrustes Implementation

The Procrustes library has an example where it demonstrates how to get the Transformation Matrix of two matrices by solving the Procrustes problem. The library seems to be old and doesn't work in Python 3. I was wondering if there's any way to use…
Cypher
  • 2,374
  • 4
  • 24
  • 36
0
votes
0 answers

Procrustes problem on tensorflow, problems for SVD calculation without projection of data onto leading POD modes

I am trying to solve the Procrustes problem min |Y - K* X| I have implemented SVD and the Procrustes problem in a code that I inherited. I have written 2 versions, First version: I am projecting X (e.g. self.__psiNNx) and Y (e.g. self.__psiNNy)…
Carlos
  • 15
  • 4
0
votes
0 answers

Plotting a function using matplotlib

I have the following function def procrustes(data1, data2, data3): ... which returns return mtx1, mtx2, mtx3_prime, disparity mxt1, mtx2 and mtx3_prime, are all arrays. disparity is a float I want to plot mtx1, mtx2, mtx3_prime…
lanny kayz
  • 59
  • 9