Nicolas ,<div><br></div><div>Interesting problem, I'll take a look at it. Would you create a ticket for this issue?</div><div><br></div><div>Best regards,</div><div><br></div><div>Tamas</div><div><br><br><div class="gmail_quote">
2012/4/18 Nicolas Garel <span dir="ltr"><<a href="mailto:ngarel@micromine.com">ngarel@micromine.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-AU" link="blue" vlink="purple"><div><p class="MsoNormal">Hi,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I’m relaunching the following thread:<u></u><u></u></p><p class="MsoNormal">
<a href="http://lists.osgeo.org/pipermail/gdal-dev/2011-June/029034.html" target="_blank">http://lists.osgeo.org/pipermail/gdal-dev/2011-June/029034.html</a><u></u><u></u></p><p class="MsoNormal">This time with GDAL 1.9 and the problem remains.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Further explanation of what I’m trying to achieve:<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I would like to write and read 3d spatial data. <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I have a 3d point(<i>-1092.7172471465196</i> <i>122.9009469897652</i> <i>-150.11442565917974)</i> and 3d line (-1363.889169860787 -5.255646621772087 -150.1144256591798,-1109.9664261063763 142.60863989620026 -17.396886979865542)<u></u><u></u></p>
<p class="MsoNormal">I’ve have the passed the DIM=3D parameter when creating the layer and made sure that the parameter is being read properly into gdal.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal" style="text-autospace:none">
It seems that writing the LINESTRING into the spatial database is correct because I’ve been debugging <span style="font-size:9.5pt;font-family:Consolas">OGRMSSQLSpatialTableLayer::CreateFeature </span>and the SQL statement seems alright for the line:<u></u><u></u></p>
<p class="MsoNormal" style="text-autospace:none">INSERT INTO [dbo].[test01] (ogr_geometry) VALUES (geometry::STGeomFromText('LINESTRING <b>(-1363.889169860787 -5.255646621772087 -150.1144256591798,-1109.9664261063763 142.60863989620026 -17.396886979865542</b>)',0)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Now if I want to read the data back <b>the Z component seems to be mixed up with other components of other points</b>. Here is the output of ogrinfo when reading <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">>ogrinfo -al "MSSQL:server=mm00786\mm00786;database=Spatial;tables=dbo.test01;trusted_connection=yes"<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">INFO: Open of `MSSQL:server=mm00786\mm00786;database=Spatial;tables=dbo.test01;trusted_connection=yes'<u></u><u></u></p><p class="MsoNormal"> using driver `MSSQLSpatial' successful.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Layer name: test01<u></u><u></u></p><p class="MsoNormal">Geometry: Unknown (any)<u></u><u></u></p><p class="MsoNormal">Feature Count: 2<u></u><u></u></p><p class="MsoNormal">
Extent: (-1363.889170, -150.114426) - (142.608640, 122.900947)<u></u><u></u></p><p class="MsoNormal">Layer SRS WKT:<u></u><u></u></p><p class="MsoNormal">(unknown)<u></u><u></u></p><p class="MsoNormal">FID Column = ogr_fid<u></u><u></u></p>
<p class="MsoNormal">Geometry Column = ogr_geometry<u></u><u></u></p><p class="MsoNormal">OGRFeature(test01):1<u></u><u></u></p><p class="MsoNormal"> POINT (-1092.7172471465196 122.9009469897652 -150.11442565917974)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">OGRFeature(test01):2<u></u><u></u></p><p class="MsoNormal"> LINESTRING <b>(-1363.889169860787 -5.255646621772087 -1109.9664261063763,142.60863989620026 -150.1144256591798 -17.396886979865542)</b><u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Coordinates of the lines are WRONG. Check above with the example given above. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">However IF I specify a different geometry format other than native the X and Y coordinates are read properly but Z is omitted. See below ogrinfo ouput:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">>ogrinfo -al "MSSQL:server=mm00786\mm00786;database=Spatial;tables=dbo.test01;trusted_connection=yes;<b>GeometryFormat=wkt;"</b><u></u><u></u></p><p class="MsoNormal">
INFO: Open of `MSSQL:server=mm00786\mm00786;database=Spatial;tables=dbo.test01;trusted_connection=yes;GeometryFormat=wkt;'<u></u><u></u></p><p class="MsoNormal"> using driver `MSSQLSpatial' successful.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Layer name: test01<u></u><u></u></p><p class="MsoNormal">Geometry: Unknown (any)<u></u><u></u></p><p class="MsoNormal">Feature Count: 2<u></u><u></u></p><p class="MsoNormal">
Extent: (-1363.889170, -5.255647) - (-1092.717247, 142.608640)<u></u><u></u></p><p class="MsoNormal">Layer SRS WKT:<u></u><u></u></p><p class="MsoNormal">(unknown)<u></u><u></u></p><p class="MsoNormal">FID Column = ogr_fid<u></u><u></u></p>
<p class="MsoNormal">Geometry Column = ogr_geometry<u></u><u></u></p><p class="MsoNormal">OGRFeature(test01):1<u></u><u></u></p><p class="MsoNormal"> POINT (-1092.7172471465196 122.9009469897652)<u></u><u></u></p><p class="MsoNormal">
<u></u> <u></u></p><p class="MsoNormal">OGRFeature(test01):2<u></u><u></u></p><p class="MsoNormal"> LINESTRING <b>(-1363.889169860787 -5.255646621772087,-1109.9664261063763 142.60863989620026</b>)<u></u><u></u></p><p class="MsoNormal">
<u></u> <u></u></p><p class="MsoNormal">Can you <b>please</b> tell me if the problem is from the writing or reading of GDAL or maybe I’m simply missing some information during writing or reading?!<u></u><u></u></p><p class="MsoNormal">
<u></u> <u></u></p><p class="MsoNormal">Please reply to me. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d">Nicolas Garel </span></b><b><span lang="EN-GB" style="font-size:10.0pt;color:#0070c0">|</span></b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d"> Junior Software Engineer<u></u><u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d">M</span></b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d">icromine Product<br><b>Email </b></span><b><span lang="EN-GB" style="font-size:10.0pt;color:#0070c0">|</span></b><b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d"> </span></b><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d"><a href="mailto:ngarel@micromine.com" target="_blank"><span style="color:blue">ngarel@micromine.com</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-GB" style="font-size:10.0pt;color:#1f497d"><a href="http://www.micromine.com/" target="_blank"><span style="color:blue">http://www.micromine.com</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p></div><br><br>
<div align="left">
<p align="left"><img style="MARGIN:0px" border="0" alt="" src="cid:mm_corp_emailsig_perth.gif@5b728f68f250493398f956fb4b364ca1"> <br> </p></div></div><br>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br></div>