ThinkGeo.com    |     Documentation    |     Premium Support

Some objects from SQL data not showing in map

#3

Any idea how I can pass the data to you? I have the date in 2 small tables and definitely can confirm the results every time. Your shapes don’t appear to be making use of the rotating symbol class which I included - not sure if that is the cause… The symbols should rotate on the VCMG field. 
  
 Quite possibly it does work with your data - it works with 1 table of mine and doesn’t work with another, even though the data originally came from the same table and is the same format - only the lat & long are slightly different.

#4

Hi David, 
  
 Please send the data to forumsupport@thinkgeo.com or ask sales@thinkgeo.com for a FTP address to upload large size data. 
  
 Thanks, 
 Peter

#5

The problem isn’t the size of the data - in fact the data is sitting in a zip file 2 messages up in a standard SQL Server export to CSV but you say you can’t use that. We are only talking about 6 records in total. What format do you want the data in? I can post it here but if you can’t take CSV then what format will you accept?

#6

Hi David, 
  
 The SQL Scripts will be better. 
  
 Thanks, 
 Peter

#7

Ok, figured out how to do that - SQL files attached.



Thanks,



Dave

dbo.tmpTest2Vessels.Table.zip (3.07 KB)

#8

Hi David, 
  
 Thanks for your share! Any question please let us know. 
  
 Regards, 
  
 Don

#9

Actually my question still stands - why do I only see 2 of the 3 features with 1 table? I am hoping you can duplicate the problem now you have the sample code & data.

#10

Hi David, 
  
 Thanks for your data and test project, I found if I get all the features then add them to other layer, I can shows all the 3 points on map. But if I render them via the SQLserver layer I missed the id 533057600 point, the problem is our MSSQL2008FeatureLayer render features by this SQL statement: 
  
 SELECT [WGS84Loc].STAsBinary() as [WGS84Loc],[ID] FROM tmpTest2Vessels WHERE (geography::STGeomFromText(‘POLYGON((100.681610107422 7.53803430125117,100.681610107422 7.3200822994113,102.952267706394 7.3200822994113,102.952267706394 7.53803430125117,100.681610107422 7.53803430125117))’,4326).STIntersects(WGS84Loc)=1); 
  
 When I tried it in MSSQL, it only return two result.  
  
 I think maybe the problem is the WGS84Loc for 533057600 have some problem, but it’s not saved as a readable format. I think maybe we can do further research about that after our MSSQL developer back office, but if you get any information from this SQL statement please let me know. 
  
 Regards, 
  
 Don

#11

Ok, thanks - I hope you can find out why it isn’t showing as SQL Server Management Studio can show the spatial object with no problems. I am not doing anything different with that point - it is just being read in from a GPS device and INSERTed into the SQL Server table with a standard POINT() command, along with a few million other points. I have found a few missing points so far but only in ThinkGeo, not in SSMS.

#12

Hi David, 
  
 Could you show me the Sal scripts that you are using in SSMS? I was just wondering what’s the difference between the one we are using. 
  
 Thanks, 
 Johnny

#13

Just 'select * from tablename" and click on the “spatial” tab.

#14

Just in case that wasn’t clear, here is the view from SQL Server Management Studio - the green circles are added by me as the dots are a bit difficult to see. But you can see both tables are showing 3 points spatially. If SSMS can show them, it is strange that MapSuite isn’t showing them.

SpatialResults.png (63 KB)

#15

Hi David,  
  
 Just like my previous reply, our result related with the SQL query statement, so I think maybe the problem is related with the  query statement, our developer will try to find the reason about this problem, any update we will let you know. 
  
 Regards, 
  
 Don

#16

Hi David,



Sorry for delay on this. The reason why it can’t display in MapSuite Map is that the missing point doesn’t intersect with the BoundingBox. The Map only draws the features intersect with the BoundingBox. In the table tmpTest2Vessels the record 533057600 doesn’t intersect with the BoundingBox so it draws 2 features. The following is the screenshot from the SSMS.







For the GetCount() method we don’t filter the data by the BoundingBox intersection and it returns the count of records in the data table. So it returns 3 as the result.



Hope it’s helpful.



Regards,

Peter

#17

I will take a better look when I get to work, but I am confused as to why it isn’t in the bounding box when the code should be zooming the map out to the extents of the bounding box and then I even zoom out another 50% just to be on the safe side. 
  
 Dim mapOutline As New ThinkGeo.MapSuite.Core.RectangleShape 
 mapOutline = VesselLayerBad.GetBoundingBox 
 MapTwo.CurrentExtent = mapOutline 
 MapTwo.ZoomOut(50) 
  
 Can you tell me what query you are using to generate the Intersects value? I guess you typed in a bounding box polygon or something like that? 
  
 Thanks, 
  
 Dave

#18

Hi David,



The first one is tracing by the SQL Server Profiler and the other one is modified based on the first one which will get the result as the screenshot I attached before.



SQL Server Profiler:

exec sp_executesql N’SELECT [VCMG],[WGS84Loc].STAsBinary() as [WGS84Loc],[ID] FROM tmpTest2Vessels WHERE (geography::STGeomFromText(@Geography,4326).STIntersects(WGS84Loc)=1);’,N’@Geography varchar(178)’,@Geography=‘POLYGON((100.681610107422 7.53803430125117,100.681610107422 7.3200822994113,102.952267706394 7.3200822994113,102.952267706394 7.53803430125117,100.681610107422 7.53803430125117))‘



SSMS:

exec sp_executesql N’SELECT [ID],[VLat],[VLon], (geography::STGeomFromText(@Geography,4326).STIntersects(WGS84Loc)) AS Intersects FROM tmpTest2Vessels;’,N’@Geography varchar(178)’,@Geography='POLYGON((100.681610107422 7.53803430125117,100.681610107422 7.3200822994113,102.952267706394 7.3200822994113,102.952267706394 7.53803430125117,100.681610107422 7.53803430125117))'



Thanks,

Peter

#19

I am in the office now and had a better look at your SSMS output. The ‘missing point’ has a Longitude of 100.7623611 which is in between the other 2 points which are 100.68 and 102.95. The Latitude of the missing point is admittedly a tiny bit south of the other 2 points, but not enough to be outside the map area. So I don’t see how it can be outside the bounding box…



I just added a small bit of code:



pointLayer.FeatureSource.Projection = proj4
Dim missingfeature As New Feature(100.7623611, 7.3200833, “Point1”)
pointLayer.InternalFeatures.Add(“Point1”, missingfeature)
pointLayer.ZoomLevelSet.ZoomLevel01.DefaultPointStyle = PointStyles.CreateSimpleCircleStyle(GeoColor.FromHtml("#ED775B"), 10, GeoColor.FromHtml("#CA5B48"))
pointLayer.ZoomLevelSet.ZoomLevel01.ApplyUntilZoomLevel = ApplyUntilZoomLevel.Level20



to the existing sample I gave you and the results show that the ‘missing point’ (or at least another point in the same place on the same map on a different layer) is in fact in the middle of the 2 points that are shown. So again, I don’t understand how it can be outside the bounding box.



Also, even if the point was outside the bounding box, when one zooms out (which is possible in the sample project) the point should be visible?



Actually after I put that orange point in I had a frame of reference and found that if I double click on the orange point 11 times (which zooms in each time) the ‘missing point’ actually does finally show up - right underneath the orange point.

So, why doesn’t it show normally? I typically have thousands of these green triangles on the map at once, overlapping very densely so I don’t think it is being trimmed because if an overlap…









#20

Hi Dave,



I guess this issue is caused by the SQL Server spatial coordinate calculation precision.  The geography is used to represent geodetic data and the geometry is used to represent planar data.  The MsSqlFeatureLayer gets the intersection features by SQL Server spatial calculation which will project planer data to spherical by the geography::STGeomFromText(@Geography,4326). As you mentioned before the ‘missing point’ is in the boundingbox, this is true in the planar coordinate but may be not true in the spherical coordinate.



The InMemoryFeatureLayer gets the Intersection features just via mathematical calculation so the ‘missing point’ will be plot on the map.



There is a workaround, we can do the buffer operation before intersects by binding the  ExecutingSqlStatement event and rebuild the SQL statement in the event, the code snippet maybe like the following:



Bind event:


AddHandler (DirectCast(VesselLayerBad.FeatureSource, MsSql2008FeatureSource)).ExecutingSqlStatement, New EventHandler(Of ExecutingSqlStatementMsSql2008FeatureSourceEventArgs)(AddressOf OnExecutingSqlStatement)



Handler:




Private Sub OnExecutingSqlStatement(ByVal sender As Object, ByVal e As ExecutingSqlStatementMsSql2008FeatureSourceEventArgs)
       If (e.ExecutingSqlStatementType = ExecutingSqlStatementType.GetFeaturesInsideBoundingBoxEx) Then
           e.SqlStatement = e.SqlStatement.Replace(".STIntersects", “.STBuffer(22).STIntersects”)
       End If
   End Sub

NOTE: In STBuffer(buffer) method above, you can set the buffer to a proper value.



Regards,

Peter




#21

Ok, that does work now. Thanks very much.



 A little concerned about the performance impact after reading:

stackoverflow.com/questions/6776721/sql-server-2008-geometry-stbuffer-really-slow

but I think most users are on 2010 or 2012 by now so hopefully that won’t be an issue.

From what I can see, by “set the buffer to a proper value” you are talking about the ‘22’ and that is metres. I played with values smaller and the point disappears so it seems you found the minimum. I guess the only effect of choosing bigger numbers will be that it would slow down the task. I will put the code into the production code and see if it still works with a larger number as there are a lot of variables such as number of points, area they cover and even map height & width so I guess there is no ‘one-size-fits-all’ number.



So, thanks for the help - that is one of my problems out of the way. Next problem is to find out why I am getting memory overflow errors with raster overlays! :-)

#22

Hi Dave, 
  
 Very glad to hear it works with you. 
  
 I did some tests based on the data you provided before to got the minimum value ‘22’. As I didn’t know what the value you need here so I mentioned “set the buffer to a proper value”. 
  
 ‘About the memory overflow error with raster overlays’ please create a new post in a new thread. 
  
  
 Hope to meet you in the other thread. 
  
 Thanks, 
 Peter