| | 145 | == Set CRS for output coverage |
| | 146 | Suppose you have a 3D coverage which has 3 axes (Lat, Long and ansi) and the compoundCRS for this coverage is: |
| | 147 | |
| | 148 | {{{ |
| | 149 | http://localhost:8080/def/crs-compound? |
| | 150 | 1=http://localhost:8080/def/crs/EPSG/0/4326& |
| | 151 | 2=http://localhost:8080/def/crs/OGC/0/AnsiDate |
| | 152 | }}} |
| | 153 | |
| | 154 | so when you want to slice on the ansi axis (i.e: subset with slicing on this axis, example: |
| | 155 | {{{ |
| | 156 | localhost:8080/rasdaman/ows?&SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&COVERAGEID=test_AverageTemperature&SUBSET=ansi(150455)&FORMAT=image/tiff |
| | 157 | }}} |
| | 158 | the output image will have EPSG:4326 as geo-referenced CRS. |
| | 159 | |
| | 160 | {{{ |
| | 161 | Coordinate System is: |
| | 162 | GEOGCS["WGS 84", |
| | 163 | DATUM["WGS_1984", |
| | 164 | SPHEROID["WGS 84",6378137,298.257223563, |
| | 165 | AUTHORITY["EPSG","7030"]], |
| | 166 | AUTHORITY["EPSG","6326"]], |
| | 167 | PRIMEM["Greenwich",0], |
| | 168 | UNIT["degree",0.0174532925199433], |
| | 169 | AUTHORITY["EPSG","4326"]] |
| | 170 | }}} |
| | 171 | |
| | 172 | However, if you want to slice on Latitude axis or Longitude axis. |
| | 173 | {{{ |
| | 174 | localhost:8080/rasdaman/ows?&SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&COVERAGEID=test_AverageTemperature&SUBSET=Lat(0)&FORMAT=image/tiff |
| | 175 | }}} |
| | 176 | |
| | 177 | the output image will have axes combined from Long&ansi or Lat&ansi (depend on which axis you want to slice above (Long or Lat)). ''There is non-existing CRS which made of a geo-referenced axis (Long or Lat) and a time axis (ansi).'' '''Therefore, the output coverage in 2D will have no CRS'''. |
| | 178 | |
| | 179 | |
| | 180 | |
| | 181 | |