<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    <br>

    <div class="moz-cite-prefix">On 12/12/2012 01:54 PM, Knut-Frode

      Dagestad wrote:<br>

    </div>

    <blockquote cite="mid:50C87E89.6080103@hotmail.com" type="cite">The

      GCPs are completely regular, every 40 pixel in each direction.

      Warping quality is very good (only slow), so the GCPs are probably

      also ok.

      <br>

      <br>

      I would appreciate very much if you could download this file and

      test the command below:

      <br>

      <a class="moz-txt-link-freetext" href="http://dl.dropbox.com/u/15885758/testgcp.tif">http://dl.dropbox.com/u/15885758/testgcp.tif</a>

      <br>

      <br>

      time gdalwarp -tps -t_srs '+proj=merc' testgcp.tif out.tif

      <br>

      <br>

      If it works fast for you, the problem should be related to my

      local GDAL build (1.10).

      <br>

      <br>

    </blockquote>

    <font face="Times New Roman, Times, serif">It's far too slow here

      too, so it's not your build.<br>

      <br>

      I see that the warp itself doesn't start but after a very long

      time, so the problem lies in the initial computation of the

      transformation matrix. It must have something to do with the

      distribution of points, perhaps too close together, so as to cause

      singularities. The target coordinates are latlon-values, very

      close numerically. I see two testing possibilities:<br>

      <br>

      1) Try to warp with subsets of points <br>

      2) Multiply all target coordinates with (e.g.) 1000. You'll get an

      impossible georeference, but if this goes much faster, you know

      where the problem lies.<br>

      <br>

      Jan<br>

    </font>

  </body>

</html>