I need to convert an hdf4 file by MODIS, the MOD13A3 v006 (NDVI Index). In particular, I'm interested in the block "h18v04", which is north-west Europe. You can download the example file I'm using in this question here (16 MB).
What I've done, taking inspiration by this already asked question, is the following python script:
import os
filename = "MOD13A3.A2019001.h18v04.006.2019038102530.hdf"
subdataset = 'HDF4_EOS:EOS_GRID:"MOD13A3.A2019001.h18v04.006.2019038102530.hdf"://MOD_Grid_monthly_1km_VI:"1 km monthly NDVI"'
f_out = "prova.tif"
s_srs = "+proj=sinu +R=6371007.181 +nadgrids=@null +wktext"
t_srs = "+proj=longlat +datum=WGS84 +no_defs"
command = 'gdalwarp -of GTiff -s_srs "' + s_srs + '" -r cubic -t_srs "' + t_srs + '" ' + subdataset + ' ' + f_out
os.system(command)
But it returns to me the following error:
Warning 1: Unknown type -1
ERROR 1: Input file HDF4_EOS:EOS_GRID:MOD13A3.A2019001.h18v04.006.2019038102530.hdf://MOD_Grid_monthly_1km_VI:1 km monthly NDVI has no raster bands.
Giving gdalinfo, I get
gdalinfo HDF4_EOS:EOS_GRID:"MOD13A3.A2019001.h18v04.006.2019038102530.hdf"://MOD_Grid_monthly_1km_VI:"1 km monthly NDVI"
Warning 1: Unknown type -1
Driver: HDF4Image/HDF4 Dataset
Files: MOD13A3.A2019001.h18v04.006.2019038102530.hdf
MOD13A3.A2019001.h18v04.006.2019038102530.hdf.aux.xml
Size is 0, 1
Metadata:
ALGORITHMPACKAGEACCEPTANCEDATE=102004
ALGORITHMPACKAGEMATURITYCODE=Normal
ALGORITHMPACKAGENAME=MOD_PR13A3
ALGORITHMPACKAGEVERSION=5
ASSOCIATEDINSTRUMENTSHORTNAME.1=MODIS
ASSOCIATEDPLATFORMSHORTNAME.1=Terra
ASSOCIATEDSENSORSHORTNAME.1=MODIS
AUTOMATICQUALITYFLAG.1=Passed
AUTOMATICQUALITYFLAG.10=Passed
AUTOMATICQUALITYFLAG.11=Passed
AUTOMATICQUALITYFLAG.2=Passed
AUTOMATICQUALITYFLAG.3=Passed
AUTOMATICQUALITYFLAG.4=Passed
AUTOMATICQUALITYFLAG.5=Passed
AUTOMATICQUALITYFLAG.6=Passed
AUTOMATICQUALITYFLAG.7=Passed
AUTOMATICQUALITYFLAG.8=Passed
AUTOMATICQUALITYFLAG.9=Passed
AUTOMATICQUALITYFLAGEXPLANATION.1=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.10=No automatic quality assessment
is performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.11=No automatic quality assessment
is performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.2=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.3=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.4=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.5=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.6=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.7=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.8=No automatic quality assessment is
performed in the PGE
AUTOMATICQUALITYFLAGEXPLANATION.9=No automatic quality assessment is
performed in the PGE
CHARACTERISTICBINANGULARSIZE=30.0
CHARACTERISTICBINSIZE=926.625433055556
COMPOSITEMONTHNAME=January
DATACOLUMNS=1200
DATAROWS=1200
DAYNIGHTFLAG=Both
DAYSOFYEAR=1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,
17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
DAYSPROCESSED=2019001, 2019002, 2019003, 2019004, 2019005, 2019006,
2019007, 2019008, 2019009, 2019010, 2019011, 2019012, 2019013,
2019014, 2019015, 2019016, 2019017, 2019018, 2019019, 2019020,
2019021, 2019022, 2019023, 2019024, 2019025, 2019026, 2019027,
2019028, 2019029, 2019030, 2019031, 2019032, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1
DESCRREVISION=6.0
EASTBOUNDINGCOORDINATE=15.5702026309755
EVI1KMMONTHQCLASSPERCENTAGE=0
EXCLUSIONGRINGFLAG.1=N
GEOANYABNORMAL=False
GEOESTMAXRMSERROR=50.0
GLOBALGRIDCOLUMNS=43200
GLOBALGRIDROWS=21600
GRINGPOINTLATITUDE.1=39.8197706210566, 50.0069689864989,
49.9989721995752, 39.8144150819193
GRINGPOINTLONGITUDE.1=0.000131401253181314, -0.00867479214979005,
15.5723926647109, 13.0379032886097
GRINGPOINTSEQUENCENO.1=1, 2, 3, 4
HDFEOSVersion=HDFEOS_V2.19
HORIZONTALTILENUMBER=18
identifier_product_doi=10.5067/MODIS/MOD13A3.006
identifier_product_doi=10.5067/MODIS/MOD13A3.006
identifier_product_doi_authority=http://dx.doi.org
identifier_product_doi_authority=http://dx.doi.org
INPUTFILENAME=MOD13A2.A2019001.h18v04.006.2019024152607.hdf,
MOD13A2.A2019017.h18v04.006.2019037160749.hdf
INPUTPOINTER=MOD13A2.A2019001.h18v04.006.2019024152607.hdf,
MOD13A2.A2019017.h18v04.006.2019037160749.hdf
INSTRUMENTNAME=Moderate-Resolution Imaging SpectroRadiometer
LOCALGRANULEID=MOD13A3.A2019001.h18v04.006.2019038102530.hdf
LOCALVERSIONID=5.2.2
LONGNAME=MODIS/Terra Vegetation Indices Monthly L3 Global 1km SIN
Grid
NDVI1KMMONTHQCLASSPERCENTAGE=0
NORTHBOUNDINGCOORDINATE=49.9999999955098
NUMBEROF16DAYPERIODS=2
NUMBEROFDAYS=31
PARAMETERNAME.1=1 km monthly NDVI
PARAMETERNAME.10=1 km monthly relative azimuth angle
PARAMETERNAME.11=1 km monthly pixel reliability
PARAMETERNAME.2=1 km monthly EVI
PARAMETERNAME.3=1 km monthly VI Quality
PARAMETERNAME.4=1 km monthly red reflectance
PARAMETERNAME.5=1 km monthly NIR reflectance
PARAMETERNAME.6=1 km monthly blue reflectance
PARAMETERNAME.7=1 km monthly MIR reflectance
PARAMETERNAME.8=1 km monthly view zenith angle
PARAMETERNAME.9=1 km monthly sun zenith angle
PERCENTLAND=73
PGEVERSION=6.0.3
PROCESSINGCENTER=MODAPS
PROCESSINGENVIRONMENT=Linux minion7032 3.10.0-957.5.1.el7.x86_64 #1
SMP Fri Feb 1 14:54:57 UTC 2019 x86_64 x86_64 x86_64 GNU/Linux
PRODUCTIONDATETIME=2019-02-07T15:25:30.000Z
QAPERCENTCLOUDCOVER.1=11
QAPERCENTCLOUDCOVER.10=11
QAPERCENTCLOUDCOVER.11=11
QAPERCENTCLOUDCOVER.2=11
QAPERCENTCLOUDCOVER.3=11
QAPERCENTCLOUDCOVER.4=11
QAPERCENTCLOUDCOVER.5=11
QAPERCENTCLOUDCOVER.6=11
QAPERCENTCLOUDCOVER.7=11
QAPERCENTCLOUDCOVER.8=11
QAPERCENTCLOUDCOVER.9=11
QAPERCENTGOODQUALITY=50
QAPERCENTINTERPOLATEDDATA.1=100
QAPERCENTINTERPOLATEDDATA.10=100
QAPERCENTINTERPOLATEDDATA.11=100
QAPERCENTINTERPOLATEDDATA.2=100
QAPERCENTINTERPOLATEDDATA.3=100
QAPERCENTINTERPOLATEDDATA.4=100
QAPERCENTINTERPOLATEDDATA.5=100
QAPERCENTINTERPOLATEDDATA.6=100
QAPERCENTINTERPOLATEDDATA.7=100
QAPERCENTINTERPOLATEDDATA.8=100
QAPERCENTINTERPOLATEDDATA.9=100
QAPERCENTMISSINGDATA.1=0
QAPERCENTMISSINGDATA.10=0
QAPERCENTMISSINGDATA.11=0
QAPERCENTMISSINGDATA.2=0
QAPERCENTMISSINGDATA.3=0
QAPERCENTMISSINGDATA.4=0
QAPERCENTMISSINGDATA.5=0
QAPERCENTMISSINGDATA.6=0
QAPERCENTMISSINGDATA.7=0
QAPERCENTMISSINGDATA.8=0
QAPERCENTMISSINGDATA.9=0
QAPERCENTNOTPRODUCEDCLOUD=0
QAPERCENTNOTPRODUCEDOTHER=0
QAPERCENTOTHERQUALITY=50
QAPERCENTOUTOFBOUNDSDATA.1=0
QAPERCENTOUTOFBOUNDSDATA.10=0
QAPERCENTOUTOFBOUNDSDATA.11=0
QAPERCENTOUTOFBOUNDSDATA.2=0
QAPERCENTOUTOFBOUNDSDATA.3=0
QAPERCENTOUTOFBOUNDSDATA.4=0
QAPERCENTOUTOFBOUNDSDATA.5=0
QAPERCENTOUTOFBOUNDSDATA.6=0
QAPERCENTOUTOFBOUNDSDATA.7=0
QAPERCENTOUTOFBOUNDSDATA.8=0
QAPERCENTOUTOFBOUNDSDATA.9=0
QAPERCENTPOORQ1KMMONTHEVI=0, 47, 15, 6, 10, 11, 5, 4, 2, 0, 0, 0, 0,
0, 0, 0
QAPERCENTPOORQ1KMMONTHNDVI=0, 47, 15, 6, 10, 11, 5, 4, 2, 0, 0, 0,
0, 0, 0, 0
QA_STRUCTURE_STYLE=C5 or later
RANGEBEGINNINGDATE=2019-01-01
RANGEBEGINNINGTIME=00:00:00
RANGEENDINGDATE=2019-01-31
RANGEENDINGTIME=23:59:59
REPROCESSINGACTUAL=reprocessed
REPROCESSINGPLANNED=further update is anticipated
SCIENCEQUALITYFLAG.1=Not Investigated
SCIENCEQUALITYFLAG.10=Not Investigated
SCIENCEQUALITYFLAG.11=Not Investigated
SCIENCEQUALITYFLAG.2=Not Investigated
SCIENCEQUALITYFLAG.3=Not Investigated
SCIENCEQUALITYFLAG.4=Not Investigated
SCIENCEQUALITYFLAG.5=Not Investigated
SCIENCEQUALITYFLAG.6=Not Investigated
SCIENCEQUALITYFLAG.7=Not Investigated
SCIENCEQUALITYFLAG.8=Not Investigated
SCIENCEQUALITYFLAG.9=Not Investigated
SCIENCEQUALITYFLAGEXPLANATION.1=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.10=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.11=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.2=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.3=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.4=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.5=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.6=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.7=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.8=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SCIENCEQUALITYFLAGEXPLANATION.9=See
http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra
for the product Science Quality status.
SEAPROCESSED=Yes
SHORTNAME=MOD13A3
SOUTHBOUNDINGCOORDINATE=39.9999999964079
SPSOPARAMETERS=2749, 4334, 2749a, 4334a
TileID=51018004
VERSIONID=6
VERTICALTILENUMBER=04
WESTBOUNDINGCOORDINATE=0.0
Can you help me?
