Latest Contributions by XanderBakker

Hi Kevin, I hate to break the news to you, but there is a lot more wrong with the script then what Jim is mentioning. I'll mention a few points: You access the environment setting "env" before they have been initialized ("from arcpy import env" is missing) You set the workspace property of the environment setting using a variable "Input_Model", which is not available at that point, since it is obtained from the first parameter later in the script You use "gp.overwrite", but gp is not set (you should use the "env" for this purpose and set it properly). The overwrite property does not exist this should be overwriteOutput You check out a Spatial Analyst extension, but the tools you're using do not require the sa extension. You load the Data Management Tools toolbox, but this is really not necessary (was done in previous versions) You set a variable wse_grid equal to a hard coded string including " + Count", when you probably want to add the count to the string As mentioned before your first argument is read after its first use "Input_Model = arcpy.GetParameterAsText(0)" should move to the top You preform an "arcpy.ListRasters". This works on the current workspace. To do this for each folder, you should set the env.workspace for each folder Next you  loop through the rasters with the statement for "wsgridp004" in rasters:  This will not work. Normally you would use something like for raster in rasters: and then access the raster object in the following steps In the arcpy.Clipmanagement statement, you specify strings, when you should specify the variables. arcpy.Clipmanagement does not exits, it should be arcpy.Clip_management "FP_Poly" is specified as hard coded string, but should be specified as variable (without double quotes) Next you add a variable to the array with "array.add", but this variable is never set correctly (the initial wrong value is assigned multiple times) In the arcpy.MosaicToNewRaster_management you again specify hard coded strings, when you should specify variables (for instance Input_Model should be without quotes. The string "array.next" should be a list of rasters separated by semicolons. Since you obtain it from different workspaces, you should include the path for each raster To get you started I've included some code. The code however has not been tested... import arcpy import os from arcpy import env env.overwriteOutput = True  # Argument 1 is the GeoRAS Directory Input_Model = arcpy.GetParameterAsText(0) # Argument 2 is the 100-Year Floodplain FP_Poly = arcpy.GetParameterAsText(1)   Count = 0 lstRas = '' rasName = 'wsgridp004'  # set workspace for list folders env.workspace = Input_Model Folders = arcpy.ListWorspaces("*", "Folder")  for Folder in Folders:     env.workspace = Folder     #Find the wsgridp004 Raster     rasters = arcpy.ListRasters(rasName , "GRID")     for raster in rasters:         Count += 1         wse_grid = Input_Model + os.sep + "{0}{1}".format(rasName ,Count)          #Clip Raster with 100-Yr Floodplain Polygon         arcpy.Clip_management(raster, "#", wse_grid, FP_Poly, "0", "ClippingGeometry", "NO_MAINTAIN_EXTENT")         print arcpy.GetMessages()         if lstRas == '':             lstRas = wse_grid         else:             lstRas = lstRas + ";" + wse_grid  #Mosaics Clipped 100-Yr Floodplain Grids print "Mosaicing Clipped Water Surface Grids" arcpy.MosaicToNewRaster_management(lstRas, Input_Model, "wse_100_Yr","NAD_1983_StatePlane_Texas_Central_FIPS_4203_Feet", "32_BIT_FLOAT", "3", "1") print arcpy.GetMessages() print "Complete" Assuming from the line where you load the toolbox that you are using ArcGIS 10.0 you should look into this topic: A quick tour of ArcPy The training.esri.com site offers some Python/arcoy courses that might be interesting: http://training.esri.com/gateway/index.cfm?fa=search.results&searchterm=arcpy&search=Search&CourseTypeID=1&ArcGISVersion=10.0 Kind regards, Xander ... View more

Hi Tom, When accessing a layer you will not have information if the layer is selected or not. A map document does not need to be open to be accessed by Python code. That's why there is no "selected" property on layers. Instead you may want to use the "visible" property of the Layer (arcpy.mapping) object. This enable you to detect whether the layer is switched on or off. Layer.visible: Controls the display of a layer. This has the same effect as checking the check box next to the layer in the table of contents in ArcMap. If set to True, the layer will draw; if set to False, the layer will not be drawn. Not all layers support the visible property (for example, restricted web service layers), so it is good practice to test for this ahead of time using the supports method. Kind regards, Xander Bakker ... View more