Subject: Marine Geospatial Ecology Tools (MGET) help
Text archives
| From: | Elizabeth Davidson <> |
|---|---|
| To: | |
| Subject: | [mget-help] Interpolating No Data |
| Date: | Wed, 27 Nov 2013 11:30:07 -0500 |
Hi all,
I'm new to MGET, and trying to get started by filling in some holes in my data. I am using MODIS chlorophyll data, and have imported them into ArcMap10.1 using SeaDAS to convert the files. Cloudy days are causing pixels with No Data in my raster, and I would like to use your interpolate tool to try to fill in the missing data.
I am receiving this error message when I try to interpolate, no matter what method I choose:
Executing: InterpolatorInpaintArcGISRaster august03_clip4 "Y:\OceanColourFoxeBasin\Ocean Colour Files\2012 Monthly Comps (L3).gdb\interpolate6" Del2a 25 false
Start Time: Wed Nov 27 11:29:01 2013
Running script InterpolatorInpaintArcGISRaster...
python.exe: Error: Invalid noDataValue: None
RuntimeError: python.exe returned exit code 1, indicating failure.
The following consequences resulted from the original error:
Execution of C:\Python27\ArcGIS10.1\python.exe failed.
Traceback (most recent call last):
File "C:\Program Files\GeoEco\ArcGISToolbox\Scripts\InterpolatorInpaintArcGISRaster.py", line 5, in <module>
ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.SpatialAnalysis.Interpolation', 'Interpolator', 'InpaintArcGISRaster')
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 "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\SpatialAnalysis\Interpolation.py", line 674, 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 1139, in _ImportDatasetsToPath
GDALDataset._ImportDatasetsToPath(tempRaster, 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 3635, in __getitem__
return getattr(self._Grid(), self._GetMethod)(key)
File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\__init__.py", line 3275, in _GetUnscaledDataAsArray
data, actualNoDataValue = self._ReadNumpyArray(reorderedSliceList)
File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Datasets\Virtual.py", line 3377, in _ReadNumpyArray
maxRunTime=None)
File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\DataManagement\Processes.py", line 575, in ExecuteProgram
process.Execute()
File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\DataManagement\Processes.py", line 498, in Execute
Logger.RaiseException(RuntimeError(_(u'%(name)s returned exit code %(code)i, indicating failure.') % {u'name' : programName, u'code' : self.ExitCode}))
File "C:\Python27\ArcGIS10.1\lib\site-packages\GeoEco\Logging.py", line 104, in RaiseException
raise exception
RuntimeError: python.exe returned exit code 1, indicating failure.
Failed to execute (InterpolatorInpaintArcGISRaster).
Failed at Wed Nov 27 11:29:03 2013 (Elapsed Time: 2.00 seconds)
Hoping someone might be able to offer some assistance.
Thanks,
Liz
Archives powered by MHonArc.