<div dir="ltr"><div class="gmail_extra"><div class="gmail_extra">I'll have to look into this more.</div><div class="gmail_extra"><br></div><div class="gmail_extra">gdal_grid --debug on -txe 440720.0 441920.0 -tye 3751320.0 3750120.0 \</div><div class="gmail_extra">> -outsize 20 20 -ot Float64 -l grid \</div><div class="gmail_extra">> -a invdist:power=2.0:smoothing=0.0:radius1=0.0:radius2=0.0:angle=0.0:max_points=0:min_points=0:nodata=0.0 data/grid.vrt out.tif</div><div class="gmail_extra">OGR: OGROpen(data/grid.vrt/0x7f8c9be00240) succeeded as VRT.</div><div class="gmail_extra">GDAL: GDALDriver::Create(GTiff,out.tif,20,20,1,Float64,0x0)</div><div class="gmail_extra">OGR: OGROpen(data/grid.csv/0x7f8c9bc24350) succeeded as CSV.</div><div class="gmail_extra">Grid data type is "Float64"</div><div class="gmail_extra">Grid size = (20 20).</div><div class="gmail_extra">Corner coordinates = (440690.000000 3750090.000000)-(441950.000000 3751350.000000).</div><div class="gmail_extra">Grid cell size = (60.000000 -60.000000).</div><div class="gmail_extra">Source point count = 400.</div><div class="gmail_extra">Algorithm name: "invdist".</div><div class="gmail_extra">Options are "power=2.000000:smoothing=0.000000:radius1=0.000000:radius2=0.000000:angle=0.000000:max_points=0:min_points=0:nodata=0.000000"</div><div class="gmail_extra"><br></div><div class="gmail_extra">GDAL_GRID: Work buffer: 20 * 20</div><div class="gmail_extra">GDAL_GRID: Using AVX optimized version</div><div class="gmail_extra">GDAL_GRID: Using 8 threads</div><div class="gmail_extra">Segmentation fault: 11</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 12, 2015 at 12:23 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Le jeudi 12 février 2015 20:54:42, Kurt Schwehr a écrit :<br>
<div><div class="h5">> Here is the full output:<br>
><br>
> <a href="https://gist.github.com/schwehr/c9e9fb595d0370e7eff7" target="_blank">https://gist.github.com/schwehr/c9e9fb595d0370e7eff7</a><br>
><br>
> fink list -i libjpeg webp<br>
> i libjpeg-bin 9a-2 Executables for<br>
> libjpeg package<br>
> i libjpeg9 9a-2 JPEG image format<br>
> handling library<br>
> i libjpeg9-shlibs 9a-2 Shared libraries<br>
> for libjpeg package<br>
> i libwebp5 0.4.0-1 Library for<br>
> manipulating WebP format<br>
> i libwebp5-shlibs 0.4.0-1 Library for<br>
> manipulating WebP format<br>
><br>
> And the details of test_gdal_grid_3:<br>
><br>
> ./test_gdal_grid.py<br>
> TEST: test_gdal_grid_1 ... success<br>
> TEST: test_gdal_grid_2 ... success<br>
> TEST: test_gdal_grid_3 ... ERROR 4: `tmp/grid_invdist.tif' not recognised<br>
> as a supported file format.<br>
><br>
> fail (blowup)<br>
> Traceback (most recent call last):<br>
> File "../pymod/gdaltest_python2.py", line 42, in run_func<br>
> result = func()<br>
> File "./test_gdal_grid.py", line 273, in test_gdal_grid_3<br>
> maxdiff = gdaltest.compare_ds(ds, ds_ref, verbose = 0)<br>
> File "../pymod/gdaltest.py", line 1089, in compare_ds<br>
> width = ds1.RasterXSize<br>
> AttributeError: 'NoneType' object has no attribute 'RasterXSize'<br>
> TEST: test_gdal_grid_4 ... success<br>
> TEST: test_gdal_grid_5 ... success<br>
> TEST: test_gdal_grid_6 ... success<br>
> TEST: test_gdal_grid_7 ... success<br>
> TEST: test_gdal_grid_8 ... success<br>
> TEST: test_gdal_grid_9 ... success<br>
> TEST: test_gdal_grid_10 ... success<br>
> TEST: test_gdal_grid_cleanup ... ERROR 4: `tmp/grid_invdist.tif' not<br>
> recognised as a supported file format.<br>
><br>
> success<br>
<br>
</div></div>Kurt,<br>
<br>
I'm wondering if it wouldn't be an issue with SSE / AVX optimizations.<br>
<br>
Could you go into autotest/utilities and run manually the following ?<br>
<br>
$ gdal_grid --debug on -txe 440720.0 441920.0 -tye 3751320.0 3750120.0 \<br>
-outsize 20 20 -ot Float64 -l grid \<br>
-a invdist:power=2.0:smoothing=0.0:radius1=0.0:radius2=0.0:angle=0.0:max_points=0:min_points=0:nodata=0.0 data/grid.vrt out.tif<br>
<br>
Here's what I get :<br>
<br>
OGR: OGROpen(data/grid.vrt/0x2482e50) succeeded as VRT.<br>
GDAL: GDALDriver::Create(GTiff,out.tif,20,20,1,Float64,(nil))<br>
OGR: OGROpen(data/grid.csv/0x24844c0) succeeded as CSV.<br>
Grid data type is "Float64"<br>
Grid size = (20 20).<br>
Corner coordinates = (440690.000000 3750090.000000)-(441950.000000 3751350.000000).<br>
Grid cell size = (60.000000 -60.000000).<br>
Source point count = 400.<br>
Algorithm name: "invdist".<br>
Options are "power=2.000000:smoothing=0.000000:radius1=0.000000:radius2=0.000000:angle=0.000000:max_points=0:min_points=0:nodata=0.000000"<br>
<br>
GDAL_GRID: Work buffer: 20 * 20<br>
GDAL_GRID: Using SSE optimized version<br>
GDAL_GRID: Using 4 threads<br>
0...10...20...30...40...50...60...70...80...90...100 - done.<br>
VRT: 400 features read on layer 'grid'.<br>
CSV: 400 features read on layer 'grid'.<br>
GDAL: GDALClose(out.tif, this=0x24baac0)<br>
<br>
If you have issues, you can try adding "--config GDAL_USE_AVX NO" and/or "--config GDAL_USE_SSE NO"<br>
You can also play with the GDAL_NUM_THREADS configuration option.<br>
<br>
Does this work with trunk ? There have been some changes by moving SSE specifc code into a specific file.<br>
<span class=""><font color="#888888"><br>
Even<br>
</font></span><div class=""><div class="h5"><br>
<br>
><br>
> Test Script: test_gdal_grid<br>
> Succeeded: 10<br>
> Failed: 1 (1 blew exceptions)<br>
> Skipped: 0<br>
> Expected fail:0<br>
> Duration: 1.09s<br>
><br>
> file tmp/grid_invdist.tif<br>
> tmp/grid_invdist.tif: empty<br>
><br>
><br>
> On Wed, Feb 11, 2015 at 12:28 AM, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>><br>
><br>
> wrote:<br>
> > Le mardi 10 février 2015 20:05:52, Kurt Schwehr a écrit :<br>
> > > I think it looks okay w.r.t. fink on mac osx 10.9. autotest failures<br>
> > > for anyone who is interested:<br>
> > > <a href="https://gist.github.com/schwehr/757bae6bd28d9a7ffc4d" target="_blank">https://gist.github.com/schwehr/757bae6bd28d9a7ffc4d</a><br>
> ><br>
> > Kurt,<br>
> ><br>
> > Some of the output traces that contain "fail" do not necessarily make the<br>
> > step<br>
> > a failure, but should skip it. The error message and information printed<br>
> > and<br>
> > line reported might be interesting.<br>
> > Most must be due to a different version of the JPEG or webp lib that must<br>
> > produce slightly different results.<br>
> > The exception stack associated with "TEST: test_gdal_grid_3 ... fail<br>
> > (blowup)"<br>
> > would be interesting<br>
> ><br>
> > Even<br>
> ><br>
> > > On Tue, Feb 10, 2015 at 4:30 AM, Even Rouault <<br>
> ><br>
> > <a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>><br>
> ><br>
> > > wrote:<br>
> > > > Hi,<br>
> > > ><br>
> > > > I have prepared a second GDAL/OGR 1.11.2 release candidate, with the<br>
> ><br>
> > fix<br>
> ><br>
> > > > for<br>
> > > > the 1.11.0 regression in the PG driver identified yesterday (#5837).<br>
> > > > Please review and test.<br>
> > > ><br>
> > > > Peek up an archive among the following ones (by ascending size):<br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdal-1.11.2RC2.tar.xz" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdal-1.11.2RC2.tar.xz</a><br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdal-1.11.2RC2.tar.gz" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdal-1.11.2RC2.tar.gz</a><br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdal1112RC2.zip" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdal1112RC2.zip</a><br>
> > > ><br>
> > > > The gdal-grass plugin archive hasn't been modified since RC1 :<br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdal-grass-1.11.2.tar.gz" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdal-grass-1.11.2.tar.gz</a><br>
> > > ><br>
> > > > A snapshot of the gdalautotest suite is also available :<br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdalautotest-1.11.2.tar.gz" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdalautotest-1.11.2.tar.gz</a><br>
> > > > <a href="http://download.osgeo.org/gdal/1.11.2/gdalautotest-1.11.2.zip" target="_blank">http://download.osgeo.org/gdal/1.11.2/gdalautotest-1.11.2.zip</a><br>
> > > ><br>
> > > > The NEWS file is here :<br>
> > > > <a href="http://trac.osgeo.org/gdal/wiki/Release/1.11.2-News" target="_blank">http://trac.osgeo.org/gdal/wiki/Release/1.11.2-News</a><br>
> > > ><br>
> > > > I'll call for a vote promoting it soon if no serious problems are<br>
> > > > reported before.<br>
> > > ><br>
> > > > Best regards,<br>
> > > ><br>
> > > > Even<br>
> > > > --<br>
> > > > Spatialys - Geospatial professional services<br>
> > > > <a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
> > > > _______________________________________________<br>
> > > > gdal-dev mailing list<br>
> > > > <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
> > > > <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
> ><br>
> > --<br>
> > Spatialys - Geospatial professional services<br>
> > <a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">--<div><a href="http://schwehr.org" target="_blank">http://schwehr.org</a></div></div>
</div></div>