Latest Contributions by DavidWasserman

Wow. This is really impressive. Is this all procedural? Where is the rule that made the bike intersection? To answer your question, it seems like your case logic is set up to do either a median or a bike lanes or something else.  Case statements in CGA are not like if statements in python/other programming languages, they are the flow of decisions down the shape hierarchy, they do not execute for every possible true statement, it takes the first true statement and goes down that path. I just want to clarify that just in case it is helpful (might not be).  To put it close to the median, you would need to maybe use some type of binary switch patter. A good example of this is how I place bus lanes in different locations in the complete street rule, with one of the configurations being adjacent to the median.  See _Transit_Lane_Width_Switch() here: Complete_Street_Rule/Complete_Street.cga at master · Holisticnature/Complete_Street_Rule · GitHub  ... View more

Wow that is pretty astounding. So it sounds like someone would have to use data driven pages to do something similar to what I am thinking based on a time field. I am not sure however that can be done with Mosaic datasets however. It might be some label hack.  I will add it to the ideas page later. Thanks Robert.  ... View more

Are you trying to create a new rule based on assets and features in the complete street rule?  https://community.esri.com/docs/DOC-6915  Can you provide more context on what you are trying to do? Do you want bike lanes closer to the median?  ... View more

Only local, though I got the tool to work if I added them a few at a time (which was painful). It starts at tool initialization, not even run time (the second the workspaces are defined). It's like one of the script initialization parameters consumes a lot of time or something. They are 32 bit kernel densities from different time periods.  ... View more

Just to add to this, I am using Pro 1.3 and ArcMap 10.4.1. 10.4.1 the workflow I am doing with mosaics works (though only when processing extent and raster mask environments are nulled), but in Pro the Add Rasters tool is stuck on a loop during the stage at which you are just selecting which datasets to pull in. The computer I am running handles decent scenes in CityEngine (2 G Nvidia, 16 GB ram, decent processor), so it is likely not the issue. OS is Microsoft 10. Just having similar performance issues. I will explore the steps to resolve below.  ... View more

I will give this a try next week when I can. Thanks for the code! ... View more

Hey Joshua, Sorry for the confusion. As some background I am familiar with working with ArcGIS server in the desktop environment, but not necessarily in python directly. I am imagining that in order to connect to an SDE that is using SQLServer, I need to use the Create DB Connection tool (using 10.3 or later correct?). The issue is I am assuming in only know the instance name, but not the name of the SDE which is also required by the tool. The goal is to run a compression tool on all SDEs encountered in the instance essentially. A script I am thinking about would look like this. # Name: CreateDatabaseConnections.py # Description: Will use a server instance and authentication to derive SDEs connections that will then have the Compress tool ran on them as # as part of on going server maintance. # Import system modules import arcpy def find_sdes_from_server_instance(server_instance):      "Function will look at all SDEs on a server instance and return their names as a list"      return ["Roads.sde","Bikes.sde","BuiltEnv.sde", SomeOther.sde] # Example return sde_name_list=find_sdes_from_server_instance("MYSQL_SERVER") for sde in sde_name_list:      arcpy.CreateDatabaseConnection_management("Database Connections", sde, "SQL_SERVER", "utah", "DATABASE_AUTH", "gdb", "gdb", "SAVE_USERNAME", "garfield", "#", "TRANSACTIONAL", "sde.DEFAULT")      .... more stuff      arcpy.Compress_management(connection_file) clean up.... ... View more

Hey Joshua, I guess I do want the files or specifically the SDE names, but I guess I want to spontaneously create the SDE connection files based on the server instance. I think were were interested in one server in particular, but I understand it is a risky proposition generally. I will give walk a try then and report back. I guess do the SDE connection files have to exist before you look for them (IE have you had to connect to the SDE before). My assumption is the script will only know the server instance to crawl. It might be what I want already exists, so I think I will give the walk function more scrutiny in this example. ... View more

Generally their subdivision function for shapes generally is one of better ones in the Land Use/GIS world which is still using fishnets to cut land  and then using eliminates to remove slivers in some cases. Recently in 2016 they also made a split area function that might be worth checking out, but I have had performance issues with CE 2016/rule interface bugs. Generally, the best example I have seen in terms of application was a lightning talk a few years ago. His python scripts are pretty simple, but effective, if above the certain size, he subdivides his parcels. Generally this function is accessed in the Python API under subdivideShapes. Example from docs: subdivideShapes(self, shapes, settings=None) Subdivides a set of shapes using a block to lot algorithm. @param shapes: The set of shapes to subdivide. @param settings: The subdivision settings. Omit to use default settings. (default  = None). @return: The new lot shapes. [sequence] @note: # subdivide shapes in scene shapes = ce.getObjectsFrom(ce.scene, ce.isShape, ce.withName("'Block 1'")) subdsettings = SubdivideShapesSettings() subdsettings.setLotAreaMax(800) ce.subdivideShapes(shapes, subdsettings) ... View more

Generally, to get coordinates of the shape, you could use the getGeoCoordinate function in CGA, put it into a report, then have python pump out the reports like you see in the tutorial (theoretically). This can be done on a shape basis. As a quick experiment I put this under the BikeLane rule. If used this way, I suggest using this function on shapes that are subdivided alone multiple shapes with scopes that are narrow and easy to understand. I have not used this function extensively, but theoretically it could be used to help guide exports for more complicated rendering programs and some people have suggested transportation modeling programs like VisSim. This function has the limitations noted in CityEngine Coordinate systems, and to the CGA function itself. BikeLane(sideflip,Width)-->   print("Bike Lane Shape Coordinates:"+(str(getGeoCoord(lat))+","+ str(getGeoCoord(lon))))   report("Coordinates:Long",float(getGeoCoord(lat)*1000000)) Console Output was:      Bike Lane Shape Coordinates:0.0003011754172411,-0.0017071481280305 Report in rule: ... View more

This is exactly what I am looking for. I will test it out when I get a chance. Thanks! ... View more

This has been my experience as well. I have had conflicts between my anaconda environment and Melinda Morang's tools before (might be a SQLite conflict). Edit: On a second computer where I installed Pro, I also found that a 2.7 Anaconda distribution can also interfere with installing 3rd Party Packages. If the details matter to anyone, the error was usually related to the numpy library not be recognized or found in conda. It would often suggest doing a search with binstar to find the library.  ... View more

Thanks, I think I will just work on it at home then. I am sure I can list them out and get something written pretty quickly. ... View more