Opened 13 years ago

Closed 13 years ago

Last modified 13 years ago

#87 closed defect (fixed)

Internal error: The index found more cells than allowed

Reported by: Alan Beccati <beccati@…> Owned by: Dimitar Misev
Priority: major Milestone:
Component: rasgeo Version:
Keywords: Cc: HerzigA@…
Complexity: Medium

Description

The use of rasgeo (rasimport) to load overlapping data into a 3d collection breaks select operation on the collection.

To reproduce install rasdaman, rasgeo and petascope
unpack the attached archive and run script.sh into it to load the images into the created "lol" collection at given raster coordinates

execute the following select to raise the error:

rasql -q "select t from lol as t" —out string
rasql: rasdaman query tool v1.0, rasdaman v8 — generated on 23.11.2011 14:34:26.
opening database RASBASE at localhost:7001...ok
Executing retrieval query…rasdaman error 2016: Internal error: The index found more cells than allowed.
aborting transaction…ok
rasql done.

Attachments (2)

utmDataRasimport.tgz (177.4 KB ) - added by Alan Beccati <beccati@…> 13 years ago.
Test dataset to reproduce the error
utmDataRasimport_alex.tgz (176.9 KB ) - added by herziga@… 13 years ago.
updated script to demonstrate usage of rasgeo for Alan's use case

Download all attachments as: .zip

Change History (6)

by Alan Beccati <beccati@…>, 13 years ago

Attachment: utmDataRasimport.tgz added

Test dataset to reproduce the error

comment:1 by Dimitar Misev, 13 years ago

Cc: HerzigA@… added

I'm not sure about this, is rasgeo meant for importing 3D data?

by herziga@…, 13 years ago

Attachment: utmDataRasimport_alex.tgz added

updated script to demonstrate usage of rasgeo for Alan's use case

in reply to:  1 comment:2 by herziga@…, 13 years ago

Replying to dmisev:

I'm not sure about this, is rasgeo meant for importing 3D data?

rasgeo supports creating and importing/updating of 2D x/y slices into/in 3D cubes. When ever rasgeo performs an update operation (i.e. -oid is specified) it automatically imports the named file into its correct 2D geospatial location as provided by its geospatial meta data info. In order to calculate the appropriate shift parameters, rasgeo needs to know about the geospatial location of the image meant to be updated. It achieves that by storing/updating meta data in the nm_meta table within the rasdaman data base after each rasgeo operation. Since in Alan's script the collection/image to be updated is not created by rasgeo, no meaningful meta data is available to calculate the proper shift vector for the file to be imported. In Alan's case that results in this spatial domain [-2147483648:1,-2147483648:1,0:10]. I'm not familiar with rasdaman internals, but obviously this causes the rasql select query to fail.

Using exclusively rasgeo to create and update the collection (see script in updated attachment) results in this spatial domain [-443:116,0:3740,10:10], and the rasql select query is executed successfully.

I'll update the README to emphasise that rasgeo depends on the geospatial meta data in nm_meta.

Alex

comment:3 by Dimitar Misev, 13 years ago

Resolution: fixed
Status: newclosed

Thanks Alex! I'll close the ticket as fixed then.

comment:4 by Dimitar Misev, 13 years ago

Component: GDAL driverrasgeo
Note: See TracTickets for help on using tickets.