issues with small extents

Yewondwossen Assefa assefa at DMSOLUTIONS.CA
Tue Dec 20 18:09:10 EST 2005


Hi There,

I am having some trouble with a wms request using a small bbx 
(&BbOx=-0.0042,-0.0024,0.0042,0.0024), lat/long values.  What happens 
from what I can see is that there seems to be invalid values  set inside 
the projection objects and the points transformation (from layer 
projection to map projection) seem to give uncorrect values.  BTW both 
layer and map are set to epsg:4326 but at one point the 
msProjectionsDiffer call returns true ( Igues this is another issue).  
Following in the code, It seems that my points get corrupted in 
msProjectPoint.


if( out && out->gt.need_geotransform )
  {
      double x_out, y_out;

      x_out = out->gt.invgeotransform[0]
          + out->gt.invgeotransform[1] * point->x
          + out->gt.invgeotransform[2] * point->y;
      y_out = out->gt.invgeotransform[3]
          + out->gt.invgeotransform[4] * point->x
          + out->gt.invgeotransform[5] * point->y;

      point->x = x_out;
      point->y = y_out;
}

 The values inside the invgeotransform array seems to be  wrong but I am 
not sure

invgeotransform    0x00ec6850    double [6]
    [0]    84.000000000000000    double
    [1]    20119.760479041917    double
    [2]    0.00000000000000000    double
    [3]    47.999999999999993    double
    [4]    0.00000000000000000    double
    [5]    20210.526315789473    double

-    geotransform    0x00ec6820    double [6]
    [0]    -0.0041749999999999999    double
    [1]    4.9702380952380948e-005    double
    [2]    0.00000000000000000    double
    [3]    -0.0023749999999999999    double
    [4]    0.00000000000000000    double
    [5]    4.9479166666666665e-005    double


  before conversion
    point->x    0.00000000000000000    double
    point->y    -1.0000000000000000    double

 after conversion
    point->x    84.000000000000000    double
    point->y    -20162.526315789473    double


 At what point the invgeotransform   are set ?  I realise this e-mail is 
far from being clear but I am more looking general ideas on how to 
quickly  debug this.


 Thanks

-- 
----------------------------------------------------------------
Assefa Yewondwossen           
Software Analyst   

Email: assefa at dmsolutions.ca    
http://www.dmsolutions.ca/

Phone: (613) 565-5056 (ext 14)
Fax:   (613) 565-0925
----------------------------------------------------------------



More information about the mapserver-dev mailing list