Thanks David, Using the correct units seems to have fixed the problem.
@Jeff : Thanks for pointer, now that I have re-read it all makes sense Cheers David On 29 Aug 2013, at 16:43, Fawcett, David (MNIT) wrote: > David, > > Take a look at the map extent and units. They should be the units used by > your output spatial reference system. > > In your current map file, you have DD for the units, for 3857, they should be > meters. You will also want to express your extent in 3857 coords. > > David. > > From: [email protected] > [mailto:[email protected]] On Behalf Of David Power > Sent: Thursday, August 29, 2013 10:01 AM > To: MapserverList OSGEO > Subject: [mapserver-users] Projections and raster files > > Hello, > > I have a number of geotiffs which are currently using epsg:4326 and I would > like to output a region using epsg:3857. > > I am new to mapserver but from what I have read the way achieve this is by > setting the projection of the map to epsg:3857 and the layer to epsg:4326. > > When I set the projections to the same value I get the output I expect > without a problem, if I set the projections to different values I do not get > an error but the output is blank. > > To make life more complicated I am using mapscript via the Java swig > interface, the version of mapserver I am using is mapserver 6.1.2.1 which I > built from source. > > If anybody has any advice I would be most grateful. > > David > > - > > gdalinfo gives the following information about one of the files: > > Driver: GTiff/GeoTIFF > Files: globcover_n10e005.tif > Size is 7200, 7200 > Coordinate System is: > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, > AUTHORITY["EPSG","7030"]], > AUTHORITY["EPSG","6326"]], > PRIMEM["Greenwich",0], > UNIT["degree",0.0174532925199433], > AUTHORITY["EPSG","4326"]] > Origin = (5.000000000000000,9.999999999999615) > Pixel Size = (0.000138888888889,-0.000138888888889) > Metadata: > DataType=Generic > AREA_OR_POINT=Area > Image Structure Metadata: > INTERLEAVE=BAND > Corner Coordinates: > Upper Left ( 5.0000000, 10.0000000) ( 5d 0'0.00"E, 10d 0'0.00"N) > Lower Left ( 5.0000000, 9.0000000) ( 5d 0'0.00"E, 9d 0'0.00"N) > Upper Right ( 6.0000000, 10.0000000) ( 6d 0'0.00"E, 10d 0'0.00"N) > Lower Right ( 6.0000000, 9.0000000) ( 6d 0'0.00"E, 9d 0'0.00"N) > Center ( 5.5000000, 9.5000000) ( 5d30'0.00"E, 9d30'0.00"N) > Band 1 Block=7200x128 Type=Byte, ColorInterp=Gray > NoData Value=0 > > - > > The map file is built programmatically, the following is the result of saving > the map object to file: > > MAP > EXTENT 5 9.9 5.1 10 > IMAGETYPE "png" > NAME "blank" > SIZE 800 800 > STATUS ON > UNITS DD > > OUTPUTFORMAT > NAME "png" > MIMETYPE "image/png" > DRIVER "AGG/PNG" > EXTENSION "png" > IMAGEMODE RGB > TRANSPARENT FALSE > END # OUTPUTFORMAT > > PROJECTION > "init=epsg:3857" > END # PROJECTION > LEGEND > KEYSIZE 20 10 > KEYSPACING 5 5 > LABEL > SIZE MEDIUM > OFFSET 0 0 > SHADOWSIZE 1 1 > TYPE BITMAP > END # LABEL > STATUS OFF > END # LEGEND > > QUERYMAP > SIZE -1 -1 > STATUS OFF > STYLE HILITE > END # QUERYMAP > > SCALEBAR > INTERVALS 4 > LABEL > SIZE MEDIUM > OFFSET 0 0 > SHADOWSIZE 1 1 > TYPE BITMAP > END # LABEL > SIZE 200 3 > STATUS OFF > UNITS MILES > END # SCALEBAR > > WEB > END # WEB > > LAYER > DATA "/home/leftworker/data/landcover/lefttiles/n10/globcover_n10e005.tif" > NAME "n10e005" > PROJECTION > "init=epsg:4326" > END # PROJECTION > STATUS ON > TILEITEM "location" > TYPE RASTER > UNITS DD > END # LAYER > > END # MAP > > > > > > ====================================== > Dr David J Power > Research Officer > Department of Computer Science > University of Oxford > tel: 01865 283670 > email: [email protected] > ====================================== > > > ====================================== Dr David J Power Research Officer Department of Computer Science University of Oxford tel: 01865 283670 email: [email protected] ======================================
_______________________________________________ mapserver-users mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/mapserver-users
