5

How can I devise an R-code that can derive a rectangular bounding box (with 90-degree angle on each corner) for a set of points with arbitrary orientation? e.g. if the 2 longer edges have a 45-degree angular direction, or 15-degree, or any angle. It seems to have been solved on the paper below, but I cannot read their mathematical notation.

http://www.cccg.ca/proceedings/2004/29.pdf

It is partially solved on the link below, however, it looks only for the "minimum area rectangle" and not the minimum bounding rectangle for a specified orientation.

Finding minimum-area-rectangle for given points?

|improve this question
  • 1
    Maybe this will help? stackoverflow.com/questions/41833933/… – GBG Oct 17 '17 at 16:56
  • No, the orientation of the bounding box on that thread is always vertical/horizontal, it means its pair of parallel sides are always 90 degrees and 0 degrees respectively. My question is, find the smallest bounding box if the orientation of one of the parallel sides is say 45 degrees, 10 degrees, etc. Thanks for your response. – Jun Folledo Oct 17 '17 at 23:43
  • do you mean minimum oriented bounding box? QGIS has an implementation but that doesn't answer the implementation in R... but it might clarify / be a starting point? – Steven Kay Oct 18 '17 at 19:06
  • @Steven Kay, I think this is what I am looking for, but I'm not familiar in Python coding in QGIS, but I'm willing to study it if needed. Can you please illustrate how to implement it? I suggest that you post your answer, but not in the comment section so you can earn the bounty. Thank you very much! – Jun Folledo Oct 18 '17 at 22:25
  • 2
    The solution is already provided by @whuber in gis.stackexchange.com/questions/22895/… The rectangle is being rotated, isnt that what you want?? – BERA Oct 19 '17 at 9:21
3
+50

I guess you need Principal Component Analysis (PCA), read this publication:

On the Bounding Boxes Obtained by Principal Component Analysis

Look at this as well:

https://www.rdocumentation.org/packages/shotGroups/versions/0.1/topics/getMinBBox

|improve this answer
  • I'm not merely looking for the minimum bounding box, I need a solution wherein I can specify first the orientation of the bounding box, then get the minimum bounding. If I cannot control the orientation of the bounding box, then it won't answer my query here. Thanks for your help though. – Jun Folledo Oct 21 '17 at 0:22
  • Then rotate (constrain: alpha parameter) the points with a rotation matrix. Use the usual boundary box: rdocumentation.org/packages/spatstat/versions/1.6-1/topics/… And rotate back the box. Or, simpler: 1. Rotate the points. 2. Create corners of boundary box (4 points) with the combinations of xmin, xmax, ymin, ymax. 3. Rotate back the corner points, and create the box with convexhull for example. – pnz Oct 21 '17 at 11:22
  • I also thought about that before but rotating the points back and forth would be cumbersome when iterated on thousands of data. The orientedMinimumBounding shared by @Steven Kay may do the trick only if I am able to implement it – Jun Folledo Oct 22 '17 at 8:00
  • You needed an R solution: orientedMinimumBoundingBox (QGIS) Returns the oriented minimum bounding box for the geometry, which is the smallest (by area) rotated rectangle which fully encompasses the geometry. getMinBBox (R) Calculates the vertices of the minimum-area, possibly oriented bounding box given a set of 2D-coordinates. - They are the same. - Look at my 2. link, that's the answer for you question. – pnz Oct 22 '17 at 10:59
  • Both algorithm return the same properties by the way: angle, width, height – pnz Oct 22 '17 at 11:20
1

Then rotate (constrain: alpha parameter) the points with a rotation matrix.

Use the usual boundary box: https://www.rdocumentation.org/packages/spatstat/versions/1.6-1/topics/bounding.box

And rotate back the box.

Or, simpler:

  1. Rotate the points.
  2. Create corners of boundary box (4 points) with the combinations of xmin, xmax, ymin, ymax.
  3. Rotate back the corner points, and create the box with convexhull for example.
|improve this answer

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.