.. _raster.georaster:

================================================================================
Oracle Spatial GeoRaster
================================================================================

.. shortname:: GeoRaster

.. build_dependencies:: Oracle client libraries

This driver supports reading and writing raster data in Oracle Spatial
GeoRaster format (10g or later). The Oracle Spatial GeoRaster driver is
optionally built as a GDAL plugin, but it requires Oracle client
libraries.

When opening a GeoRaster, its name should be specified in the form:

| georaster:<user>{,/}<pwd>{,@}[db],[schema.][table],[column],[where]
| georaster:<user>{,/}<pwd>{,@}[db],<rdt>,<rid>

Where:

| user   = Oracle server user's name login
| pwd    = user password
| db     = Oracle server identification (database name)
| schema = name of a schema
| table  = name of a GeoRaster table (table that contains GeoRaster
  columns)
| column = name of a column data type MDSYS.SDO_GEORASTER
| where  = a simple where clause to identify one or multiples
  GeoRaster(s)
| rdt    = name of a raster data table
| rid    = numeric identification of one GeoRaster

Examples:

| geor:scott,tiger,demodb,table,column,id=1
| geor:scott,tiger,demodb,table,column,"id = 1"
| "georaster:scott/tiger@demodb,table,column,gain>10"
| "georaster:scott/tiger@demodb,table,column,city='Brasilia'"
| georaster:scott,tiger,,rdt_10$,10
| geor:scott/tiger,,rdt_10$,10

Note: do note use space around the field values and the commas.

Note: like in the last two examples, the database name field could be
left empty (",,") and the TNSNAME will be used.

Note: If  the query results in more than one GeoRaster it will be
treated as a GDAL metadata's list of sub-datasets (see below)

Driver capabilities
-------------------

.. supports_createcopy::

.. supports_create::

.. supports_georeferencing::

.. supports_virtualio::

Browsing the database for GeoRasters
------------------------------------

By providing some basic information the GeoRaster driver is capable of
listing the existing rasters stored on the server:

To list all the GeoRaster table on the server that belongs to that user
name and database:

.. code-block:: bash

   gdalinfo georaster:scott/tiger@db1

To list all the GeoRaster type columns that exist in that table:

.. code-block:: bash

   gdalinfo georaster:scott/tiger@db1,table_name

That will list all the GeoRaster objects stored in that table.

.. code-block:: bash

   gdalinfo georaster:scott/tiger@db1,table_name,georaster_column

That will list all the GeoRaster existing on that table according to a
Where clause.

.. code-block:: bash

   gdalinfo georaster:scott/tiger@db1,table_name,georaster_column,city='Brasilia'


Note that the result of those queries are returned as GDAL metadata
sub-datasets, e.g.:

.. code-block:: bash

  gdalinfo georaster:scott/tiger
  # Driver: GeoRaster/Oracle Spatial GeoRaster
  # Subdatasets:
  # SUBDATASET_1_NAME=georaster:scott,tiger,,LANDSAT
  # SUBDATASET_1_DESC=Table:LANDSAT
  # SUBDATASET_2_NAME=georaster:scott,tiger,,GDAL_IMPORT
  # SUBDATASET_2_DESC=Table:GDAL_IMPORT

Creation Options
----------------

-  .. co:: BLOCKXSIZE

      The number of pixel columns per raster block.

-  .. co:: BLOCKYSIZE

      The number of pixel rows per raster block.

-  .. co:: BLOCKBSIZE

      The number of bands per raster block.

-  .. co:: BLOCKING

      Decline the use of blocking (NO) or request an
      automatic blocking size (OPTIMUM).

-  .. co:: SRID

      Assign a specific EPSG projection/reference system
      identification to the GeoRaster.

-  .. co:: INTERLEAVE
      :choices: BAND, LINE, PIXEL

      Band interleaving mode, BAND, LINE, PIXEL (or BSQ,
      BIL, BIP) for band sequential, Line or Pixel interleaving.
      Starting with GDAL 3.5, when copying from a source dataset with multiple bands
      which advertises a INTERLEAVE metadata item, if the INTERLEAVE creation option
      is not specified, the source dataset INTERLEAVE will be automatically taken
      into account, unless the :co:`COMPRESS` creation option is specified.

-  .. co:: DESCRIPTION

      A simple description of a newly created table in SQL
      syntax. If the table already exist, this create option will be
      ignored, e.g.:

      .. code-block:: bash

         gdal_translate -of georaster landsat_823.tif geor:scott/tiger@orcl,landsat,raster \
           -co DESCRIPTION="(ID NUMBER, NAME VARCHAR2(40), RASTER MDSYS.SDO_GEORASTER)" \
           -co INSERT="VALUES (1,'Scene 823',SDO_GEOR.INIT())"

-  .. co:: INSERT

      A simple SQL insert/values clause to inform the driver
      what values to fill up when inserting a new row on the table, e.g.:

      .. code-block:: bash

         gdal_translate -of georaster landsat_825.tif geor:scott/tiger@orcl,landsat,raster \
           -co INSERT="(ID, RASTER) VALUES (2,SDO_GEOR.INIT())"

-  .. co:: COMPRESS
      :choices: JPEG-F, JP2-F, DEFLATE, NONE

      Compression options.
      The JPEG-F options is lossy, meaning that the original pixel values
      are changed. The JP2-F compression is lossless if :co:`JP2_QUALITY=100`.

-  .. co:: GENPYRAMID

      Generate pyramid after a GeoRaster object have been
      loaded to the database. The content of that parameter must be the
      resampling method of choice NN (nearest neighbor) , BILINEAR,
      BIQUADRATIC, CUBIC, AVERAGE4 or AVERAGE16. If :co:`GENPYRLEVELS` is not
      informed the PL/SQL function sdo_geor.generatePyramid will calculate
      the number of levels to generate.

-  .. co:: GENPYRLEVELS

      Define the number of pyramid levels to be
      generated. If :co:`GENPYRAMID` is not informed the resample method NN
      (nearest neighbor) will apply.

-  .. co:: GENSTATS
      :choices: TRUE, FALSE
      :default: FALSE

      To generate statistics from the given rasters, set this value to TRUE.
      Default value is FALSE.
      This option must be present in order to generate the stats, otherwise,
      other GENSTATS options are ignored.


-  .. co:: GENSTATS_SAMPLINGFACTOR
      :choices: <integer>
      :default: 1

      Defines the number of cells skipped in both row and column dimensions when 
      the statistics are computed.
      For example, when setting this value to 4, one-sixteenth of the cells are sampled.
      The higher the value, the less accurate the statistics are likely to be, 
      but the more quickly they will be computed.
      Defaults to 1, which means that all cells are sampled.

-  .. co:: GENSTATS_SAMPLINGWINDOW
      :choices: <integer\,integer\,integer\,integer>

      This parameter identifies the upper-left (row, column) and lower-right 
      (row, column) coordinates of a rectangular window, and raster space is assumed.
      The intersection of the MBR (minimum bounding rectangle) of the
      GeoRaster object in raster space is used for computing statistics.
      When this value is not specified, statistics are computed for the entire raster.

-  .. co:: GENSTATS_HISTOGRAM
      :choices: TRUE, FALSE
      :default: FALSE

      When this value is set to TRUE, a histogram will be computed and stored.
      Defaults to FALSE, so a histogram won't be generated.

-  .. co:: GENSTATS_LAYERNUMBERS
      :choices: <integer\,integer\,...>, <integer>-<integer>

      Defines the numbers of the layers for which to compute the statistics.
      This can include numbers, ranges (indicated by hyphens), and commas 
      to separate any combination of those.
      For example, '1,3-5,7', '1,3,6', '1-6'.
      If this value is not specified, statistics will be computed for all
      layers.

-  .. co:: GENSTATS_USEBIN
      :choices: TRUE, FALSE
      :default: TRUE

      Defaults to TRUE.
      Specifies whether or not to use a provided bin function 
      (specified in :co:`GENSTATS_BINFUNCTION`).
      When this value is set to TRUE, the bin function to be
      used follows the following sequence: (1) the bin function 
      specified in :co:`GENSTATS_BINFUNCTION`. (2) the bin
      function specified by the <binFunction> element in the
      GeoRaster XML metadata. (3) a dynamically generated bin
      function generated as follows: 
      Min and max are the actual min and max values of the raster
      Numbins is defined as:
      * celldepth = 1, numbins = 2.
      * cellDepth = 2, numbins = 4.
      * cellDepth = 4, numbins = 8.
      * cellDepth >= 8, numbins = 256. 

      When this value is set to FALSE, then the bin function
      to use is the one generated dynamically as previously
      described.

-  .. co:: GENSTATS_BINFUNCTION
      :choices: <integer\,integer\,integer\,integer\,integer>
      
      An array whose element specify the bin type, total
      number of bins, first bin number, minimum cell value,
      and maximum cell value. Bin type must be linear (0).
      When this value is not specified, and :co:`GENSTATS_USEBIN`
      is TRUE, then the bin function to use is as follows: 

      1. A valid function defined in the GeoRaster metadata.
      2. The same bin function generated when :co:`GENSTATS_USEBIN` is FALSE.

-  .. co:: GENSTATS_NODATA
      :choices: TRUE, FALSE
      :default: FALSE

      Specifies whether or not to compare each cell values
      with NODATA values defined in the metadata when computing
      statistics. When set to TRUE, all pixels with a NODATA
      value won't be considered. When set to FALSE, NODATA
      pixels will be considered as regular pixels.

      A NODATA value is used for cells whose values are either not known
      or meaningless

-  .. co:: QUALITY
      :choices: 0-100
      :default: 75

      Quality compression option for JPEG ranging from 0 to 100.

-  .. co:: JP2_QUALITY
      :choices: <float_value\,float_value\,...>

      Only if :co:`COMPRESS=JP2-f`.
      Percentage between 0 and 100. A value of 50 means the file will be
      half-size in comparison to uncompressed data, 33 means 1/3, etc..
      Defaults to 25 (unless the dataset is made of a single band with
      color table, in which case the default quality is 100).

-  .. co:: JP2_REVERSIBLE
      :choices: YES, NO

      Only if :co:`COMPRESS=JP2-f`. YES means use of
      reversible 5x3 integer-only filter, NO use of the irreversible DWT
      9-7. Defaults to NO (unless the dataset is made of a single band with
      color table, in which case reversible filter is used).

-  .. co:: JP2_RESOLUTIONS
      :choices: <integer>

      Only if :co:`COMPRESS=JP2-f`. Number of
      resolution levels. Default value is selected such the smallest
      overview of a tile is no bigger than 128x128.

-  .. co:: JP2_BLOCKXSIZE
      :choices: <integer>
      :default: 1024

      Only if :co:`COMPRESS=JP2-f`. Tile width.

-  .. co:: JP2_BLOCKYSIZE
      :choices: <integer>
      :default: 1024

      Only if :co:`COMPRESS=JP2-f`. Tile height.

-  .. co:: JP2_PROGRESSION
      :choices: LRCP, RLCP, RPCL, PCRL, CPRL
      :default: LRCP

      Only if :co:`COMPRESS=JP2-f`.
      Progression order.

-  .. co:: NBITS
      :choices: 1, 2, 4

      Sub byte data type.

-  .. co:: SPATIALEXTENT
      :choices: TRUE, FALSE
      :default: TRUE

      Generate Spatial Extents. The default value
      is TRUE, which means that this option only need to be set
      to force the Spatial Extent to remain as NULL. If :co:`EXTENTSRID` is not
      set the Spatial Extent geometry will be generated with the same
      SRID as the GeoRaster object.

-  .. co:: EXTENTSRID

      SRID code to be used on the Spatial Extent geometry.
      If the table/column has already a spatial index, the value specified
      should be the same as the SRID on the Spatial Extents of the other
      existing GeoRaster objects, on which the spatial index is built.

-  .. co:: OBJECTTABLE
      :choices: TRUE, FALSE
      :default: FALSE

      To create RDT as SDO_RASTER object set to TRUE.
      Otherwise, the RDT will be created as
      regular relational tables. That does not apply for Oracle version
      older than 11.

Importing GeoRaster
-------------------

During the process of importing raster into a GeoRaster object it is
possible to give the driver a simple SQL table definition and also a SQL
insert/values clause to inform the driver about the table to be created
and the values to be added to the newly created row. The following
example does that:

.. code-block:: bash

    gdal_translate -of georaster Newpor.tif georaster:scott/tiger,,landsat,scene \
      -co "DESCRIPTION=(ID NUMBER, SITE VARCHAR2(45), SCENE MDSYS.SDO_GEORASTER)" \
      -co "INSERT=VALUES(1,'West fields', SDO_GEOR.INIT())" \
      -co "BLOCKXSIZE=512" -co "BLOCKYSIZE=512" -co "BLOCKBSIZE=3" \
      -co "INTERLEAVE=PIXEL" -co "COMPRESS=JPEG-F"

Note that the create option DESCRIPTION requires to inform table name
(in bold). And column name (underlined) should match the description:

.. code-block:: bash

    gdal_translate -of georaster landsat_1.tif georaster:scott/tiger,,landsat,scene \
      -co "DESCRIPTION=(ID NUMBER, SITE VARCHAR2(45), SCENE MDSYS.SDO_GEORASTER)" \
      -co "INSERT=VALUES(1,'West fields', SDO_GEOR.INIT())"

If the table "landsat" exist, the option "DESCRIPTION" is ignored. The
driver can only update one GeoRaster column per run of
gdal_translate. Oracle create default names and values for RDT and RID
during the initialization of the SDO_GEORASTER object but user are also
able to specify a name and value of their choice.

.. code-block:: bash

   gdal_translate -of georaster landsat_1.tif georaster:scott/tiger,,landsat,scene \
     -co "INSERT=VALUES(10,'Main building', SDO_GEOR.INIT('RDT', 10))"

If no information is given about where to store the raster the driver
will create (if doesn't exist already) a default table named GDAL_IMPORT
with just one GeoRaster column named RASTER and a table GDAL_RDT as the
RDT, the RID will be given automatically by the server, example:

.. code-block:: bash

   gdal_translate -of georaster input.tif “geor:scott/tiger@dbdemo”

Exporting GeoRaster
-------------------

A GeoRaster can be identified by a Where clause or by a pair of RDT & RID:

.. code-block:: bash

   gdal_translate -of gtiff geor:scott/tiger@dbdemo,landsat,scene,id=54 output.tif
   gdal_translate -of gtiff geor:scott/tiger@dbdemo,st_rdt_1,130 output.tif

Cross schema access
-------------------

As long as the user was granted full access the GeoRaster table and
the Raster Data Table, e.g.:

::

    % sqlplus scott/tiger
    SQL> grant select,insert,update,delete on gdal_import to spock;
    SQL> grant select,insert,update,delete on gdal_rdt to spock;

It is possible to an user access to extract and load GeoRaster from
another user/schema by informing the schema name as showed here:

Browsing:

.. code-block:: bash

   gdalinfo geor:spock/lion@orcl,scott.
   gdalinfo geor:spock/lion@orcl,scott.gdal_import,raster,"t.raster.rasterid > 100"
   gdalinfo geor:spock/lion@orcl,scott.gdal_import,raster,t.raster.rasterid=101

Extracting:

.. code-block:: bash

   gdal_translate geor:spock/lion@orcl,scott.gdal_import,raster,t.raster.rasterid=101out.tif
   gdal_translate geor:spock/lion@orcl,gdal_rdt,101 out.tif

Note: On the above example that accessing by RDT/RID doesn't need
schame name as long as the users is granted full access to both
tables.

Loading:

.. code-block:: bash

    gdal_translate -of georaster input.tif geor:spock/lion@orcl,scott.
    gdal_translate -of georaster input.tif geor:spock/lion@orcl,scott.cities,image \
      -co INSERT="(1,'Rio de Janeiro',sdo_geor.init('cities_rdt'))"

General use of GeoRaster
------------------------

GeoRaster can be used in any GDAL command line tool with all the available options.
Like a image subset extraction or re-project:

.. code-block:: bash

    % gdal_translate -of gtiff geor:scott/tiger@dbdemo,landsat,scene,id=54 output.tif \
    -srcwin 0 0 800 600
    % gdalwarp -of png geor:scott/tiger@dbdemo,st_rdt_1,130 output.png
   -t_srs EPSG:9000913

Two different GeoRaster can be used as input and output on the same operation:

.. code-block:: bash

    % gdal_translate -of georaster geor:scott/tiger@dbdemo,landsat,scene,id=54 \
    geor:scott/tiger@proj1,projview,image -co INSERT="VALUES (102, SDO_GEOR.INIT())"

Applications that use GDAL can theoretically read and write from GeoRaster just like
any other format but most of then are more inclined to try to access files on the file
system so one alternative is to create VRT to represent the GeoRaster description, e.g.:

.. code-block:: bash

    % gdal_translate -of VRT geor:scott/tiger@dbdemo,landsat,scene,id=54 view_54.vrt
    % openenv view_54.vrt
