Latest Contributions by JamesCrandall

Hello, I'm looking to create a legend for a raster-catalog (stored in a personal geodatabase). The catalog contains 50 raster grid data sets. With a Time Layer Animation I visualize the runoff after a rain event in time. The raster catalog displays no default symbology thus when you create a legend it just display's the word "Legend". Has anyone an idea what I have to do? Best Regards Sophie In our case, all of the individual rasters in the catalog were to be symbolized with the same classified renderer properties (ie, the same minimum/maximum range).  So, we determined the min/max for all rasters and set the desired symbology of ONE of the rasters.  Then we exported the class breaks to an .xml and just applied this .xml to the raster catalog loaded in the TOC.  Then it was just a matter of setting the desired color scheme. To have this show up in the legend we just used the single raster we symbolized in the inital step but renamed it to something appropriate. ... View more

i have tried using ADO.NET Entity Data Model in my ArcDesktop project using dockable windows. but every time i start debugging im getting a error. COMException was unhandled by user code Catastrophic failure (Exception from HRESULT: 0x8000FFFF (E_UNEXPECTED)) Private Sub SetupDockableWindow()         If m_dockableWindow Is Nothing Then             Dim dockWindowManager As IDockableWindowManager             dockWindowManager = CType(m_application, IDockableWindowManager)             If Not dockWindowManager Is Nothing Then                 Dim windowID As UID = New UIDClass                 windowID.Value = regKey.GetValue("dwUID") 'DockableWindowGuid                 m_dockableWindow = dockWindowManager.GetDockableWindow(windowID) ' <-----on this line             End If         End If End Sub Add a Try...Catch with some error checking.  Or at minimum, step through and determine the exact line it is failing on. ... View more

Glen -- something that completely slipped my thought process for this thread: What if your Database repository is NOT ArcSDE? This is where and why cx_Oracle is needed for our implementation.  That is, we are integrating non-spatial Oracle attirbute repositories (ancillary to spatial databases), joining this data in the in_memory space to "make it" spatial data, then processing it from there. Hope that makes sense.  Again, if your attributes are stored in SDE then your ArcSDESQLExecute is probably a good choice --- I really have no benchmark to comment on performace differences between this and cx_Oracle.  I'd gather if I moved my SQL into procedural packages on the database server I could really boost things up (but even that may be negligable gains). Anyway -- good luck! ... View more

Not sure if this will help or hurt, but instead of Createfeatureclass_management like you are doing, why not create a table (CreateTable_management) and use your sdeReturn to populate it?  Then use that table in the arcpy.MakeXYEventLayer_management, which a table is the parameter it is asking for. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00170000006z000000 Edit: forgot --- this would mean you don't need cx_Oracle at all.  How is the performance of ArcSDESQLExecute? ... View more

Glen, I have several Geoprocessor implementations that use cx_Oracle to pull in attribute data for processing and it works very well for us.  The datasets that must be processed are large (15 minute sample collections over several years), so since we can issue SQL thru the cx_Oracle we only return the results that are needed over each iteration (date).  That is, instead of having to chug through an entire huge dataset, we simply return what is needed, process it, then move on to the next result that is needed. The trickiest part to get right was making the cx_Oracle cursor (the resulting rows from the SQL) compatible with ESRI-related data.  Specifically, this means we must build an empty data container (table) on the GIS side of things and then populate it with the results in the cx_Oracle cursor.  This of course means it is a field mapping issue, but I think all-in-all it really isn't that difficult to get right.  Once you get the result into an ESRI table format of your choosing, you can then process your XY Event theme as normal. Here is an example of doing this field mapping and populating the in_memory table:          gp.createtable_management("IN_MEMORY", "TempDT", "", "")                     cxRows = cursor.fetchall()  #the cursor here was already filled.          rowcount = cursor.rowcount          for row in range(1, rowcount - 1):   tables = gp.ListTables()   for tbl in tables:        if tbl=="TempDT":                                ### add the fields         for i in range(0, len(cursor.description)):    val1 = str(cursor.description[0])    val2 = str(cursor.description[1])    val3 = str(cursor.description[2])    if val2=="<type 'cx_Oracle.STRING'>":       fldType = "Text"       val3 = cursor.description[2]       gp.AddField(tbl, str(cursor.description[0]), fldType, val3)    if val2=="<type 'cx_Oracle.NATIVE_FLOAT'>":       fldType = "Float"       gp.AddField(tbl, str(cursor.description[0]), fldType)    if val2=="<type 'cx_Oracle.DATETIME'>":       fldType = "Date"       gp.AddField(tbl, str(cursor.description[0]), fldType)                          ### now populate the table      insRows = gp.InsertCursor(tblNew)      for cxRow in cxRows:         insRow = insRows.newRow()         for i in range(0, len(cursor.description)):      insRow.setvalue(str(cursor.description[0]), cxRow)         insRows.insertRow(insRow)          cursor.close()          db.close() ... View more

Faster approach: Copy your FC's to the IN_MEMORY space then peform the join and querying there. http://gis.stackexchange.com/questions/31699/ways-to-speed-up-python-scripts-running-as-arcgis-tools ... View more

AML script merge two SAS cvs. files to an Arc cover file. Aren't coverages going away soon?  You will have to be a bit more descriptive in what was actually being processed from the csv.  What type of features are you generating?  What is in the csv? etc, etc, etc... (typically people post up specific code that they are having difficulty getting working, which as a community we can all understand and contribute ideas on how to work through problems.  Wide, open-ended "how can I" questions most often go unanswered because they require a ton of investigation to resolve.  Be very specificl, show your data, existing code, etc...) ... View more

Hi Is there any possibility to convert AML script to ArcPy script? Tom Anders Engebakken yes.  What does the AML do? ... View more

Oh, this is my bad.  It has been a while since I implemented this, but I incorporated a 3rd party Python library to setup a table with a csv source -- this maybe a bit much for your requirements.  How about if you invoke Table To Table to bring in the csv into a .gdb (like in your scratch workspace or something), then use that table in the MakeXYEventLayer method. Import .csv into gdb: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001200000027000000 MakeXYEventLayer: http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000006z000000 Sorry I don't have a specific example for you, maybe someone else can post up their solution. James ... View more

You can start here, but I am not sure if you will need additional tweaks for the .csv format (this page shows ex with .dbf): http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000006z000000 ... View more

Thanks, Duncan -- I will take a look. Take Care, j ... View more

I just create my own .log file and write lines where needed. #oLoc is just a network path string variable where I want to save this #strD1 and strD2 are just variables of dates set previously logpath = oLoc + "\\" + strD1 + "_to_" + strD2 + ".log" logfile = open(logpath, 'w') logfile.write("Log begin: " + str(log_startdate) + "\n") #here is an example of how I am writting errors to the log file # see if spatial analyst extension is available for use availability = gp.CheckExtension("Spatial") if availability=="Available":    gp.CheckOutExtension("Spatial") else:    logfile.write("ERROR: Spatial Analyst extension not available" + "\n")    return logfile.close() ... View more

Thanks, but my question was really about how to specify the table join. ive tried putting the tables in the input parameter, but it doesnt make sense anyway as you need to tell arcpy how they are related. ive also tried to put sql in the expression parameter but i dont think i got it right. ideally id be able to find an example of where this is already done. Thanks anyway. Post up the section of code you are having an error with. ... View more

I have not tried this myself, so this is just an observation.  But...  Did you attempt to add all of the tables into the in_table parameter? arcpy.MakeQueryTable_management (["rain","geol","tab3","tab4","tab5"], (that's what appears to be allowable from the help reference) http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00170000006r000000 "...The name of the table or tables to be used in the query" Wish I had a difinitive answer for you -- hopefully someone will follow up soon. ... View more