> I've also tried using the proj4 strings offered
> for that EPSG, but the results are the same.
Hmm. It it looks to me that this should work:
gdal_translate -a_srs "+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1
+x_0=0 +y_0=0 +a=6378273 +b=6356889.449 +units=m +no_defs" -a_ullr ...
See [1] for gdalinfo's output from the resulting GeoTIFF.
Your case seems to be related to the recently opened ticket
http://trac.osgeo.org/gdal/ticket/3220
Hermann
[1]
$ gdalinfo asi-n6250-20020619-v5i_proj4.tif
Driver: GTiff/GeoTIFF
Files: asi-n6250-20020619-v5i_proj4.tif
Size is 1216, 1792
Coordinate System is:
PROJCS["unnamed",
GEOGCS["GCS Name = unnamed ellipse|Datum = unknown|Ellipsoid =
unnamed|Primem = Greenwich|",
DATUM["unknown",
SPHEROID["unnamed",6378273,298.279411123064]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Polar_Stereographic"],
PARAMETER["latitude_of_origin",70],
PARAMETER["central_meridian",-45],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-3800000.000000000000000,5600000.000000000000000)
Pixel Size = (6250.000000000000000,-6250.000000000000000)
Metadata:
(...)
Seb wrote:
Hi,
I originally posted to the GRASS list, but it was suggested that it
would be more appropriate here, given that GRASS uses GDAL for this
particular problem.
I'm in the process of importing sea ice concentration data from:
ftp://ftp-projects.zmaw.de/seaice/AMSR-E_ASI_IceConc
I have chosen to use their HDF product because their Geotiff uses a
WGS84 datum that results in up to 500 m error at the borders, according
to their info files at the root of each subdirectory.
The grids are not georeferenced, but they use the same grid that NSIDC
does with the specifications given here:
http://nsidc.org/data/docs/daac/nsidc0002_ssmi_seaice.gd.html
Fortunately, there's an EPSG code (3411, for the northern hemisphere)
for that projection.
gdalinfo on the HDF file:
---<--------------------cut here---------------start------------------->---
$ gdalinfo asi-n6250-20020619-v5i.hdf
Driver: HDF4Image/HDF4 Dataset
Files: asi-n6250-20020619-v5i.hdf
Size is 1216, 1792
Coordinate System is `'
Metadata:
valid_range=0, 1
long_name=ASI Ice Concentration, 20020619, res: 6.25000, AMSR-E, ASI Version:
5.5i, missing data interpolated.
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 1792.0)
Upper Right ( 1216.0, 0.0)
Lower Right ( 1216.0, 1792.0)
Center ( 608.0, 896.0)
Band 1 Block=1216x82 Type=Float32, ColorInterp=Gray
---<--------------------cut here---------------end--------------------->---
Converting to GTiff for importing into GRASS, using the EPSG code and
bounds in the NSIDC's info:
---<--------------------cut here---------------start------------------->---
$ gdal_translate -of GTiff -a_srs EPSG:3411 -a_ullr -3800000 5600000 3800000 -5600000 asi-n6250-20020619-v5i.hdf asi-n6250-20020619-v5i.tif
Warning 6: A dataset opened by GDALOpenShared should have the same filename (asi-n6250-20020619-v5i.hdf) and description (HDF4_SDS:UNKNOWN:"asi-n6250-20020619-v5i.hdf":0)
Input file size is 1216, 1792
0...10...20...30...40...50...60...70...80...90...100 - done.
$ gdalinfo asi-n6250-20020619-v5i.tif
Driver: GTiff/GeoTIFF
Files: asi-n6250-20020619-v5i.tif
Size is 1216, 1792
Coordinate System is:
PROJCS["NSIDC Sea Ice Polar Stereographic North",
GEOGCS["Unspecified datum based upon the Hughes 1980 ellipsoid",
DATUM["Not_specified_based_on_Hughes_1980_ellipsoid",
SPHEROID["Hughes 1980",6378273,298.279411123061,
AUTHORITY["EPSG","7058"]],
AUTHORITY["EPSG","6054"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4054"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3411"]]
Origin = (-3800000.000000000000000,5600000.000000000000000)
Pixel Size = (6250.000000000000000,-6250.000000000000000)
Metadata:
AREA_OR_POINT=Area
valid_range=0, 1
long_name=ASI Ice Concentration, 20020619, res: 6.25000, AMSR-E, ASI Version:
5.5i, missing data interpolated.
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-3800000.000, 5600000.000)
Lower Left (-3800000.000,-5600000.000)
Upper Right ( 3800000.000, 5600000.000)
Lower Right ( 3800000.000,-5600000.000)
Center ( 0.0000000, 0.0000000)
Band 1 Block=1216x1 Type=Float32, ColorInterp=Gray
---<--------------------cut here---------------end--------------------->---
So far so good, but when importing into a GRASS location:
---<--------------------cut here---------------start------------------->---
r.in.gdal in=/home/sluque/tmp/asi-n6250-20020619-v5i.tif out=test location=Test
WARNING: No projection name! Projection parameters likely to be
meaningless.
WARNING: Datum <Not_specified_based_on_Hughes_1980_ellipsoid> not
recognised by GRASS and no parameters found
Location <Test> created
100%
r.in.gdal complete. Raster map <test> created.
---<--------------------cut here---------------end--------------------->---
What is missing here? I've also tried using the proj4 strings offered
for that EPSG, but the results are the same. Any tips would be appreciated.
Cheers,
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
