Latest Contributions by PeterWilson

I've writen a Python script that resamples ECW images and at the same time converts the tiles to TiFF images. When I try to run the following script I receive the following error message: Processing COCT_W34A_23_8cm_2015 raster Unhandled Exception: System.IO.FileNotFoundException: Could not load file or assembly 'GpMetadataFunctions, Version=10.3.0.0, Culture=neutral, PublicKeyToken=8fc3cc631e44ad86' or one of its dependencies. The system cannot find the file specified. File name: 'GpMetadataFunctions, Version=10.3.0.0, Culture=neutral, PublicKeyToken=8fc3cc631e44ad86'    at System.Reflection.Assembly._nLoad(AssemblyName fileName, String codeBase, Evidence assemblySecurity, Assembly locationHint, StackCrawlMark& stackMark, Boolean throwOnFileNotFound, Boolean forIntrospection)    at System.Reflection.Assembly.InternalLoad(AssemblyName assemblyRef, Evidence assemblySecurity, StackCrawlMark& stackMark, Boolean forIntrospection)    at System.Reflection.Assembly.InternalLoad(String assemblyString, Evidence assemblySecurity, StackCrawlMark& stackMark, Boolean forIntrospection)    at System.Reflection.Assembly.Load(String assemblyString)    at GetManagedType(Char* _assembly, Char* _path, Char* _className) WRN: Assembly binding logging is turned OFF. To enable assembly bind failure logging, set the registry value [HKLM\Software\Microsoft\Fusion!EnableLog] (DWORD) to 1. Note: There is some performance penalty associated with assembly bind failure logging. To turn this feature off, remove the registry value [HKLM\Software\Microsoft\Fusion!EnableLog]. Error Message ''' Created on 05 Aug 2016 Resample and Convert ECW images to TiFF @author: PeterW ''' # import site packages and modules import arcpy from pathlib import Path # set arguments index_layer = r"E:\Projects\2016\G112502\lyr\SubcatchA.lyr" output_folder = r"E:\Projects\2016\G112502\CoCT_Imagery\SubCatchA\Tiles" # set environment variables arcpy.env.overwriteOutput = True # create single image from tiles based on index selection def tiles_to_mosaic(index_layer, output_folder):     with arcpy.da.SearchCursor(index_layer, ["RasterPath"]) as scur:  # @UndefinedVariable         for row in scur:             raster_name = Path(row[0]).stem             arcpy.AddMessage("Processing {0} raster".format(raster_name))             out_raster = "{0}\\{1}.tif".format(output_folder, raster_name)             arcpy.Resample_management(row[0], out_raster, cell_size = 0.32, resampling_type = "NEAREST") tiles_to_mosaic(index_layer, output_folder) Python  Script I've even go so far as to set the input arguments to a single ECW image and the output to the output TiFF image, yet I still receive the same error message. I'm able to run the Resample Tool from ArcMap or ArcCatalog with no error messages. I've tested the following on TiFF images and the script runs successfully. Any advice in how to resolve the following will be appreciated. ... View more

I've added a new field "CatchAreaKm2" (Double) to a File Geodatabase Table. I then create a list of fields of the table that will be used in an UpdateCursor. The problem that I have is that the newly added field is listed within the list of fields yet when I iterate through a sliced version of the fields (landuse) it doesn't recognize the newly added "CatchAreaKm2" field and I can't understand for the life of me why not. In order the get a list of the Landuse fields the index slice should be row[1:-1] yet if I use the following it skips the last landuse field ignoring the newly added field. in order to get the following to work I need to using the index slice row[1:] Any help solving why the newly field is not being recognized will be appreciated. [u'HydroID', u'Commercial_Forestry', u'Cultivated', u'Indigenous', u'Mines', u'Natural_Vegetation_Forest', u'Urban', u'Waterbodies', u'Wetlands', u'CatchAreaKm2'] Python: landuse_fields # calculate percentage of landuse per watershed def landuse_percentage(watershed, output_pivot):     watershed_fields = ["HydroID", "Shape_Area"]     valuedict = {r[0]: (r[1:]) for r in arcpy.da.SearchCursor(watershed, watershed_fields)}  # @UndefinedVariable     arcpy.AddField_management(output_pivot, "CatchAreaKm2", "DOUBLE")     landuse_fields = [f.name for f in arcpy.ListFields(output_pivot)[1:]]     print(landuse_fields)     with arcpy.da.UpdateCursor(output_pivot, landuse_fields) as upcur:  # @UndefinedVariable         for row in upcur:             keyvalue = row[0]             if keyvalue in valuedict:                 arcpy.AddMessage("Processing Watershed {0}".format(keyvalue))                 # add catchment area km² to newly added field                 row[-1] = valuedict[keyvalue][0]/1000000                 for i, landuse in enumerate(row[1:]): # index slice not picking up newly added field                     row = (row*(28*28))/(valuedict[keyvalue][0])*100                     upcur.updateRow(row) landuse_percentage(watershed, output_pivot) Python: Code (Currently working) # calculate percentage of landuse per watershed def landuse_percentage(watershed, output_pivot):     watershed_fields = ["HydroID", "Shape_Area"]     valuedict = {r[0]: (r[1:]) for r in arcpy.da.SearchCursor(watershed, watershed_fields)}  # @UndefinedVariable     arcpy.AddField_management(output_pivot, "CatchAreaKm2", "DOUBLE")     landuse_fields = [f.name for f in arcpy.ListFields(output_pivot)[1:]]     print(landuse_fields)     with arcpy.da.UpdateCursor(output_pivot, landuse_fields) as upcur:  # @UndefinedVariable         for row in upcur:             keyvalue = row[0]             if keyvalue in valuedict:                 arcpy.AddMessage("Processing Watershed {0}".format(keyvalue))                 # add catchment area km² to newly added field                 row[-1] = valuedict[keyvalue][0]/1000000                 for i, landuse in enumerate(row[1:-1]): # index slice should be                     row = (row*(28*28))/(valuedict[keyvalue][0])*100                     upcur.updateRow(row) landuse_percentage(watershed, output_pivot) Python: Code (Not working) ... View more