Opened 18 years ago
Closed 18 years ago
#1509 closed defect (fixed)
Incomplete extent of raster after reprojection - corner missing
Reported by: | Owned by: | warmerdam | |
---|---|---|---|
Priority: | high | Milestone: | |
Component: | GDAL Support | Version: | 4.6 |
Severity: | normal | Keywords: | |
Cc: | jacob.delfos@… |
Description
When reprojecting a raster on-the-fly, a corner of the image disappears. It seems that the image gets cropped too much prior to reprojecting, such that the distortion by the reprojection causes it not to cover the entire map image anymore. The original raster is in an LCC projection. The problem occurs when projecting to UTM or Geographic LL. Image format is ecw, but it also occurs as TIFF. The problem exists in mapserver 4.4 and 4.6, using gdal 1.2.6 and 1.3.1. The distortion after reprojection is not unreasonably large, and should be acceptable, because the extent falls within the zone it is projected to, which falls within the geographic extent of the projected image (Australia). I'm attaching a mapfile, with a downsampled version of the image. The problem is reproducable with these files, using the shp2img command in the file "SHP2IMG_CALL.TXT".
Attachments (1)
Change History (4)
by , 18 years ago
comment:1 by , 18 years ago
Cc: | added |
---|
comment:3 by , 18 years ago
Resolution: | → fixed |
---|---|
Status: | assigned → closed |
I have determined that the msTransformMapSource() function in mapresample.c was not including the bottom right corner in it's calculations and was therefor not loading the right extent when there was a non-trivial rotation as part of the rotation. I have corrected this in CVS. 4.8b1 should include the fix when released later this week. Thanks!
Note:
See TracTickets
for help on using tickets.
Test files and shp2img command