<div dir="ltr"><div>Interesting. One difference is that I am running this on an input of type float32 with NODATA. I am on Widows 7. The data I am working with is from</div><div><br></div><div><a href="http://gis3.oit.ohio.gov/geodatadownload/osip.aspx">http://gis3.oit.ohio.gov/geodatadownload/osip.aspx</a></div><div><br></div><div>I am using the "ESRI GRID DEM Mosaic". Here is gdalinfo.</div><div><br></div><div>gdalinfo per_dem<br>Driver: AIG/Arc/Info Binary Grid<br>Files: per_dem<br>       per_dem.aux<br>       per_dem.rrd<br>       per_dem.aux.xml<br>       per_dem\dblbnd.adf<br>       per_dem\hdr.adf<br>       per_dem\metadata.xml<br>       per_dem\prj.adf<br>       per_dem\sta.adf<br>       per_dem\w001000.adf<br>       per_dem\w001000x.adf<br>       per_dem\w001001.adf<br>       per_dem\w001001x.adf<br>       per_dem\z001001.adf<br>       per_dem\z001001x.adf<br>       per_dem\z001002.adf<br>       per_dem\z001002x.adf<br>       per_dem\z001003.adf<br>       per_dem\z001003x.adf<br>       per_dem\z001004.adf<br>       per_dem\z001004x.adf<br>       per_dem\z001005.adf<br>       per_dem\z001005x.adf<br>Size is 52000, 58000<br>Coordinate System is:<br>PROJCS["unnamed",<br>    GEOGCS["Unknown datum based upon the GRS 1980 ellipsoid",<br>        DATUM["Not_specified_based_on_GRS_1980_ellipsoid",<br>            SPHEROID["GRS 1980",6378137,298.257222101,<br>                AUTHORITY["EPSG","7019"]],<br>            AUTHORITY["EPSG","6019"]],<br>        PRIMEM["Greenwich",0,<br>            AUTHORITY["EPSG","8901"]],<br>        UNIT["degree",0.0174532925199433,<br>            AUTHORITY["EPSG","9122"]],<br>        AUTHORITY["EPSG","4019"]],<br>    PROJECTION["Lambert_Conformal_Conic_2SP"],<br>    PARAMETER["standard_parallel_1",38.73333333333333],<br>    PARAMETER["standard_parallel_2",40.03333333333333],<br>    PARAMETER["latitude_of_origin",38],<br>    PARAMETER["central_meridian",-82.5],<br>    PARAMETER["false_easting",1968500],<br>    PARAMETER["false_northing",0],<br>    UNIT["Foot_US",0.3048006096012192]]<br>Origin = (1974999.999999996507540,705000.000000005704351)<br>Pixel Size = (2.499999999999996,-2.499999999999996)<br>Corner Coordinates:<br>Upper Left  ( 1975000.000,  705000.000) ( 82d28'36.55"W, 39d56' 8.22"N)<br>Lower Left  ( 1975000.000,  560000.000) ( 82d28'37.03"W, 39d32'15.15"N)<br>Upper Right ( 2105000.000,  705000.000) ( 82d 0'47.65"W, 39d56' 4.59"N)<br>Lower Right ( 2105000.000,  560000.000) ( 82d 0'57.63"W, 39d32'11.54"N)<br>Center      ( 2040000.000,  632500.000) ( 82d14'44.72"W, 39d44'10.72"N)<br>Band 1 Block=512x4 Type=Float32, ColorInterp=Undefined<br>  Description = per_dem<br>  Min=618.830 Max=1199.330<br>  NoData Value=-3.4028234663852886e+38<br>  Overviews: 13000x14500, 6500x7250, 3250x3625, 1625x1813, 813x907, 407x454, 204<br>x227, 102x114, 51x57<br>  Metadata:<br>    LAYER_TYPE=athematic<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 9, 2018 at 12:53 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Calvin,<br>
<br>
The terrain analysis plugin in 2.18 is QGIS specific code.<br>
<br>
Your observations don't match mine. I tried on a 60,000x80,000 input TIFF of <br>
type Byte, without nodata. Ubuntu 16.04 / gcc 5.5<br>
<br>
The timings I get:<br>
- debug GDAL -dev & QGIS 2.18 builds:<br>
    * gdaldem hillshade: 11m10s<br>
    * terrain analysis plugin: 14m10s<br>
<br>
- optimized GDAL -dev & QGIS 2.18 builds:<br>
     * gdaldem hillshade: 2m35s<br>
     * terrain analysis plugin: 12m00s<br>
<br>
I've compared source codes between GDAL and QGIS. The overall structure is <br>
identical, process line by line, same type of algorithms, but the GDAL code <br>
has a few optimized paths compared to QGIS one (much less trigonometric <br>
computations in critical path), which explains the difference of timings in <br>
optimized builds. But not in the direction you mention.<br>
<br>
can you post the output of gdalinfo on your input DEM ?<br>
<br>
Even<br>
<div class="HOEnZb"><div class="h5"><br>
> I was using the Hillshade from the Terrain Analysis plugin. I don't know if<br>
> it uses GDAL for the processing or not, but it does not give a GDAL command<br>
> line. The Terrain Analysis plugin seems to be dropped for QGIS 3, but if<br>
> the GDAL hillshade is going to be as slow as it was, that is going to be a<br>
> serious problem for me.<br>
> <br>
> Thanks,<br>
> <br>
> Calvin<br>
> <br>
> On Sat, Jun 9, 2018 at 3:01 AM, Nyall Dawson <<a href="mailto:nyall.dawson@gmail.com">nyall.dawson@gmail.com</a>> wrote:<br>
> > On Sat, 9 Jun 2018 at 06:25, C Hamilton <<a href="mailto:adenaculture@gmail.com">adenaculture@gmail.com</a>> wrote:<br>
> > > OS: Windows 7<br>
> > > QGIS 3.0.3<br>
> > > <br>
> > > I have encountered several issues with the Hillshade algorithm in QGIS<br>
> > > 3.<br>
> > > <br>
> > > 1. It is much slower that the QGIS 2 hillshade perhaps as much as 10x<br>
> > <br>
> > slower. This is a problem because I process very large images. 60,000 x<br>
> > 80,000 pixels or greater. I eventually killed the process and went back to<br>
> > QGIS 2 to complete the job.<br>
> > <br>
> > > 2. I tried using Multidirectional shading and received an error. This is<br>
> > <br>
> > basically the command issued.<br>
> > <br>
> > > gdaldem hillshade input.tif output.tif -of GTiff -b 1 -z 1.0 -s 1.0 -az<br>
> > <br>
> > 315.0 -alt 45.0 -multidirectional<br>
> > <br>
> > Can you post the command created in QGIS 2.18 for comparison?<br>
> > <br>
> > > I get the gdal error "ERROR 6: -multidirectional and -az cannot be used<br>
> > <br>
> > together". If multidirectional is enabled then you cannot add the -az<br>
> > parameter. If have tried running the same command in the command line and<br>
> > removed the -az and it worked fine.<br>
> > <br>
> > Confirmed, thanks - fixed by <a href="https://github.com/qgis/QGIS/pull/7209" rel="noreferrer" target="_blank">https://github.com/qgis/QGIS/<wbr>pull/7209</a><br>
> > and unit tests added.<br>
> > <br>
> > Nyall<br>
<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>