Latest Contributions by MaryBandziukas

OK, thanks.  I am completely new to numpy and pandas.  I'm trying to move away from excel.  This is my first question posted to GeoNet, so thanks for showing me what a great resource it is.  Looking forward to learning what numpy can do (and maybe Pandas) ... View more

so many choices ExtendTable joins an array to a table, I need to join a table (or array) to a feature class.  Would it still work?   Also, the feature class does not have unique values in its join field; if I provide only unique values in the join table, I will have a N:1 join, which should work fine.  e.g. if I have 3 records with the same GenericID (the join field) in the feature class, they will each be matched to the same record in the join table, which has two fields - GID (genericID - the join field) and countGID, which in this case would have the value of 3, and which tells me that this particular segment joined to 3 "valid intersection nodes" as a result of a spatial join.  I noticed that the inputs into arcpy.JoinField_management are supposed to be table views, and I didn't convert my files, just input them as feature class and table.  Perhaps that also needs attention.  I did it because I wanted to see if it would work (and then forgot that part).  Maybe not  Thanks for your continuing attention!  The whole purpose of the script is to avoid having to export to and process in excel and then re-import the spreadsheet back into the gdb.  I hope I can get it; it will make file management a bit easier ... View more

Thanks, Dan, for your reply.  You may have hit it on the head; my path does indeed have spaces and a period, because whoever set up the directories numbered them, such as:  3. Geometric Design\Official Drawings ...  I have no way around this.  The join does work fine in ArcMap, which is how I program - do it first in ArcMap, then write the script.   However, the Join Validation Report gives me a cardinality warning that there are non-distinct matching records, which I am aware of and which doesn't cause errors in the join when done in ArcMap, but maybe in Python I need to present a join table with only unique records.   It is probably pretty obvious that I don't write many programs - I am just discovering that there are nuanced differences between how ArcMap and Python handle things.    Interestingly, the whole beginning of the script works OK using this wksp - spaces, period and all - i.e. read the feature class, place ID numbers into a python list, use .count to find out how many of each ID number there are and put the counts into a second list, zip the lists into a list of tuples, convert the tuples into a nympy structured array and convert the array to a table.  Python will confirm that the table exists, continuing to use this same workspace.  It's just the last step, the join, that hangs up.  Which leads me to believe the non-unique records might be more of an issue than the path.  I noticed another interesting oddity too:  when I first converted the array to a table, I started the table name with upper case letters.  When I did an arcpy.ListTables(), it found other tables with lower case names, but not my new one, or another one I had previously created in ArcMap also beginning with upper case letters.  So I changed the name and took out the upper case characters, and then arcpy found it.  I have a habit of prefacing my gdb tables with GDBT_ (geodatabase table) so I can differentiate them from imported spreadsheets with the same name.   And finally, my feature class filename begins with BK_ (uppercase), and arcpy.ListFeatureClasses() does find it.  So it's not simply a matter of upper/lower case.   In any case, I think I'll try paring my join table to unique records first, and if necessary I may try moving everything to another temporary path without spaces etc, to see if it works there.  Thanks so much for your help ... View more

 Hi, I am Mary Bandziukas, writing from New York City.  I was recently hired by the NYC Dept.of Transportation, where I am on the asset management GIS team.  We perform record-keeping and research functions.  We currently are refining a model to prioritize refurbishment of pavement markings.  I was with this team part time for almost 2 years, during which I began to develop a pedestrian trip model for the agency.  One purpose of the model will be to calculate pedestrian risk of being involved in a ped-vehicle accident.  Using documented trip generation rates, I am calculating trips per lot based upon land use, and I will need to decide how to calculate modal split.  The agency has short-term (2- or 3-hour) ped counts and is beginning to collect some long-term counts which may be useful in validating the model.   The main innovative GIS successes for me over the past few years have been:  - learning about and using the SatScan tool to generate center lines for sidewalk and crosswalk polygons for the purpose of developing a pedestrian network - investigating potential applications of the agency's DEMs generated from aerial LiDAR.  We were interested in the potential to identify road imperfections such as potholes.  I implemented ESRI's Bluespots model in this work.  We will have access to mobil LiDAR in the future, which will likely give much better results. - automating mapping of work orders and other documents with a Python program If I have to identify areas of expertise, I would say I have had an eclectic career so far, with experience in urban planning and higher ed, as well as recently-developed skills in research and with respect specifically to GIS, in map design and cartography.  Long-term, my "dream" would be to implement service-learning projects between students and community groups.  Shorter-term, I want to provide GIS services to community groups on a volunteer basis.  I expect to learn a great deal about community need and technical aspects of GIS in this way. Thank you for implementing this group.  I am looking forward to seeing all the interesting things people are doing around the world.   ... View more