<p>Hi Frank,<br />
<br />
I&nbsp;tested the program below but it does not out put the .tab or
.mif/.mid with the coordinate system.&nbsp; It is always non-earth.&nbsp;
Could you please check it against your gdal/ogr&nbsp;installation to see
if something wrong&nbsp;happened during my compilation of gdal 1.6.2?<br
/>
<br />
Thanks a lot,<br />
<br />
Reinaldo<br />
<br />
*******************************<br />
#include &quot;gdal.h&quot;<br />
#include &quot;ogr_spatialref.h&quot;<br />
#include &quot;ogrsf_frmts.h&quot;<br />
#include &quot;stdio.h&quot;</p>
<p>double Lat[3], Lon[3];<br />
int n=0;&nbsp;char *pszWKT = NULL;</p>
<p>void main()<br />
{<br />
&nbsp;OGRSpatialReference oSRS;<br />
//&nbsp;oSRS.SetWellKnownGeogCS( &quot;EPSG:4326&quot; );<br />
&nbsp;oSRS.importFromMICoordSys( &quot;CoordSys Earth Projection 1,
104&quot; );</p>
<p>&nbsp;const char *pszDriverName = &quot;MapInfo File&quot; ;</p>
<p>&nbsp;OGRSFDriver *poDriver;<br />
&nbsp;&nbsp;&nbsp; OGRRegisterAll();</p>
<p>&nbsp;&nbsp;&nbsp; poDriver =
OGRSFDriverRegistrar::GetRegistrar()-&gt;GetDriverByName(<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
pszDriverName );<br />
&nbsp;&nbsp;&nbsp; if( poDriver == NULL )<br />
&nbsp;&nbsp;&nbsp; {<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf( &quot;%s driver not
available.\n&quot;, pszDriverName );<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit( 1 );<br />
&nbsp;&nbsp;&nbsp; }</p>
<p>&nbsp;OGRDataSource *poDS;<br />
&nbsp;poDS = poDriver-&gt;CreateDataSource( &quot;c:/test01.mif&quot;,
NULL );<br />
//&nbsp;poDS = poDriver-&gt;CreateDataSource( &quot;c:/test01.tab&quot;,
NULL );<br />
&nbsp;&nbsp;&nbsp; if( poDS == NULL )<br />
&nbsp;&nbsp;&nbsp; {<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf( &quot;Creation of
output file failed.\n&quot; );<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit( 1 );<br />
&nbsp;&nbsp;&nbsp; }</p>
<p>&nbsp;&nbsp;&nbsp; OGRLayer *poLayer;<br />
&nbsp;&nbsp;&nbsp; poLayer = poDS-&gt;CreateLayer( &quot;test01&quot;,
NULL, wkbPoint, NULL );<br />
&nbsp;&nbsp;&nbsp; if( poLayer == NULL )<br />
&nbsp;&nbsp;&nbsp; {<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf( &quot;Layer creation
failed.\n&quot; );<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit( 1 );<br />
&nbsp;&nbsp;&nbsp; }</p>
<p>&nbsp;OGRFieldDefn oField0( &quot;NAME1&quot;, OFTString );<br />
&nbsp;OGRFieldDefn oField1( &quot;NAME2&quot;, OFTString );</p>
<p>&nbsp;oField0.SetWidth(32);<br />
&nbsp;&nbsp;&nbsp; oField1.SetWidth(32);</p>
<p>&nbsp;if( poLayer-&gt;CreateField( &amp;oField0 ) != OGRERR_NONE )<br />
&nbsp;&nbsp;&nbsp; {<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf( &quot;Creating Name
field failed.\n&quot; );<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit( 1 );<br />
&nbsp;&nbsp;&nbsp; }</p>
<p>&nbsp;if( poLayer-&gt;CreateField( &amp;oField1 ) != OGRERR_NONE )<br />
&nbsp;&nbsp;&nbsp; {<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf( &quot;Creating Name
field failed.\n&quot; );<br />
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit( 1 );<br />
&nbsp;&nbsp;&nbsp; }<br />
<br />
*******************************<br />
.mif output:<br />
<br />
Version 300<br />
Charset &quot;Neutral&quot;<br />
Delimiter &quot;,&quot;<br />
Columns 2<br />
&nbsp; NAME1 Char(32)<br />
&nbsp; NAME2 Char(32)<br />
Data</p>
<p>Point -45.13975272 -22.69426671<br />
&nbsp;&nbsp;&nbsp; Symbol (35,0,12)<br />
Point -45.86066928 -23.21117775<br />
&nbsp;&nbsp;&nbsp; Symbol (35,0,12)<br />
Point -44.33023872 -22.68505859375<br />
&nbsp;&nbsp;&nbsp; Symbol (35,0,12)<br />
<br />
************************************<br />
.mid output:<br />
<br />
&quot;TEST01&quot;,&quot;TEST02&quot;<br />
&quot;TEST01&quot;,&quot;TEST02&quot;<br />
&quot;TEST01&quot;,&quot;TEST02&quot;<br />
<br />
&gt; Reinaldo escada Chohfi wrote: <br />
&gt;&gt; Hello List, <br />
&gt;&gt; <br />
&gt;&gt; I am still trying to get a x,y (lon/Lat) point file in MapInfo
format <br />
&gt;&gt; with <br />
&gt;&gt; coordinate information. <br />
&gt;&gt; <br />
&gt;&gt; I get the correct MapInfo files out (.tab or .mif) with the x,y
points <br />
&gt;&gt; and <br />
&gt;&gt; attributes, but without the CoordSys information. <br />
&gt;&gt; The output file is always Non-Earth. <br />
&gt;&gt; <br />
&gt;&gt; I've tried everything: <br />
&gt;&gt; <br />
&gt;&gt; OGRSpatialReference oSRS; <br />
&gt;&gt; oSRS.SetWellKnownGeogCS(&quot;EPSG:4326&quot;); <br />
&gt;&gt; oSRS.importFromMICoordSys( &quot;CoordSys Earth Projection 1,
104&quot; ); <br />
&gt;&gt; <br />
&gt;&gt; <br />
&gt;&gt; Could someone PLEASE let me know what to do? <br />
&gt; <br />
&gt; Reinaldo, <br />
&gt; <br />
&gt; I just tried: <br />
&gt; <br />
&gt; ogr2ogr -skipfailures -f &quot;MapInfo File&quot; out.mif
GB4X0000.000 TUNNEL <br />
&gt; <br />
&gt; and this produced a mapinfo mid/mif fileset with &quot;CoordSys Earth
<br />
&gt; Projection 1, 104&quot; as the coord sys line. I also tested a small
<br />
&gt; program: <br />
&gt; <br />
&gt; <br />
&gt; #include &quot;gdal.h&quot; <br />
&gt; #include &quot;ogr_spatialref.h&quot; <br />
&gt; <br />
&gt; int main() <br />
&gt; <br />
&gt; { <br />
&gt; OGRSpatialReference oSRS; <br />
&gt; <br />
&gt; oSRS.importFromMICoordSys( &quot;CoordSys Earth Projection 1,
104&quot; ); <br />
&gt; } <br />
&gt; <br />
&gt; and it appeared to correctly parse the Coordsys into a WGS84 srs. <br />
&gt; <br />
&gt; I even tried: <br />
&gt; <br />
&gt; ogr2ogr -a_srs wgs84.wkt -skipfailures -f &quot;MapInfo File&quot;
out.mif <br />
&gt; data/s57/GB4X0000.000 TUNNEL <br />
&gt; <br />
&gt; with wgs84.wkt holding: <br />
&gt; <br />
&gt;
GEOGCS[&quot;unnamed&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS
<br />
&gt;
84&quot;,6378137,298.257223563],TOWGS84[0,0,0,-0,-0,-0,0]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433]]
<br />
&gt; <br />
&gt; and things seemed to work. <br />
&gt; <br />
&gt; Can you provide a more complete example that is failing for you? <br />
&gt; <br />
&gt; Best regards, <br />
&gt; -- <br />
&gt;
---------------------------------------+--------------------------------------
<br />
&gt; I set the clouds in motion - turn up | Frank Warmerdam, <br />
&gt; warmerdam@pobox.com <br />
&gt; light and sound - activate the windows | http://pobox.com/~warmerdam
<br />
&gt; and watch the world go round - Rush | Geospatial Programmer for Rent
<br />
&gt; <br />
&gt; <br />
&nbsp;</p>