NETCDF export broken - file has incorrect orientation

[olcl@…]

when exporting data in the NetCDF format the resulting image is rotated incorrectly. However the same data when exported in either PNG of GeoTIFF is correctly orientated.

the query below is used to generate the image and example output files are included.

you can run this query at ​http://earthserver.pml.ac.uk/petascope/wcps, by replacing the tiff with netcdf or png you can test the other format images.

for c in (myocean_test)
return
encode (
c[t(29)], "tiff"
)

We are using the latest RPM's, version 8.3.1-0.el6

[olcl@…]

screen shot showing difference between png and netcdf output

[Dimitar Misev]

When I look at the file produced by this query now with ncview, it looks pretty messed up:

for c in (myocean_test)
return
encode (
c[t(29)], "netcdf"
)

[olcl@…]

The latest version of petascope from the new EOX repository fixes this bug and NetCDF output is correctly orientated.

For reference the version with the fix is :

Installed Packages
rasdaman.x86_64 8.3.1-3.el6 @eox-testing
rasdaman-debuginfo.x86_64 8.3.1-3.el6 @eox-testing
rasdaman-devel.x86_64 8.3.1-3.el6 @eox-testing
rasdaman-docs.noarch 8.3.1-3.el6 @eox-testing-noarch
rasdaman-examples.noarch 8.3.1-3.el6 @eox-testing-noarch
rasdaman-petascope.noarch 8.3.1-3.el6 @eox-testing-noarch
rasdaman-rasdaview.x86_64 8.3.1-3.el6 @eox-testing
rasdaman-rasgeo.x86_64 8.3.1-3.el6 @eox-testing

[Dimitar Misev]

In the latest petascope it's fixed because we are now using the rasql encode function; the netcdf function still needs to be fixed.

[abeccati]

To reporter, please verify and let us know if this can be marked as resolved, thanks.

[olcl@…]

Hi Alan,

as far as i am concerned this is fixed with the introduction of the Encode function. however Dimitar has commented above saying it still needs fixing in the netcdf function. so i will leave it up to him to decide if its closed.

[abeccati]

Thank you Olly.

Dimitar, please check whether to close, reassign or open another ticket. I've lowered priority since problem has been solved for the user anyway.

[Dimitar Misev]

It's not fixed yet, can be reproduced by

rasql -q 'select netcdf(c[29,*:*,*:*]) from myocean_test as c'

[mrusu]

Most likely caused by the precedence of latitude over longitude in netcdf dimensions. By the way, even if ncview shows a rotated image, the values on the axes are correct, the lower left corner being denoted as (0, 0) and the higher right as (max_x, max_y). As far as I know, netcdf() will be deprecated in favor of encode(), also on the server side.

[Dimitar Misev]

Oki, let's close.