Thanks for clarifying!<div><br></div><div>Cheers, Mike<br><div dir="auto"><br><div class="gmail_quote"><div dir="ltr">On Sat, Dec 1, 2018, 02:04 Andre Joost <<a href="mailto:andre%2Bjoost@nurfuerspam.de">andre+joost@nurfuerspam.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Michael,<br>
<br>
running eosdump on the file, I get:<br>
<br>
File Size: 94926773<br>
Number of grids: 1<br>
Number of swaths: 0<br>
Number of points: 0<br>
Grid: GLASS02B06<br>
        projection: SNSOID<br>
        griddataclash:<br>
<br>
<br>
   GRID GLASS02B06 {{{<br>
     xdim: 7200<br>
     ydim: 3600<br>
     upleft: -2.00151e+07, 1.00076e+07<br>
     lowright: -1.89032e+07, 8.8956e+06<br>
     PROJECTION {<br>
       projection: 16<br>
       zone: -1<br>
       sphere: -1<br>
       param: ...<br>
         6.37101e+06 0 0 0 0 0 0 0<br>
         0 0 0 0 0 3.14511e+35 4.21472e-153 1.72352e+45<br>
       pix: 0<br>
       origin: 0<br>
   }<br>
<br>
So they defined the wrong projection, and gdal can't do much about it.<br>
<br>
You can sanitize it with gdal_translate -a_ullr -180 90 180 -90 -a_srs <br>
EPSG:4326 src dst<br>
<br>
HTH,<br>
Andre Joost<br>
<br>
Am 29.11.18 um 21:31 schrieb Michael Sumner:<br>
> Hello, these data in HDF4 format are interpreted as being in sinusoidal<br>
> projection, I guess some conflation with the MODIS logic about some<br>
> products (?).<br>
><br>
> <a href="ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/" rel="noreferrer" target="_blank">ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/</a><br>
><br>
> Particular file is here, but note there's also an auxiliarly .hdf.xml with<br>
> it:<br>
> <a href="ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/GLASS02B06.V04.A2013169.2017128.hdf" rel="noreferrer" target="_blank">ftp://ftp.glcf.umd.edu/glcf/GLASS/ABD/MODIS/0.05D/2013/GLASS02B06.V04.A2013169.2017128.hdf</a><br>
><br>
> A plot of the first subdataset shows that it is trivially in regular<br>
> longlat (or eqc), the gdalinfo output is below.<br>
><br>
> I'd appreciate any insights into why this is interpreted this way, and<br>
> whether the GDAL detection logic could be improved - or whether the data<br>
> providers need to update the internal metadata.<br>
><br>
> Cheers, Mike.<br>
><br>
> gdalinfo --version GDAL 2.3.2, released 2018/09/21<br>
><br>
> gdalinfo<br>
> HDF4_EOS:EOS_GRID:"GLASS02B06.V04.A2013169.2017128.hdf":GLASS02B06:ABD_BSA_VIS<br>
> Driver: HDF4Image/HDF4 Dataset<br>
> Files: GLASS02B06.V04.A2013169.2017128.hdf<br>
> Size is 7200, 3600<br>
> Coordinate System is:<br>
> PROJCS["unnamed",<br>
>      GEOGCS["Unknown datum based upon the custom spheroid",<br>
>          DATUM["Not specified (based on custom spheroid)",<br>
>              SPHEROID["Custom spheroid",6371007.181,0]],<br>
>          PRIMEM["Greenwich",0],<br>
>          UNIT["degree",0.0174532925199433]],<br>
>      PROJECTION["Sinusoidal"],<br>
>      PARAMETER["longitude_of_center",0],<br>
>      PARAMETER["false_easting",0],<br>
>      PARAMETER["false_northing",0],<br>
>      UNIT["Meter",1]]<br>
> Origin = (-20015109.353999998420477,10007554.676999999210238)<br>
> Pixel Size = (154.437572175972150,-308.875144351944300)<br>
> Metadata:<br>
>    add_offset=0<br>
>    add_offset_err=0<br>
>    ASSOCIATEDINSTRUMENTSHORTNAME.1=AVHRR<br>
>    ASSOCIATEDPLATFORMSHORTNAME.1=NOAA<br>
>    ASSOCIATEDSENSORSHORTNAME.1=AVHRR<br>
>    BANDDEFINITION=Surface<br>
>    CHARACTERISTICBINSIZE=926.6<br>
>    DATACOLUMNS=7200<br>
>    DATAROWS=3600<br>
>    DATASETNAME=GLASS02B06<br>
>    DAYNIGHTFLAG=Both<br>
>    EASTHBOUNDINGCOORDINATE=180<br>
>    EXCLUSIONGRINGFLAG=N<br>
>    GLOBALGRIDCOLUMNS=7200<br>
>    GLOBALGRIDROWS=3600<br>
>    GRIDTYPE=Geographic Lat/Lon<br>
>    GRINGPOINTLATITUDE.1=90,90,-90,-90<br>
>    GRINGPOINTLONGITUDE.1=-180,180,180,-180<br>
>    GRINGPOINTSEQUENCENO.1=1,2,3,4<br>
>    HDFEOSVersion=HDFEOS_V2.16<br>
>    HORIZONTALTILENUMBER=00<br>
>    HORIZONTALTILENUMBER=00<br>
>    INPUTFILESNAME=Multi-Source Data<br>
>    INPUTPOINTER=Multi-Source Data<br>
>    INSTITUTENAME=BEIJING<br>
>    LOCALGRANULEID=GLASS02B06.V04.A2013169.2017128.hdf<br>
>    LOCALVERSIONID=V04<br>
>    long_name=Black-sky Albedo in Visible band<br>
>    NORTHBOUNDINGCOORDINATE=90<br>
>    PARAMETERNAME.1=GLASS02B06<br>
>    PGEVERSION=V04<br>
>    PROCESSINGENVIRONMENT=High Performance Computing System<br>
>    PRODUCTDATEANDTIME=2017-05-08 16:00:12<br>
>    PRODUCTIONDATETIME=2017-05-08 16:00:12<br>
>    PRODUCTNAME=Surface<br>
>    PRODUCTQUALITY=CLOUDREMOVED<br>
>    REPROCESSINGACTUAL=reprocessed<br>
>    REPROCESSINGPLANNED=further update is anticipated<br>
>    scale_factor=0.0001<br>
>    scale_factor_err=0<br>
>    SHORTNAME=GLASS02B06<br>
>    SOUTHBOUNDINGCOORDINATE=-90<br>
>    TILENUMBER=H00V00<br>
>    units=Albedo,no units<br>
>    valid_range=0, 10000<br>
>    VERSIONID=1<br>
>    VERTICALTILENUMBER=00<br>
>    VERTICALTILENUMBER=00<br>
>    WEBSITE=<a href="http://glass-product.bnu.edu.cn" rel="noreferrer" target="_blank">http://glass-product.bnu.edu.cn</a><br>
>    WESTHBOUNDINGCOORDINATE=-180<br>
>    _FillValue=-1<br>
> Corner Coordinates:<br>
> Upper Left  (-20015109.354,10007554.677) (124d17'29.21"W, 90d 0' 0.00"N)<br>
> Lower Left  (-20015109.354, 8895604.157) ( 43d25'16.73"E, 80d 0' 0.00"N)<br>
> Upper Right (-18903158.834,10007554.677) ( 77d21'53.56"W, 90d 0' 0.00"N)<br>
> Lower Right (-18903158.834, 8895604.157) (101d 0'32.47"E, 80d 0' 0.00"N)<br>
> Center      (-19459134.094, 9451579.417) (152d 6' 0.67"E, 85d 0' 0.00"N)<br>
> Band 1 Block=7200x3600 Type=Int16, ColorInterp=Gray<br>
>    Description = Black-sky Albedo in Visible band<br>
>    NoData Value=-1<br>
>    Unit Type: Albedo,no units<br>
>    Offset: 0,   Scale:<br>
><br>
><br>
><br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
><br>
<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div></div></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span>Dr. Michael Sumner</span><br><span>Software and Database Engineer</span><br><span>Australian Antarctic Division</span><br><span>203 Channel Highway</span><br><span>Kingston Tasmania 7050 Australia</span><br><br></div></div>