search

MapServer / MapServer-import

3 stars 2 forks source link

[PATCH] Fix reprojection bug in mapwcs.c #2180

Closed tbonfort closed 12 years ago

tbonfort commented 12 years ago

Reporter: rouault Date: 2007/07/26 - 19:26 Producing non-square pixels or reprojected results with MapServer WCS does not work with tileindexed layers. The issue is described in the email:

http://lists.umn.edu/cgi-bin/wa?A2=ind0702&L=mapserver-users&T=0&F=&S=&P=852)

With the patch, I can make GetCoverage requests on a tileindex of DTED files with resolutions (or dimensions) not necesseraly a multiple of the original resolution, or into another projection system.

The fix consists of adding the following lines after msMapComputeGeotransform, as it's done in mapdraw.c


    /* Do we need to fake out stuff for rotated support? */
    if( map->gt.need_geotransform )
        msMapSetFakedExtent( map );
tbonfort commented 12 years ago

Author: warmerdam Date: 2007/07/26 - 21:27 I've made a first attempt to reproduce this problem by requesting an odd aspect ratio in a WCS request but things seemed to work fine. I'll try again tomorrow with reprojection.

tbonfort commented 12 years ago

Author: rouault Date: 2007/07/26 - 22:06 The problem with requesting arbitrary resolutions only appears when the LAYER is a TILEINDEX, even if the TILEINDEX only contains one DTED file. If your layer is made of just one single DTED file, you can query arbitrary resolutions.

tbonfort commented 12 years ago

Author: warmerdam Date: 2007/07/27 - 19:37 The patch has been applied in the 4.10 branch (for 4.10.3) and in trunk (for 5.0 beta2).

I have also incorporated a two extra tests in msautotest/wcs/wcs_simple.map for this issue, and changed that map to use a one-file tileindex so the problem is triggered.