[GIS] Problem in importing raster data from PostGIS to Geoserver

geoserverpostgisraster

I am using Opengeo Suite 3.0.2, PostGIS 2.0.3 and PostgreSQl 9.2.

I am trying to import raster data from PostGIS to Geoserver but facing an error.

I follow this tutorial link to import the raster data from PostGIS.

http://dev.horizon.opengeo.org/opengeo-docs/geoserver/tutorials/imagemosaic-jdbc/imagemosaic-jdbc_tutorial.html

All steps are going good.

I used FWtools 2.2.4 to make tiles.
gdal_retile works fine.

These are my XML configuration files…

rasterdata.pgraster.xml

<?xml version="1.0" encoding="UTF-8" standalone="no"?> 
<!DOCTYPE ImageMosaicJDBCConfig [ <!ENTITY mapping PUBLIC "mapping" "mapping.pgraster.xml.inc"> 
<!ENTITY connect PUBLIC "connect" "connect.pgraster.xml.inc"> ]> 
<config version="1.0"> 
<coverageName name="slp"/> 
<coordsys name="EPSG:4326"/> 
<scaleop interpolation="1"/> &mapping; &connect; 
</config>

connect.pgraster.xml.inc

<connect>
<dstype value="DBCP"/> 
<username value="postgres" /> 
<password value="postgre" /> 
<jdbcUrl value="jdbc:postgresql://localhost:5432/spatial_db" /> 
<driverClassName value="org.postgresql.Driver"/> 
<maxActive value="10"/> 
<maxIdle value="0"/> 
</connect>

mapping.pgraster.xml.inc

<spatialExtension name="postgis"/>
<mapping> 
<masterTable name="rasterdata_master2" > 
  <coverageNameAttribute name="slp"/> 
  <maxXAttribute name="maxx"/> 
  <maxYAttribute name="maxy"/> 
  <minXAttribute name="minx"/> 
  <minYAttribute name="miny"/> 
  <resXAttribute name="resx"/> 
  <resYAttribute name="resy"/> 
  <tileTableNameAtribute name="TileTable" /> 
  <spatialTableNameAtribute name="SpatialTable" />
</masterTable> 
<tileTable> 
  <blobAttributeName name="rast" /> 
  <keyAttributeName name="rid" /> 
</tileTable> 
<spatialTable>
  <keyAttributeName name="rid" />
  <geomAttributeName name="geom" />
  <tileMaxXAttribute name="maxX"/>
  <tileMaxYAttribute name="maxY"/>
  <tileMinXAttribute name="minX"/>
  <tileMinYAttribute name="minY"/>
</spatialTable>
</mapping>

After that I used the DDL generation utility which creates 4 Sql Scripts.

I execute the required script which creates the metadata tables.

now I used java utility to import the image data. That also works fine.

But when I try to use the image mosaic jdbc plugin to import the raster data through jdbc database, it gives an error, i.e. "Could not list layers for this store, an error occurred retrieving them: Argument "value" should not be null."

I am not able to insert the screenshot right now, so can't show you the screenshot of error.

Thanks

Akhil

Best Answer

I recommend the same as in How to publish PostGIS Rasters in Geoserver?: Following the configuration instructions for the PostGIS raster plugin, you should try replacing

<spatialExtension name="postgis"/>

with

<spatialExtension name="pgraster"/>

in configration file mapping.postgis.xml.inc.

Also, maybe <password value="postgre" /> in connect.pgraster.xml.inc is a typo for <password value="postgres" />.

Related Solutions

Implementing Search Functionality in OpenLayers and GeoServer

You can pass query terms in the URL to an SQL View in GeoServer 2.1 and above. Create the view with an ilike clause and you should get the basic functionality. Have OpenLayers fetch the view with the supplied name as a parameter. Then, tell OpenLayers zoom to the extent of the result.

If you don't have GeoServer 2.1:

Get the name of the feature with Javascript/OpenLayers,
pass the name asynchronously to a php script on the server,
the php script will query the database directly for the st_envelope or st_box2d of the feature with that name (select st_envelope(the_geom) from table where name ilike passed_parameter),
the php script returns the extent/bounding box of the feature which is used in Javascript/OpenLayers to zoom the map.

You'll have to parse the st_envelope result into something that OpenLayers understands. You also might want to add a little buffer to the zoom extent. You can add this manually or use st_expand instead of st_envelope.

[GIS] Publishing PostGIS Rasters in GeoServer

I worked through several issues trying to use a postgis raster in geoserver similar to the OP. The below is how I successfully loaded my raster data.

Preconfiguration: A BNG geotiff layer, generated from FME, has been imported into postgis using raster2pgsql. The command used was:

raster2pgsql -s 27700 -I -M -C {geotiffName}.tif -t 100x100 -l 2,4,6 {schema}.{tableName} | psql -U {user} -d {database}

This created the raster table along with its pyramid tables.

Currently my server has an instance of Geoserver 2.6.2. To use Postgis rasters it requires the extension JDBC Image Mosaic which can be downloaded from http://geoserver.org/release/2.6.2/

The extension .jar was then extracted into WEB-INF/lib. This will then enable geoserver to use rasters from postgis. However to access the data several settings need configuring.

Following the instructions in: http://docs.geotools.org/latest/userguide/library/coverage/pgraster.html but slightly tweaked:

create table {schema}.MOSAIC (NAME varchar(254) not null,
     TileTable varchar(254)not null,
     minX FLOAT8,minY FLOAT8, maxX FLOAT8, maxY FLOAT8,resX FLOAT8, resY FLOAT8, primary key (NAME,TileTable))

Now insert records for a "coverage", named oek in the example, these are your raster tables: the master one and pyramid ones.

insert into {schema}.MOSAIC (name, tiletable) values ('oek', {Main raster table name});
insert into {schema}.MOSAIC (name, tiletable) values ('oek', {raster pyramid 2 table name});
insert into {schema}.MOSAIC (name, tiletable) values ('oek', {raster pyramid 4 table name});
insert into {schema}.MOSAIC (name, tiletable) values ('oek', {raster pyramid 6 table name});

Also according to the docs, for Postgis JDBC drivers < version 9 it is necessary to set a property in the database:

ALTER DATABASE dbname SET bytea_output TO 'escape'

Now Geoserver configuration:

Within your Geoserver data directory create a directory /data/coverages, this needs to contain 3 files: connect.pgraster.xml.inc, mapping.pgraster.xml.inc, oek.pgraster.xml

Follow http://docs.geotools.org/latest/userguide/library/coverage/pgraster.html for their naming and contents. Obviously adjust the database connection string and username/password settings in connect.pgraster.xml.inc

mapping.pgraster.xml.inc is where we need to reference the raster tables correctly.

<spatialExtension name="pgraster"/>
 <mapping>
  <masterTable name="MOSAIC" >
          <coverageNameAttribute name="name"/>
          <maxXAttribute name="maxX"/>
          <maxYAttribute name="maxY"/>
          <minXAttribute name="minX"/>
          <minYAttribute name="minY"/>
          <resXAttribute name="resX"/>
          <resYAttribute name="resY"/>
          <tileTableNameAtribute  name="tiletable" />
  </masterTable>
  <tileTable>
          <blobAttributeName name="rast" />
  </tileTable>
</mapping>

So this will get the values we added to the MOSAIC table, the name values are the columns of MOSAIC, which link to our raster table and pyramid tables. rast is the name of the raster column in the raster tables.

In oek.pgraster.xml make sure your SRS is correctly set, 27700 BNG in this case.

<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE ImageMosaicJDBCConfig [
  <!ENTITY mapping PUBLIC "mapping"  "mapping.pgraster.xml.inc">
  <!ENTITY connect PUBLIC "connect"  "connect.pgraster.xml.inc">
]>

<config version="1.0">
  <coverageName name="oek"/>
  <coordsys name="EPSG:27700"/>
  <!-- interpolation 1 = nearest neighbour, 2 = bipolar, 3 = bicubic -->
  <scaleop  interpolation="1"/>
<axisOrder ignore="false"/>
  &mapping;
  &connect;
</config

In Geoserver we can now add a new datastore.

Go to new Add new store.

Select ImageMosaicJDBC

Give it a name and the set the URL to file:coverages/oek.pgraster.xml

Save this, and fingers crossed no errors are returned. If they are, check geoserver logs which hopefully will indicate what the issue was, hopefully just a typo or naming issue.

Go to add a new layer and select the new raster data store.

Select your new layer/coverage named oek and publish it.

Then in the layer settings adjust your coordinate system and bounding boxes appropriately.

To style the raster follow docs here: http://docs.geoserver.org/latest/en/user/styling/sld-reference/rastersymbolizer.html

I added a style:

<RasterSymbolizer>
 <Opacity>1.0</Opacity>
  <ColorMap>
   <ColorMapEntry color="#ffffff" quantity="1" label="label1" opacity="1"/>
   <ColorMapEntry color="#0000ff" quantity="200" label="label2" opacity="1"/>
  </ColorMap>
</RasterSymbolizer>

This colored my rasters blue where they had a value.

You should then have a usable layer.

Disclaimer: I am no GIS expert, and a lot of what i worked through was trial and error. I have no doubt that there are better ways to do it, but for the time being I'm only working on a proof of concept to confirm if we have the capability to use raster data in postgis. Finer details like raster values and styling are things I still need to iron out. If I make any further progress I will update my answer.

Best Answer

Related Solutions

Implementing Search Functionality in OpenLayers and GeoServer

[GIS] Publishing PostGIS Rasters in GeoServer

Related Question