py_tools_ds issueshttps://git.gfz-potsdam.de/danschef/py_tools_ds/-/issues2018-09-17T14:08:33+02:00https://git.gfz-potsdam.de/danschef/py_tools_ds/-/issues/4py_tools_ds.geo.raster.warp_ndarray() freezes / crashes in case of very large...2018-09-17T14:08:33+02:00Daniel Schefflerpy_tools_ds.geo.raster.warp_ndarray() freezes / crashes in case of very large input rastersCrashes at line 430:
```python
res_ds = gdal_Warp(
'', in_ds, format='MEM',
dstSRS=get_SRS(out_prj),
outputType=get_GDT(out_dtype) if out_dtype else get_GDT(in_dtype_np),
xRes=out_gsd[0],
yRes=out...Crashes at line 430:
```python
res_ds = gdal_Warp(
'', in_ds, format='MEM',
dstSRS=get_SRS(out_prj),
outputType=get_GDT(out_dtype) if out_dtype else get_GDT(in_dtype_np),
xRes=out_gsd[0],
yRes=out_gsd[1],
outputBounds=out_bounds,
outputBoundsSRS=get_SRS(out_bounds_prj),
width=out_XYdims[0],
height=out_XYdims[1],
resampleAlg=rspAlg,
srcNodata=in_nodata,
dstNodata=out_nodata,
srcAlpha=in_alpha,
dstAlpha=out_alpha,
options=options if options else [],
warpOptions=warpOptions or [],
transformerOptions=transformerOptions or [],
targetAlignedPixels=targetAlignedPixels,
tps=True if gcpList else False,
polynomialOrder=polynomialOrder,
warpMemoryLimit=warpMemoryLimit,
callback=ProgressBar(prefix='Warping progress ', timeout=None) if progress and not q else None,
callback_data=[0],
errorThreshold=0.125, # this is needed to get exactly the same output like the console version of GDAL warp
)
```Daniel SchefflerDaniel Scheffler