Latest Contributions by JakeSkinner

Looks like this may be a bug for the replaceDataSource method and raster datasets.  However, this works with Mosaic Datasets.  Have you considered migrating your large raster mosaics to Mosaic Datasets?  This is a much more efficient way to manage your imagery at ArcGIS 10.  The mosaic dataset makes it easy for you to manage, search, and discover imagery in collections of any size. The mosaic dataset catalogs image and raster data sources, stores detailed metadata, and defines how imagery should be transformed into different products on the fly without the need to preprocess. The mosaic dataset includes: �?� Dynamic mosaicking for efficient handling of overlapping imagery �?� On-the-fly processing to create multiple products from a single source �?� Management and accessibility to metadata �?� Integration for multiple disparate image sensors and formats �?� Streamlined data maintenance and updates of new imagery The mosaic dataset also saves time and storage space.  The mosaic dataset simply reads the images from their native file locations and format.  Since it reads the data from their native file locations it is extremely fast to create mosaic datasets.  Loading large rasters into SDE that took hours will only take seconds now and you will no longer have to allocated the storage space in your SDE database to accomodate these large rasters. ... View more

You can update the layer name by searching for a keyword and then replacing the text.  I was working on a similar script earlier this morning.  I added this functionality to it. import arcpy, os from arcpy import env from arcpy import mapping env.workspace = workspace = arcpy.GetParameterAsText(0) env.overwriteOutput = True for (path, dirs, files) in os.walk(workspace):     for file in files:         if ".mxd" in file.lower():             mxd = os.path.join(path, file)             print mxd + " is being processed"             mxd = mapping.MapDocument(mxd)             for df in mapping.ListDataFrames(mxd, "*"):                 for lyr in mapping.ListLayers(mxd, "*", df):                     lyr.replaceDataSource(arcpy.GetParameterAsText(1), "SDE_WORKSPACE", "")                     print "Successfully updated data sources"                     if "vector" in lyr.name:                         lyr.name = "RASTER.DATALOADER." + lyr.name.split(".")[2]                         print "Succesfully changed layer name"             mxd.save() del mxd The previous database name was "vector", so the feature classes were listed as "vector.dataloader.parcels".  I created and IF statement to find all layer names containing vector and renamed the layers to the new database raster: if "vector" in lyr.name:       lyr.name = "raster.dataloader." + lyr.name.split(".")[2] ... View more

Are you trying to Sum the two lists together based on the index of each list? Ex: listA = [1,2,3,4,3,2,1] listB = [7,8,9,5,4,3,1] index[0] = 7 + 1 = 8 index[1] = 8 + 2 = 10 index[2] = 9 + 3 = 12 If you are, you can use the following code to do this: import os listA = [1,2,3,4,3,2,1] listB = [7,8,9,5,4,3,1] listC = [] x = 0 while x < len(listA):     listC.append(listB + listA )     x += 1                  print listC ... View more

Thanks, Ruth!  The 'replaceDataSource' method was exactly what I needed.  I was able to get this working with the following code: import arcpy from arcpy import env from arcpy import mapping env.workspace = workspace = arcpy.GetParameterAsText(0) env.overwriteOutput = True mxdList = arcpy.ListFiles("*.mxd") for mxd in mxdList:     mxd = workspace + "\\" + mxd     print mxd + " is being processed"     mxd = mapping.MapDocument(mxd)     for df in mapping.ListDataFrames(mxd, "*"):         for lyr in mapping.ListLayers(mxd, "*", df):             lyr.replaceDataSource(arcpy.GetParameterAsText(1), "SDE_WORKSPACE", "")             print "Successfully updated data sources"     mxd.save() del mxd ... View more

You can use the 'Create File Geodatabase' or 'Create Personal Geodatabase' tool to create a 9.3 geodatabase and then copy the data from your 10 geodatabase to the 9.3 geodatabase using ArcGIS 10. ... View more

It looks you are missing a ',' between OBJECTID and SDE for the -C option.  Here is an example on registering a table: SQL> CREATE TABLE building (OBJECTID INTEGER, shape SDE.ST_GEOMETRY); SQL> INSERT INTO building VALUES (1, SDE.ST_Polygon ('polygon ((0 0, 0 10, 10 10, 10 0, 0 0))', 2)); SQL> INSERT INTO building VALUES (2, SDE.ST_Polygon ('polygon ((20 0, 30 20, 40 0, 20 0))', 2)); SQL> INSERT INTO building VALUES (3, SDE.ST_Polygon ('polygon ((20 30, 25 35, 30 30, 20 30))', 2)); SQL> create index bf_indx on building(shape) indextype is SDE.ST_spatial_index parameters = ('SDE.ST_grids=1,0,0 SDE.ST_srid=2'); $ sdelayer -o regiSDE.STer -l building,shape -e nac+ -C OBJECTID,SDE -P HIGH -t SDE.ST_GEOMETRY -i 5151 -u vector -p vector ... View more

Try removing the 417th file and then re-execute the script to see if it works.  The file may be corrupted. ... View more

Below is an example on how you can do this using the Field Calculator.  Be sure to check 'Python' at the top and 'Show Codeblock'.  Under 'Pre-logic Script Code' enter the following: def update():   import arcpy   list = []   rows = arcpy.SearchCursor(r"C:\data\Parcels")   for row in rows:     list.append(row.AREA)   S = sum(list)   del row, rows   rows = arcpy.UpdateCursor(r"C:\data\Parcels")   for row in rows:     row.AREA2 = S     val = row.AREA2   del row, rows   return val Then, under the second block enter: update() ... View more

Below is an example on how to do this.  The code adds all the shapefiles first two characters to a list, then removes the duplicates.  It will then merge the related shapefiles together to one shapefile. lstFCs = arcpy.ListFeatureClasses("*") list1 = [] # Append first two characters of each shapefile to list for fc in lstFCs:     list1.append(fc[0:2]) # Remove duplicates from list list1 = dict.fromkeys(list1) list1 = list1.keys() # Merge feature classes for n in list1:     m = str(n) + "*"     lstFCs = arcpy.ListFeatureClasses(m)     arcpy.Merge_management(lstFCs, n + "_merge") ... View more

In order to update the symbology for layers within MXD using another MXD as the source, you will need to use the 'arcpy.mapping.ListLayers' function rather the the 'arcpy.mapping.Layer'.  Here is an example below on how you can accomplish this: import arcpy mxd = arcpy.mapping.MapDocument(r"C:\temp\Final_Project.mxd") df = arcpy.mapping.ListDataFrames(mxd, "*")[0] lstLayers = arcpy.mapping.ListLayers(mxd, "Rivers*", df) mxd2 = arcpy.mapping.MapDocument(r"C:\temp\Source.mxd") df2 = arcpy.mapping.ListDataFrames(mxd2, "*")[0] for layer in lstLayers:     updateLayer = layer     sourceLayer = arcpy.mapping.ListLayers(mxd2, "River", df2)[0]     arcpy.mapping.UpdateLayer(df, updateLayer, sourceLayer, True)       mxd.save() print "Updated symbology successfully" del mxd You can copy/paste the syntax above to update the Roads and Building layers as well. ... View more

If you have an ArcEditor license for ArcGIS 10 you can use the 'arcpy.FlipLine_edit' function to reverse the direction of your polyline data. ... View more

I would recommend using the arcpy.AddFieldDelimiters function to select the field your using.  That way you can be sure you have the expression correct. Ex: fc = "Roads" Segment_ID = arcpy.AddFieldDelimiters(fc, "Segment_ID") Expression = "Segment_ID = 'C1'" Give that a go and see if you receive the correct results. ... View more

The append should work inside or outside an edit session.  If you perform the append outside of an edit session, make sure the 'Connection Details' within the connection properties of the geodatabase are set to the version you wish to update. ... View more

You can use the arcpy.Append_management command to do this. ... View more