Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

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 = <your file geodatabase>    
arcpy.env.outputCoordinateSystem = arcpy.SpatialReference("North America Albers Equal Area Conic")
arcpy.env.OverwriteOutput = True
arcpy.env.pyramid = "PYRAMIDS -1 NEAREST DEFAULT"
#arcpy.env.snapRaster = <your snap raster> # I snap to an NLCD layer, not sure its pertinent to the example, though
arcpy.env.cellSize = 30

inFolder = <the folder which contains the extracted zip file contents>
inList = []
arcpy.env.workspace= inFolder
for 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: enter image description here

And as I zoom in incrementally, the gaps fill: enter image description here enter image description here

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.
share|improve this question
add comment

1 Answer 1

Check to see if build overviews or pyramids is on when you run this process. The missing display indicates that the overview or pyramid for the "white" area did not generate correctly.

share|improve this answer
    
I initially suspected this also, but I've rebuilt overviews and pyramids and also recalculated statistics - none of which have provided improvement. –  metasequoia May 23 '13 at 3:24
    
Hmmm....I wonder if generating the pyramids is picking up white space cell values and with nearest neighbor they are being promoted to larger bands at the pyramid levels. Is your file geodatabase local on your computer, and the standard file outputs it local? Another maybe is network bandwidth. –  user18424 May 23 '13 at 10:25
    
All files are local. I'm curious if you've been able to replicate the issue. –  metasequoia May 23 '13 at 12:13
    
For the record I just downloaded the data and ran it through 10.2 and used the Mosaic To New Raster tool storing the final raster in a file geodatabase. I did not get the white lines as shown in your images. –  Hornbydd Nov 21 '13 at 16:45
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.