<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html; charset=windows-1252"
 http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Hello list,<br>
I posted a question 3 month ago and lost track of it and found that
Daniel Morissette answered practically the same day. So first of all
thank you Daniel, your help is greatly approciated, sorry for the
silence: second of all I pick up were I left off:<br>
<br>
PROBLEM:<br>
I have a mapfile with layers which use two projections:<br>
<b>main project projection</b><br>
  PROJECTION<br>
    "proj=utm"<br>
    "zone=32"<br>
    "ellps=WGS84"<br>
    "datum=WGS84"<br>
    "units=m"<br>
    "no_defs"<br>
  END<br>
<b>layer projection</b><br>
    PROJECTION<br>
      "proj=tmerc"<br>
      "lat_0=0"<br>
      "lon_0=9"<br>
      "k=0.9996"<br>
      "x_0=1500000"<br>
      "y_0=0"<br>
      "ellps=intl"<br>
      "units=m"<br>
      "no_defs"<br>
      "nadgrids=lombardia_tot.gsb,null"<br>
    END<br>
<br>
Both do <b>not </b>use the epsg file. the layer is reprojected
"perfectly" (sorry for the non-carthographic emotional term) in the
correct project coordinates. Next step was to zoom to an object in the
layer so I calculated the minum bounding box of the object in the layer
and reprojected it to project coordinates so that I could zoom to the
map correctly. This works but does not consider the "nadgrids"
parameter and thus has quite an error (30 m). Daniel suggested that the
PROJ_LIB configure parameter might be different in the map and in the
script. I checked with a mapobj->configureoptions  and got back an
empty array in the script and in the mapfile. In my environment
variables PROJ_LIB is set to the correct folder contiaining the proj
epsg & co. files. <br>
I will keep trying and testing and let you know. It is quite strange,
almost like mapscript does not consider the "nadgrids" parameter in the
projection object. Will keep testing.<br>
Cheers,<br>
<div class="moz-signature">-- <br>
<font face="Arial, Helvetica, sans-serif" size="-1"><b>Francesco Pirotti</b>
<br>
Dep. TESAF <br>
CIRGEO – Interdepartmental Research Center on Cartography
Photogrammetry <br>
Remote Sensing and G.I.S. <br>
University of Padova <br>
Web: <a href="http://www.cirgeo.unipd.it/cirgeo/francescopirotti.htm">http://www.cirgeo.unipd.it/cirgeo/francescopirotti.htm</a>
<br>
Email: <a href="mailto:francesco.pirotti@unipd.it">francesco.pirotti@unipd.it</a>
<br>
Phone: +39 049 827 2710 <br>
Phone: +39 349 55 39 261
<p><br>
</p>
</font></div>
</body>
</html>