<div dir="ltr">Thanks for your reply.<br><br>I'm using gdal 1.10.1.  I'm using cubic resampling, but i have noticed the dark line artifact with bilinear as well.  I unfortunately do not reproduce the problem with a basic gdalwarp with cubic resampling.<br><br>I have written this code to do some registration. It's hard to tell in the example, but my transformation actually pushed pixels up a bit. The transformation I use handles interior pixels nicely and extends beyond the image, so there is nothing 'weird' going on with that near the boundary that I can think of.<br><br>I was hoping there is a warp option i have not yet thought of.<br><br>Thanks again,<br>Mike<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 16, 2015 at 11:10 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">Le jeudi 16 avril 2015 19:01:57, Michael Aschenbeck a écrit :<br>
> Hello,<br>
><br>
> I'm using ChunkAndWarpMulti to warp an image.  The warping is working<br>
> nicely, however, at the boundary I seem to be getting some artifacts.<br>
><br>
> The first artifact i see is a DARK boundary of pixels in some location.  My<br>
> guess is that the interpolator is interpolating with blackfill (intensity<br>
> zero pixels).  Note that I am using the following setup:<br>
> psWarpOptions->padfSrcNoDataReal = (double *)<br>
> CPLMalloc(psWarpOptions->nBandCount*sizeof(double));<br>
> for (int ii = 0; ii < psWarpOptions->nBandCount; ii++)<br>
> {<br>
> psWarpOptions->padfSrcNoDataReal[ii] = 0;<br>
> }<br>
> which i thought was supposed to treat zeros as nodata.  It doesn't seem to<br>
> be doing what I think it should.<br>
><br>
> In some cases, I am also noticing a BRIGHT band strip of pixels adjacent to<br>
> the dark boundary pixels.  I don't have any thoughts on where this artifact<br>
> is coming from.<br>
><br>
> Below you can find links to an example.  The orange you see is just the<br>
> background color of my viewer.  Zero pixels are set to transparent so you<br>
> can see the dark boundary artifact.  (Note that all of the orange region is<br>
> covered with zero-intensity pixels, so we haven't reached the boundary of<br>
> the file, just the boundary of the non-zero pixels.)  Sorry if that's<br>
> confusing.<br>
> The before image: <a href="https://www.dropbox.com/s/3kmyi8wu0qybsq9/before.JPG?dl=0" target="_blank">https://www.dropbox.com/s/3kmyi8wu0qybsq9/before.JPG?dl=0</a><br>
> The after image with the artifacts:<br>
> <a href="https://www.dropbox.com/s/fo0m8q95b26s61m/after.JPG?dl=0" target="_blank">https://www.dropbox.com/s/fo0m8q95b26s61m/after.JPG?dl=0</a><br>
><br>
<br>
</div></div>Mike,<br>
<br>
Which gdal version ? Which resampling method ?<br>
Do you reproduce with gdalwarp ?<br>
If so, providing input file + full gdalwarp command line would help.<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>