Bug#690884: qgis fails to access HDF-EOS datasets directly

2013-09-29 Thread Sebastiaan Couwenberg
Hi Ivan,

I've tried to reproduce your problem with QGIS 2.0.1, but I'm not able
to download the HDF file.

ftp://e4ftl01.cr.usgs.gov/MODIS_Composites/MOTA/MCD12Q1.005/2008.01.01/MCD12Q1.A2008001.h22v03.005.2011215183430.hdf

The FTP server claims to allow anonymous access with an email address as
password, but refuses login when done so.

Is there some HDF file publicly available with which I can reproduce
your problem?

If you don't mind building QGIS yourself, you can try building QGIS
2.0.1 from the pkg-grass git to see if it fixes your problem:

http://anonscm.debian.org/gitweb/?p=pkg-grass/qgis.git

Kind Regards,

Bas

-- 
GnuPG: 0xE88D4AF1 (new) / 0x77A975AD (old)

___
Pkg-grass-devel mailing list
Pkg-grass-devel@lists.alioth.debian.org
http://lists.alioth.debian.org/cgi-bin/mailman/listinfo/pkg-grass-devel


Bug#690884: qgis fails to access HDF-EOS datasets directly

2012-10-18 Thread Ivan Shmakov
Package: qgis
Version: 1.7.4+1.7.5~20120320-1.1+b1

Somehow, qgis fails to access HDF-EOS (MODIS L3, as per [1])
datasets “directly.”

E. g., trying to “Add Raster Layer” against, say,
MCD12Q1.A2008001.h22v03.005.2011215183430.hdf [2], results in
a warning (“Cannot get GDAL raster band:”), and also in the
following output at the terminal:

Warning: QgsRasterLayer::setDataProvider: Data provider is invalid.
ERROR 10: Pointer 'hObject' is NULL in 'GDALGetMetadata'.

Naturally, the gdal-bin tools accept such a file (consider,
e. g., the gdalinfo(1) output MIME'd.)

However, should I create a .vrt dataset out of particular
HDF-EOS “subdataset” (which a single HDF-EOS file can contain an
arbitrary number of) with gdal_translate(1) (see below), the
resulting dataset is accepted by qgis.

$ gdal_translate -of VRT \
  
HDF4_EOS:EOS_GRID:/.../files/e4ftl01u.ecs.nasa.gov/MODIS_Composites/MOTA/MCD12Q1.005/2008.01.01/MCD12Q1.A2008001.h22v03.005.2011215183430.hdf:MOD12Q1:Land_Cover_Type_1
 \
  foo.vrt 

ISTR, that there was some older version of QGIS (though perhaps
not in Debian) that was presenting a list of “subdatasets”
contained within the HDF-EOS file selected.  Thus, this bug may
be a regression.

[1] 
http://en.wikipedia.org/wiki/Moderate-Resolution_Imaging_Spectroradiometer#MODIS_Level_3_datasets
[2] 
ftp://e4ftl01.cr.usgs.gov/MODIS_Composites/MOTA/MCD12Q1.005/2008.01.01/MCD12Q1.A2008001.h22v03.005.2011215183430.hdf

-- 
FSF associate member #7257
Driver: HDF4/Hierarchical Data Format Release 4
Files: 
files/e4ftl01u.ecs.nasa.gov/MODIS_Composites/MOTA/MCD12Q1.005/2008.01.01/MCD12Q1.A2008001.h22v03.005.2011215183430.hdf
Size is 512, 512
Coordinate System is `'
Metadata:
  ALGORITHMPACKAGEACCEPTANCEDATE=1998-01-01
  ALGORITHMPACKAGEMATURITYCODE=LAUNCH
  ALGORITHMPACKAGENAME=MOD12Q1
  ALGORITHMPACKAGEVERSION=V2.0
  ASSOCIATEDINSTRUMENTSHORTNAME=MODIS
  ASSOCIATEDINSTRUMENTSHORTNAME=MODIS
  ASSOCIATEDPLATFORMSHORTNAME=Aqua
  ASSOCIATEDPLATFORMSHORTNAME=Terra
  ASSOCIATEDSENSORSHORTNAME=MODIS
  ASSOCIATEDSENSORSHORTNAME=MODIS
  AUTOMATICQUALITYFLAG=Passed
  AUTOMATICQUALITYFLAGEXPLANATION=To be set as 'Passed' or 'Failed' to indicate 
failure of PGE test.
  CHARACTERISTICBINANGULARSIZE=15.0
  CHARACTERISTICBINSIZE=463.3127165
  COVERAGECALCULATIONMETHOD=volume
  DATACOLUMNS=2400
  DATAROWS=2400
  DAYNIGHTFLAG=Day
  DESCRREVISION=5.0
  EASTBOUNDINGCOORDINATE=100.016667
  EXCLUSIONGRINGFLAG=N
  GEOANYABNORMAL=False
  GEOESTMAXRMSERROR=50.0
  GLOBALGRIDCOLUMNS=86400
  GLOBALGRIDROWS=43200
  GRINGPOINTLATITUDE=49.762291, 60.09218, 59.929919, 49.676752
  GRINGPOINTLONGITUDE=61.938978, 80.073996, 100.833098, 78.016106
  GRINGPOINTSEQUENCENO=1, 2, 3, 4
  HDFEOSVersion=HDFEOS_V2.9
  HORIZONTALTILENUMBER=22
  INPUTPOINTER=MOD12Q1.A2004001.h22v03.004.2006117174145.hdf, 
MOD12Q1_TEMP.A2008001.h22v03.002.2009278182148.hdf, LC1.A2008001.h22v03.hdf, 
LC2.A2008001.h22v03.hdf, LC3.A2008001.h22v03.hdf, LC4.A2008001.h22v03.hdf, 
LC5.A2008001.h22v03.hdf, LC1A.A2008001.h22v03.hdf, LC2A.A2008001.h22v03.hdf, 
LC3A.A2008001.h22v03.hdf, LC4A.A2008001.h22v03.hdf, LC5A.A2008001.h22v03.hdf, 
LCS.A2008001.h22v03.hdf, LCS_P.A2008001.h22v03.hdf, LCP1.A2008001.h22v03.hdf, 
LCP2.A2008001.h22v03.hdf, LCP3.A2008001.h22v03.hdf
  LOCALGRANULEID=MCD12Q1.A2008001.h22v03.005.2011215183430.hdf
  LOCALVERSIONID=2.2.1
  LONGNAME=MODIS/Terra+Aqua Land Cover Type Yearly L3 Global 500m SIN Grid
  MAXIMUMOBSERVATIONS=2
  NADIRDATARESOLUTION=500m
  NORTHBOUNDINGCOORDINATE=60.0
  NUMBEROFGRANULES=5
  PARAMETERNAME=Land_Cover_Type_1
  PGEVERSION=5.1.7
  PROCESSINGCENTER=MODAPS
  PROCESSINGENVIRONMENT=Linux minion5557 2.6.18-238.19.1.el5PAE #1 SMP Fri Jul 
15 08:15:44 EDT 2011 i686 i686 i386 GNU/Linux
  PRODUCTIONDATETIME=2011-08-03T18:34:30.000Z
  QAPERCENTGOODQUALITY=50
  QAPERCENTINTERPOLATEDDATA=14
  QAPERCENTMISSINGDATA=0
  QAPERCENTNOTPRODUCEDCLOUD=0
  QAPERCENTNOTPRODUCEDOTHER=1
  QAPERCENTOTHERQUALITY=49
  QAPERCENTOUTOFBOUNDSDATA=0
  RANGEBEGINNINGDATE=2008-01-01
  RANGEBEGINNINGTIME=04:45:00.00
  RANGEENDINGDATE=2008-12-31
  RANGEENDINGTIME=06:40:00.00
  REPROCESSINGACTUAL=reprocessed
  REPROCESSINGPLANNED=further update is anticipated
  SCIENCEQUALITYFLAG=Not Investigated
  SCIENCEQUALITYFLAGEXPLANATION=See 
http://landweb.nascom/nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=aqua the 
product Science Quality status.
  SHORTNAME=MCD12Q1
  SOUTHBOUNDINGCOORDINATE=50.0
  SPSOPARAMETERS=2669
  TileID=51022003
  VERSIONID=5
  VERTICALTILENUMBER=03
  WESTBOUNDINGCOORDINATE=62.228953
Subdatasets:
  
SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:files/e4ftl01u.ecs.nasa.gov/MODIS_Composites/MOTA/MCD12Q1.005/2008.01.01/MCD12Q1.A2008001.h22v03.005.2011215183430.hdf:MOD12Q1:Land_Cover_Type_1
  SUBDATASET_1_DESC=[2400x2400] Land_Cover_Type_1 MOD12Q1 (8-bit unsigned 
integer)