Hello, I'm trying to reproject on-the-fly some geotiffs from custom coordinate system, defined in "epsg.properties" file to the EPSG:4326. I can do that for single image, defined in a GeoTIFF store, however the same image, imported by ImageMosaic store can't be displayed properly.
Apparently, ImageMosaic doesn't make right reprojection - image is displayed in original coordinates (or throws a "Dimensions (width=87581 height=87595) are too large" error, depending on the original CRS - rotated pole degrees or meters). More setup details: CRS settings: "Force declared" (from EPSG:404000 to custom, eg. EPSG:100102) 100102=PROJCS["WRF Lambert",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenw ich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic longitude",EAST],AXIS["Geodetic latitude",NORTH],AUTHORITY["EPSG","4326"]],PROJECTION["Lambert_Confor mal_Conic_2SP"],PARAMETER["central_meridian",18.20199966430664],PARAMETER["latitude_of_origin",54.13700866699219],PARAMETER["standard_parallel_1",30.0],PARAMETER["false_easting",0.0],PARA METER["false_northing",0.0],PARAMETER["standard_parallel_2",60.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH]] ImageMosaic boundaries handled by PostGis store. Below I'm attaching appropriate parts of the log file - same request, two different storage. Regards Maciej ************** GeoTIFF store: ************** 2015-06-09 12:27:19,247 INFO [org.geoserver.wms] - Request: getMap Tiled = false Width = 2048 Height = 2048 Layers = [org.geoserver.wms.MapLayerInfo@aeb606aa] Crs = GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]] Bbox = SRSEnvelope[-16.50941800460898 : 62.99640914526753, 31.87980262617571 : 64.40683143260988] RemoteOwsType = null RemoteOwsURL = null Env = {} FormatOptions = {AUTOFIT=true, KMATTR=true, KMPLACEMARK=false, KMSCORE=40, MODE=refresh, SUPEROVERLAY=false} Angle = 0.0 CQLFilter = null Elevation = [] FeatureId = null StartIndex = null ViewParams = null MaxFeatures = null Palette = null FeatureVersion = null SRS = EPSG:4326 Styles = [StyleImpl[ name=temperatura_barn1]] BgColor = java.awt.Color[r=255,g=255,b=255] Transparent = false StyleFormat = sld SldBody = null Interpolations = [] ValidateSchema = false Sld = null StyleUrl = null TilesOrigin = null Exceptions = SE_XML StyleBody = null SldVersion = null StyleVersion = null ScaleMethod = null Format = application/vnd.google-earth.kml xml Time = [] Filter = null Buffer = 0 Filters = null Get = true RawKvp = {TRANSPARENT=false, BBOX=-16.50941800460898,31.87980262617571,62.99640914526753,64.40683143260988, VERSION=1.1.1, FORMAT=application/vnd.google-earth.kml xml, FORMAT_OPTIONS=AUTOFIT:true;KMATTR:true;KMPLACEMARK:false;KMSCORE:40;MODE:refresh;SUPEROVERLAY:false, SERVICE=wms, HEIGHT=2048, REQUEST=GetMap, LAYERS=um:wrf_d02_T2_201504040600_005, STYLES=temperatura_barn1, WIDTH=2048, SRS=EPSG:4326} BaseUrl = http://geo.nebula.grid.icm.edu.pl:8080/geoserver/ RequestCharset = UTF-8 Version = 1.1.1 Request = GetMap 2015-06-09 12:27:19,307 INFO [org.geoserver.wms] - Request: getServiceInfo 2015-06-09 12:27:19,308 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:27:19,308 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only. 2015-06-09 12:27:19,308 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:27:19,308 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only. 2015-06-09 12:27:19,310 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory) 2015-06-09 12:27:19,310 TRACE [org.geotools.factory] - RETURN (GridCoverageFactory): found implementation GridCoverageFactory. 2015-06-09 12:27:19,311 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:27:19,311 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:27:19,311 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:27:19,313 DEBUG [org.geotools.data.coverage.grid] - Coverage read: width = 519 height = 795 2015-06-09 12:27:19,313 TRACE [org.geotools.factory] - ENTRY (CoordinateOperationFactory, COORDINATE_OPERATION_FACTORY) 2015-06-09 12:27:19,314 TRACE [org.geotools.factory] - RETURN (CoordinateOperationFactory, COORDINATE_OPERATION_FACTORY): found implementation BufferedCoordinateOperationFactory. 2015-06-09 12:27:19,327 DEBUG [org.geotools.referencing.operation.transform] - Optimizing the warp into an grid warp 67 x 33 2015-06-09 12:27:19,328 DEBUG [org.geotools.coverage.processing] - Resampled coverage "wrf_d02_T2_201504040600_005" from coordinate system "WRF Lambert" (for an image of size 519×795) to coordinate system "WGS 84" (image size 519×795). JAI operation is "Warp" with "Nearest" interpolation on packed pixels values. Background value is (255, 255, 255). 2015-06-09 12:27:19,328 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Reprojecting to crs GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]] 2015-06-09 12:27:19,328 TRACE [org.geotools.factory] - ENTRY (MathTransformFactory, MATH_TRANSFORM_FACTORY) 2015-06-09 12:27:19,328 TRACE [org.geotools.factory] - RETURN (MathTransformFactory, MATH_TRANSFORM_FACTORY): found implementation DefaultMathTransformFactory. 2015-06-09 12:27:19,341 DEBUG [org.geotools.referencing.operation.transform] - Optimizing the warp into an grid warp 65 x 32 2015-06-09 12:27:19,341 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Applying Raster Symbolizer ***************** ImageMosaic store ***************** 2015-06-09 12:29:16,559 INFO [org.geoserver.wms] - Request: getMap Tiled = false Width = 2048 Height = 2048 Layers = [org.geoserver.wms.MapLayerInfo@87d5335b] Crs = GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]] Bbox = SRSEnvelope[-16.50941800460898 : 62.99640914526753, 31.87980262617571 : 64.40683143260988] RemoteOwsType = null RemoteOwsURL = null Env = {} FormatOptions = {AUTOFIT=true, KMATTR=true, KMPLACEMARK=false, KMSCORE=40, MODE=refresh, SUPEROVERLAY=false} Angle = 0.0 CQLFilter = null Elevation = [] FeatureId = null StartIndex = null ViewParams = null MaxFeatures = null Palette = null FeatureVersion = null SRS = EPSG:4326 Styles = [StyleImpl[ name=temperatura_barn1]] BgColor = java.awt.Color[r=255,g=255,b=255] Transparent = false StyleFormat = sld SldBody = null Interpolations = [] ValidateSchema = false Sld = null StyleUrl = null TilesOrigin = null Exceptions = SE_XML StyleBody = null SldVersion = null StyleVersion = null ScaleMethod = null Format = application/vnd.google-earth.kml xml Time = [] Filter = null Buffer = 0 Filters = null Get = true RawKvp = {TRANSPARENT=false, BBOX=-16.50941800460898,31.87980262617571,62.99640914526753,64.40683143260988, VERSION=1.1.1, FORMAT=application/vnd.google-earth.kml xml, FORMAT_OPTIONS=AUTOFIT:true;KMATTR:true;KMPLACEMARK:false;KMSCORE:40;MODE:refresh;SUPEROVERLAY:false, SERVICE=wms, HEIGHT=2048, REQUEST=GetMap, LAYERS=wrf:WRF2_T2, STYLES=temperatura_barn1, WIDTH=2048, SRS=EPSG:4326} BaseUrl = http://geo.nebula.grid.icm.edu.pl:8080/geoserver/ RequestCharset = UTF-8 Version = 1.1.1 Request = GetMap 2015-06-09 12:29:16,606 INFO [org.geoserver.wms] - Request: getServiceInfo 2015-06-09 12:29:16,607 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,607 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only. 2015-06-09 12:29:16,607 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,607 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only. 2015-06-09 12:29:16,609 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,609 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,609 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,609 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,609 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory) 2015-06-09 12:29:16,609 TRACE [org.geotools.factory] - RETURN (GridCoverageFactory): found implementation GridCoverageFactory. 2015-06-09 12:29:16,609 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:29:16,609 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:29:16,609 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@9704ee 2015-06-09 12:29:16,610 DEBUG [org.geotools.gce.imagemosaic] - Reading mosaic from file:/var/lib/geoserver/data_dir/data/wrf/WRF2_T2/datastore.properties 2015-06-09 12:29:16,610 DEBUG [org.geotools.gce.imagemosaic] - Highest res 3400.0 3400.0 2015-06-09 12:29:16,610 TRACE [org.geotools.gce.imagemosaic] - Adjusted Requested Envelope = ReferencedEnvelope[-17.24702089320647 : 63.811654443191074, 31.06468615250612 : 65.14431776592998] Requested raster dimension = java.awt.Rectangle[x=-10,y=-10,width=2088,height=878] Corresponding raster source region = java.awt.Rectangle[x=-10,y=-10,width=2088,height=878] Corresponding source Envelope = ReferencedEnvelope[-17.24702089320647 : 63.811654443191074, 31.06468615250612 : 65.14431776592998] 2015-06-09 12:29:16,610 DEBUG [org.geotools.gce.imagemosaic] - Loading level 0 with subsampling factors 1 1 2015-06-09 12:29:16,610 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,610 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,610 DEBUG [org.geotools.gce.imagemosaic] - Prepping to manage SortBy Clause 2015-06-09 12:29:16,610 DEBUG [org.geotools.gce.imagemosaic.catalog] - Index Loaded 2015-06-09 12:29:16,610 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,610 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,611 DEBUG [org.geotools.jdbc] - CREATE CONNECTION 2015-06-09 12:29:16,611 INFO [org.geotools.jdbc] - Unable to find aggregate function matching visitor: class org.geotools.gce.imagemosaic.catalog.CachingDataStoreGranuleCatalog$1 2015-06-09 12:29:16,611 DEBUG [org.geotools.jdbc] - CLOSE CONNECTION 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,611 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - ENTRY (GeometryFactory, JTS_GEOMETRY_FACTORY) 2015-06-09 12:29:16,612 TRACE [org.geotools.factory] - RETURN (GeometryFactory, JTS_GEOMETRY_FACTORY): found implementation GeometryFactory. 2015-06-09 12:29:16,612 DEBUG [org.geotools.jdbc] - CREATE CONNECTION 2015-06-09 12:29:16,612 DEBUG [org.geotools.filter] - exporting PropertyName 2015-06-09 12:29:16,612 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only. 2015-06-09 12:29:16,612 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only. 2015-06-09 12:29:16,612 DEBUG [org.geotools.jdbc] - SELECT "fid",encode(ST_AsEWKB("the_geom"),'base64') as "the_geom","location","time","forecast" FROM "public"."WRF2_T2" WHERE "the_geom" && ST_GeomFromWKB(?, 404000) 2015-06-09 12:29:16,612 DEBUG [org.geotools.jdbc] - 1 = POLYGON ((-17.24702089320647 31.06468615250612, -17.24702089320647 65.14431776592998, 63.811654443191074 65.14431776592998, 63.811654443191074 31.06468615250612, -17.24702089320647 31.06468615250612)) 2015-06-09 12:29:16,614 DEBUG [org.geotools.gce.imagemosaic] - Thread:http-apr-8080-exec-3 Loading raster data for granuleDescriptor Description of a granuleDescriptor BBOX: ReferencedEnvelope[-941667.6 : 822932.4, -1106746.8 : 1596253.2] file: file:/var/lib/geoserver/data_dir/data/wrf/WRF2_T2/wrf_d02_T2_201504040600_005.tiff gridToWorld: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3], PARAMETER["elt_0_0", 3400.0], PARAMETER["elt_0_2", -939967.6], PARAMETER["elt_1_1", -3400.0], PARAMETER["elt_1_2", 1594553.2]] Description of level 2 Description of a granuleDescriptor level width: 519 height: 795 scaleX: 1.0 scaleY: 1.0 baseToLevelTransform: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3]] gridToWorldTransform: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3], PARAMETER["elt_0_0", 3400.0], PARAMETER["elt_0_2", -941667.6], PARAMETER["elt_1_1", -3400.0], PARAMETER["elt_1_2", 1596253.2]] 2015-06-09 12:29:16,614 DEBUG [org.geotools.gce.imagemosaic] - Loading level 0 with source region: java.awt.Rectangle[x=274,y=467,width=5,height=5] subsampling: 1,1 for granule:file:/var/lib/geoserver/data_dir/data/wrf/WRF2_T2/wrf_d02_T2_201504040600_005.tiff 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - We added the granule Description of a granuleDescriptor BBOX: ReferencedEnvelope[-941667.6 : 822932.4, -1106746.8 : 1596253.2] file: file:/var/lib/geoserver/data_dir/data/wrf/WRF2_T2/wrf_d02_T2_201504040600_005.tiff gridToWorld: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3], PARAMETER["elt_0_0", 3400.0], PARAMETER["elt_0_2", -939967.6], PARAMETER["elt_1_1", -3400.0], PARAMETER["elt_1_2", 1594553.2]] Description of level 2 Description of a granuleDescriptor level width: 519 height: 795 scaleX: 1.0 scaleY: 1.0 baseToLevelTransform: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3]] gridToWorldTransform: PARAM_MT["Affine", PARAMETER["num_row", 3], PARAMETER["num_col", 3], PARAMETER["elt_0_0", 3400.0], PARAMETER["elt_0_2", -941667.6], PARAMETER["elt_1_1", -3400.0], PARAMETER["elt_1_2", 1596253.2]] 2015-06-09 12:29:16,615 DEBUG [org.geotools.jdbc] - CLOSE CONNECTION 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - Producing the final mosaic, step 1, loop through granule collectors 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - Using collector with filter:Filter.INCLUDE 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - Adding to mosaic granule file:/var/lib/geoserver/data_dir/data/wrf/WRF2_T2/wrf_d02_T2_201504040600_005.tiff 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - Producing the final mosaic, step 2, final mosaicking 2015-06-09 12:29:16,615 DEBUG [org.geotools.gce.imagemosaic] - Loaded bbox ReferencedEnvelope[-17.24702089320647 : 63.811654443191074, 31.06468615250612 : 65.14431776592998] while crop bbox ReferencedEnvelope[-17.24702089320647 : 63.811654443191074, 31.06468615250612 : 65.14431776592998] 2015-06-09 12:29:16,617 INFO [org.geotools.coverage.processing] - Apply operation "CoverageCrop" on coverage "WRF2_T2" with interpolation "Nearest". 2015-06-09 12:29:16,617 TRACE [org.geotools.factory] - ENTRY (CoordinateOperationFactory, COORDINATE_OPERATION_FACTORY) 2015-06-09 12:29:16,617 TRACE [org.geotools.factory] - RETURN (CoordinateOperationFactory, COORDINATE_OPERATION_FACTORY): found implementation BufferedCoordinateOperationFactory. 2015-06-09 12:29:16,617 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Reprojecting to crs GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]] 2015-06-09 12:29:16,621 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Applying Raster Symbolizer 2015-06-09 12:29:16,621 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory) 2015-06-09 12:29:16,621 TRACE [org.geotools.factory] - RETURN (GridCoverageFactory): found implementation GridCoverageFactory. 2015-06-09 12:29:16,621 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory) 2015-06-09 12:29:16,621 TRACE [org.geotools.factory] - RETURN (GridCoverageFactory): found implementation GridCoverageFactory. 2015-06-09 12:29:16,621 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory) ------------------------------------------------------------------------------ _______________________________________________ Geoserver-users mailing list Geoserver-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/geoserver-users