<div>Hi</div>
<div> </div>
<div>bingo! I tested wms using qgis 2.0.1 and there is a bug regarding revert axis order.</div>
<div>Than you Jukka very much :)</div>
<div> </div>
<div>Bartek</div>
<div> </div>
<div>W dniu 2014-09-02 16:43:35 użytkownik Rahkonen Jukka (Tike) <jukka.rahkonen@mmmtike.fi> napisał:</div>
<blockquote style="margin-left: 7px; border-left: 2px solid orange; padding-left: 8px;">
<div>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US">Hi, </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US"> </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US">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" target="_blank"> 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> </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US"> </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US">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. </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US"> </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US">-Jukka Rahkonen- </span></p>
<p><span style="font-size: 11.0pt; font-family: 'Calibri','sans-serif'; color: #1f497d;" lang="EN-US"> </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><span style="font-size: 10.0pt; font-family: 'Tahoma','sans-serif'; color: windowtext;" lang="EN-US">Bartlomiej Burkot wrote: </span></p>
</div>
</div>
<p><span lang="EN-US"> </span></p>
<p><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" target="_blank">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</p>
</div>
</div>
</blockquote>
<div> </div>