ThinkGeo.com    |     Documentation    |     Premium Support

Projection problem?

#1

Hi,


I have a layer containing points with the internal projection as Extended Lambert II (ESRI 102582). The external projection is Spherical Mercator. I use also an OpenStreetMap layer.


As you can see on the picture, there's a shift.


The red triangle is the point shown by thinkgeo. The blue circle is the point in WGS84 (conversion by an external tool). The green star is the point in WGS84 obtained by geocoding using google api. The point represent a drugstore "Pharmacie de la croix bleue".



I've attached a project that reproduce the shift.


Thanks,


P.OM



#2

 Puthirak,


 
This issue is caused by projection mismatch. You can adjust your ExternalPorjectionParametersString to meet your requirement. The code snippet below is for your reference.
 
    ManagedProj4Projection proj = new ManagedProj4Projection();
// NTF Lambert II Etendu
// 
proj.InternalProjectionParametersString = ManagedProj4Projection.GetEsriParametersString(102582);
            //proj.ExternalProjectionParametersString = ManagedProj4Projection.GetGoogleMapParametersString();
            proj.ExternalProjectionParametersString = "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"; 
 
For detail please see another thread below:
gis.thinkgeo.com/Support/DiscussionForums/tabid/143/aff/39/aft/9175/afv/topic/Default.aspx
 
Thanks,
 
James
 

#3

 Puthirak,


 
This issue is caused by projection mismatch. You can adjust your ExternalPorjectionParametersString to meet your requirement. The code snippet below is for your reference.
 
    ManagedProj4Projection proj = new ManagedProj4Projection();
// NTF Lambert II Etendu
proj.InternalProjectionParametersString = ManagedProj4Projection.GetEsriParametersString(102582);
            //proj.ExternalProjectionParametersString = ManagedProj4Projection.GetGoogleMapParametersString();
            proj.ExternalProjectionParametersString = "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"; 
 
 
For detail please see another thread below:
gis.thinkgeo.com/Support/DiscussionForums/tabid/143/aff/39/aft/9175/afv/topic/Default.aspx
 
Thanks,
 
James

#4

Hi James,


I was aware of this projection mismatch and I already tried it. But I get the same result...


Did you try it with my attached sample ?


Thanks,


P.OM



#5

 Puthirak,



 


Sorry about that, my machine couldn’t display the background OpenStreetMapOverlay, I think maybe my network service provider block the remote server of OSM or network speed is slow to connect the server, so that I couldn’t compare the distance of point. 


I don’t have good solution so far, you know, our projection system rely on a third part open source library Proj4, so probably we couldn’t fix it soon, you can find more information for it at trac.osgeo.org/proj/


 


Thanks,


James



#6

Hi James,


I find that the Extended Lambert II could also be represented by EPSG:27582. This, combined with the projection mismatch correction, give me the correct location :)


 


Thanks,


P.OM



#7

Puthirak, 
  
 I am glad you got the correct location, sorry about that I didn’t give you much help, but if you have any more problem related with our map suite, just feel free to let me know. 
  
 Thanks, 
  
 James