Petascope ingestion of non-geo-referenced images

[Alex Dumitru]

When importing this GML (​http://pastebin.com/RfhmM4T1) I get a coverage that starts at (1,1), although the origin is set at (0.5, 0.5). Is this correct behavior? I am assuming some rounding is performed somewhere for Index2D CRSes.

[Vlad Merticariu]

Can you attach the coverage description after ingestion?

[Bang Pham Huu]

I think the reason is based on this function and also want to know why need to do this? Please answer if you have time.

  def get_origin(self):
	        """
	        Returns the origin of the dataset. This is calculated using the origin from the dataset + 0.5 of an offset
	                vector, as petascope requires it so.
	        :rtype: list[str]
	        """
	        geo = self.gdal_dataset.GetGeoTransform()
	        # Add 0.5 * size of pixel to make sure there's no error in petascope
	        # Should be removed once petascope handles numeric computations correctly
	        return str(geo[0] + 0.5 * self.gdal_dataset.GetGeoTransform()[1]), \
	               str(geo[3] + 0.5 * self.gdal_dataset.GetGeoTransform()[5])

[Alex Dumitru]

@Bang: The origin in petascope is considered to be in the middle of the geopixel, that is top-left corner + 0.5 offset vector. The whole geo-model is described here ​http://www.rasdaman.org/wiki/PetascopeUserGuide
The comment delimited by # is not correct, only the one between """ """ is.

This is a bug in petascope.wcst, not in WCSTImport, the created coverage is correct.

[Bang Pham Huu]

Thanks Alex for your time I will read the article, I also would like to add the file as Vlad requested.

[Alex Dumitru]

Please leave it assigned to Vlad. This is more of a standards decision than actual implementation.
@Vlad Should the origin be in the middle of the pixel for non-georef grids as well?