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.