Show Posts
This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.
chunk.crs = PhotoScan.CoordinateSystem('LOCAL_CS["Local CS",LOCAL_DATUM["Local Datum",0],UNIT["millimetre",1]]')to the script.# Project to automate Agisoft Process
# Automation of the floor alignment and setting the correct bounding box.
# Initially created by: Richard Garsthagen - www.pi3dscan.com
# SCRIPT DERIVED FROM: http://www.pi3dscan.com/index.php/instructions/item/agisoft-how-to-process-a-scan-with-projection
# Changed by Vincent
#
import PhotoScan
import math
doc = PhotoScan.app.document
chunk = doc.chunk # supposed to point to the active chunk
#chunk.crs = PhotoScan.CoordinateSystem("EPSG::1025") # set our coordinate-sys to "mm"
chunk.crs = PhotoScan.CoordinateSystem('LOCAL_CS["Local CS",LOCAL_DATUM["Local Datum",0],UNIT["millimetre",1]]')
# py
#
# |
# |
# |
#
# p0 ------------ px
p0 = "target 2"
px = "target 3"
py = "target 1"
distancepx = 0.15
distancepy = 0.2
regionSizeX = 0.15
regionSizeY = 0.35
regionSizeZ = 0.1
mp0 = 0
mpy = 0
mpx = 0
fp0 = 0
fpy = 0
fpx = 0
#setting for Z up, Y forward
vector0 = PhotoScan.Vector((0,0,0))
vectorY = PhotoScan.Vector((0,distancepy,0)) # Specify Y Distance
vectorX = PhotoScan.Vector((distancepx,0,0)) # Specify X Distance
#doing some marker magic, finding/assigning them, I guess...
c = 0
for m in chunk.markers:
if m.label == p0:
mp0 = c
fp0 = 1
m.reference.location = vector0
m.reference.enabled = 1
print ("Found center point: ", vector0)
#print("chunk.markers[0].position", chunk.markers[0].position)
if m.label == py:
mpy = c
fpy = 1
m.reference.location = vectorY
m.reference.enabled = 1
print ("found Y point: ", vectorY)
#print("chunk.markers[y].position", chunk.markers[1].position)
if m.label == px:
mpx = c
fpx = 1
m.reference.location = vectorX
m.reference.enabled = 1
print ("found X point: ", vectorX)
#print("chunk.markers[x].position", chunk.markers[2].position)
c = c + 1
if fp0 and fpx and fpy:
print ("Found all markers")
chunk.updateTransform()
else:
print ("Error: not all markers found")
# creating a new region from the existing one and rotating it
newregion = chunk.region
T = chunk.transform.matrix
v_t = T * PhotoScan.Vector( [0,0,0,1] )
m = PhotoScan.Matrix().diag([1,1,1,1])
m = m * T
s = math.sqrt(m[0,0] ** 2 + m[0,1] ** 2 + m[0,2] ** 2) #scale factor
R = PhotoScan.Matrix( [[m[0,0],m[0,1],m[0,2]], [m[1,0],m[1,1],m[1,2]], [m[2,0],m[2,1],m[2,2]]])
R = R * (1. / s)
newregion.rot = R.t()
#scaling the new region
dist = chunk.markers[mp0].position - chunk.markers[mpy].position
dist = dist.norm()
ratio = dist / distancepy
newregion.size = PhotoScan.Vector( [regionSizeX*ratio, regionSizeY*ratio, regionSizeZ*ratio] )
# moving the new region to a different location
posReg_world = PhotoScan.Vector( [0.1 , (0.5*regionSizeY) , (0.5*regionSizeZ+0.002) ] )
posReg_int = chunk.transform.matrix.inv().mulp(posReg_world)
newregion.center = posReg_int
chunk.region = newregion
chunk.updateTransform()
print ("Bounding box should be aligned now")
Never mind. I managed to do it
). But as for FirePro cards they should not be very fast in terms of PhotoScan performance.
