Latest Contributions by JúlioLuta

Edit: It worked. At first I thought it didn't but after I reloaded the file the AngleMin column had the correct values. Thanks! ... View more

Hi, I wanted to share my progress in case it's useful to someone in the future. I also have a question because I can't understand why one step is not working. I calculated the angle of the two feature classes using the following code I found online and which uses the !shape! and the field calculator: # Pre-Loci Script Code import math def GetGeographicalDegrees(shape):   radian = math.atan2(shape.lastpoint.y - shape.firstpoint.y,                       shape.lastpoint.x - shape.firstpoint.x)   radian = radian - (math.pi /2 ) # turn minus 90°   if (radian > 0):      degrees = 360 - ( radian  *  360) / ( 2 * math.pi  )   else:      degrees = 360 - ((2* math.pi + radian  ) * 360) / ( 2 * math.pi  )   return degrees # angle = GetGeographicalDegrees( !SHAPE! )‍‍‍‍‍‍‍‍‍‍‍‍‍‍ Now both line features have an angle field and after the spatial join (join one to many – within a distance of XXm) I created a new field that has the difference between angles (AngleDif): abs(!AngleArcs! - !AngleSegments!)‍ Next using the python console I used code I found as an answer to a simular problem, which needed to find for every field A with the same value the minimum value of field B: import arcpy fc = "C:\\Dropbox\\NewThesis\\databases\\Final.gdb\\ArcosSolSombra_VGA2M_Segments" data = [row for row in arcpy.da.SearchCursor(fc, ["TARGET_FID", "AngleDif"])] import collections ddict = collections.defaultdict(list) for x,y in data:    ddict[x].append(y) minvals = dict((key, min(ddict[key])) for key in ddict.keys()) with arcpy.da.UpdateCursor(fc, ["TARGET_FID", "AngleMin"]) as rows:     for row in rows:         row[1] = minvals[row[0]] ‍‍‍‍‍‍‍‍‍‍‍ The problem I had here is after getting the minvals list the updatecursor is not writing anything in the AngleMin field.  To finish I'll create and calculate a field to identify when AngleMin is equal to AngleDif. I know this solution isn't very elegant but as I mentioned before I didn't practice enough after I learnt a bit of Python and now I barely remember anything. ... View more

Thanks Clinton and Dan! As soon as I try the functions I'll come back with more detailed feedback. Happy Holidays!  ... View more

I'm bumping this thread just to add that I think this result might be caused by the "Drainage enforcement process". At this time I don't have the possibility and time to redo the model and test it. It was unfortunate that I din't notice this detail/option before or, if I did, I didn't correlate it with this "issue". This information might be useful if somebody is noticing the same thing. ... View more

Thanks for the reply. I was getting the warning followed by another saying it couldn't find a route. I thought they were connected but then I discovered that it was a restriction that had a wrong value and was prohibiting several arcs. But it was your reply that made me double check all restrictions. Thankfully everything seems to be working now... phew! Cheers ... View more

This methodology didn't have the results I wanted but with additional steps it might work. With the Intersect 3D Line With Multipatch I can get the lines that touch the Shadow Volume surface (it doesn't output the ones that are inside the Shadow Volume but don't touch the surface. I can cut those lines at the output intersection points. Then I can use Inside 3D to these lines and I'll get the lines inside the Shadow Volume and touch the surface (already cut) and Inside 3D to the original lines to get the ones inside that don't touch the surface. Unfortunately the Create Sun Shadow Volume can't have as secondary input the terrain. I tried to cut the Sun Shadow Volume at the terrain by intersecting it with a extrusion between the terrain and a plane above the buildings to then use multipatch footprint... but I always get errors with the Intersection 3D, always. I even tried with smaller volumes. I could later use the resulting polygons instead of all the lines I'm creating with the above process. ... View more

I just read about the Intersect 3D Line With Multipatch so in the next days I'll try to test the workflow. Creating buildings in AutoCAD -> Sketch Up -> ArcGIS -> Create Sun Shadow Volume -> Intercept with the polylines I have. ... View more

Thanks for the reply. In the Sun Shadow Volume: - To confirm if I understood correctly, I need to use the 3d buildings and the terrain right? - The terrain can be in a raster layer or should I model it in a different feature? To simplify and dividing it into steps, what I need to do: 1a- Create sun polygons and shade polygons from the sun shadow raster/multipatch. or 1b- Intercept polylines with the shade part of the raster/multipatch and do a second intercept with the sun part. and 2- Store in the polylines I already have in which polygon(or polyline group) they are the most. If they are >50% inside the sun polygon I'd store 1, otherwise 0. ... View more

I edited my post adding more screenshots. With that kind of modelling I don't know how much can I trust TopoToRaster. I have counted at least four instances where this "buggy" values appear. ... View more