It looks like my knack for stumbling into weird edge cases is still alive and well! Whee.
I'm using gdal from Qgis v3.16 on Windows. I was trying the (mode,max,min,med) resampling methods and get the error reliably with `mode`. I thought I got it on a 2nd or 3rd method also but didn't repeat that this morning. I've been playing with GDAL_CACHEMAX and -wm numbers so that might be a factor. Working directory is a local SSD drive and there is 32 GB physical memory, generally 50% or more free. Free drive space is on the low end, 10 GB on C: and 7 GB on D:, where the working folder is. ~~~ gdalwarp ^ --config GDAL_CACHEMAX 4096 ^ -wm 512 ^ -multi ^ -co TILED=YES ^ -co COMPRESS=ZSTD ^ -t_srs EPSG:3579 ^ -r mode ^ -co predictor=2 ^ tri_mergedDEM_clipped_16bit.tif ^ tri_mergedDEM_clipped_16bit_ytalb-rs-mode.tif ~~~ Console: ~~~ $ gdalwarp --config GDAL_CACHEMAX 4096 -wm 512 -multi -co TILED=YES -co COMPRESS=ZSTD -t_srs EPSG:3579 -r mode -co predictor=2 tri_mergedDEM_clipped_16bit.tif tri_mergedDEM_clipped_16bit_ytalb-rs-mode.tif Creating output file that is 89772P x 73148L. Processing tri_mergedDEM_clipped_16bit.tif [1/1] : 0Using internal nodata values (e.g. -99) for image tri_mergedDEM_clipped_16bit.tif. Copying nodata values from source tri_mergedDEM_clipped_16bit.tif to destination tri_mergedDEM_clipped_16bit_ytalb-rs-mode.tif. ..ERROR 1: Failed to acquire WarpMutex in WarpRegion(). . ~~~ The info report for source image is attached, which can be downloaded from https://sft.gov.yk.ca/link/gV2K9F4PyH1PTJLQH8fy99 until April 14. (500mb). -Matt From: Even Rouault <[email protected]> Sent: March 30, 2021 11:39 PM To: Matt.Wilkie <[email protected]>; [email protected] Subject: Re: [gdal-dev] More info on gdalwarp resampling methods? Matt, (leaving out a bit your general question) The below looks is an error unrelated to whatever resampling method you apply. This is completely unfamiliar though. Do you build GDAL yourself ? Did you disable threading at build time ? Or maybe you're enabling threading with a huge number of threads ? Does that happen all the time ? Even Background: I'm getting error messages when trying to use some of the methods with 16bit imagery and I don't know if I'm seeing a gdal bug, a data bug, or just "you can't use that method with this raster type". ERROR 1: Failed to acquire WarpMutex in WarpRegion(). Thanks! Matt Wilkie Geomatics Analyst Environment | Technology, Innovation and Mapping T 867-667-8133 | Yukon.ca<https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2f%2fyukon.ca&umid=AC1EC954-BECF-5E05-B6CB-AE995C90C7E2&auth=c132af8ee7c9d1278d61a701569070a095ce962e-90d4f63b9d26941043d059b2a373176c1c6f200f> Hours: 08:30-16:30, Tue-Wed: Office, Thu-Fri: Remote. _______________________________________________ gdal-dev mailing list [email protected]<mailto:[email protected]> https://lists.osgeo.org/mailman/listinfo/gdal-dev -- http://www.spatialys.com<https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2f%2fwww.spatialys.com&umid=AC1EC954-BECF-5E05-B6CB-AE995C90C7E2&auth=c132af8ee7c9d1278d61a701569070a095ce962e-76fdcd2f1d6f12f56ce5a5d39e4c63f2993393ad> My software is free, but my time generally not.
tri_mergedDEM_clipped_16bit.tif.info
Description: tri_mergedDEM_clipped_16bit.tif.info
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
