<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 01/02/2013 04:51 PM, Andre Joost
      wrote:<br>
    </div>
    <blockquote cite="mid:kc1l1r$72i$1@ger.gmane.org" type="cite">
      <pre wrap="">

</pre>
      <pre wrap="">
Simpliest way would be to play with the +towgs84 parameters of the 
current EPSG definition.
The first three are shifts of the center of the earth in all three 
dimensions, the next three rotations around these axes, and the last a 
scale factor.

So I would change the values four to six one after the other and look 
for the result. Quantum GIS can handle self-defined proj strings, 
create vector grids in all projections, and show them togehter with 
other grids.


</pre>
    </blockquote>
    <font face="Times New Roman, Times, serif">Thanks Andre, it's in the
      right direction, but it is a bit more complex. The towgs84
      parameters are around the center of the earth. The rotations of
      the historical maps, and the coefficients I computed, are around
      the center of the projection of epsg:28992 (church tower of
      Amersfoort).  <br>
      <br>
      However, I can compute wgs84 latlon values for every point on the
      historical map, by first converting the point to epsg:28992 and
      then to epsg:4326 (wgs84). Do you know of any software that can
      compute the towgs84 parameters based on a set of input coordinates
      and the corresponding wgs84 coordinates?<br>
      <br>
      Jan  </font><br>
  </body>
</html>