<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7654.12">
<TITLE>RE: reprojection issue, using grass6.2</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/rtf format -->

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">Update, I replicated the correct projection by appending</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">no_uoff: defined</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">rot_conv: defined</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">to the PROJ_INFO file</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">as per advice here</FONT></SPAN>

<BR><SPAN LANG="en-ca"></SPAN><A HREF="http://osgeo-org.1803224.n2.nabble.com/GRASSLIST-8028-Oblique-Mercator-Projection-in-GRASS-td1874326.html"><SPAN LANG="en-ca"><U></U><U><FONT COLOR="#0000FF" FACE="Arial">http://osgeo-org.1803224.n2.nabble.com/GRASSLIST-8028-Oblique-Mercator-Projection-in-GRASS-td1874326.html</FONT></U></SPAN></A><SPAN LANG="en-ca"></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">maybe my version (6.2.1 (2006)) of g.proj is old (6.2.1 (2006)) but it for sure didn't work as it should in this instance.</FONT></SPAN></P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">Garret</FONT></SPAN>
</P>
<BR>
<BR>

<P><B><SPAN LANG="en-us"><FONT FACE="Corbel">Dr. Garret Duffy</FONT></SPAN></B>

<BR><B><SPAN LANG="en-us"><FONT FACE="Corbel">______________________</FONT></SPAN></B>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Department of Earth and Ocean Sciences</FONT></I></SPAN>
</P>

<P><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">National University of Ireland, Galway</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">University Road</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Galway</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">IRELAND</FONT></I></SPAN>
</P>

<P><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Office:&nbsp; +353 (0)91 492169</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Fax:&nbsp; +353 (0)91 750533</FONT></I></SPAN><I><SPAN LANG="en-ca"></SPAN></I>
</P>
<BR>
<BR>

<P><SPAN LANG="en-ca"><FONT SIZE=1 FACE="Tahoma">_____________________________________________ </FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">From: &nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">Duffy, Garret&nbsp; </FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">Sent:&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">12-Oct-11 2:46 PM</FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">To:&nbsp;&nbsp;&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">'grass-user@lists.osgeo.org'</FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">Subject:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">RE: reprojection issue, using grass6.2</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">Hi,</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">After doing some more investigating, it appears g.proj&nbsp; is not reading all the parameters in :</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><B><FONT FACE="Arial">proj=omerc +ellps=WGS84 +k_0=1.0 +lat_0=53.25 +no_uoff +rot_conv +alpha=-88 +lonc=-5.7</FONT></B></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">For example, the above command projects (-5.7d, 53.25d) correctly to (X,Y)=(-5901478, 7901021).&nbsp; But after g.proj, the same point projects to (X,Y)=(0,0).</FONT></SPAN></P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Can someone help with this?</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Thanks,</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Garret</FONT></SPAN>
</P>
<BR>

<P><B><SPAN LANG="en-us"><FONT FACE="Corbel">Dr. Garret Duffy</FONT></SPAN></B>

<BR><B><SPAN LANG="en-us"><FONT FACE="Corbel">______________________</FONT></SPAN></B>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Department of Earth and Ocean Sciences</FONT></I></SPAN>
</P>

<P><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">National University of Ireland, Galway</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">University Road</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Galway</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">IRELAND</FONT></I></SPAN>
</P>

<P><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Office:&nbsp; +353 (0)91 492169</FONT></I></SPAN>

<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Fax:&nbsp; +353 (0)91 750533</FONT></I></SPAN><I><SPAN LANG="en-ca"></SPAN></I>
</P>
<BR>
<BR>

<P><SPAN LANG="en-ca"><FONT SIZE=1 FACE="Tahoma">_____________________________________________ </FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">From: &nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">Duffy, Garret&nbsp; </FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">Sent:&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">12-Oct-11 12:29 PM</FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">To:&nbsp;&nbsp;&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">'grass-user@lists.osgeo.org'</FONT></SPAN>

<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">Subject:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT></B> <FONT SIZE=1 FACE="Tahoma">reprojection issue, using grass6.2</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Hi,</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">I have three locations defined, one in Geographic 'projection', one UTM and the other described below.&nbsp; When I try to reproject data, e.g. a vector rectangle from either Geographic or UTM to the Oblique Mercator projection defined below, the rectangle isn't rotated as I expect.&nbsp; The origin of the projected coordinates looks to be in the correct location but there is no rotation applied (geographic north is supposed to ~45 deg from the y axis).&nbsp; BTW, is it even possible to display a geodetic grid for a non-rectified projection?</FONT></SPAN></P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">output from g.proj -p:</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">-PROJ_INFO-------------------------------------------------</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">name&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : Oblique Mercator</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">proj&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : omerc</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">a&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 6378137</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">es&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 0.006694379990141316</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">lat_0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 53.25</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">lonc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : -5.7</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">alpha&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : -88</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">k&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 1</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">x_0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 0</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">y_0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 0</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">towgs84&nbsp;&nbsp;&nbsp; : 0,0,0,0,0,0,0</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">no_defs&nbsp;&nbsp;&nbsp; : defined</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">-PROJ_UNITS------------------------------------------------</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">unit&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : meter</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">units&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : meters</FONT></SPAN>

<BR><SPAN LANG="en-ca"><FONT FACE="Arial">meters&nbsp;&nbsp;&nbsp;&nbsp; : 1</FONT></SPAN>
</P>
<BR>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">The command I used to specify the oblique projection is:</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><B><FONT FACE="Arial">proj +proj=omerc +ellps=WGS84 +k_0=1.0 +lat_0=53.25 +no_uoff +rot_conv +alpha=-88 +lonc=-5.7</FONT></B></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">I must emphasize that I have checked the proj.4 parameters using the command-line proj.4 and the *</FONT><B><FONT FACE="Arial">parameters are correct*.</FONT></B><FONT FACE="Arial">&nbsp; I applied these parameters in GlobalMapper to get the final map to the state that I would like eventually in GRASS.</FONT></SPAN></P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Any advice would be *greatly* appreciated, I have spent a lot of time on this.</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Regards,</FONT></SPAN>
</P>

<P><SPAN LANG="en-ca"><FONT FACE="Arial">Garret</FONT><B></B><B><I></I></B></SPAN>
</P>

</BODY>
</HTML>