[gdal-dev] HDF data mis-applied sinusoidal

[Michael Sumner]

Hello, these data in HDF4 format are interpreted as being in sinusoidal
projection, I guess some conflation with the MODIS logic about some
products (?).

ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/

Particular file is here, but note there's also an auxiliarly .hdf.xml with
it:
ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/GLASS02B06.V04.A2013169.2017128.hdf

A plot of the first subdataset shows that it is trivially in regular
longlat (or eqc), the gdalinfo output is below.

I'd appreciate any insights into why this is interpreted this way, and
whether the GDAL detection logic could be improved - or whether the data
providers need to update the internal metadata.

Cheers, Mike.

gdalinfo --version GDAL 2.3.2, released 2018/09/21

gdalinfo
HDF4_EOS:EOS_GRID:"GLASS02B06.V04.A2013169.2017128.hdf":GLASS02B06:ABD_BSA_VIS
Driver: HDF4Image/HDF4 Dataset
Files: GLASS02B06.V04.A2013169.2017128.hdf
Size is 7200, 3600
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Unknown datum based upon the custom spheroid",
        DATUM["Not specified (based on custom spheroid)",
            SPHEROID["Custom spheroid",6371007.181,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Sinusoidal"],
    PARAMETER["longitude_of_center",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
Origin = (-20015109.353999998420477,10007554.676999999210238)
Pixel Size = (154.437572175972150,-308.875144351944300)
Metadata:
  add_offset=0
  add_offset_err=0
  ASSOCIATEDINSTRUMENTSHORTNAME.1=AVHRR
  ASSOCIATEDPLATFORMSHORTNAME.1=NOAA
  ASSOCIATEDSENSORSHORTNAME.1=AVHRR
  BANDDEFINITION=Surface
  CHARACTERISTICBINSIZE=926.6
  DATACOLUMNS=7200
  DATAROWS=3600
  DATASETNAME=GLASS02B06
  DAYNIGHTFLAG=Both
  EASTHBOUNDINGCOORDINATE=180
  EXCLUSIONGRINGFLAG=N
  GLOBALGRIDCOLUMNS=7200
  GLOBALGRIDROWS=3600
  GRIDTYPE=Geographic Lat/Lon
  GRINGPOINTLATITUDE.1=90,90,-90,-90
  GRINGPOINTLONGITUDE.1=-180,180,180,-180
  GRINGPOINTSEQUENCENO.1=1,2,3,4
  HDFEOSVersion=HDFEOS_V2.16
  HORIZONTALTILENUMBER=00
  HORIZONTALTILENUMBER=00
  INPUTFILESNAME=Multi-Source Data
  INPUTPOINTER=Multi-Source Data
  INSTITUTENAME=BEIJING
  LOCALGRANULEID=GLASS02B06.V04.A2013169.2017128.hdf
  LOCALVERSIONID=V04
  long_name=Black-sky Albedo in Visible band
  NORTHBOUNDINGCOORDINATE=90
  PARAMETERNAME.1=GLASS02B06
  PGEVERSION=V04
  PROCESSINGENVIRONMENT=High Performance Computing System
  PRODUCTDATEANDTIME=2017-05-08 16:00:12
  PRODUCTIONDATETIME=2017-05-08 16:00:12
  PRODUCTNAME=Surface
  PRODUCTQUALITY=CLOUDREMOVED
  REPROCESSINGACTUAL=reprocessed
  REPROCESSINGPLANNED=further update is anticipated
  scale_factor=0.0001
  scale_factor_err=0
  SHORTNAME=GLASS02B06
  SOUTHBOUNDINGCOORDINATE=-90
  TILENUMBER=H00V00
  units=Albedo,no units
  valid_range=0, 10000
  VERSIONID=1
  VERTICALTILENUMBER=00
  VERTICALTILENUMBER=00
  WEBSITE=http://glass-product.bnu.edu.cn
  WESTHBOUNDINGCOORDINATE=-180
  _FillValue=-1
Corner Coordinates:
Upper Left  (-20015109.354,10007554.677) (124d17'29.21"W, 90d 0' 0.00"N)
Lower Left  (-20015109.354, 8895604.157) ( 43d25'16.73"E, 80d 0' 0.00"N)
Upper Right (-18903158.834,10007554.677) ( 77d21'53.56"W, 90d 0' 0.00"N)
Lower Right (-18903158.834, 8895604.157) (101d 0'32.47"E, 80d 0' 0.00"N)
Center      (-19459134.094, 9451579.417) (152d 6' 0.67"E, 85d 0' 0.00"N)
Band 1 Block=7200x3600 Type=Int16, ColorInterp=Gray
  Description = Black-sky Albedo in Visible band
  NoData Value=-1
  Unit Type: Albedo,no units
  Offset: 0,   Scale:
-- 
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20181130/b158979f/attachment.html>