Get center coordinates

[r0xx]

How can I get the center coordinates from the project and set them to the same coordinates but with another bounding box size?

[Alexey Pasumansky]

Hello r0xx,

Could you please provide additional info on this matter? I just don't get what task you are trying to solve.

[r0xx]

I am trying to automatically export mesh files in different sizes.
So for example I have a grid with center coordinates (in coordinate system 21781) and generate the dense cloud with a bounding box size of 1 km.
How could I get the center coordinates of the current bounding box, set the size of the bounding box to 500 meter, generate a new mesh, export this mesh and do the same thing again with 250 meters?

Basically I have all the parts except how I can get the center coordinates from the current bounding box.

[Alexey Pasumansky]

Hello r0xx,

You can use the following script as a reference:
http://wiki.agisoft.com/wiki/Split_in_chunks.py

In splitChunks() function copies of the original chunk are created with the different bounding box.

Note that chunk.region.center is in the internal coordinate system, so if you wish to use geographic coordinates you need to convert them to the internal system:

chunk.transform.matrix.inv().mulp(chunk.crs.unproject(point)), where point is PhotoScan.Vector with geographic coordinate system values.

[r0xx]

This is what I'm using at the moment but this doesn't work:

doc = PhotoScan.app.document

  # Get the current chunk
  chunk  = doc.chunk
 
  c = chunk.region.center
  centerCoords = chunk.transform.matrix.inv().mulp(chunk.crs.unproject(c))

  # Get center coordinates
  print("Center Coordinates: " + str(centerCoords))

This does output some other coordinates... What is wrong with the code?

[Alexey Pasumansky]

Hello r0xx,

chunk.region.center coordinates are already in internal system, they are not geographic.

[r0xx]

Ok thank you. And how can I print the geographic coordinates now?

[Alexey Pasumansky]

Hello r0xx,

To print geographic coordinates of the region center you can use the following:

Code: [Select]

chunk.crs.project(chunk.transform.matrix.mulp(chunk.region.center))

[r0xx]

Awesome thanks alot!