tiff produces invalid file

[Dimitar Misev]

Reported on the users mailing list.

Data can be downloaded ​here, and the recipe for ingestion is below:

{
  "config": {
    "service_url": "http://localhost:8080/rasdaman/ows",
    "tmp_directory": "/tmp/importmp/",
    "crs_resolver": "http://localhost:8080/def/",
    "default_crs": "http://www.opengis.net/def/crs/EPSG/0/4258",
    "mock": false,
    "automated": false,
    "subset_correction": true
  },
  "input": {
    "coverage_id": "wcs_sst",
    "paths": [
      "19_15_26_30/*.nc"
    ]
  },
  "recipe": {
    "name": "time_series_irregular",
    "options": {
      "time_parameter": {
        "metadata_tag": {
          "tag_name": "NC_GLOBAL#NC_GLOBAL.start_date"
        },
        "datetime_format": "YYYY-MM-DD"
      },
      "time_crs": "http://opengis.net/def/crs/OGC/0/AnsiDate",
      "tiling": "ALIGNED [0:10, 0:1000, 0:500]"
    }
  }
}

This query produces a 91x91 tiff file, which is definitely not correct.

rasql -q 'select tiff(c[-274:90,273:364,4328]) from wcs_sst AS c' --out file 

@Bang - could you try it when you don't have any other tasks?

[Bang Pham Huu]

Sorry Dimitar for late replying.

I could not use your recipe with error in ANSI time but I think the time in here is not important. Yes, the query result correct is '365 x 92'

here is my recipe

{
  "config": {
    "service_url": "http://localhost:8080/rasdaman/ows",
    "tmp_directory": "/tmp/",
    "crs_resolver": "http://localhost:8080/def/",
    "default_crs": "http://localhost:8080/def/crs/EPSG/0/4258",
    "mock": false,
    "automated": false,
    "subset_correction": true
  },
  "input": {
    "coverage_id": "wcs_sst2",
    "paths": [
      "/home/rasdaman/test_dimitar/*.nc"
    ]
  },
  "recipe": {
    "name": "time_series_regular",
    "options": {
       "time_start": "2012-12-02T20:12:02",
       "time_format": "auto",
       "time_crs": "http://localhost:8080/def/crs/OGC/0/AnsiDate",
       "time_step": "2 days 10 minutes 3 seconds",
       "tiling": "ALIGNED [0:1000, 0:1000, 0:2]"
    }
  }
}

please view attachment file for the detail

[Bang Pham Huu]

Describe coverage if you need

<?xml version="1.0" encoding="UTF-8"?><wcs:CoverageDescriptions xmlns:wcs="http://www.opengis.net/wcs/2.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:crs="http://www.opengis.net/wcs/service-extension/crs/1.0" xmlns:ows="http://www.opengis.net/ows/2.0" xmlns:gml="http://www.opengis.net/gml/3.2" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/wcs/2.0 http://schemas.opengis.net/wcs/2.0/wcsAll.xsd">
  <wcs:CoverageDescription xmlns="http://www.opengis.net/gml/3.2" xmlns:gmlcov="http://www.opengis.net/gmlcov/1.0" xmlns:swe="http://www.opengis.net/swe/2.0" gml:id="wcs_sst2">
    <boundedBy>
      <Envelope srsName="http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/4258&amp;2=http://localhost:8080/def/crs/OGC/0/AnsiDate" axisLabels="Lat Long ansi" uomLabels="degree degree d" srsDimension="3">
        <lowerCorner>26.00460616027095 -18.99539384107375 150450.838200231478435077667236328125</lowerCorner>
        <upperCorner>30.00369271366015 -14.99630633138375 150468.901012731508435077667236328125</upperCorner>
      </Envelope>
    </boundedBy>
    <wcs:CoverageId>wcs_sst2</wcs:CoverageId>
    <gmlcov:metadata/>
    <domainSet>
      <RectifiedGrid dimension="3" gml:id="wcs_sst2-grid">
        <limits>
          <GridEnvelope>
            <low>0 0 0</low>
            <high>363 363 8</high>
          </GridEnvelope>
        </limits>
        <axisLabels>Long Lat ansi</axisLabels>
        <origin>
          <Point gml:id="wcs_sst2-origin" srsName="http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/4258&amp;2=http://localhost:8080/def/crs/OGC/0/AnsiDate" axisLabels="Lat Long ansi" uomLabels="degree degree d" srsDimension="3">
            <pos>29.9981994629 -18.989900589 150451.841689814835</pos>
          </Point>
        </origin>
        <offsetVector srsName="http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/4258&amp;2=http://localhost:8080/def/crs/OGC/0/AnsiDate" axisLabels="Lat Long ansi" uomLabels="degree degree d" srsDimension="3">0 0.0109865041475 0</offsetVector>
        <offsetVector srsName="http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/4258&amp;2=http://localhost:8080/def/crs/OGC/0/AnsiDate" axisLabels="Lat Long ansi" uomLabels="degree degree d" srsDimension="3">-0.0109865015203 0 0</offsetVector>
        <offsetVector srsName="http://localhost:8080/def/crs-compound?1=http://localhost:8080/def/crs/EPSG/0/4258&amp;2=http://localhost:8080/def/crs/OGC/0/AnsiDate" axisLabels="Lat Long ansi" uomLabels="degree degree d" srsDimension="3">0 0 2.00697916667</offsetVector>
      </RectifiedGrid>
    </domainSet>
    <gmlcov:rangeType>
      <swe:DataRecord>
        <swe:field name="field_1">
          <swe:Quantity>
            <swe:nilValues>
              <swe:NilValues/>
            </swe:nilValues>
            <swe:uom code="10^0"/>
          </swe:Quantity>
        </swe:field>
      </swe:DataRecord>
    </gmlcov:rangeType>
    <wcs:ServiceParameters>
      <wcs:CoverageSubtype>RectifiedGridCoverage</wcs:CoverageSubtype>
      <CoverageSubtypeParent xmlns="http://www.opengis.net/wcs/2.0">
        <CoverageSubtype>AbstractDiscreteCoverage</CoverageSubtype>
        <CoverageSubtypeParent>
          <CoverageSubtype>AbstractCoverage</CoverageSubtype>
        </CoverageSubtypeParent>
      </CoverageSubtypeParent>
      <wcs:nativeFormat>application/octet-stream</wcs:nativeFormat>
    </wcs:ServiceParameters>
  </wcs:CoverageDescription>
</wcs:CoverageDescriptions>

[Bang Pham Huu]

Sorry yesterday I can import with Dimitar's recipe, the result is the same with size 365x92. All the detail, please view here

​http://pastebin.com/CTnmFfTC

[Dimitar Misev]

Ok I'll have a look directly on Ricardo's server tomorrow.

[Bang Pham Huu]

Can I do anything else for you Dimitar? Because I don't know clearly the reason also result. If I can do something else or it is invalid, I may want to know as you assigned.

[Dimitar Misev]

Continuing on ticket #1049