Skip to Content.

mget-help - [mget-help] Bug with Hycom data

Please Wait...

Subject: Marine Geospatial Ecology Tools (MGET) help

Text archives


From: Grant Humphries <>
To:
Subject: [mget-help] Bug with Hycom data
Date: Thu, 04 Jun 2015 20:26:08 +0000
Hi there,

I'm currently trying to use the Interpolate HYCOM GLBu0.08 4D Variables at Points tool to get HYCOM data for the california current ecosystem.  I'm getting the following error:

ServerError: 'Server error 3: "Invalid Parameter Exception: DArrayDimension.setProjection: Bad Projection Request: stop >= size: 386:386"'
Failed script HYCOMGLBu0084DInterpolateAtArcGISPoints...

Traceback (most recent call last):
  File "C:\Program Files\GeoEco\ArcGISToolbox\Scripts\HYCOMGLBu0084DInterpolateAtArcGISPoints.py", line 5, in <module>
    ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.DataProducts.HYCOM', 'HYCOMGLBu0084D', 'InterpolateAtArcGISPoints')
  File "C:\Python27\ArcGIS10.3\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.3\lib\site-packages\GeoEco\DataProducts\HYCOM.py", line 3020, in InterpolateAtArcGISPoints
    Interpolator.InterpolateGridsValuesForTableOfPoints(grids, ArcGISTable(points), valueFields, zField=zField, tField=tField, zValue=zValue, where=where, orderBy=orderBy, method=method, noDataValue=noDataValue, useAbsZ=True, seafloorZValue=20000., numBlocksToCacheInMemory=numBlocksToCacheInMemory, xBlockSize=xBlockSize, yBlockSize=yBlockSize, zBlockSize=zBlockSize, tBlockSize=tBlockSize)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\SpatialAnalysis\Interpolation.py", line 568, in InterpolateGridsValuesForTableOfPoints
    value = cls._InterpolatePointOnGrid_NearestNeighbor(grids[i], transformedCoords, noDataValues[i], debugLoggingEnabled)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\SpatialAnalysis\Interpolation.py", line 42, in _InterpolatePointOnGrid_NearestNeighbor
    nearestIndices = grid.GetIndicesForCoords(coords)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2779, in GetIndicesForCoords
    index = bisect.bisect_right(self.MaxCoords[d], coord)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 3670, in __getitem__
    return getattr(self._Grid(), self._GetMethod)(key)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2845, in _GetMaxCoords
    return self._GetCoordsForOffset(key, 0.5)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2881, in _GetCoordsForOffset
    self._CachedTCoords[fixedIncrementOffset] = self._GetTCoordsList(fixedIncrementOffset, self.Shape[coordNum])
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2625, in _GetShape
    return self.GetLazyPropertyValue('Shape')
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\Virtual.py", line 1691, in _GetLazyPropertyPhysicalValue
    return self._Grid.GetLazyPropertyValue(name)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\Virtual.py", line 1446, in _GetLazyPropertyPhysicalValue
    return self._Grids[0].GetLazyPropertyValue(name)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\DataProducts\HYCOM.py", line 2651, in _GetLazyPropertyPhysicalValue
    endDay = grid.CenterCoords['t', -1]
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 3670, in __getitem__
    return getattr(self._Grid(), self._GetMethod)(key)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2842, in _GetCenterCoords
    return self._GetCoordsForOffset(key, 0.0)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\Datasets\__init__.py", line 2857, in _GetCoordsForOffset
    return self._GetCoords(coord, coordNum, slices, sliceDims, fixedIncrementOffset)
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\DataProducts\HYCOM.py", line 2529, in _GetCoords
    time_ = list(self.ParentCollection._PydapDataset['time'][:])
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\AggregatedModules\pydap\model.py", line 192, in __getitem__
    return self.data[index]
  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\AggregatedModules\pydap\proxy.py", line 112, in __getitem__
    resp, data = "request(url)  File "C:\Python27\ArcGIS10.3\lib\site-packages\GeoEco\AggregatedModules\pydap\util\http.py", line 57, in request
    raise ServerError(msg)
ServerError: 'Server error 3: "Invalid Parameter Exception: DArrayDimension.setProjection: Bad Projection Request: stop >= size: 386:386"'

Failed to execute (HYCOMGLBu0084DInterpolateAtArcGISPoints).


Projection of my dataset is WGS 1984 - and have been able to use with the other data products (i.e. University of Maine ROMS) ,

Any help on this error would be greatly appreciated

Cheers
Grant Humphries


Archives powered by MHonArc.
Top of Page