10

USGS provides different resolution orthoimage quadrangles in the WorldImage format (*.tif with an accompanying *.tfw world file). Some of these quadrangles can be used by GeoServer (2.1.0 - 2.1.2) by creating a WorldImage store and publishing the layer (some of them throw an exception when trying to publish the layer). When publishing a GeoTiff layer that has been converted from a currently working WorldImage layer, I am getting the following exception:

2012-01-10 08:38:33,228 WARN [referencing.factory] - Axis elements found in a wkt definition, the force longitude first axis order hint might not be respected:
PROJCS["WGS84 / Simple Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295]], PROJECTION["Mercator_1SP_Google"], PARAMETER["latitude_of_origin", 0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH], AUTHORITY["EPSG","54004"]]
2012-01-10 08:38:39,674 INFO [geoserver.web] - Error saving layer
java.lang.NullPointerException
    at org.geoserver.catalog.impl.CoverageStoreInfoImpl.getFormat(CoverageStoreInfoImpl.java:42)
...

I tried using gdal_translate in the following way (substitute 1234 with real values): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tif The conversion succeeds, and I can successfully create a GeoTiff store in GeoServer, but cannot publish the layer because of the error described above. This leads me to believe that I am trying to publish layers (in both cases) that are not properly formatted. Thus, what is the correct way to convert a [working] USGS orthoimage in WorldImage format to a tiled GeoTiff for GeoServer 2.1.2?

Improve this question
9
  • I think we need to know which epsg code you are actually using as it is probably the important bit
    – Ian Turton ♦
    Jan 10, 2012 at 14:35
  • Edited to show an EPSG I've tried with a 2007 data set retrieved from EarthExplorer.
    – Matt
    Jan 10, 2012 at 14:37
  • Do you know what the input tiff's SRS is?
    – Ian Turton ♦
    Jan 10, 2012 at 14:52
  • 1
    It looks like you've not got to the root cause because Parameter message can't be null is an Apache error thrown when a FeedbackMessage doesn't have a message string (see here: bit.ly/zBUjy3 for the gory details). I suggest you post the whole stack trace so we might be able to identify what message GeoServer is trying to tell us.
    – MerseyViking
    Jan 10, 2012 at 15:51
  • 3
    gdal_translate -a_srs doesn't reproject, it merely assigns (or overwrites) the image's projection. If the image has no bounding information, it'll assume 1 US foot (for EPSG:2239) per pixel and place the image at the coordinate system's origin. If the bounding box is in a unit other than US feet, then your image could be anywhere or at any size depending on the units and origin of whatever projection it was digitised in.
    – MerseyViking
    Jan 10, 2012 at 16:14

2 Answers 2

Reset to default
2

If you're getting geopdf's from USGS and you want to convert to geotiff - and you have a Linux box with perl or mod_perl capability and also GDAL/OGR libraries - you could batch run them with geopdf2gtiff.pl script as shown here: https://github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl

Improve this answer
1

Why not just download in geotiff? USGS used to have that option. The site USGS list for this format are The National Map and US Topo.

Improve this answer

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.