Latest Contributions by TomGeo

Hi, I want to convert some date strings into correct dates, to be recognized by Python and ArcGIS. Setup:      ArcGIS 10.2.1      Python 2.7.5 All data are located within a FGDB... I set the workspace to be the FGDB, created a new Field of type data in my table and went on with the following things within the IDLE Shell import datetime with arcpy.da.UpdateCursor(myTable, [field01, field02]) as UCursor:      for row in UCursor:           origDate = datetime.datetime.strptime(row[0], '%d%b%Y')           row[1] = datetime.date.strftime(origDate, '%Y-%m-%d')           UCursor.updateRow(row) Things work fine until the input string is '11MAY2004', where I get the following error message: Traceback (most recent call last):   File "<pyshell#9>", line 4, in <module>     origDate = datetime.datetime.strptime(row[0], '%d%b%Y')   File "C:\Python27\ArcGISx6410.2\lib\_strptime.py", line 325, in _strptime     (data_string, format)) ValueError: time data '11MAY2004' does not match format '%d%b%Y' I did try the whole thing in plain Python (same shell) before, where I simply did: orig = '11MAY2004' x = datetime.datetime.strptime(orig, '%d%b%Y') and this works without any problem!!! Can somebody enlighten me where I go wrong? Cheers, Thomas ... View more

Hi, is there any way to avoid getting a error 210 when running one instance of ArcMap or ArcCatalog, doing some geo-processing and starting another instance of ArcMap or ArcCatalog? It is very frustrating when your software is opening the Add Data dialog or a Save Data dialog while some model or Python script is working in that folder/database. The result seems obvious...  the model or the script crash, unable to write any longer to the location! Bests Thomas ... View more

Hi, I want to do a routing for an incinerator and have different companies in the supply chain. Of course the amount of waste offered by a single company is exceeding the capacity of the truck. The result of that is that non of the locations with pickup quantities above the truck load is considered for the routing. In reality, it would be necessary that in a fleet of vehicles the waste supplied by one point is shared, or a single vehicle is routed until all the waste is picked up. I there a way to do this? Bests, Thomas ... View more

Hi everybody, while working with Python scripts or with ModelBuilder, at some point I usually end up that I have to write some data to my HDD instead of keeping it in memory. The result might be that I create a couple thousand featureclasses where every single one might hold tens of thousands features. Since those featureclasses still represent only temporary products they need to be deleted in a later step. When I try this in a script I would use 'ListFeatureClasses' with some wildcard and a for-loop to delete the items. That will finally not work, for whatever reason and the script will crash. If I try the same thing in ArcCatalog by selecting the featureclasses and go for the delete button then ArcCatalog will simply get unresponsive and finally the system tells me that the application stopped working. The only way I found so far is to use the Python window from within ArcCatalog or ArcMap, list all the classes of interest, do the for-loop and delete the files listed. --> That's terrible slow!!! I have three questions: How can it be that the same code works in the Python window, but not in an external script? Is there a way to speed up the process of deleting the files? Since ArcSDE is invoked in ArcGIS Desktop, is there a way to use SQL to access the FGDB and delete the data that way? Bests Thomas ... View more

Hi, I have some speed issues with the datasets I am using to find the shortest route. I have a road network, the points of origin and the destination points in one file geodatabase. Points of origin I have 2500 and destination points I have a bit more than 74.000! The script I wrote works but process is really slow and during the process my hard drive is getting filled up with some temporary data. In the following there is the script I'm using to process the data and it would be great if some of you could have a look at it and point out how to speed up process. import arcpy, gc arcpy.env.workspace = 'D:/Projects/gis/processed/Network/Distance.gdb/' arcpy.env.overwriteOutput = True arcpy.CheckOutExtension('Network') network = 'DM/DMnetwork_ND' origin = 'data/Start_N2500' destination = arcpy.MakeFeatureLayer_management('data/rIntersect_0Tolerance', 'in_memory') fail = [] okay = [] OD_Cost_Matrix = arcpy.MakeODCostMatrixLayer_na(network, 'OD Cost Matrix', 'Length', '' , '', '',                                                 'ALLOW_UTURNS', '', 'NO_HIERARCHY', '', 'STRAIGHT_LINES', '') arcpy.na.AddLocations(OD_Cost_Matrix, 'Destinations', destination, 'Name Areas #', '5000 Meters', '',                       'DMnetwork SHAPE;DMnetwork_ND_Junctions NONE', 'MATCH_TO_CLOSEST', 'APPEND',                       'NO_SNAP', '5 Meters', 'INCLUDE', 'DMnetwork #;DMnetwork_ND_Junctions #') arcpy.Delete_management(destination) for row in arcpy.da.SearchCursor(origin, ['Respondent']):     expr = 'Respondent = ' + str(int(row[0]))     arcpy.na.AddLocations(OD_Cost_Matrix, 'Origins', origin, 'Name Respondent #', '5000 Meters', '',                           'DMnetwork SHAPE;DMnetwork_ND_Junctions NONE', 'MATCH_TO_CLOSEST', 'CLEAR',                           'NO_SNAP', '5 Meters', 'INCLUDE', 'DMnetwork #;DMnetwork_ND_Junctions #')         try:         arcpy.Solve_na(OD_Cost_Matrix, 'SKIP', 'TERMINATE', '')                 output = 'Min_Output_' + str(int(row[0]))                 tmpLyr = arcpy.MakeFeatureLayer_management('OD Cost Matrix/Lines', str(int(row[0])), '', '', 'ObjectID ObjectID HIDDEN NONE;Shape Shape HIDDEN NONE;Name Name VISIBLE NONE;OriginID OriginID VISIBLE NONE;DestinationID DestinationID VISIBLE NONE;DestinationRank DestinationRank HIDDEN NONE;Total_Length Total_Length VISIBLE NONE')                 arcpy.Statistics_analysis(tmpLyr, output, 'Total_Length MIN', 'OriginID;DestinationID')                 arcpy.Delete_management(tmpLyr)         okay.append(output)         print str(int(row[0])) + ' processed'     except:         fail.append(int(row[0])) arcpy.Merge_management(okay, 'finaloutput') arcpy.Delete_management(origin) arcpy.Delete_management(tmpLyr) for item in okay:     arcpy.Delete_management(item) gc.collect() As I mentioned my hard drive is getting flooded and it all happens in 'C:\Users\TomGeo\AppData\Local\Temp'. I'm talking about more than 100GB that are written to the directory and I wonder why such a data amount is created. Thanks in advance for your comments! Bests Thomas ... View more

Dear all, I'm trying to install Geoportal on a Linux machine following manual. It is written at page six " Create the user like you created the geoportal user, but do not make it a member of the wheel or sudo group." However, when trying to set up Geoportal user and Geoportal schema in PostgreSQL (page 10) then I'm asked to give the password for postgres... [sudo] password for postgres . The password I assigned to the user postgres when creating it as written at page six will not work at this point and I'm a bit confused about the "[sudo]" when being asked for the password. Do I have to make postgres member of the sudo group or how do I overcome that problem? Also, can you confirm that I have to change some lines in grants_linuxpg.sh and create_schema_linuxpg.sh since they otherwise come back with some error about and unexpected "("? I'm working with Linux Mint 14 (Debian related in the context to the manual linked above). Cheers Thomas ... View more