ThinkGeo.com    |     Documentation    |     Premium Support

How to use GetFeaturesWithinDistanceOf Properly?

#1

 Hello,


I would like to understand how to use the GetFeaturesWithinDistanceOf properly. I am now either getting no features or the whole dataset. I am using the 4.0 version.


I have a map using the OpenStreetMapOverlay and I am loading some shapefiles (defined in WGS:84 / Esri 4326) on it.


I am trying to allow the user to select the shapes and I am using the GetFeaturesWithinDistanceOf method.


 



wgsToMercator = new ManagedProj4Projection(ManagedProj4Projection.GetEsriParameters(4326),
ManagedProj4Projection.GetSphericalMercatorParameter());
mercatorToWgs = new ManagedProj4Projection(ManagedProj4Projection.GetSphericalMercatorParameter(),
ManagedProj4Projection.GetEsriParameters(4326));
roadLayer.Open();
mercatorToWgs.Open();
var shape = mercatorToWgs.ConvertToExternalProjection(clickWorldPosition);
mercatorToWgs.Close();
var featuresWithinDistance = roadLayer.QueryTools.GetFeaturesWithinDistanceOf(shape,GeographyUnit.DecimalDegree, DistanceUnit.Meter, 300, ReturningColumnsType.AllColumns);

The shape point is coming out as a sensible lat/long value that is actually on spot to where I click. 


If I set the Geographic unit to meter I get all the shapes in the feature source. If I set the Geographic unit to degrees I get none. The shape file itself is in decimal degrees. What am I doing wrong  ? Am I wrong in the assumption (and since it isn't documented) that the distance parameter of the method relates to distance unit ? What unit should the distance be ?


Thanks for any pointers.


EDIT


I decided to fiddle with the distance parameter, and I am finding the results rather odd.


Selecting the same point:


With GeographyUnit.Meter, DistanceUnit.Meter 



    Distance vs features found
    
        
            1
            0
        
        
            10
            0
        
        
            50
            0
        
        
            50.5
            0
        
        
            51
            230
        
        
            52
            51207
        
        
            53
            89497
        
        
            56
            121143 (complete dataset)
        
        
            62
            121143
        
        
            75
            121143
        
    

With GeographyUnit.DecimalDegreer, DistanceUnit.Meter 



    
        
            100
            0
        
        
            1000
            0
        
        
            5 000 000
            0
        
        
            5 100 000 
            0
        
        
            5 300 000
            0
        
        
            5 500 000
            0 (long search)
        
        
            5 550 000 
            0
        
        
            5 600 000
            4134
        
        
            5 750 000
            78741
        
        
            6 000 000
            121143
        
    

I cannot make of these values. All I want is to get the feature I clicked on, I undesrtand that lines and so on make this bit more complex and that distance comes in play, but it seems the distance parameter is waay to sensitive and I do not know to what real units can I relate the distance to. The searches take very very long as well. I must be doing something wrong.


 



#2

I found my problem. 
 If you assign a projection to a layer than you must NOT convert the point you get from the world click position using the projecting. This is done implicitly under the hood. Once I sent the clickworldPosition directly to the query my results made sense.  
  
 I kind of wish this was documented somewhere that when you assign projection to a layer anything you do with that layer convert the incoming points based on that projection. 
  
 Tomas

#3

 We actually have a Code Community that addresses that issue of doing Distance Query on a projected layer, Distance Query on Projected Layers  wiki.thinkgeo.com/wiki/Map_Suite_De...ted_Layers


 But there still seems to be a bug for the case of projected layers from Geodetic to Spherical Mercator where the behavior of that distance query is unexpected. I raised this issue to the developement team and they should provide a fix soon. Thank you for pointing that out.



#4

Thanks for the kind words.

#5

Tomas, 
   
  The development team told me that this bug have been fixed for while but the fix was only available in the Development Branch. Now the latest Release Branch has been updated. So, if you get the latest Bug Daily Release (4.5.0.86 or later), you should have the fix. Thank you.

#6

Great. Thanks.

#7

Tomas, 
  
  You are welcome. If you have any other concerns, let us know.