Morning all, <br><br>I am struggling to get a coordinate transformation
to work in a python script. Below is the code snippet that is causing
the following error:<br><br>"ERROR 6: No translation for 'Azimuthal_Equidistant' to PROJ.4 format is known."<br>
<br>I am pulling the input spatial reference from the input shapefile
(Mars_2000), defining the custom projection (custom in that it is
centered on a point), and then trying to transform to the custom
projection.<br><br><blockquote style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;" class="gmail_quote">
<font size="1">Create a new shapefile and copy, one geometry at a time, from the input to the new output<br>...<br>#Define the new projection to ortho centered on the crater<br> azimuthal = "PROJCS['Azimuthal_Equidistant',GEOGCS['GCS_Mars_2000',DATUM['D_Mars_2000',SPHEROID['Mars_2000_IAU_IAG',3396190.0,169.8944472236118]],PRIMEM['Reference_Meridian',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Azimuthal_Equidistant'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',%f],PARAMETER['Latitude_Of_Origin',%f],UNIT['Meter',1.0]]" %(x,y)<br>
<br> intermediate_SR.ImportFromWkt(azimuthal)<br> <br> #Transform the point from the spatial input reference to an azimuthal projection centered on the point<br> <span style="color: rgb(255, 0, 0);">outputGeom.TransformTo(intermediate_SR)</span><br>
...<br></font></blockquote>Do I need to define the coordinate system differently or am I making an error somewhere else?<br><br>Thanks!<br>Jay<div style="visibility: hidden; display: inline;" id="avg_ls_inline_popup"></div>
<style type="text/css">#avg_ls_inline_popup { position:absolute; z-index:9999; padding: 0px 0px; margin-left: 0px; margin-top: 0px; width: 240px; overflow: hidden; word-wrap: break-word; color: black; font-size: 10px; text-align: left; line-height: 13px;}</style>