Latest Contributions by JohnWall

I have several thousand TIFF files that I need to convert from 32-bit to 16-bit, but the Copy Raster geoprocessing tool appears to be providing erroneous results. As an example, I have a TIF in 32-bit format with values ranging 0.0241 to 0.0674. When I use Copy Rater and set the "NoData Value" to blank, "Pixel Type" to 16-bit unsigned, set the "Scale Pixel Value" to true, I get a TIFF with values ranging 0 to 32768, but the image appears blank no matter what histogram stretch I use. I have tried Calculating Statistics within ArcCatalog with no results. When I set "Scale Pixel Value" to false, I get values ranging 0 to 1 with an image that appears blank. I am using ArcGIS 10.3.1. Is there something I am missing? ... View more

I am trying to run Train ISO Cluster Classifier using Python. Can someone help me identify why I am getting Error 999999? My code looks like: #Libraries import arcpy, os from arcpy import env from arcpy.sa import * inRaster = "C:/gis/oberlin/unsupervised/inrasternorm" path = os.path.dirname(inRaster) #Get directory name env.workspace = path #Set workspace env.overwriteOutput = True #Checkout extensions arcpy.CheckOutExtension("Spatial") #Build directory structure segmentFolder = path + "/segments" #Stores the segmented images try:    os.stat(segmentFolder) except:    os.mkdir(segmentFolder)    classifierFolder = path + "/classifier" #Stores the ECD classifier rules try:    os.stat(classifierFolder) except:    os.mkdir(classifierFolder)    classifiedFolder = path + "/classified" #Stores the classified images try:    os.stat(classifiedFolder) except:    os.mkdir(classifiedFolder) #Begin with the Segment Mean Shift variables that need to change for x in range(5, 16, 1):       spectral_detail = x    for y in range (5, 16, 1):    spatial_detail = y    min_segment_size = 100    seg_raster = SegmentMeanShift(inRaster, spectral_detail, spatial_detail, min_segment_size)    outName = segmentFolder + '/sms{:02d}_{:02d}.tif'.format(spectral_detail, spatial_detail)    seg_raster.save(outName)    TrainIsoClusterClassifier("C:/gis/oberlin/unsupervised/segments/sms05_06.tif", 5,    "C:/gis/oberlin/unsupervised/classifier/sms05_06.ecd", "C:/gis/oberlin/unsupervised/inrasternorm", 50, 100, "",    "MEAN")    #Classify Raster    classifiedraster = ClassifyRaster(seg_raster, out_definition, in_additional_raster)    outName2 = classifiedFolder + '/clsfid{:02d}_{:02d}.tif'.format(spectral_detail, spatial_detail)    classifiedraster.save(outName2) print "Done!" However, I get the following error: Traceback (most recent call last): File "C:\Python27\ArcGIS10.3\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 326, in RunScript exec codeObject in __main__.__dict__ File "C:\Users\John\Desktop\unspervised_classifier.py", line 68, in <module> TrainIsoClusterClassifier("C:/gis/oberlin/unsupervised/segments/sms05_06.tif", 5, "C:/gis/oberlin/unsupervised/classifier/sms05_06.ecd", "C:/gis/oberlin/unsupervised/inrasternorm", 50, 100, "", "MEAN") File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7741, in TrainIsoClusterClassifier used_attributes) File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py", line 53, in swapper result = wrapper(*args, **kwargs) File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7731, in Wrapper used_attributes) File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\geoprocessing\_base.py", line 504, in <lambda> return lambda *args: val(*gp_fixargs(args, True)) ExecuteError: ERROR 999999: Error executing function. Failed to execute (TrainIsoClusterClassifier). ... View more

I am new to image classification and would like to know about optimizing image segmentation within ArcGIS. From what I have read, image segmentation is possible with the Segment Mean Shift tool in ArcGIS 10.x, but I have not been able to find anything mentioning how one determines optimal segmentation has been reached. One article I found developed a method for identifying optimal segmentation by minimizing local variance [1], but it is not within ArcGIS (or any other software I have access to). Assuming minimizing variance within a segment is considered optimal segmentation, is this possible within ArcGIS? [1] Drǎguţ, L., Tiede, D., & Levick, S. R. (2010). ESP: a tool to estimate scale parameter for multiresolution image segmentation of remotely sensed data. International Journal of Geographical Information Science, 24(6), 859-871. ... View more

I have the following code: import arcpy, os from arcpy import env from arcpy.sa import * arcpy.env.workspace = "C:/gis/enforcer" #Dummy path arcpy.env.overwriteOutput = True dem = "C:/gis/enforcer/dem" #Dummy path rmseDEM = 0.09 #Dummy value roads = "C:/gis/enforcer/roads.shp" #Dummy path rails = "C:/gis/enforcer/rail.shp" #Dummy path fgdb = "enforcer.gdb" #Dummy path arcpy.env.cellSize = dem arcpy.env.extent = dem arcpy.env.sanpRaster = dem desc = arcpy.Describe(dem) demResolution = desc.meanCellHeight #Assumes a square dirName = os.path.dirname(dem) fgdbPath = dirName + "/" + fgdb arcpy.CheckOutExtension("Spatial") #-----------------------------------Identify Low Points---------------------------------- analysisSinks = fgdbPath + "/" + "sinkanalysis" outName = fgdbPath + "/" + "zonemin01" outZonalStats = ZonalStatistics(analysisSinks, "OBJECTID", dem, "MINIMUM", "NODATA") outZonalStats.save(outName) outName02 = "minpnt01" #fgdbPath + "/" outCon = Con(dem == outName, 1) outCon.save(outName02) However, Con in the next to last line seems to crash PythonWin every chance it gets or throws ERROR 010240. I've tried everything from saving it as a TIF to switching to a different directory to no avail. Any recommendations? ... View more