Merge branch 'bugfix/fix_stretching' into 'master'

Bugfix/fix stretching

See merge request !6

geoarray/baseclasses.py

@@ -1054,7 +1054,10 @@ class GeoArray(object):
 
         transOpt = ['SRC_METHOD=NO_GEOTRANSFORM'] if tuple(gt) == (0, 1, 0, 0, 0, -1) else None
         xdim, ydim = None, None
-        nodataVal = nodataVal if nodataVal is not None else self.nodata
+        in_nodata = nodataVal if nodataVal is not None else self.nodata
+        out_nodata = in_nodata if in_nodata is not None else -9999
+        if not np.can_cast(out_nodata, image2plot.dtype):
+            image2plot = image2plot.astype(np.int32)
 
         # rotated images always have to be resampled for plotting
         if gt[2] or gt[4]:
@@ -1070,8 +1073,12 @@ class GeoArray(object):
         if xdim or ydim or out_prj:
             from py_tools_ds.geo.raster.reproject import warp_ndarray
             image2plot, gt, prj = warp_ndarray(image2plot, self.geotransform, self.projection,
-                                               out_XYdims=(xdim, ydim), in_nodata=nodataVal, out_nodata=nodataVal,
-                                               transformerOptions=transOpt, out_prj=out_prj, q=True)
+                                               out_XYdims=(xdim, ydim),
+                                               in_nodata=in_nodata,
+                                               out_nodata=out_nodata,
+                                               transformerOptions=transOpt,
+                                               out_prj=out_prj,
+                                               q=True)
             if transOpt and 'NO_GEOTRANSFORM' in ','.join(transOpt):
                 image2plot = np.flipud(image2plot)
                 gt = list(gt)
@@ -1112,7 +1119,7 @@ class GeoArray(object):
         band = (band if band is not None else 0) if not interactive else band
 
         # get image to plot
-        nodataVal = nodataVal if nodataVal is not None else self.nodata
+        nodataVal = nodataVal if nodataVal is not None else self.nodata if self.nodata is not None else -9999
         image2plot, gt, prj = \
             self._get_plottable_image(xlim, ylim, band, boundsMap=boundsMap, boundsMapPrj=boundsMapPrj,
                                       res_factor=res_factor, nodataVal=nodataVal)

geoarray/version.py

-__version__ = '0.8.15'
-__versionalias__ = '20190429.01'
+__version__ = '0.8.16'
+__versionalias__ = '20190429.02'