diff --git a/speckle_toolbox/esri/toolboxes/speckle/speckle/converter/features/feature_conversions.py b/speckle_toolbox/esri/toolboxes/speckle/speckle/converter/features/feature_conversions.py index 0c8bf3b..2dd0ec9 100644 --- a/speckle_toolbox/esri/toolboxes/speckle/speckle/converter/features/feature_conversions.py +++ b/speckle_toolbox/esri/toolboxes/speckle/speckle/converter/features/feature_conversions.py @@ -193,6 +193,7 @@ def rasterFeatureToSpeckle( rasterBandCount = my_raster.bandCount rasterBandNames = my_raster.bandNames rasterDimensions = [my_raster.width, my_raster.height] + print(rasterDimensions) # ds = gdal.Open(selectedLayer.source(), gdal.GA_ReadOnly) extent = my_raster.extent @@ -241,17 +242,20 @@ def rasterFeatureToSpeckle( rb = my_raster.getRasterBands(item) print(rb) - print(np.shape(rb.read())) + size = np.shape(rb.read()) + print(size) + rasterDimensions = [size[1], size[0]] + valMin = rb.minimum valMax = rb.maximum bandVals = np.swapaxes(rb.read(), 1, 2).flatten() # .tolist() np.flip( , 0) - bandValsFlat = [] - bandValsFlat.extend(bandVals.tolist()) - # look at mesh chunking + bandValsFlat = bandVals.tolist() + # print(bandValsFlat) const = float(-1 * math.pow(10, 30)) defaultNoData = rb.noDataValue + print(defaultNoData) # check whether NA value is too small or raster has too small values # assign min value of an actual list; re-assign NA val; replace list items to new NA val @@ -298,7 +302,7 @@ def rasterFeatureToSpeckle( # if rasterBandNoDataVal[len(rasterBandNoDataVal) - 1] is None: # rasterBandNoDataVal[len(rasterBandNoDataVal) - 1] = np.nan - + print(len(bandValsFlat)) rasterBandVals.append(bandValsFlat) rasterBandMinVal.append(valMin) rasterBandMaxVal.append(valMax)