EPSG:3821

GEOGCS["TWD67",DATUM["D_Taiwan_Datum_1967",SPHEROID["GRS_1967_Truncated",6378160,298.25]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]

+proj=longlat +ellps=aust_SA +no_defs 

proj4.defs("EPSG:3821","+proj=longlat +ellps=aust_SA +no_defs");

<?xml version="1.0" encoding="UTF-8"?>
 <gml:GeodeticCRS xmlns:epsg="urn:x-ogp:spec:schema-xsd:EPSG:1.0:dataset" xmlns:gml="http://www.opengis.net/gml/3.2" xmlns:xlink="http://www.w3.org/1999/xlink" gml:id="epsg-crs-3821">
   <gml:metaDataProperty>
      <epsg:CommonMetaData>
         <epsg:type>geographic 2D</epsg:type>
         <epsg:informationSource>National Land Surveying and Mapping Center (NLSC), http://www.nlsc.gov.tw</epsg:informationSource>
         <epsg:revisionDate>2008-08-11</epsg:revisionDate>
         <epsg:show>true</epsg:show>
         <epsg:isDeprecated>false</epsg:isDeprecated>
      </epsg:CommonMetaData>
   </gml:metaDataProperty>
   <gml:identifier codeSpace="OGP">urn:ogc:def:crs:EPSG::3821</gml:identifier>
   <gml:name>TWD67</gml:name>
   <gml:remarks>Shares the same origin point with the earlier Hu Tzu Shan system (CRS code 4236) but away from this point coordinates differ. Do not confuse!</gml:remarks>
   <gml:domainOfValidity xlink:href="urn:ogc:def:area:EPSG::3315" />
   <gml:scope>Geodetic survey.</gml:scope>
   <gml:ellipsoidalCS xlink:href="urn:ogc:def:cs:EPSG::6422" />
   <gml:geodeticDatum xlink:href="urn:ogc:def:datum:EPSG::1025" />
</gml:GeodeticCRS>

<?xml version="1.0" encoding="UTF-8"?>
 <gml:GeographicCRS gml:id="ogrcrs38544" xmlns:gml="http://www.opengis.net/gml/3.2">
  <gml:srsName>TWD67</gml:srsName>
  <gml:srsID>
    <gml:name gml:codeSpace="urn:ogc:def:crs:EPSG::">3821</gml:name>
  </gml:srsID>
  <gml:usesEllipsoidalCS>
    <gml:EllipsoidalCS gml:id="ogrcrs38545">
      <gml:csName>ellipsoidal</gml:csName>
      <gml:csID>
        <gml:name gml:codeSpace="urn:ogc:def:cs:EPSG::">6402</gml:name>
      </gml:csID>
      <gml:usesAxis>
        <gml:CoordinateSystemAxis gml:id="ogrcrs38546" gml:uom="urn:ogc:def:uom:EPSG::9102">
          <gml:name>Geodetic latitude</gml:name>
          <gml:axisID>
            <gml:name gml:codeSpace="urn:ogc:def:axis:EPSG::">9901</gml:name>
          </gml:axisID>
          <gml:axisAbbrev>Lat</gml:axisAbbrev>
          <gml:axisDirection>north</gml:axisDirection>
        </gml:CoordinateSystemAxis>
      </gml:usesAxis>
      <gml:usesAxis>
        <gml:CoordinateSystemAxis gml:id="ogrcrs38547" gml:uom="urn:ogc:def:uom:EPSG::9102">
          <gml:name>Geodetic longitude</gml:name>
          <gml:axisID>
            <gml:name gml:codeSpace="urn:ogc:def:axis:EPSG::">9902</gml:name>
          </gml:axisID>
          <gml:axisAbbrev>Lon</gml:axisAbbrev>
          <gml:axisDirection>east</gml:axisDirection>
        </gml:CoordinateSystemAxis>
      </gml:usesAxis>
    </gml:EllipsoidalCS>
  </gml:usesEllipsoidalCS>
  <gml:usesGeodeticDatum>
    <gml:GeodeticDatum gml:id="ogrcrs38548">
      <gml:datumName>Taiwan_Datum_1967</gml:datumName>
      <gml:datumID>
        <gml:name gml:codeSpace="urn:ogc:def:datum:EPSG::">1025</gml:name>
      </gml:datumID>
      <gml:usesPrimeMeridian>
        <gml:PrimeMeridian gml:id="ogrcrs38549">
          <gml:meridianName>Greenwich</gml:meridianName>
          <gml:meridianID>
            <gml:name gml:codeSpace="urn:ogc:def:meridian:EPSG::">8901</gml:name>
          </gml:meridianID>
          <gml:greenwichLongitude>
            <gml:angle gml:uom="urn:ogc:def:uom:EPSG::9102">0</gml:angle>
          </gml:greenwichLongitude>
        </gml:PrimeMeridian>
      </gml:usesPrimeMeridian>
      <gml:usesEllipsoid>
        <gml:Ellipsoid gml:id="ogrcrs38550">
          <gml:ellipsoidName>GRS 1967 Modified</gml:ellipsoidName>
          <gml:ellipsoidID>
            <gml:name gml:codeSpace="urn:ogc:def:ellipsoid:EPSG::">7050</gml:name>
          </gml:ellipsoidID>
          <gml:semiMajorAxis gml:uom="urn:ogc:def:uom:EPSG::9001">6378160</gml:semiMajorAxis>
          <gml:secondDefiningParameter>
            <gml:inverseFlattening gml:uom="urn:ogc:def:uom:EPSG::9201">298.25</gml:inverseFlattening>
          </gml:secondDefiningParameter>
        </gml:Ellipsoid>
      </gml:usesEllipsoid>
    </gml:GeodeticDatum>
  </gml:usesGeodeticDatum>
</gml:GeographicCRS>

3821=GEOGCS["TWD67",DATUM["Taiwan_Datum_1967",SPHEROID["GRS 1967 Modified",6378160,298.25,AUTHORITY["EPSG","7050"]],AUTHORITY["EPSG","1025"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","3821"]]

PROJECTION
	"proj=longlat"
	"ellps=aust_SA"
	"no_defs"
END

wkt = '''GEOGCS["TWD67",DATUM["Taiwan_Datum_1967",SPHEROID["GRS 1967 Modified",6378160,298.25,AUTHORITY["EPSG","7050"]],AUTHORITY["EPSG","1025"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","3821"]]'''
m = mapObj('')
m.setWKTProjection(ref.ExportToWkt())
lyr = layerObj(m)
lyr.setWKTProjection(ref.ExportToWkt())

<?xml version="1.0" encoding="utf-8"?>
<Map srs="+proj=longlat +ellps=aust_SA +no_defs">
	<Layer srs="+proj=longlat +ellps=aust_SA +no_defs">
	</Layer>
</Map>

proj4 = '+proj=longlat +ellps=aust_SA +no_defs'
m = Map(256,256,proj4)
lyr = Layer('Name',proj4)

INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 3821, 'EPSG', 3821, '+proj=longlat +ellps=aust_SA +no_defs ', 'GEOGCS["TWD67",DATUM["Taiwan_Datum_1967",SPHEROID["GRS 1967 Modified",6378160,298.25,AUTHORITY["EPSG","7050"]],AUTHORITY["EPSG","1025"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","3821"]]');

GEOGCS["TWD67",DATUM["Taiwan_Datum_1967",SPHEROID["GRS 1967 Modified",6378160,298.25,AUTHORITY["EPSG","7050"]],AUTHORITY["EPSG","1025"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","3821"]]