Since upgrading to ArcGIS 10.1 (SP 1) from 10.0, I've been encountering an issue when running the arcpy.MosaicToNewRaster_management() tool with output to a file geodatabase. I'm principally working with python, but I've duplicated the issue in ArcMap.
Consider a sample dataset, available here, which includes two states (MA and CT) of gSSURGO soil rasters. Create a file geodatabase for outputs, then:
outGDB = arcpy.env.outputCoordinateSystem = arcpy.SpatialReference("North America Albers Equal Area Conic")arcpy.env.OverwriteOutput = Truearcpy.env.pyramid = "PYRAMIDS -1 NEAREST DEFAULT"#arcpy.env.snapRaster = # I snap to an NLCD layer, not sure its pertinent to the example, thougharcpy.env.cellSize = 30inFolder = inList = []arcpy.env.workspace= inFolderfor f in arcpy.ListFiles(): inList.append(inFolder+f+"/"+f)arcpy.MosaicToNewRaster_management(ify(inList,";"),outGDB,"soils",albers,"32_BIT_UNSIGNED","#","1",'FIRST')The operation successfully produces a new raster, but for me it has some unusual properties. Namely, I find large sections of the file do not display when zoomed out:
And as I zoom in incrementally, the gaps fill:
A few notes:
أكثر...
Consider a sample dataset, available here, which includes two states (MA and CT) of gSSURGO soil rasters. Create a file geodatabase for outputs, then:
outGDB = arcpy.env.outputCoordinateSystem = arcpy.SpatialReference("North America Albers Equal Area Conic")arcpy.env.OverwriteOutput = Truearcpy.env.pyramid = "PYRAMIDS -1 NEAREST DEFAULT"#arcpy.env.snapRaster = # I snap to an NLCD layer, not sure its pertinent to the example, thougharcpy.env.cellSize = 30inFolder = inList = []arcpy.env.workspace= inFolderfor f in arcpy.ListFiles(): inList.append(inFolder+f+"/"+f)arcpy.MosaicToNewRaster_management(ify(inList,";"),outGDB,"soils",albers,"32_BIT_UNSIGNED","#","1",'FIRST')The operation successfully produces a new raster, but for me it has some unusual properties. Namely, I find large sections of the file do not display when zoomed out:
And as I zoom in incrementally, the gaps fill:
A few notes:
- The "Identify" tool tip shows that there is data present even when not depicted. The same issue occurs in both ArcMap and ArcCatalog.
- The issue does not occur when I perform arcpy.MosaicToNewRaster_management() with output to a regular system folder.
- Scale range is set to "Show layers at all scales".
- As I zoom out, the whole extent disappears eventually.
أكثر...