Anyone know how I can use Rotations.jl to rotate a 3D point around another? I'm struggling with the README.md:
function rotateby(pt::Point3D, about::Point3D, angleX, angleY, angleZ)
...
return pt1
end
It seems like they have XYZ-axis rotations. That's totally not my field, but I would guess that given points x and y and a rotation R around the origin, a rotation of x around y would simply decompose to R*(x-y) + y
(?)
I think you'll need to better define exactly what you mean by angleX, angleY, angleZ.. It sounds from your comment in the issue you opened that you might need to get familiar with the different ways rotations can be represented -- probably what you want corresponds to one of them but from what you've said so far I can't tell which.
There are lots of resources I could point you to but it'd be useful to know where you're starting from. E.g. https://rotations.berkeley.edu/contents/ is very comprehensive but pretty mathy and likely overkill.
So, maybe you can say a bit more about where angleX/Y/Z come from and the context of how you intend to use rotateby
.
Hi Patrick! Thanks for helping me out, I appreciate it.
So let's say I'm thinking of a point A in 3D space, at Point3D(1, 1, 0)
. I want to rotate this by π radians around the z axis. So I get B, Point3D(1, -1, 0)
. Then I want to rotate this around the y axis by π radians, so I get C, Point3D(-1.0, -1.0, 0)
. I could also get to C by rotating A around the x-axis.
Screenshot-2021-01-25-at-09.38.42.png
In this case obviously rotations can be either clockwise or anticlockwise (I used π for simplicity ), but any convention is OK.
This was the kind of thing I was thinking for "3D rotations made easy". I thought it could probably be done in a few lines... Perhaps I'm wrong.
I think what you need here are relative poses. You can define/think of a point as a reference frame, and since you want the axes to be parallel, you can have all games at 0 relative angle to the origin. Then convert one point (or frame) to another point's frame of reference (simply subtract coordinates) and then use a rotx/roty/rotz. Then compose the result with the reference point's translation. In all it'll look like this:
rotz(p2-p1, angle) + p1
The same was also mentioned by Florian above
Yes, that's what I meant. I think what @Patrick Bouffard is hinting it is, that simply saying you want tp rotate degrees a, b, c around the x, y, z axes doesn't get you too far, because the order of application is not irrelevant for rotations! Imagine you start at a point on your x-axis, say (1., 0., 0.) and you just want to rotate 90° around the x and y axes. If you start by rotating around x, you don't change your point on the x-axis, so the rotation around y will bring you to (0., 0., ±1.). But if you start by rotating around y, you go to (0., 0., ±1.) initially and suddenly, your rotation around x is now actually doing something and you end up at (0., ±1., 0.). So you either have to do one rotation at a time or specify which convention you use when you say a° around x and b° around y.
So that's why I've been having trouble! :) Perhaps I shouldn't be trying to do three rotations at once, just provide a rotateX() rotateY() and rotate(Z) functions instead...
That said, I found the package you are probably looking for:
https://github.com/JuliaGeometry/CoordinateTransformations.jl
Which was linked from Rotations.jl. It has a recenter method:
"""
recenter(trans::Union{AbstractMatrix,Transformation}, origin::AbstractVector) -> ctrans
Return a new transformation ctrans
such that point origin
servesas the origin-of-coordinates for trans
. Translation by ±origin
occurs both before and after applying trans
,...
"""
Yup, or there are also matrices for rotation along any axis. So if you can define one axis around which you have one rotation corresponding to 3 rotations from the origin, then that might be what you want. But only limited transformations will fit this criteria
Yes I'm sure these packages can easily do this sort of thing without breaking into a sweat. It's just sometimes hard for me to find out how to do the simple things. I suppose they're not really intended for that though.
Did angleX/Y/Z come from another source and you are trying to apply them correctly? Or are you saying that you need to rotate points and are trying to work out a way to do so? If it's the former then it's important you interpret the 3 values correctly--there are a number of ways you could do so. Some of them are more common than others depending on if you're doing aerospace, or physics, or whatever. But if it's the latter the problem is reversed: there are a lot of different ways you can parametrize rotations and they all have pros and cons that depend on the application.
If you give us a bit more context that will help give you better advice!
Well as Florian said I'm probably getting confused by the fact that rotations have to applied in a particular order, so perhaps it isn't a brilliant idea to try to supply all rotations in one function call.
Screenshot-2021-01-25-at-17.36.33.png
I originally thought that there might be an easy alternative to all these rotation functions by just switching to calls to Rotations.jl or CoordinateTransformations.jl functions, but I wasn't able to get much from the READMEs, since they're obviously aimed at a different level of user (PhD rather than GCSE :)). I'll probably stick to what I have at the moment, just switching to separate functions...
Which library are you using to display text in 3D plots? I'd like to do this for a package too
I use a small package called Thebes.jl, but it's fairly simplistic, and I want to improve the rotations code... :joy:
Haha, fair point about the docs, though don't let that intimidate you! Assuming you don't just have this one thing to do and then never touch rotations in 3D again I think it would be good for you to find a resource that meets you where you're at mathematically. You definitely don't need a PhD to grasp the basics or even the advanced stuff.
But you will need to be comfortable at least with vectors, matrix multiplication, and basic trigonometry (just sin & cos). If you are then it'd be easy to help you, but again also knowing more about what you're trying to do here will really be useful.
Thanks Patrick! I like to encounter challenges sequentially, but multiple things arrive very quickly in JuliaGeometry!:joy:
Fair enough. Neat package! Can you say more specifically what you want to improve? It looks like you have something that works.
Well either I don't understand my code or my code doesn't understand me, I'm not sure which it is. So I thought I'd investigate using someone else's code, which - I usually find - is better... :joy: Perhaps I just have to watch more videos...
Thanks for all the help! The answer is:
using Rotations
function rotateby(point::Point3D, about::Point3D, angleX, angleY, angleZ)
RotXYZ(angleX, angleY, angleZ) * (point - about) + about
end
(Chris Foster helped me out!)
Last updated: Nov 06 2024 at 04:40 UTC