Subject: Marine Geospatial Ecology Tools (MGET) help
Text archives
| From: | "Jason Roberts" <> |
|---|---|
| To: | "'Elizabeth Davidson'" <> |
| Cc: | <> |
| Subject: | RE: [mget-help] Interpolating No Data |
| Date: | Wed, 27 Nov 2013 11:39:03 -0500 |
Hi Liz, I should be able to help you fix this. First question: did you install the MATLAB MCR (from here)? The tool you are using requires it. The tool is supposed to check for it and report a more understandable error message if it is missing, so perhaps that is not the problem, but I want to check with you first before going on… Best, Jaso From: Elizabeth Davidson [mailto:] 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 |