resize bounding box to sparse point cloud

[Seboon]

Thanks Alexey,

I was suspecting it was something like that.
Maybe you have in your magic luggage an already existing "pythonesc"  formula for this kind of calculation : Axe of the mesh compared to the North (0°), resulting in  a rotation degrees  to apply on the Bbox?

As a workaround, i use a mix of your both scipts"BBox to CRS" and "CRS to BBox", see attached script.

Thanks for any suggestion for improving it, escpecillay with the magic formula :-)

Best regards!

[Alexey Pasumansky]

Hello Seboon,

Do you mean manual (custom) rotation from the Northing direction?

[Seboon]

Hello Alexey,

Yeah I suppose it's this kind of idea.
Also drawing a temporary line shape with two points could indicate the difference azimuth to apply.
Thanks in any case!

[chitty1989]

Hello Susann,

For the unreferenced chunks it should be shorter:

Code: [Select]

import PhotoScan, math

doc = PhotoScan.app.document
chunk = doc.chunk
region = chunk.region
T = chunk.transform.matrix

m = PhotoScan.Vector([10E+10, 10E+10, 10E+10])
M = -m

for point in chunk.point_cloud.points:
if not point.valid:
continue
coord = T * point.coord
for i in range(3):
m[i] = min(m[i], coord[i])
M[i] = max(M[i], coord[i])

center = (M + m) / 2
size = M - m
region.center = T.inv().mulp(center)
region.size = size * (1 / T.scale())

region.rot = T.rotation().t()

chunk.region = region
print("Script finished.")

Hey Alexey,

Is there any way I can run this script as part of a batch?  At the moment the batch stops at the first chunk and doesnt run the script through the rest of the chunks!