One thing you could try is to use the middle position as the center of your bounding box and rotate it to use your own covariance matrix vectors as your axes. See, for example, the diagram at http://en.wikipedia.org/wiki/Principal_component_analysis . This does not guarantee that you will receive the maximum possible bounding box - you can see this if you notice that all points, including inside the convex hull, that should not affect the minimum possible bounding box, will affect the eigenvectors - but it can be a decent approximation for some kinds of data.
source share