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: mapserver-users-boun...@lists.osgeo.org
[mailto:mapserver-users-boun...@lists.osgeo.org] 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.000,9.615)
Pixel Size = (0.0001389,-0.0001389)
Metadata:
DataType=Generic
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 5.000, 10.000) ( 5d 0'0.00E, 10d 0'0.00N)
Lower Left ( 5.000, 9.000) ( 5d 0'0.00E, 9d 0'0.00N)
Upper Right ( 6.000, 10.000) ( 6d 0'0.00E, 10d 0'0.00N)
Lower Right ( 6.000, 9.000) ( 6d 0'0.00E, 9d 0'0.00N)
Center ( 5.500, 9.500) ( 5d30'0.00E, 9d30'0.00N)
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: david.po...@cs.ox.ac.uk
==
==
Dr David J Power
Research Officer
Department of Computer Science
University of Oxford
tel: 01865 283670
email: david.po...@cs.ox.ac.uk
==
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users