topic Re: CalculateStatistics on Mosaic Dataset in Python Questions

CalculateStatistics on Mosaic Dataset

JamesCrandall — Mon, 16 Dec 2013 12:28:41 GMT

I am attempting to calculate statistics on mosaic datasets in a FGDB. It appears to work correctly however it runs and I add the mosaic to the TOC and access the symbology properties of the Mosaic dataset, it does not allow for Classified type (only Stretched and Discrete color are available).

I am expecting to be able to apply a classified renderer on the mosaic, which it works as expected if I manually run the Calculate Statistics from the "Enhance" menu (by right clicking on the mosaic in the catalog tree and accessing this menu option there).

Is there any reason why I can access the classified symbol property after running calculate stat from the Enhance menu, but not from the code below?

Thanks,
j


def calcStatsMOSAIC(t):

  try:
      
    arcpy.env.workspace = ws_output

    bndy_mosaic = str(ws_output) + "\\rasext" + str(t)
    mosaic_name = "ConcMosaic" + str(t)
    ws_mosaic = str(ws_output) + "\\" + str(mosaic_name)
    tr = "*Mosaic*"
    rasters = arcpy.ListDatasets(tr) 
    for r in rasters:
        
       ##CALC STATISTICS
       arcpy.CalculateStatistics_management(ws_mosaic, 1, 1, "#", "OVERWRITE", bndy_mosaic)

Re: CalculateStatistics on Mosaic Dataset

JamesCrandall — Thu, 02 Jan 2014 14:05:01 GMT

bump to the top

Re: CalculateStatistics on Mosaic Dataset

JakeSkinner — Sat, 11 Dec 2021 23:36:57 GMT

Hi James,

I've noticed after a mosaic dataset is created and statistics are calculated using a script, they will not show until you restart ArcMap. You may want to follow up with Tech Support on this. Also, you can simplify your script:

def calcStatsMOSAIC(t):
    arcpy.env.workspace = ws_output
    mosaic_name = "ConcMosaic" + str(t)
    arcpy.CalculateStatistics_management(mosaic_name, 1, 1)

Try the following:

1. Create your mosaic dataset
2. Run the script
3. Restart ArcMap
4. Add the mosaic dataset to ArcMap and see if you 'Classified' for a symbology type

Re: CalculateStatistics on Mosaic Dataset

JamesCrandall — Sat, 11 Dec 2021 23:37:00 GMT

Hi James,

I've noticed after a mosaic dataset is created and statistics are calculated using a script, they will not show until you restart ArcMap. You may want to follow up with Tech Support on this. Also, you can simplify your script:
def calcStatsMOSAIC(t):
    arcpy.env.workspace = ws_output
    mosaic_name = "ConcMosaic" + str(t)
    arcpy.CalculateStatistics_management(mosaic_name, 1, 1)
Try the following:

1. Create your mosaic dataset
2. Run the script
3. Restart ArcMap
4. Add the mosaic dataset to ArcMap and see if you 'Classified' for a symbology type

Thanks for your input, Jake,

I ran the script/.py from PythonWin and from a Toolbox Script, both have the same result even after restarting ArcMap.

Also, I appreciate the simplification suggestion. I need to iterate over raster datasets as shown because there are both mosaic's and individual rasters, which I just need to acquire the mosaic to calc the statistics.

Re: CalculateStatistics on Mosaic Dataset

JamesCrandall — Fri, 03 Jan 2014 13:50:21 GMT

I think I have narrowed down the problem to the Area of Interest parameter.

1. In my script, I dynamically create a polygon feature class that is the same extent of the footprint/boundary of the mosaic. This FC does reside in the same FGDB as the mosaic, but is its own feature class.

2. Running Calculate Statistics toolbox script (which successfully adds classified symbology property) auto-populates the Area of Interest parameter with: CalculateStatistics::area_of_interest and adds an empty polygon fc to the Table of Contents.

So, there is something going on there with the tool adding CalculateStatistics::area_of_interest that I am not understanding. Any additional input is appreciated.

j