[mapserver-users] Projections and raster files

[David Power]

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-bounces at lists.osgeo.org [mailto:mapserver-users-bounces at 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.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: David.Power at cs.ox.ac.uk
> ======================================
>  
>  
>  

======================================
Dr David J Power
Research Officer
Department of Computer Science
University of Oxford
tel:   01865 283670
email: David.Power at cs.ox.ac.uk
======================================



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20130829/b03f486b/attachment.htm>