Changes between Version 8 and Version 9 of Ticket #1318
- Timestamp:
- May 24, 2016, 10:07:59 AM (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #1318
- Property Owner set to
- Property Status new → accepted
-
Ticket #1318 – Description
v8 v9 23 23 ---- 24 24 25 26 The coverage is 3D and slicing in Time axis to 2D and native CRS is 4326, it is valid to project to output 3857. However, in handle class, such as GeotiffFormatExtension.java. What it was done is: 25 E.g: WCS on a coverage which has compoundCRS 27 26 28 27 29 28 {{{ 30 // First, transform possible non-native CRS subsets 31 CRSExtension crsExtension = (CRSExtension) ExtensionsRegistry.getExtension(ExtensionsRegistry.CRS_IDENTIFIER); 32 crsExtension.handle(request, m, meta); 33 34 //Handle the range subset feature 35 RangeSubsettingExtension rsubExt = (RangeSubsettingExtension) ExtensionsRegistry.getExtension(ExtensionsRegistry.RANGE_SUBSETTING_IDENTIFIER); 36 rsubExt.handle(request, m); 37 29 service=WCS&version=2.0.1&request=GetCoverage&coverageId=eobstest& 30 subset=Lat(45,60)&subset=Long(25,30)&subset=t(0)&format=image/tiff& 31 outputCrs=http://localhost:8080/def/crs/EPSG/0/4269 38 32 }}} 39 33 40 So it needs the 2D CRS (4326) from the coverage before the coverage is sliced from 3D to 2D. Then in CrsExtension.java, it will throw error due to the compound CRS like this is not valid 34 will return a Rasql with error in CRS of projection() 41 35 42 36 43 37 {{{ 44 http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/3857&2=http://localhost:8080/def/crs/OGC/0/AnsiDate" 38 service=WCS&version=2.0.1&request=GetCoverage&coverageId=eobstest 39 &subset=Lat(45,60)&subset=Long(25,30)&subset=t(0)&format=image/tiff 40 &outputCrs=http://localhost:8080/def/crs/EPSG/0/4269 45 41 }}} 46 47 for GeoTools.JTS library to transform (this library works well, before I had a mistake with parameters then it is transformed to wrong coordiniates).48 49 So, what need to change here:50 + Get the right grid/geo CRS (4326) in compound CRS (it is tricky as the order in compoundCRS is not consistent, e.g: 1=time&2=4326 or 1=4326&2=time)51 or52 + Handle the coverage first to make it 2D and in GetCoverageMetadata (somehow it is handled correctly to get only 2D axis (4326)). Then it is input for the CRSExtension class.53 42 54 43 55 44 56 45 46 47