<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=us-ascii">
<META content="MSHTML 6.00.2900.3243" name=GENERATOR></HEAD>
<BODY>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008>Hi,</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN class=562272910-21012008>I&nbsp;deal 
with&nbsp;climate model data and have regional model data defined on a long/lat 
grid with the poles rotated away from their true position (in order to keep the 
pole singularities away from the region of interest).&nbsp; Typically this grid 
would be defined as follows:</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008>grid_north_pole_latitude:&nbsp; 39.2 deg (true latitude 
of&nbsp;rotated north pole)</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008>grid_north_pole_longitude: 198.0 deg (true longitude of 
rotated north pole)</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008>north_pole_grid_longitude: 0 deg (longitude of the 
prime meridian in the rotated coordinate system)</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN class=562272910-21012008>How would I 
represent this coordinate system with proj4 parameters?&nbsp; I have been trying 
to use +towgs84:</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN class=562272910-21012008>+proj=longlat 
+towgs84=0,0,0,-grid_north_pole_longitude*60*60,90-grid_north_pole_latitude*60*60,north_pole_grid_longitude,0</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN class=562272910-21012008>So far this doesn't 
seem to work (I'm using gdal to transform an image to EPSG:4326).&nbsp; Am I on 
the right lines?&nbsp; Standard GIS terminology such as datum, spheroid etc. are 
unfortunately not part of the standard vocabulary of climate 
modellers.</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV><FONT face=Arial size=2><SPAN class=562272910-21012008>Many 
thanks,</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008>Stephen.</SPAN></FONT></DIV>
<DIV><FONT face=Arial size=2><SPAN 
class=562272910-21012008></SPAN></FONT>&nbsp;</DIV>
<DIV align=left><FONT face=Arial size=2>---</FONT></DIV>
<DIV align=left><FONT face=Arial size=2>Stephen Pascoe&nbsp; +44 (0)1235 
445980</FONT></DIV>
<DIV align=left><FONT face=Arial size=2>British Atmospheric Data 
Centre</FONT></DIV>
<DIV align=left><FONT face=Arial size=2>Rutherford Appleton 
Laboratory</FONT></DIV>
<DIV>&nbsp;</DIV></BODY></HTML>