<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><meta name="ProgId" content="Word.Document"><meta name="Generator" content="Microsoft Word 11"><meta name="Originator" content="Microsoft Word 11"><link rel="File-List" href="file:///C:%5CDOCUME%7E1%5CADMINI%7E1%5CLOCALS%7E1%5CTemp%5Cmsohtml1%5C01%5Cclip_filelist.xml"><style>
<!--
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{mso-style-parent:"";
margin:0mm;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman";
mso-fareast-font-family:"Times New Roman";
mso-bidi-font-family:"Times New Roman";
mso-bidi-language:AR-SA;}
@page Section1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;
mso-header-margin:36.0pt;
mso-footer-margin:36.0pt;
mso-paper-source:0;}
div.Section1
{page:Section1;}
-->
</style>
<p class="MsoNormal">Hi,</p><p class="MsoNormal">Thank you for the support</p>
<p class="MsoNormal">I couldn’t solve the issue yet</p>
<p class="MsoNormal">Please note that I have <span style=""> </span>redone the work with updated SL99 information.
Still I have the same issue.</p>
<p class="MsoNormal">Earlier I have used projection parameters in the ogr2ogr
command and this time</p>
<p class="MsoNormal">the WKT based on a projection file which I am using regularly
is used.</p>
<p class="MsoNormal">Detailed information is as follows</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">The content of the file :</p>
<p class="MsoNormal">SL99 WKT</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">PROJCS["SL_grid_99",</p>
<p class="MsoNormal">GEOGCS["GCS_Kandawala",</p>
<p class="MsoNormal">DATUM["D_Kandawala",</p>
<p class="MsoNormal">SPHEROID["Everest_1830",6377276.345,300.8017]],</p>
<p class="MsoNormal">PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],</p>
<p class="MsoNormal">PROJECTION["Transverse_Mercator"],</p>
<p class="MsoNormal">PARAMETER["False_Easting",200000.000],</p>
<p class="MsoNormal">PARAMETER["False_Northing",200000.00],</p>
<p class="MsoNormal">PARAMETER["Central_Meridian",80.77171],</p>
<p class="MsoNormal">PARAMETER["Scale_Factor",0.9999238418],</p>
<p class="MsoNormal">PARAMETER["Latitude_Of_Origin",7.00048],</p>
<p class="MsoNormal">UNIT["Meter",1]]</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Following command is used to project the vector data</p>
<p class="MsoNormal">ogr2ogr -a_srs slgrid99.wkt F:\UDAdata\test\cmc_b_sl99.shp F:\UDAdata\test\cmc_bound.shp</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Raster data’s original projection information is as follows
( from gdalinfo)</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">PROJCS["UTM Zone 44, Northern Hemisphere",</p>
<p class="MsoNormal"><span style=""> </span>GEOGCS["WGS_1984",</p>
<p class="MsoNormal"><span style=""> </span>DATUM["WGS_1984",</p>
<p class="MsoNormal"><span style=""> </span>SPHEROID["WGS
84",6378137,298.2572235630016],</p>
<p class="MsoNormal"><span style=""> </span>TOWGS84[0,0,0,0,0,0,0]],</p>
<p class="MsoNormal"><span style=""> </span>PRIMEM["Greenwich",0],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["degree",0.0174532925199433],</p>
<p class="MsoNormal"><span style=""> </span>AUTHORITY["EPSG","4326"]],</p>
<p class="MsoNormal"><span style=""> </span>PROJECTION["Transverse_Mercator"],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["latitude_of_origin",0],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["central_meridian",81],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["scale_factor",0.9996],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["false_easting",500000],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["false_northing",0],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["Meter",1],</p>
<p class="MsoNormal"><span style=""> </span>AUTHORITY["EPSG","32644"]]</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Then I have used the following command to re-project raster
data.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">gdalwarp -t_srs slgrid99.wkt -s_srs "+proj=utm +zone=44
+ellps=WGS84 +datum=WGS84 +units=m" F:\UDAdata\po_310491_rgb_0000000_Clip.img
F:\</p>
<p class="MsoNormal">UDAdata\test\raster_sl99.tiff</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Then , when over layering I use the following custom
projection in QGIS. It worked fine with different vector data in the local
context. This is the first time I’m using raster data in my work.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">+proj=tmerc +lat_0=7.00047 +lon_0=80.77171 +x_0=200000
+y_0=200000 +ellps=evrst30 +to_meter=0.999924 +units=m</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">The shift is around 200m.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">WKT of the re projected files are as below</p>
<p class="MsoNormal">Vector</p>
<p class="MsoNormal">PROJCS["SL_grid_99",</p>
<p class="MsoNormal"><span style=""> </span>GEOGCS["GCS_Kandawala",</p>
<p class="MsoNormal"><span style=""> </span>DATUM["Kandawala",</p>
<p class="MsoNormal"><span style=""> </span>SPHEROID["Everest_1830",6377276.345,300.8017]],</p>
<p class="MsoNormal"><span style=""> </span>PRIMEM["Greenwich",0],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["Degree",0.017453292519943295]],</p>
<p class="MsoNormal"><span style=""> </span>PROJECTION["Transverse_Mercator"],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["False_Easting",200000.000],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["False_Northing",200000.00],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["Central_Meridian",80.77171],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["Scale_Factor",0.9999238418],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["Latitude_Of_Origin",7.00048],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["Meter",1]]</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Raster</p>
<p class="MsoNormal">PROJCS["SL_grid_99",</p>
<p class="MsoNormal"><span style=""> </span>GEOGCS["GCS_Kandawala",</p>
<p class="MsoNormal"><span style=""> </span>DATUM["unknown",</p>
<p class="MsoNormal"><span style=""> </span>SPHEROID["unnamed",6377276.345,300.8017000000115]],</p>
<p class="MsoNormal"><span style=""> </span>PRIMEM["Greenwich",0],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["degree",0.0174532925199433]],</p>
<p class="MsoNormal"><span style=""> </span><span style=""> </span>PROJECTION["Transverse_Mercator"],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["latitude_of_origin",7.00048],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["central_meridian",80.77171],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["scale_factor",0.9999238418],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["false_easting",200000],</p>
<p class="MsoNormal"><span style=""> </span>PARAMETER["false_northing",200000],</p>
<p class="MsoNormal"><span style=""> </span>UNIT["metre",1,</p>
<p class="MsoNormal"><span style=""> </span>AUTHORITY["EPSG","9001"]]]</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Any clarification would be much appreciated.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Thanking you</p>
<p class="MsoNormal">nimalika</p>