Skip to Content.

mget-help - [mget-help] Problem using Find ArcGIS Rasters and Interpolate No Data Cells tool

Please Wait...

Subject: Marine Geospatial Ecology Tools (MGET) help

Text archives


From: Bryan Costa - NOAA Affiliate <>
To:
Subject: [mget-help] Problem using Find ArcGIS Rasters and Interpolate No Data Cells tool
Date: Fri, 26 Sep 2014 12:44:07 -0400
Hi,
I'm writing because I keep getting an error when I run MGET's Find ArcGIS Rasters and Interpolate No Data Cells tool. The error is listed below.  I tried reducing the size of the rasters (because it says the array is too big), but still no luck.  

I look forward to hearing from you at your convenience.

Regards,
Bryan C

---------------------------------------------
Error

Executing: InterpolatorFindAndInpaintArcGISRasters F:\Model_Predictors\Resolution_1p2x1p2km F:\Model_Predictors\Fill_Holidays Del2a # false * false # os.path.join(outputWorkspace, inputRaster[len(workspaceToSearch)+1:]) os.path false
Start Time: Fri Sep 26 12:39:42 2014
Running script InterpolatorFindAndInpaintArcGISRasters...
Finding ArcGIS rasters and inserting rows into table "work": workspace="F:\Model_Predictors\Resolution_1p2x1p2km", wildcard="*", searchTree=False, rasterType="None"
Finished inserting rows: 0:00:00 elapsed, 48 rows inserted, 0:00:00.002750 per row.
Updating field "outputRaster" in table "work"...
Update complete: 0:00:00 elapsed, 48 rows updated, 0 deleted, 0 unchanged, 0:00:00.000270 per row.
Querying the table "work" to build lists of inputs and outputs...
Query complete: 0:00:00 elapsed, 48 rows retrieved, 0:00:00.000166 per row.
Processing 48 inputs...
RuntimeError: Failed to retrieve a block of data of 37912 columns by 35946 rows at offsets x=0, y=0 from band 1 of band 1 of GDAL dataset "F:\Model_Predictors\Resolution_1p2x1p2km\Chlorophylla_01Jan_Mean_2003to2013_1p2x1p2km.tif" with the Geospatial Data Abstraction Library (GDAL). Verify that the dataset exists, is accessible, and has the expected dimensions. Detailed error information: band.ReadAsArray(0, 0, 37912, 35946) reported ValueError: array is too big..

Traceback (most recent call last):
  File "C:\Program Files\GeoEco\ArcGISToolbox\Scripts\InterpolatorFindAndInpaintArcGISRasters.py", line 5, in <module>
    ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.SpatialAnalysis.Interpolation', 'Interpolator', 'FindAndInpaintArcGISRasters')
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\ArcGISScripts.py", line 210, in ExecuteMethodFromCommandLineAsArcGISTool
    exec sourceCode in globals(), locals()
  File "<string>", line 1, in <module>
  File "<string>", line 7, in FindAndInpaintArcGISRasters
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\BatchProcessing.py", line 2161, in ExecuteFindAndProcessMethod
    boundProcessTableMethod(*tuple(argList))
  File "<string>", line 5, in InpaintArcGISRasterTable
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\BatchProcessing.py", line 1295, in ExecuteProcessTableMethod
    results = boundProcessListMethod(*tuple(argList))
  File "<string>", line 5, in InpaintArcGISRasterList
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\BatchProcessing.py", line 912, in ExecuteProcessListMethod
    result = boundMethod(*tuple(argList))
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\SpatialAnalysis\Interpolation.py", line 726, in InpaintArcGISRaster
    ArcGISWorkspace(os.path.dirname(outputRaster), ArcGISRaster, pathCreationExpressions=[os.path.basename(outputRaster)]).ImportDatasets([InpaintedGrid(ArcGISRasterBand.ConstructFromArcGISPath(inputRaster), method, maxHoleSize, xEdgesWrap)], {False: u'Add', True: u'Replace'}[overwriteExisting], reportProgress=False)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 816, in ImportDatasets
    self._ImportDatasets(datasets, mode.lower(), reportProgress, options)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\Collections.py", line 694, in _ImportDatasets
    self._ImportDatasetsToPath(pathComponentsForPath[path], datasetsForPath[path], mode, progressReporter, options)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\ArcGIS.py", line 607, in _ImportDatasetsToPath
    self.DatasetType._ImportDatasetsToPath(os.path.join(self.Path, *pathComponents), sourceDatasets, mode, progressReporter, options)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\ArcGIS.py", line 1119, in _ImportDatasetsToPath
    GDALDataset._ImportDatasetsToPath(path, sourceDatasets, mode, None, {'useArcGISSpatialReference': True, 'useUnscaledData': useUnscaledData, 'calculateStatistics': False, 'blockSize': blockSize})
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\GDAL.py", line 1073, in _ImportDatasetsToPath
    data = "sourceDatasets[i].Data[rowsCopied:rowsCopied+rowsToCopy," :]
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3666, in __getitem__
    return getattr(self._Grid(), self._GetMethod)(key)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3306, in _GetUnscaledDataAsArray
    data, actualNoDataValue = self._ReadNumpyArray(reorderedSliceList)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\Virtual.py", line 3337, in _ReadNumpyArray
    grid.Data[:].tofile(holesFile)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3666, in __getitem__
    return getattr(self._Grid(), self._GetMethod)(key)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3306, in _GetUnscaledDataAsArray
    data, actualNoDataValue = self._ReadNumpyArray(reorderedSliceList)
  File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\GDAL.py", line 1374, in _ReadNumpyArray
    raise RuntimeError(_(u'Failed to retrieve a block of data of %(win_xsize)i columns by %(win_ysize)i rows at offsets x=%(xoff)i, y=%(yoff)i from band %(band)i of %(dn)s with the Geospatial Data Abstraction Library (GDAL). Verify that the dataset exists, is accessible, and has the expected dimensions. Detailed error information: band.ReadAsArray(%(xoff)i, %(yoff)i, %(win_xsize)i, %(win_ysize)i) reported %(e)s: %(msg)s.') % {u'band': self.Band, u'dn': self.DisplayName, u'xoff': xoff, u'yoff': yoff, u'win_xsize': win_xsize, u'win_ysize': win_ysize, u'e': e.__class__.__name__, u'msg': self._Unicode(e)})
RuntimeError: Failed to retrieve a block of data of 37912 columns by 35946 rows at offsets x=0, y=0 from band 1 of band 1 of GDAL dataset "F:\Model_Predictors\Resolution_1p2x1p2km\Chlorophylla_01Jan_Mean_2003to2013_1p2x1p2km.tif" with the Geospatial Data Abstraction Library (GDAL). Verify that the dataset exists, is accessible, and has the expected dimensions. Detailed error information: band.ReadAsArray(0, 0, 37912, 35946) reported ValueError: array is too big..

Failed to execute (InterpolatorFindAndInpaintArcGISRasters).
Failed at Fri Sep 26 12:39:50 2014 (Elapsed Time: 8.00 seconds)


--
Bryan Costa
Geospatial Scientist
Ecosystem Modeler
CSS-Dynamac 

NOAA|CCMA|Biogeography Branch
1305 East West Highway
N-SCI-1, SSMC 4, 9th Floor, #9232
Silver Spring, MD 20910

Phone: (301) 713-3028 x146
Fax: (301) 713-4384
Email:

The contents of this message are mine and do not necessarily reflect any position of NOAA
Archives powered by MHonArc.
Top of Page