<!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">Hi,</FONT></SPAN>
</P>
<P><SPAN LANG="en-ca"><FONT COLOR="#0000FF" FACE="Arial">After doing some more investigating, it appears g.proj 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). 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: +353 (0)91 492169</FONT></I></SPAN>
<BR><SPAN LANG="en-us"><I><FONT SIZE=2 FACE="Corbel">Fax: +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: </FONT></B> <FONT SIZE=1 FACE="Tahoma">Duffy, Garret </FONT></SPAN>
<BR><SPAN LANG="en-ca"><B><FONT SIZE=1 FACE="Tahoma">Sent: </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: </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: </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. 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. 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). 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 : Oblique Mercator</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">proj : omerc</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">a : 6378137</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">es : 0.006694379990141316</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">lat_0 : 53.25</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">lonc : -5.7</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">alpha : -88</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">k : 1</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">x_0 : 0</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">y_0 : 0</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">towgs84 : 0,0,0,0,0,0,0</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">no_defs : defined</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">-PROJ_UNITS------------------------------------------------</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">unit : meter</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">units : meters</FONT></SPAN>
<BR><SPAN LANG="en-ca"><FONT FACE="Arial">meters : 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"> 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>