Here is another script: This one uses neither intersectRayEx nor RayMeshGridIntersect.
It does not convert the object either. The intersection is calculated directly in the script.
There is probably a means to optimize that. In particular we could add a test to know if the ray crosses the bounding box of the object. Currently if the ray not intersect the object the routine test all faces (and only a little number if the ray intersects the object)
struct intersectionStruct ( type=0, meshFace=0, polyFace=0, edge=0, vertex=0, point=0 )
fn isFacePoint pA pB pC pD = (
local vAB=normalize (pB-pA)
local vAC=normalize (pC-pA)
local vAD=normalize (pD-pA)
local angleBC=dot vAB vAC
local angleBD=dot vAB vAD
if angleBD<angleBC do return false
local angleCD=dot vAC vAD
if angleCD<angleBC do return false
local vBA=-vAB
local vBC=normalize (pC-pB)
local vBD=normalize (pD-pB)
local angleAC=dot vBA vBC
local angleAD=dot vBA vBD
if angleAD<angleAC do return false
local angleCD=dot vBC vBD
if angleCD<angleAC do return false
true
)
fn planelineintersect planePoint planeNormal vectorPoint vectorDir =
(
local planeNormal=normalize planeNormal
local vectorDir=normalize vectorDir
local tol=.0000000754
local vectorPoint2=vectorPoint+vectorDir
local d1=dot (planePoint-vectorPoint) planeNormal
local d2=dot vectorDir planeNormal
local absd1=abs(d1), absd2=abs(d2)
if absd2<tol and absd1>tol do return 0
if absd1<tol and absd2<tol do return -1
(vectorPoint+d1/d2*(vectorPoint2-vectorPoint))
)
fn findFirstSmaller theArray theValue = (
local theArrayCount=theArray.count
for idx=1 to theArrayCount do (if theArray[idx]<theValue do return idx)
0
)--fn
fn getIntersection obj theFaces vRay vRayPos = (
for f in theFaces do (
local faceVerts=(meshop.getVertsUsingFace obj.mesh #{f})as array
pA=polyOp.getVert obj faceVerts[1]
pB=polyOp.getVert obj faceVerts[2]
pC=polyOp.getVert obj faceVerts[3]
local planeNormal=cross (pB-pA) (pC-pA)
local pD=planelineintersect pA planeNormal vRayPos vRay
if (isFacePoint pA pB pC pD) do return (intersectionStruct point:pD meshFace:f)
)
(intersectionStruct())
)--fn
fn getPolyfaceUsingMeshface obj meshFace = (
local mpolyface=meshop.getPolysUsingFace obj.mesh #{meshFace} ignoreVisEdges:false threshhold:360.0
local mpolyfaceVerts=meshop.getVertsUsingFace obj.mesh mpolyface
local polyfaces=polyOp.getFacesUsingVert obj mpolyfaceVerts
local retFace
for face in polyfaces do (
local polyfaceVerts=polyOp.getVertsUsingFace obj face
if ((polyfaceVerts*mpolyfaceVerts).numberset == mpolyfaceVerts.numberset) do retFace=face
)
retFace
)
fn pointLineDist2 pA pB pC = (
local vAB=pB-pA
local vAC=pC-pA
(length (cross vAB vAC))/(length vAB)
)
fn isEdgeVert obj theIntersection theEdge tol = (
local pos=theIntersection.point
local currentVerts=polyOp.getEdgeVerts obj theEdge
local pA=polyOp.getVert obj currentVerts[1]
local pB=polyOp.getVert obj currentVerts[2]
local d=pointLineDist2 pA pB pos
if d<=tol do (
local L1=distance pA pB
local L2=distance pA pos
local L3=distance pB pos
if not ( L2>L1 or L3>L1 ) do (
if not (L2<=tol or L3<=tol)
then ( theIntersection.type=2; theIntersection.edge=theEdge; return theIntersection )
else (
if L2<=tol do ( theIntersection.type=1; theIntersection.vertex=currentVerts[1]; return theIntersection )
if L3<=tol do ( theIntersection.type=1; theIntersection.vertex=currentVerts[2]; return theIntersection )
)
)
)
theIntersection
)
fn findEdgeVert obj theIntersection theEdges tol = (
local pos=theIntersection.point
for thisEdge in theEdges do (
isEdgeVert obj theIntersection thisEdge tol
if (theIntersection.type==1 or theIntersection.type==2) do return theIntersection
)
)
fn intersectPolyRay obj theRay tol:0.000001 = (
local tolDot=0.000001
local vRay=-(normalize theRay.dir)
local vRayPos=theRay.pos
local theFaces=#{1..polyOp.getNumFaces obj}
local theVerts=#{}
-- face visibility test
for f in theFaces do (
local d=dot vRay (polyOp.getFaceNormal obj f)
if d>tolDot do (
local currentVerts=polyOp.getVertsUsingFace obj #{f}
for v in currentVerts do theVerts[v]=true
)
)
-- proxy vertices array
local vertsIdxArray=#()
local vertsDotArray=#()
for v in theVerts do (
local d=dot vRay ( normalize (vRayPos-(polyOp.getVert obj v)) )
local newIdx=findFirstSmaller vertsDotArray d
if newIdx==0
then (append vertsDotArray d ; append vertsIdxArray v)
else (insertItem d vertsDotArray newIdx ; insertItem v vertsIdxArray newIdx)
)
-- test the ray on meshFaces
local testedFaces=#{}
local vertexCount=1
local vertsIdxArrayCount=vertsIdxArray.count
local theIntersection=intersectionStruct()
while (theIntersection.type==0 and vertexCount<=vertsIdxArrayCount) do (
local currentFaces=meshOp.getFacesUsingVert obj.mesh #{vertsIdxArray[vertexCount]}
currentFaces-=testedFaces
theIntersection=getIntersection obj currentFaces vRay vRayPos
if theIntersection.meshFace!=0 do (
theIntersection.type=3
theIntersection.polyFace=getPolyfaceUsingMeshface obj theIntersection.meshFace
local theFaceEdges=polyOp.getEdgesUsingFace obj #{theIntersection.polyFace}
findEdgeVert obj theIntersection theFaceEdges tol
)
testedFaces+=currentFaces
vertexCount+=1
)
theIntersection
)--fn
if (isValidNode $Plane01 and isValidNode $Point01)
then (
obj=$Plane01
if classof obj==Editable_Poly then (
r=ray $Point01.pos [0,0,-1]
local theIntersection=intersectPolyRay obj r tol:0.25
format "theIntersection=% \n" theIntersection
case theIntersection.type of
(
1: (select obj; subobjectlevel=1; polyOp.setVertSelection obj #{theIntersection.vertex})
2: (select obj; subobjectlevel=2; polyOp.setEdgeSelection obj #{theIntersection.edge})
3: (select obj; subobjectlevel=4; polyOp.setFaceSelection obj #{theIntersection.polyFace})
)
)
else messageBox("not an editable poly")
)
else messageBox("rename yours objects")
