[gdal-dev] Unknown file type
Greg Ederer
greg at ergonosis.com
Fri Jan 9 13:22:22 EST 2009
Hi Frank,
Thanks you very much for your help. I think I'm close. I have
created a .vrt file with:
<VRTDataset rasterXSize="751" rasterYSize="801">
<SRS>
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
</SRS>
<VRTRasterBand dataType="Float32" band="1"
subClass="VRTRawRasterBand">
<SourceFilename relativetoVRT="1">10day_precip.bin.1999121</
SourceFilename>
<ImageOffset>0</ImageOffset>
<PixelOffset>4</PixelOffset>
<LineOffset>3004</LineOffset>
<ByteOrder>LSB</ByteOrder>
<NoDataValue>-999.0</NoDataValue>
</VRTRasterBand>
</VRTDataset>
I then run:
gdal_translate -of GTiff 10day_precip.bin.1999121.vrt
10day_precip.bin.1999121.tif
and this produces a GeoTIFF, as expected. Running:
gdalinfo 10day_precip.bin.1999121.tif
gives:
Driver: GTiff/GeoTIFF
Files: 10day_precip.bin.1999121.tif
Size is 751, 801
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 801.0)
Upper Right ( 751.0, 0.0)
Lower Right ( 751.0, 801.0)
Center ( 375.5, 400.5)
Band 1 Block=751x2 Type=Float32, ColorInterp=Gray
NoData Value=-999
However, when I attempt to load this GeoTIFF into GeoServer, my
ultimate goal, I get:
class org.vfny.geoserver.global.ConfigurationException: Could not
obtain a reader for the CoverageDataSet. Please check the
CoverageDataSet configuration!
I have successfully loaded other GeoTIFFs into GeoServer that were
prepared by another member of our team from the same raw binary data
using IDL So, I suspect that there is something amiss with my GeoTIFF.
Here's the gdalinfo output for a GeoTIFF that loads properly:
Driver: GTiff/GeoTIFF
Files: /Users/gregederer/tmp/starspan_tests/rfe_2005_07_mm.tif
Size is 751, 801
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-20.050000000000001,40.049999999999997)
Pixel Size = (0.100000000000000,-0.100000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -20.0500000, 40.0500000) ( 20d 3'0.00"W, 40d 3'0.00"N)
Lower Left ( -20.0500000, -40.0500000) ( 20d 3'0.00"W, 40d 3'0.00"S)
Upper Right ( 55.0500000, 40.0500000) ( 55d 3'0.00"E, 40d 3'0.00"N)
Lower Right ( 55.0500000, -40.0500000) ( 55d 3'0.00"E, 40d 3'0.00"S)
Center ( 17.5000000, -0.0000000) ( 17d30'0.00"E, 0d 0'0.00"S)
Band 1 Block=751x2 Type=Float32, ColorInterp=Gray
Any further suggestions?
Thanks, again!
Greg
On Jan 8, 2009, at 11:56 AM, Frank Warmerdam wrote:
> Greg Ederer wrote:
>> Hi,
>> I am a GDAL newbie. I just got handed a group of files with the
>> following description:
>> description of the data
>> 1) content:
>> rain == 10 day precipitation analysis by
>> merging GTS gauge observations and 3 kinds
>> of satellite estimates
>> (GPI,SSM/I, and AMSU)
>> 2) coverage:
>> -40.00S -- 40.00N Northward;
>> 20.00W -- 55.00E Eastward 751 grid points in east -
>> west direction
>> 801 grid points in south - north direction
>> 3) resolution:
>> 0.1 deg lat x 0.1 deg lon
>> dekad
>> 4) missing values
>> -999.0
>> I'd like to read and manipulate these files using GDAL. When I run
>> gdalinfo, I get something like the following:
>> Macintosh-4:RFE_dekedal gregederer$ gdalinfo 10day_precip.bin.1999121
>> ERROR 4: `10day_precip.bin.1999121' not recognised as a supported
>> file format.
>> gdalinfo failed - unable to open '10day_precip.bin.1999121'.
>> I suspect that GDAL requires some additional metadata in order to
>> read these files. Any suggestions?
>
> Greg,
>
> Assuming this is a raw binary file you can write a VRT "raw"
> description
> for it based on:
>
>
> http://www.gdal.org/gdal_vrttut.html
>
> Read down in particular to the ".vrt Descriptions for Raw Files"
> section.
>
> There are also several other methods for accomplishing something
> similar.
>
> In your case I'd assume the file is 751x801 x 1band x Float32 with
> ImageOffset=0, PixelOffset=4, LineOffset = 3004 and ByteOrder LSB.
>
> Best regards,
> --
> ---------------------------------------
> +--------------------------------------
> I set the clouds in motion - turn up | Frank Warmerdam, warmerdam at pobox.com
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush | Geospatial Programmer for
> Rent
>
More information about the gdal-dev
mailing list