<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";
color:black;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.Shkpostityyli17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 2.0cm 70.85pt 2.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body bgcolor="white" lang="FI" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Your system seems to have Northing-Easting axis order
<a href="http://epsg-registry.org/report.htm?type=selection&entity=urn:ogc:def:crs:EPSG::31255&reportDetail=short&style=urn:uuid:report-style:default-with-code&style_name=OGP%20Default%20With%20Code&title=EPSG:31255">
http://epsg-registry.org/report.htm?type=selection&entity=urn:ogc:def:crs:EPSG::31255&reportDetail=short&style=urn:uuid:report-style:default-with-code&style_name=OGP%20Default%20With%20Code&title=EPSG:31255</a><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Check if your BBOX has north coordinates first. If not, flip the coordinates or alternatively change your WMS version into 1.1.1 and CRS into SRS
and keep your BBOX.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">-Jukka Rahkonen-<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext">Bartlomiej Burkot wrote:<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Hi<br>
<br>
I spent much time searching for reason why my raster does not appear in WMS layer.<br>
My raster is in EPSG:31255 projection<br>
<br>
</span>gdalinfo samle_raster.tif<br>
<span style="font-size:7.5pt">Driver: GTiff/GeoTIFF<br>
Files: samle_raster.tif<br>
Size is 5000, 4000<br>
Coordinate System is:<br>
PROJCS["MGI / Austria GK Central",<br>
GEOGCS["MGI",<br>
DATUM["Militar_Geographische_Institute",<br>
SPHEROID["Bessel 1841",6377397.155,299.1528128000008,<br>
AUTHORITY["EPSG","7004"]],<br>
TOWGS84[577.326,90.129,463.919,5.137,1.474,5.297,2.4232],<br>
AUTHORITY["EPSG","6312"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4312"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],<br>
PARAMETER["central_meridian",13.33333333333333],<br>
PARAMETER["scale_factor",1],<br>
PARAMETER["false_easting",0],<br>
PARAMETER["false_northing",-5000000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","31255"]]<br>
Origin = (-112500.000000057683792,274999.999999893829226)<br>
Pixel Size = (0.249999999999589,-0.249999999999589)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)<br>
TIFFTAG_XRESOLUTION=72<br>
TIFFTAG_YRESOLUTION=72<br>
Image Structure Metadata:<br>
INTERLEAVE=PIXEL<br>
Corner Coordinates:<br>
Upper Left ( -112500.000, 275000.000) ( 11d50'13.22"E, 47d36'16.17"N)<br>
Lower Left ( -112500.000, 274000.000) ( 11d50'14.14"E, 47d35'43.80"N)<br>
Upper Right ( -111250.000, 275000.000) ( 11d51'13.05"E, 47d36'16.95"N)<br>
Lower Right ( -111250.000, 274000.000) ( 11d51'13.96"E, 47d35'44.57"N)<br>
Center ( -111875.000, 274500.000) ( 11d50'43.59"E, 47d36' 0.37"N)<br>
Band 1 Block=5000x1 Type=Byte, ColorInterp=Red<br>
Description = Band_1<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Band 2 Block=5000x1 Type=Byte, ColorInterp=Green<br>
Description = Band_2<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Band 3 Block=5000x1 Type=Byte, ColorInterp=Blue<br>
Description = Band_3<br>
Metadata:<br>
LAYER_TYPE=athematic</span><br>
<br>
I can show it in qgis and reproject to any coordinate system for example epsg:3857.<br>
<br>
My mapfile is:<br>
<span style="font-size:7.5pt">MAP<br>
IMAGETYPE PNG<br>
EXTENT -120000 270000 -110000 275000<br>
SIZE 800 600<br>
RESOLUTION 72<br>
CONFIG "MS_ERRORFILE" "../log/ms_error_standard.log"<br>
DEBUG 5<br>
NAME myname<br>
CONFIG "CPL_DEBUG" "ON" <br>
CONFIG "PROJ_DEBUG" "ON"<br>
CONFIG "PROJ_LIB" "/usr/share/proj/"<br>
<br>
PROJECTION<br>
"init=epsg:31255"<br>
END<br>
WEB<br>
METADATA<br>
wms_title "standard data in wgw"<br>
"wms_enable_request" "GetMap GetFeatureInfo GetCapabilities"<br>
"wms_srs" "EPSG:31255 EPSG:4326 EPSG:3857 EPSG:2180 EPSG:31287"<br>
END<br>
END<br>
<br>
LAYER<br>
NAME "oneRaster"<br>
DATA "samle_raster.tif"<br>
TYPE RASTER<br>
STATUS ON<br>
PROJECTION<br>
"init=epsg:31255"<br>
END<br>
METADATA<br>
wms_title "Sample raster"<br>
wms_srs "EPSG:31255"<br>
END<br>
<br>
END<br>
END #MAP</span><br>
<br>
I'm using mapserver 6.4.1 in cgi mode on the Ubuntu 14.04 linux installed with apt-get install cgi-mapserver:<br>
<span style="font-size:7.5pt">/usr/lib/cgi-bin/mapserv -v<br>
MapServer version 6.4.1 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG OUTPUT=KML SUPPORTS=PROJ SUPPORTS=GD SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO SUPPORTS=SVG_SYMBOLS SUPPORTS=RSVG SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER
SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE</span><br>
<br>
When Im requesting a map:<br>
<span style="font-size:7.5pt"><a href="http://myhost/cgi-bin/wmsstandard?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-112516.67630140134133399,273973.2257158225402236,-111233.66681915351364296,275020.36552757304161787&CRS=EPSG:31255&WIDTH=1038&HEIGHT=846&LAYERS=oneRaster&STYLES=&FORMAT=image/jpeg&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96">http://myhost/cgi-bin/wmsstandard?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-112516.67630140134133399,273973.2257158225402236,-111233.66681915351364296,275020.36552757304161787&CRS=EPSG:31255&WIDTH=1038&HEIGHT=846&LAYERS=oneRaster&STYLES=&FORMAT=image/jpeg&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96</a></span><br>
<br>
I got a message in /var/log/apache2/error.log<br>
<span style="font-size:7.5pt">[Tue Sep 02 15:01:30.386845 2014] [cgi:error] [pid 2057] [client myip:65403] AH01215: GDAL: GDALOpen(/pathto/sample.tif, this=0x819310) succeeds as GTiff.<br>
[Tue Sep 02 15:01:30.396083 2014] [cgi:error] [pid 2057] [client myip:65403] AH01215: GDAL: GDALClose(/pathto/sample.tif,</span>,<span style="font-size:7.5pt"> this=0x819310)<br>
<br>
<br>
</span>ms_error_standard.log:<br>
<span style="font-size:7.5pt">[Tue Sep 2 15:04:45 2014].609914 CGI Request 1 on process 5909<br>
[Tue Sep 2 15:04:45 2014].610243 msWMSLoadGetMapParams(): enabling non-square pixels.<br>
[Tue Sep 2 15:04:45 2014].614988 msDrawMap(): kicking into non-square pixel preserving mode.<br>
[Tue Sep 2 15:04:45 2014].615034 msDrawMap(): rendering using outputformat named AGG_JPEG (AGG/JPEG).<br>
[Tue Sep 2 15:04:45 2014].615064 msDrawMap(): WMS/WFS set-up and query, 0.000s<br>
[Tue Sep 2 15:04:45 2014].617587 msDrawRasterLayerLow(oneRaster): entering.<br>
[Tue Sep 2 15:04:45 2014].619764 msDrawMap(): Layer 1 (oneRaster), 0.005s<br>
[Tue Sep 2 15:04:45 2014].619805 msDrawMap(): Drawing Label Cache, 0.000s<br>
[Tue Sep 2 15:04:45 2014].619816 msDrawMap() total time: 0.009s<br>
[Tue Sep 2 15:04:45 2014].627898 msSaveImage(stdout) total time: 0.008s<br>
[Tue Sep 2 15:04:45 2014].628129 mapserv request processing time (msLoadMap not incl.): 0.018s<br>
[Tue Sep 2 15:04:45 2014].628152 msFreeMap(): freeing map at 0xd0c570.</span><br>
<br>
The requested image is white/empty<br>
<br>
The proj4 epsg file does include my projection:<br>
<span style="font-size:7.5pt">more /usr/share/proj/epsg | grep 31255<br>
#<31255> +proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +datum=hermannskogel +units=m +no_defs <><br>
<31255> +proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs <><br>
</span><br>
And pcs.csv as well:<br>
<span style="font-size:7.5pt"> more /usr/share/gdal/1.10/pcs.csv | grep 31255<br>
31255,"MGI / Austria GK Central",9001,4312,18005,9807,1,0,4530,8801,0,9110,8802,13.2,9110,8805,1,9201,8806,0,9001,8807,-5000000,9001,,,,,,</span><br>
<br>
When I convert the rasterfile to epsg:3857 using gdalwarp and request Layer in the same crs<span style="font-size:7.5pt"> &CRS=EPSG:31255&</span> I get the correct getMap request. If I switch to
<span style="font-size:7.5pt">&CRS=EPSG:31255&</span> to force mapserver to reproject I got the same white/empty image.<br>
<br>
Could somebody help me?<br>
I can send you a sample of raster on email address.<br>
<br>
Bartek<o:p></o:p></p>
</div>
</div>
</body>
</html>