Thank you Andrea and Jukka!
Changing ExpandToRGB to true, did not help here. All the images should have
same palette, but it is difficult to be 100% sure.
Making the imagemosaic again (for may-be 5th time) helped, but it was kind of
playing bingo. It likely is the same with single files, sometimes the same file
will show up, sometimes not.
But changing the .prj file seems to be the trick :) We have tried now with 3
different datasets.
I am wondering why this strange .prj file is generated by GeoServer. Is it
possible, that GeoServer ImageMosaic has problems with EPSG:2393?
In GeoServer UI, the layer SRS is:
PROJCS["KKJ / Finland Uniform Coordinate System",
GEOGCS["KKJ",
DATUM["Kartastokoordinaattijarjestelma (1966)",
SPHEROID["International 1924", 6378388.0, 297.0, AUTHORITY["EPSG","7022"]],
TOWGS84[-96.062, -82.428, -121.753, 4.801, -0.345, 1.376, 1.496],
AUTHORITY["EPSG","6123"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH],
AUTHORITY["EPSG","4123"]],
PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]],
PARAMETER["central_meridian", 27.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 3500000.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
AUTHORITY["EPSG","2393"]]
And for geotiffs there are external .prj files for each .tiff file (copied from
ArcGIS, so may-be something else would be better to use?)
PROJCS["KKJ / Finland Uniform Coordinate
System",GEOGCS["KKJ",DATUM["Kartastokoordinaattijarjestelma_1966",SPHEROID["International
1924",6378388,297,AUTHORITY["EPSG","7022"]],TOWGS84[-96.062,-82.428,-121.753,4.801,0.345,-1.376,1.496],AUTHORITY["EPSG","6123"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4123"]],
UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",27],PARAMETER["scale_factor",1],
PARAMETER["false_easting",3500000],PARAMETER["false_northing",0],
AUTHORITY["EPSG","2393"],AXIS["X",NORTH],AXIS["Y",EAST]]
None of them has " Geocentric Y ", so I am wondering from where it comes from?
We have same GeoServer version as Jukka, but we are running it on Linux. May
that give different .prj files?
Kylli
PS. We have the same files in use with MapServer with no problems..
From: "Rahkonen Jukka (MML)" <[email protected]>
To: "Andrea Aime" <[email protected]>, "Kylli Ek" <[email protected]>
Cc: "geoserver-users" <[email protected]>
Sent: Wednesday, 23 September, 2015 14:04:20
Subject: Re: [Geoserver-users] Paletted geotiff displayed as white
Hi,
I found what makes failure with the 2005.shp index file. It is the .prj
PROJCS["KKJ / Finland Uniform Coordinate System", GEOGCS["KKJ",
DATUM["Kartastokoordinaattijarjestelma_1966", SPHEROID["International 1924",
6378388.0, 297.0, AUTHORITY["EPSG","7022"]], TOWGS84[-96.062, -82.428,
-121.753, 4.801, 0.345, -1.376, 1.496], AUTHORITY["EPSG","6123"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree",
0.017453292519943295], AXIS["Longitude", EAST], AXIS["Latitude", NORTH],
AUTHORITY["EPSG","4123"]], PROJECTION["Transverse_Mercator"],
PARAMETER["central_meridian", 27.0], PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 3500000.0],
PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["X", NORTH],
AXIS["Geocentric Y", EAST], AUTHORITY["EPSG","2393"]]
When I create a new mosaic with Geoserver it write .prj like this:
PROJCS["KKJ / Finland Uniform Coordinate System", GEOGCS["KKJ",
DATUM["Kartastokoordinaattijarjestelma (1966)", SPHEROID["International 1924",
6378388.0, 297.0, AUTHORITY["EPSG","7022"]], TOWGS84[-96.062, -82.428,
-121.753, 4.801, -0.345, 1.376, 1.496], AUTHORITY["EPSG","6123"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree",
0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic
latitude", NORTH], AUTHORITY["EPSG","4123"]],
PROJECTION["Transverse_Mercator"], PARAMETER["central_meridian", 27.0],
PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 3500000.0], PARAMETER["false_northing", 0.0],
UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH],
AUTHORITY["EPSG","2393"]]
By using the latter .prj instead of the original 2005.prj the mosaic shows.
Perhaps these odd axis definitions lead to failure:
AXIS["X", NORTH], AXIS["Geocentric Y", EAST]
I have been thinking why the image mosaic system of GeoServer is making so much
trouble while tileindex in Mapserver never fails. Any ideas about how to make
image mosaic more robust or at least easier to debug? I had DEBUG level
logfiles and I studied them for a long time but it did not help me. I found the
issue with the odd .prj axis definitions only by looking at the contents with a
text editor.
-Jukka Rahkonen-
Rahkonen Jukka wrote:
Hi,
The case is not so simple. I got white mosaic as well (GS 2.7.1 on Windows),
but when I deleted all 2005.* files and indexer.properties file which were
included in the zip and re-created the mosaic again from the scratch and voilà,
it works. However, the new “kylli.properties” file has also
“ExpandToRGB=false”. The whole properties file below:
#-Automagically created from GeoTools-
#Wed Sep 23 10:57:36 EEST 2015
Levels=100.06253908692933,100.06253908692933
Heterogeneous=false
AbsolutePath=false
Name=kylli
TypeName=kylli
Caching=false
ExpandToRGB=false
LocationAttribute=location
SuggestedSPI=it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReaderSpi
CheckAuxiliaryMetadata=false
LevelsNum=1
I suggest to repeat the test by deleting the existing mosaic shapefile and
letting your Geoserver to create a new one.
File 32_33op.tif that was supposed to be while looks OK for me.
-Jukka Rahkonen-
Andrea Aime wrote:
Hi,
made a quick test locally with a 2.8.x I had handy, checked it,
and indeed it's white.
Then I've looked into its property config file of the mosaic and found this:
ExpandToRGB=false
Switched it to ExpandToRGB=true and reloaded the geoserver configuration,
and voilà, it works (well, almost, I still see some vertical black lines where
the tiles match up, not sure if it's the data or some other issue).
Could you open a bug report at https://osgeo-org.atlassian.net/projects/GEOS ,
it might be that we are wrongly detecting the palettes to be all the same,
or there might be some issue in the management of the paletted images
processing chain in the mosaic
Cheers
Andrea
On Wed, Sep 23, 2015 at 8:47 AM, Kylli Ek < [email protected] > wrote:
Hi,
The data is available here:
https://filesender.funet.fi/?vid=01ebcfb1-a5e0-9149-7a24-00000f6ba256
For example mapsheets 32_33op and 30_31mn are white, but 24_25op is ok.
Some test layers are visible here:
http://avoin-test.csc.fi/geoserver/test/wms ?
Kylli
From: "Rahkonen Jukka (MML)" < [email protected] >
To: "Kylli Ek" < [email protected] >, "geoserver-users" <
[email protected] >
Sent: Tuesday, 22 September, 2015 18:05:26
Subject: Re: Paletted geotiff displayed as white
Hi Kylli,
Could it be possible to put a sample image somewhere awailable?
-Jukka Rahkonen-
Kylli Ek wrote:
Hi,
We are trying to publish to GeoServer some older paletted tiff-images but have
a strange problem - GeoServer returns a totally white image, no error messages
(checked the logs).
First we tried with ImageMosaic, then also single images. Some of the single
images were showing up normally, some are also white, although with gdalinfo I
do not see any differences for these images and they are from same vendor from
same delivery (example below).
The files are paletted geotiffes and are about 10 years old, but for example in
QGIS they display without any problems. According to this post (
http://osgeo-org.1560.x6.nabble.com/Publishing-GeoTiffs-ColorPalette-with-Alpha-td5128601.html
), I could understand, that paletted geotiffs with alpha channel can not be
published with GeoServer? But some of the images do show up??
* I tried also changing the style, so that certain values would have certain
color, but no change. The images showing up have the default raster style. .
Clicking on the map preview does not give any results, could think that if for
some reason whole image is just pained white, it still would give values for
GetFeatureInfo request?
* Tried also setting Coverage Band Details minRange and MaxRange, but no
change. By default both were 0 for imagemosaic, but -indefinite to indefinite
for single image.
* The granule grid for ImageMosaic looks good. And the coordinates for Preview
seem ok.
* GeoServer version 2.7.1
gdalinfo gives:
gdalinfo 28_29op.tif
Driver: GTiff/GeoTIFF
Files: 28_29op.tif
28_29op.aux
28_29op.rrd
28_29op.tfw
Size is 1600, 1600
Coordinate System is:
PROJCS["KKJ / Finland Uniform Coordinate System",
GEOGCS["KKJ",
DATUM["Kartastokoordinaattijarjestelma_1966",
SPHEROID["International 1924",6378388,297.0000000000014,
AUTHORITY["EPSG","7022"]],
AUTHORITY["EPSG","6123"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4123"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",27],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",3500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","2393"]]
Origin = (3339950.000000000000000,7519950.000000000000000)
Pixel Size = (100.062539086929334,-100.062539086929334)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
TIFFTAG_SOFTWARE=Arc/Info
TIFFTAG_XRESOLUTION=0.0099999998
TIFFTAG_YRESOLUTION=0.0099999998
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 3339950.000, 7519950.000) ( 23d12'59.77"E, 67d43'18.78"N)
Lower Left ( 3339950.000, 7359849.937) ( 23d25'59.61"E, 66d17'21.51"N)
Upper Right ( 3500050.063, 7519950.000) ( 27d 0' 4.27"E, 67d45'56.68"N)
Lower Right ( 3500050.063, 7359849.937) ( 27d 0' 4.02"E, 66d19'48.82"N)
Center ( 3420000.031, 7439899.969) ( 25d 9'46.89"E, 67d 2'14.75"N)
Band 1 Block=1600x6 Type=Byte, ColorInterp=Palette
Description = Band_1
Overviews: 400x400, 200x200, 100x100, 50x50
Metadata:
LAYER_TYPE=thematic
Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 241,148,241,255
2: 223,255,202,255
3: 255,0,0,255
4: 186,223,234,255
etc until 255
Best regards,
Kylli Ek
GIS-coordinator
CSC - IT Center for Science Ltd
PL 405, 02101 Espoo
Tel: +358 50 38 12 838
[email protected]
------------------------------------------------------------------------------
Monitor Your Dynamic Infrastructure at Any Scale With Datadog!
Get real-time metrics from all of your servers, apps and tools
in one place.
SourceForge users - Click here to start your Free Trial of Datadog now!
http://pubads.g.doubleclick.net/gampad/clk?id=241902991&iu=/4140
_______________________________________________
Geoserver-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geoserver-users
--
==
GeoServer Professional Services from the experts! Visit
http://goo.gl/it488V for more information.
==
Ing. Andrea Aime
@geowolf
Technical Lead
GeoSolutions S.A.S.
Via Poggio alle Viti 1187
55054 Massarosa (LU)
Italy
phone: +39 0584 962313
fax: +39 0584 1660272
mob: +39 339 8844549
http://www.geo-solutions.it
http://twitter.com/geosolutions_it
AVVERTENZE AI SENSI DEL D.Lgs. 196/2003
Le informazioni contenute in questo messaggio di posta elettronica e/o nel/i
file/s allegato/i sono da considerarsi strettamente riservate. Il loro utilizzo
è consentito esclusivamente al destinatario del messaggio, per le finalità
indicate nel messaggio stesso. Qualora riceviate questo messaggio senza esserne
il destinatario, Vi preghiamo cortesemente di darcene notizia via e-mail e di
procedere alla distruzione del messaggio stesso, cancellandolo dal Vostro
sistema. Conservare il messaggio stesso, divulgarlo anche in parte,
distribuirlo ad altri soggetti, copiarlo, od utilizzarlo per finalità diverse,
costituisce comportamento contrario ai principi dettati dal D.Lgs. 196/2003.
The information in this message and/or attachments, is intended solely for the
attention and use of the named addressee(s) and may be confidential or
proprietary in nature or covered by the provisions of privacy act (Legislative
Decree June, 30 2003, no.196 - Italy's New Data Protection Code).Any use not in
accord with its purpose, any disclosure, reproduction, copying, distribution,
or either dissemination, either whole or partial, is strictly forbidden except
previous formal approval of the named addressee(s). If you are not the intended
recipient, please contact immediately the sender by telephone, fax or e-mail
and delete the information in this message that has been received in error. The
sender does not give any warranty or accept liability as the content, accuracy
or completeness of sent messages and accepts no responsibility for changes made
after they were sent or for other risks which arise as a result of e-mail
transmission, viruses, etc.
-------------------------------------------------------
------------------------------------------------------------------------------
Monitor Your Dynamic Infrastructure at Any Scale With Datadog!
Get real-time metrics from all of your servers, apps and tools
in one place.
SourceForge users - Click here to start your Free Trial of Datadog now!
http://pubads.g.doubleclick.net/gampad/clk?id=241902991&iu=/4140
_______________________________________________
Geoserver-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geoserver-users
