ray collider backface result

On 29/04/2014 at 01:20, xxxxxxxx wrote:

Using the ray collider, I get some results that I don't understand.
The order of face/backface changes when the ray is running through 4 polygons in comparison to 2 polygons.
Any ideas?

Thanks,
M

  
import c4d  
  
def SetGlobalPosition(obj, pos) :  
  m = obj.GetMg()  
  m.off = pos  
  obj.SetMg(m)  
  
def main() :  
  if op==None: return None  
  doc.StartUndo()  
  doc.AddUndo(c4d.UNDOTYPE_CHANGE, op)  
    
  HitList=[]  
    
  boundingbox = op.GetRad()*2  
  dim=boundingbox.x/10  
    
  collider = c4d.utils.GeRayCollider()   
  collider.Init(op,True)     
    
  offsetZ=op.GetMp().z-boundingbox.z/2-0.1   
  y= -50      
  inter = collider.Intersect(c4d.Vector(0,y,offsetZ),c4d.Vector(0,0,1), boundingbox.z+0.2)  
  print y  
  if inter ==True:  
      count= collider.GetIntersectionCount()  
      c=0  
        
      while c<count:  
            
            
          collpos = collider.GetIntersection(c)["hitpos"]  
            
            
          if collpos.z not in HitList:  
                                
              back = collider.GetIntersection(c)["backface"]   
              print back  
                
              target = c4d.BaseObject(c4d.Onull)  
              SetGlobalPosition(target,collpos)  
              target.InsertUnder(op)   
              target.Message(c4d.MSG_UPDATE)    
              target.SetBit(c4d.BIT_ACTIVE)        
              HitList.append(collpos.z)  
              print c,"HitCount"  
          c+=1  
            
  op.DelBit(c4d.BIT_ACTIVE)     
  c4d.EventAdd()  
if __name__=='__main__':  
  main()  

On 30/04/2014 at 08:31, xxxxxxxx wrote:

Hi all

It's me again. Reading other posts I have seen that most questions are answered basically on the spot, what makes me wonder about mine.
Maybe I just don't see the wood for the trees and my described problem is just no problem? Or is it the other way around, and it is really a hard nut to crack?
Or, reading my text once again, I kind of barged in on you without even saying "Hello" to everybody? If so, sorry and...
Hello to all of you
M

On 30/04/2014 at 09:11, xxxxxxxx wrote:

Don't read much into it, I think it just depends on who's online and has time to check it out.

In the meantime, I don't know, ray colliders always seem to act strangely around me, too. oO

On 30/04/2014 at 10:09, xxxxxxxx wrote:

I personally answer questions based two things:
1. How much I have in my notes about it.
2. How clearly the person asks the question.

In this case. I have no idea what you're asking about.
I have no idea what 2 polygons vs. 4 polygons means.

If you don't get a reply soon. Then I would provide more details about what it is you're doing.
And remember that we can't read your mind. And we are all very busy doing our own projects too. So try to be very specific and spell out exactly what your doing. As if you are explaining it to your mother who knows nothing about C4D.
If you do that. You will probably get faster replies.

Providing a scene file is also sometimes helpful too.

-ScottA

On 01/05/2014 at 03:14, xxxxxxxx wrote:

Hi,

to ScottA:
picking up your idea about talking to my mum about the issue, I have to say that she would probably interrupt my wonderful speech and question why I have never studied something decent

Anyway, while actually trying to explain and illustrate the problem in a simple way, I stumbled across the following:

https://www.dropbox.com/sh/ity9gya86a5tht3/QgBALH5amP

As you can see in the picture attached, I send rays along the z axis through an object.
Using the While Loop I'd like to receive unique z positions. Cause if that ray hits a point or an edge directly, I would get 4 respectively 2 identical z positions.
So far I figured out that you have to decide between 3 different cases so as to get the right hit order.

unique hit count = uhc

1 uhc is odd->volume is not closed or we hit the max object expansion in Y or X direction( I didn´t go further with this yet)
2 uhc is even and uhc/2 is even-> face/backface result is in the right order
3 uhc is even and uhc/2 is odd->you have to reverse face/backface results

I have absolutely no clue whatsoever why it behaves like that.
I hope this helps to explain that weird ray behaviour in different issues.
I´ll keep you updated

M

On 01/05/2014 at 09:56, xxxxxxxx wrote:

That helps me understand what you're doing better.
I'm not a big ray user. So I could be wrong about this.
But AFAIK. The GetIntersectionCount() function returns hits for a polygon's front & back facing faces. Which can be a bit tricky when looping through the hit polygons.
You will get double the expected results if you don't skip over the backfacing polygons in your loop counter code.

This is a simpler example that just sends a ray down the Z-axis in the scene. And looks for only the front facing polygon collisions.
It's returning the expected result for me. Even if I add/delete more polygons to the object that the ray is colliding with:

import c4d  
def main() :  
  
  #Make sure there is an active object  
  if op==None: return None  
    
  collider = c4d.utils.GeRayCollider()   
  collider.Init(op,True)  
  
  rayStartPos = c4d.Vector(0, 150, -500) #The origin of the ray  
  rayDirection = c4d.Vector(0,0,1)            #Point the ray down the Z axis  
  rayLength = 10000  
  
  inter = collider.Intersect(rayStartPos, rayDirection, rayLength, False)  
  
  if inter == True:  
        
      #Get the number of times the ray hits polygons  
      #IMPORATNT: Each polygon is counted twice here(frontFace and BackFace)  
      count = collider.GetIntersectionCount()  
  
      collisions = 0          
      while collisions<count:  
            
          polygonHit = collider.GetIntersection(collisions)["hitpos"]  
            
          if polygonHit:  
              print "frontFacing: ", polygonHit  
                
              #Skip counting the backfacing polygons in our counter              
              collisions +=2  
                
          print "HitCount ", collisions/2  
  
  c4d.EventAdd()  
  
if __name__=='__main__':  
  main()

Not sure if that helps you or not.

-ScottA

On 01/05/2014 at 11:02, xxxxxxxx wrote:

Thanks a lot for your help ScottA,

unfortunately that isn´t a solution for the aim to get unique collision points on a surface and preserve the right face/backface order of collisions.

If you test your script with a surface and the ray hits a point/vertex of the mesh directly, you will usually receive 8 collision vectors (cause this vertex is shared by 4 polygons) and with your script you get 4, cause you´re only dividing hits by 2.
As I mentioned before, if the ray hits an edge (cause this edge is shared by 2 polygons), you will usually receive 4 collision vectors and with your script you get 2.

.....???, a nut
M

On 01/05/2014 at 13:12, xxxxxxxx wrote:

Ouch! I see what you mean.
If the ray hits polygon faces it works as expected.
But If the ray hits vertices instead of the faces. It seems to do weird things.

I changed part of my example to add little cubes where the rays collide with the polygons to make it easier to see:

            if polygonHit:  
              print "frontFacing: ", polygonHit  
              cube = c4d.BaseObject(c4d.Ocube)  
              cube.SetAbsScale(c4d.Vector(0.1,0.1,0.1))  
              cube.SetAbsPos(polygonHit)  
              doc.InsertObject(cube)  

I see it. But I don't have an answer for handling it.
All I can think of is maybe getting the mid point of the polygons. Then check if the ray is a certain distance from that mid point.
And if it isn't. Then you do some sort of "don't use the ray values in this case" in your code.

Someone else might know a better (easier) way to handle it.

-ScottA

On 02/05/2014 at 04:47, xxxxxxxx wrote:

Maybe i didn´t get  exactly what you are trying to achieve with a polygon mid point calculation, but if the ray hits a vertex directly and all attached polygons have the same size, and therefore the same distance, there is no way to decide.

I have updated the link above.
Now you will find a c4d file with two polygon objects to test ray results.
If you do so, the Y-coordinates for your

ayStartPos = c4d.Vector(0, 150, -500) #The origin of the ray

could be -400,-200,0,200,400 to get the interesting results.

You will recognize some rounding errors at the 400 Y-coordinate at the ray check_Y_expansion object
So there is also a problem with the accuracy of ray collision module.

Still trying to find a solution...
M

On 02/05/2014 at 10:31, xxxxxxxx wrote:

Hi all,

almost a solution
I´m not satisfied with it, but it's a step forward.

It´s  n^2 writing the collision positions to a list and sort it afterwards, but anyway.
There is still a problem caused by the way computers calculate real values, I guess??
That´s why it's possible that if a ray hits a single vertex it can calculate 2 position points for it,
for example 500.0000001 and 499.9999999.
I tried to minimize the errors by rounding the Hitpositions for comparison, but store the exact values in a position list.
Now, only if the hitted vertex has an exact value off  "value".5 you will get wrong results.

The other errors were eliminated.

Still searching for a solution....
Is there a way to sort a list and delete double items in one go??
Does anyone have another idea ???
Or can confirm my results???

Thanks in advance
martin

  
import c4d  
  
def main() :  
    
  #Make sure there is an active object  
  if op==None: return None  
    
  #Add an undo sequence  
  doc.StartUndo()  
  doc.AddUndo(c4d.UNDOTYPE_CHANGE, op)  
    
  HitList=[]  
  PosList=[]  
    
  #Initialize the collider  
  collider = c4d.utils.GeRayCollider()   
  collider.Init(op,True)   
    
  #Set the origin, direction and expansion of the ray that the whole bounding box is hit by the ray  
  #Be aware that the ray collider works in local space, therefore reset your axis to the bouningbox center  
  boundingbox = op.GetRad()*2  
  rayZposition=op.GetMp().z-boundingbox.z/2-0.1   
  rayXposition=0  
  rayYposition=400  
    
  rayposition=c4d.Vector(rayXposition,rayYposition,rayZposition)  
  rayzdirection=c4d.Vector(0,0,1)  
  raylength=boundingbox.z+0.2  
    
  #Start the ray  
  inter = collider.Intersect(rayposition,rayzdirection,raylength)  
  if inter ==True:  
      count= collider.GetIntersectionCount()  
      collision=0  
        
      #Loop through all collisions   
      while collision<count:  
            
          hitposition = collider.GetIntersection(collision)["hitpos"]  
            
          #Limit the accuracy of z positions in the HitList due to the rounding errors with real type calculation **  
           **#Therefore we store the rounded z position to be compared in the HitList  
          #but preserve the real collisionposition and store it in the PosList  
            
          if round(hitposition.z) not in HitList:  
                                
              backface = collider.GetIntersection(collision)["backface"]   
                      
              HitList.append(round(hitposition.z))  
              PosList.append([hitposition,backface])  
                
          collision+=1  
            
            
  #Sort the PosList in z direction  
  PosList=sorted(PosList,key=lambda x: x[0][2])    
  print PosList   
    
  #Place cubes on every unique collision position     
  for i,position in enumerate(PosList) :  
      cube = c4d.BaseObject(c4d.Ocube)  
      cube.SetAbsScale(c4d.Vector(0.1,0.1,0.1))  
      cube.SetAbsPos(position[0])  
        
      cube.SetName(str(i)+"___"+str(position[1]))  
      doc.InsertObject(cube)  
      doc.AddUndo(c4d.UNDOTYPE_NEW, cube)  
       
  c4d.EventAdd()  
  doc.EndUndo()  
if __name__=='__main__':  
  main()  

On 06/05/2014 at 00:55, xxxxxxxx wrote:

Originally posted by xxxxxxxx

Is there a way to sort a list and delete double items in one go??

Python lists have a built-in sort() method that modifies the list in-place and a sorted() built-in function that builds a new sorted list from an iterable.

The common approach to get a unique collection of items is to use a set. Sets are  unordered collections of  distinct  objects. To create a set from any iterable, you can simply pass it to the built-in set() function. If you later need a real list again, you can similarly pass the set to the list()function.
The following example should cover whatever you are trying to do:

>>> t = [1, 2, 3, 1, 2, 5, 6, 7, 8]
>>> t
[1, 2, 3, 1, 2, 5, 6, 7, 8]
>>> list(set(t))
[1, 2, 3, 5, 6, 7, 8]

On 06/05/2014 at 01:16, xxxxxxxx wrote:

Thanks for your reply,

I like your blog.
I knew this one and I should have been more specific.
I wanted to sort a vector list and vectors are not hashable, therefore set() won´t work in this case, I guess.

martin

On 13/05/2014 at 06:31, xxxxxxxx wrote:

finally....

to sort a list with tuples of mixed values, including float values, this is a solution:

  
import c4d  
from c4d import utils  
  
def SortFloatUniq(inList) :  
  
  outList= []  
  counter= 0  
  while counter< len(inList) :  
        
      if c4d.utils.CompareFloatTolerant(inList[counter][0][2], inList[counter-1][0][2])== False:  
        
          outList.append(inList[counter])  
      counter+= 1    
  outList= sorted(outList,key= lambda x: x[0][2])    
  return outList  
  
def main() :  
    
  Listtest= [(c4d.Vector(1.0,1.0,8.0),False),(c4d.Vector(1.0,1.0,5.0),False),(c4d.Vector(1.0,1.0,5.0),False),(c4d.Vector(1.0,1.0,6.0),False),(c4d.Vector(1.0,1.0,1.0),False)]  
  Listtest= SortFloatUniq(Listtest)  
  print Listtest  
if __name__=='__main__':  
  main()  

Using the ray collider to get unique collision points on a surface and preserve the right face/backface order of collisions,
you have to limit the precision of the collider like this: (code line 48)

  
import c4d  
  
def main() :  
    
  #Make sure there is an active object  
  if op==None: return None  
    
  #add an undo sequence  
  doc.StartUndo()  
  doc.AddUndo(c4d.UNDOTYPE_CHANGE, op)  
    
  HitList=[]  
  PosList=[]  
  matr=op.GetMg()  
    
  #initialize the collider  
  collider = c4d.utils.GeRayCollider()   
  collider.Init(op,True)   
    
  #set the origin, direction and expansion of the ray that the whole bounding box is hit by the ray  
  #be aware that the ray collider works in local space so reset your axis to bouningbox center  
  boundingbox = op.GetRad()*2  
  rayZposition=op.GetMp().z-boundingbox.z/2-0.1   
  rayXposition=0  
  rayYposition=400  
    
  rayposition=c4d.Vector(rayXposition,rayYposition,rayZposition)  
  rayzdirection=c4d.Vector(0,0,1)  
  raylength=boundingbox.z+0.2  
    
  #start the ray  
  precision = 6  
  inter = collider.Intersect(rayposition,rayzdirection,raylength)  
  if inter ==True:  
      count= collider.GetIntersectionCount()  
      collision=0  
        
      #loop through all collisions   
      while collision<count:  
            
          hitposition = collider.GetIntersection(collision)["hitpos"]  
          print hitposition  
            
          #limit the accuracy of zpositions in the HitList with a precision value because of the rounding errors with real type calculation  
          #store the rounded zposition to compare in the HitList  
          #but preserve the real collisionposition and store it in the PosList  
            
          if round(hitposition.z,precision) not in HitList:  
                                
              backface = collider.GetIntersection(collision)["backface"]   
              PosList.append([hitposition,backface])        
              HitList.append(round(hitposition.z,precision))  
                
                
          collision+=1  
            
            
  #sort the PosList in Zdirection  
  PosList=sorted(PosList,key=lambda x: x[0][2])    
  print PosList   
    
  #place cubes on every unique collision position     
  for i,position in enumerate(PosList) :  
      cube = c4d.BaseObject(c4d.Ocube)  
      cube.SetAbsScale(c4d.Vector(0.1,0.1,0.1))  
      cube.SetAbsPos(position[0]*matr)  
        
      cube.SetName(str(i)+"___"+str(position[1]))  
      doc.InsertObject(cube)  
      doc.AddUndo(c4d.UNDOTYPE_NEW, cube)  
       
  c4d.EventAdd()  
  doc.EndUndo()  
if __name__=='__main__':  
  main()  
  

Cheers,
Martin