<div dir="ltr"><div class="gmail_extra"><div class="gmail_extra" style="font-size:13px">Thank you for your response Even.  </div><div class="gmail_extra" style="font-size:13px"><br></div><div class="gmail_extra" style="font-size:13px">Unfortunately when I split the chart down the 180 degree longitude line, GDALSuggestedWarpOutput() still produces smaller chart dimensions than the dimensions of the inputs.  I think this would rule out any wrapping problem because the split pieces exhibited the same behavior despite never crossing the dateline.  </div><div class="gmail_extra" style="font-size:13px"><br></div><div class="gmail_extra" style="font-size:13px">I am going to try and step through the GDALSuggestedWarpOutput() function and examine it's logic more closely.</div><div class="gmail_extra" style="font-size:13px"><br></div><div class="gmail_extra" style="font-size:13px">-- Joe</div><div><br></div></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 10, 2014 at 5:01 AM, 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"><div id=":4bv" class="" style="overflow:hidden">"joem",<br>
<br>
It seems that your input dataset is very close to the dateline, and very<br>
likely crossing it. I'm not that surprised that GDALSuggestedWarpOutput() has<br>
issues at that discontinuity. But I think it might have indeed operated<br>
correctly. If the extent of the dataset crosses the dateline the resulting<br>
eastings in Web Mercator coordinates will go from about -20 million to about<br>
20 million. A huge width. GDALSuggestedWarpOutput() tries to scale the<br>
computed dimensions so that the result has roughly the same number of pixels<br>
as the input. So the huge width becomes a reasonable width, and the reasonable<br>
height becomes a small height.<br>
In that situation, you could try using modified input and output proj.4 strings<br>
by adding +over to the normal proj.4 definition. See<br>
<a href="https://trac.osgeo.org/proj/wiki/GenParms" target="_blank">https://trac.osgeo.org/proj/wiki/GenParms</a><br>
This should produce an output datasets with eastings that go beyond the 20<br>
million limit. You might want to cut it into 2 parts afterwards.<br>
<br>
Best regards,<br>
<br>
Even<br>
<br>
><br>
><br>
> I am working on a tool to project FAA charts from Lambert Conformal Conic<br>
> to Web Mercator (EPSG 900913).  It seems to be producing correct output<br>
> except for a FAA chart of the East Aleutian Islands.    When I traced the<br>
> problem in my code, it appears that GDALSuggestedWarpOutput() is giving me<br>
> a width and height far too small given the input tiff had larger<br>
> dimensions.  The source GeoTiff has dimensions of 15555x5786 but<br>
> GDALSuggestedWarpOutput() gives a suggested width of 17050  and a height<br>
> of 178.  A height of 178 pixels surely has to be wrong.    I have added<br>
> the well-known text and proj.4  output from gdalinfo.exe of the input tiff<br>
> below.   Although I am not entirely sure if I should, I have been using<br>
> the values supplied by GDALSuggestedWarpOutput() for creating my output<br>
> Geotiff as if the output will always be valid.<br>
><br>
> * Why is GDALSuggestedWarpOutput() giving bad values for the width and<br>
> height of the output tiff?<br>
> * Is this incorrect usage of GDALSuggestedWarpOutput() or is it a bug?<br>
><br>
><br>
> PROJ.4 : '+proj=lcc +lat_1=54.66666666666666 +lat_2=49.33333333333334<br>
> +lat_0=52.16666666666666 +lon_0=-177.5 +x_0=0 +y_0=0 +datum=NAD83 +units=m<br>
> +no_defs '<br>
><br>
> OGC WKT :<br>
> PROJCS["unnamed",<br>
>     GEOGCS["NAD83",<br>
>         DATUM["North_American_Datum_1983",<br>
>             SPHEROID["GRS 1980",6378137,298.2572221010002,<br>
>                 AUTHORITY["EPSG","7019"]],<br>
>             AUTHORITY["EPSG","6269"]],<br>
>         PRIMEM["Greenwich",0],<br>
>         UNIT["degree",0.0174532925199433],<br>
>         AUTHORITY["EPSG","4269"]],<br>
>     PROJECTION["Lambert_Conformal_Conic_2SP"],<br>
>     PARAMETER["standard_parallel_1",54.66666666666666],<br>
>     PARAMETER["standard_parallel_2",49.33333333333334],<br>
>     PARAMETER["latitude_of_origin",52.16666666666666],<br>
>     PARAMETER["central_meridian",-177.5],<br>
>     PARAMETER["false_easting",0],<br>
>     PARAMETER["false_northing",0],<br>
>     UNIT["metre",1,<br>
>         AUTHORITY["EPSG","9001"]]]<br>
><br>
><br>
><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></div></blockquote></div><br><br></div></div>