<div dir="ltr">Hi,<div><br></div><div>I want to change my Dataset's projection but I have a trouble with it.</div><div><br></div><div>I have custom projection like below:</div><div><br></div><div><div>PROJCS["UTM3 Degree ITRF96 DOM39", </div><div>  GEOGCS["GCS_ITRF_1996", </div><div>    DATUM["D_ITRF_1996", </div><div>      SPHEROID["GRS_1980", 6378137.0, 298.257222101]], </div><div>    PRIMEM["Greenwich", 0.0], </div><div>    UNIT["degree", 0.017453292519943295], </div><div>    AXIS["Longitude", EAST], </div><div>    AXIS["Latitude", NORTH]], </div><div>  PROJECTION["Transverse_Mercator"], </div><div>  PARAMETER["central_meridian", 39.0], </div><div>  PARAMETER["latitude_of_origin", 0.0], </div><div>  PARAMETER["scale_factor", 1.0], </div><div>  PARAMETER["false_easting", 500000.0], </div><div>  PARAMETER["false_northing", 0.0], </div><div>  UNIT["m", 1.0], </div><div>  AXIS["x", EAST], </div><div>  AXIS["y", NORTH], </div><div>  AUTHORITY["EPSG","10103390"]]</div></div><div><br></div><div>and this is how I changed it:</div><div><br></div><div>String test = "PROJCS[\"UTM3 Degree ITRF96 DOM39\",   GEOGCS[\"GCS_ITRF_1996\",     DATUM[\"D_ITRF_1996\",       SPHEROID[\"GRS_1980\", 6378137.0, 298.257222101]],     PRIMEM[\"Greenwich\", 0.0],    UNIT[\"degree\", 0.017453292519943295],    AXIS[\"Longitude\", EAST],    AXIS[\"Latitude\", NORTH]],  PROJECTION[\"Transverse_Mercator\"],  PARAMETER[\"central_meridian\", 39.0],  PARAMETER[\"latitude_of_origin\", 0.0],  PARAMETER[\"scale_factor\", 1.0],  PARAMETER[\"false_easting\", 500000.0],  PARAMETER[\"false_northing\", 0.0],  UNIT[\"m\", 1.0],  AXIS[\"x\", EAST],  AXIS[\"y\", NORTH],  AUTHORITY[\"EPSG\",\"10103390\"]]";<br></div><div><br></div><div>SpatialReference oSRS = new SpatialReference();<br></div><div><div>oSRS.ImportFromWkt(test);</div><div>inputDataset.SetProjection(oSRS.ExportToWkt());</div></div><div><br></div><div>and it returns</div><div><br></div><div><div>PROJCS["unnamed", </div><div>  GEOGCS["GCS Name = GCS_ITRF_1996", </div><div>    DATUM["Datum = D_ITRF_1996", </div><div>      SPHEROID["Ellipsoid = GRS_1980", 6378137.0, 298.257222101]], </div><div>    PRIMEM["Greenwich", 0.0], </div><div>    UNIT["degree", 0.017453292519943295], </div><div>    AXIS["Geodetic longitude", EAST], </div><div>    AXIS["Geodetic latitude", NORTH]], </div><div>  PROJECTION["Transverse_Mercator"], </div><div>  PARAMETER["central_meridian", 39.0], </div><div>  PARAMETER["latitude_of_origin", 0.0], </div><div>  PARAMETER["scale_factor", 1.0], </div><div>  PARAMETER["false_easting", 500000.0], </div><div>  PARAMETER["false_northing", 0.0], </div><div>  UNIT["m", 1.0], </div><div>  AXIS["Easting", EAST], </div><div>  AXIS["Northing", NORTH]]</div></div><div><br></div><div>I have created custom.wkt file and include it in epsg.wkt at "gdal_data" folder but this is also not working for me.</div><div><br></div><div>custom.wkt content:</div><div><br></div><div>10103390,PROJCS["UTM3 Derece ITRF96 DOM39",GEOGCS["GCS_ITRF_1996",DATUM["D_ITRF_1996",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",39.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0], AUTHORITY["EPSG","10103390"]]<br></div><div><br></div><div>Also I tried to set PROJCS like this</div><div><br></div><div>oSRS.SetProjCS("My Proj"); but didn't work.</div><div><br></div><div>Any idea with this?</div><div><br></div><div>Thanks</div><div>Rebo</div></div>