Christopher, I have no problem reproducing the problem. Please file a ticket referencing this particular file. I'm going to look into it briefly now but I don't want the issue lost track of so it is important to have a ticket filed.
Best regards, Frank On Fri, May 18, 2012 at 9:57 AM, Christopher Mutel <[email protected]> wrote: > Frank- > > Here is the dataset in the given images: > > ftp://ftp.ntsg.umt.edu/pub/MODIS/Mirror/MOD16/MOD16A2_MONTHLY.MERRA_GMAO_1kmALB/Y2000/M01/MOD16A2.A2000M01.h12v12.105.2010357170459.hdf > > However, I see the same behaviour with all the files in the directory: > > ftp://ftp.ntsg.umt.edu/pub/MODIS/Mirror/MOD16/MOD16A2_MONTHLY.MERRA_GMAO_1kmALB/Y2000/M01/ > > Yours, Chris > > On Fri, May 18, 2012 at 6:34 PM, Frank Warmerdam <[email protected]> wrote: >> Christopher, >> >> Can you provide a pointer to one of the hdf4 images giving problems? >> >> I don't see how this could relate to disabling netcdf support. >> >> Best regards, >> >> On Fri, May 18, 2012 at 8:54 AM, Christopher Mutel <[email protected]> wrote: >>> Dear all- >>> >>> I am having trouble reading a HDF4 dataset using GDAL 1.9.0. The >>> datasets I am want to process come from the MODIS evapotranspiration >>> project (http://www.ntsg.umt.edu/project/mod16), and can be understood >>> correctly in e.g. HDFlook >>> (http://www-loa.univ-lille1.fr/Hdflook/hdflook_gb.html). They look >>> something like this: >>> >>> http://imgur.com/nw7ts >>> >>> The problem is that when GDAL opens a subdataset (any of the 4 >>> possible ones), the resulting image data >>> >>> http://imgur.com/fILW2 >>> >>> This effect happens using gdal_translate and the python bindings to >>> write the raw array data directly to a file. Basically, the lower >>> third of the image is fine; the next half of the image is masked >>> (given the nodata value); and the upper fifth, more peculiarly, has a >>> another section of the image repeated several times. >>> >>> I am using GDAL 1.9.0 installed through MacPorts. The supported >>> formats are listed below. The gdalinfo output is also below. >>> >>> My questions are: >>> >>> 1. Has anyone else experienced something like this? Is this fixable by >>> overriding some metadata in the gdal commands, or doing some other >>> transform? >>> >>> 2. Could this be related to the fact that I have NetCDF support enabled? >>> See e.g. >>> http://trac.osgeo.org/gdal/wiki/NetCDF >>> http://trac.osgeo.org/gdal/ticket/3957 >>> >>> Any help would be greatly appreciated, >>> Yours, >>> Chris >>> >>> =============== >>> Supported Formats: >>> =============== >>> GRASS (ro): GRASS Database Rasters (5.7+) >>> VRT (rw+v): Virtual Raster >>> GTiff (rw+v): GeoTIFF >>> NITF (rw+v): National Imagery Transmission Format >>> RPFTOC (rov): Raster Product Format TOC format >>> ECRGTOC (rov): ECRG TOC format >>> HFA (rw+v): Erdas Imagine Images (.img) >>> SAR_CEOS (rov): CEOS SAR Image >>> CEOS (rov): CEOS Image >>> JAXAPALSAR (rov): JAXA PALSAR Product Reader (Level 1.1/1.5) >>> GFF (rov): Ground-based SAR Applications Testbed File Format (.gff) >>> ELAS (rw+v): ELAS >>> AIG (rov): Arc/Info Binary Grid >>> AAIGrid (rwv): Arc/Info ASCII Grid >>> GRASSASCIIGrid (rov): GRASS ASCII Grid >>> SDTS (rov): SDTS Raster >>> DTED (rwv): DTED Elevation Raster >>> PNG (rwv): Portable Network Graphics >>> JPEG (rwv): JPEG JFIF >>> MEM (rw+): In Memory Raster >>> JDEM (rov): Japanese DEM (.mem) >>> GIF (rwv): Graphics Interchange Format (.gif) >>> BIGGIF (rov): Graphics Interchange Format (.gif) >>> ESAT (ro): Envisat Image Format >>> XPM (rwv): X11 PixMap Format >>> BMP (rw+v): MS Windows Device Independent Bitmap >>> DIMAP (rov): SPOT DIMAP >>> AirSAR (ro): AirSAR Polarimetric Image >>> RS2 (ro): RadarSat 2 XML Product >>> PCIDSK (rw+v): PCIDSK Database File >>> ILWIS (rw+v): ILWIS Raster Map >>> SGI (rw+): SGI Image File Format 1.0 >>> SRTMHGT (rwv): SRTMHGT File Format >>> Leveller (rw+): Leveller heightfield >>> Terragen (rw+): Terragen heightfield >>> GMT (rw): GMT NetCDF Grid Format >>> netCDF (rw+): Network Common Data Format >>> HDF4 (ro): Hierarchical Data Format Release 4 >>> HDF4Image (rw+): HDF4 Dataset >>> ISIS3 (rov): USGS Astrogeology ISIS cube (Version 3) >>> ISIS2 (rw+v): USGS Astrogeology ISIS cube (Version 2) >>> PDS (rov): NASA Planetary Data System >>> TIL (rov): EarthWatch .TIL >>> ERS (rw+v): ERMapper .ers Labelled >>> L1B (rov): NOAA Polar Orbiter Level 1b Data Set >>> FIT (rwv): FIT Image >>> RMF (rw+v): Raster Matrix Format >>> WCS (rov): OGC Web Coverage Service >>> WMS (rwv): OGC Web Map Service >>> MSGN (ro): EUMETSAT Archive native (.nat) >>> RST (rw+v): Idrisi Raster A.1 >>> INGR (rw+v): Intergraph Raster >>> GSAG (rwv): Golden Software ASCII Grid (.grd) >>> GSBG (rw+v): Golden Software Binary Grid (.grd) >>> GS7BG (rov): Golden Software 7 Binary Grid (.grd) >>> COSAR (ro): COSAR Annotated Binary Matrix (TerraSAR-X) >>> TSX (rov): TerraSAR-X Product >>> COASP (ro): DRDC COASP SAR Processor Raster >>> R (rwv): R Object Data Store >>> PNM (rw+v): Portable Pixmap Format (netpbm) >>> DOQ1 (rov): USGS DOQ (Old Style) >>> DOQ2 (rov): USGS DOQ (New Style) >>> ENVI (rw+v): ENVI .hdr Labelled >>> EHdr (rw+v): ESRI .hdr Labelled >>> GenBin (rov): Generic Binary (.hdr Labelled) >>> PAux (rw+): PCI .aux Labelled >>> MFF (rw+): Vexcel MFF Raster >>> MFF2 (rw+): Vexcel MFF2 (HKV) Raster >>> FujiBAS (ro): Fuji BAS Scanner Image >>> GSC (rov): GSC Geogrid >>> FAST (rov): EOSAT FAST Format >>> BT (rw+v): VTP .bt (Binary Terrain) 1.3 Format >>> LAN (rw+v): Erdas .LAN/.GIS >>> CPG (ro): Convair PolGASP >>> IDA (rw+): Image Data and Analysis >>> NDF (rov): NLAPS Data Format >>> EIR (rov): Erdas Imagine Raw >>> DIPEx (rov): DIPEx >>> LCP (rov): FARSITE v.4 Landscape File (.lcp) >>> GTX (rw+v): NOAA Vertical Datum .GTX >>> LOSLAS (rov): NADCON .los/.las Datum Grid Shift >>> NTv2 (rw+v): NTv2 Datum Grid Shift >>> ACE2 (rov): ACE2 >>> SNODAS (rov): Snow Data Assimilation System >>> RIK (ro): Swedish Grid RIK (.rik) >>> USGSDEM (rwv): USGS Optional ASCII DEM (and CDED) >>> GXF (ro): GeoSoft Grid Exchange Format >>> HTTP (ro): HTTP Fetching Wrapper >>> BAG (ro): Bathymetry Attributed Grid >>> HDF5 (ro): Hierarchical Data Format Release 5 >>> HDF5Image (ro): HDF5 Dataset >>> NWT_GRD (rov): Northwood Numeric Grid Format .grd/.tab >>> NWT_GRC (rov): Northwood Classified Grid Format .grc/.tab >>> ADRG (rw+v): ARC Digitized Raster Graphics >>> SRP (rov): Standard Raster Product (ASRP/USRP) >>> BLX (rw): Magellan topo (.blx) >>> PostGISRaster (ro): PostGIS Raster driver >>> SAGA (rw+v): SAGA GIS Binary Grid (.sdat) >>> KMLSUPEROVERLAY (rwv): Kml Super Overlay >>> XYZ (rwv): ASCII Gridded XYZ >>> HF2 (rwv): HF2/HFZ heightfield raster >>> OZI (rov): OziExplorer Image File >>> CTG (rov): USGS LULC Composite Theme Grid >>> E00GRID (rov): Arc/Info Export E00 GRID >>> ZMap (rwv): ZMap Plus Grid >>> NGSGEOID (rov): NOAA NGS Geoid Height Grids >>> >>> ======= >>> gdalinfo >>> ======= >>> >>> Driver: HDF4Image/HDF4 Dataset >>> Files: none associated >>> Size is 1200, 1200 >>> 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 = (-6671703.117999999783933,-3335851.558999999891967) >>> Pixel Size = (926.625433055833014,-926.625433055833355) >>> Metadata: >>> _FillValue=32767 >>> add_offset=0 >>> add_offset_err=0 >>> Affi.=NTSG, School of Forestry, U. of MT >>> ALGORITHMPACKAGEACCEPTANCEDATE=10-01-2004 >>> ALGORITHMPACKAGEMATURITYCODE=Normal >>> ALGORITHMPACKAGENAME=MODPR_15A2 >>> ALGORITHMPACKAGEVERSION=5 >>> ASSOCIATEDINSTRUMENTSHORTNAME=MODIS >>> ASSOCIATEDPLATFORMSHORTNAME=Terra >>> ASSOCIATEDSENSORSHORTNAME=MODIS >>> AUTOMATICQUALITYFLAG=Passed >>> AUTOMATICQUALITYFLAGEXPLANATION=No automatic quality assessment is >>> performed in the PGE >>> calibrated_nt=22 >>> CHARACTERISTICBINANGULARSIZE=30.0 >>> CHARACTERISTICBINSIZE=926.625433055556 >>> DATACOLUMNS=1200 >>> DATAROWS=1200 >>> DAYNIGHTFLAG=Day >>> DESCRREVISION=5.0 >>> EASTBOUNDINGCOORDINATE=-57.7254044077251 >>> [email protected] >>> EXCLUSIONGRINGFLAG=N >>> GEOANYABNORMAL=False >>> GEOESTMAXRMSERROR=50.0 >>> GLOBALGRIDCOLUMNS=43200 >>> GLOBALGRIDROWS=21600 >>> GRANULEBEGINNINGDATETIME=2006-12-11T09:38:17.000Z >>> GRANULEDAYNIGHTFLAG=Day >>> GRANULEENDINGDATETIME=2006-12-11T09:38:17.000Z >>> GRINGPOINTLATITUDE=-39.9999999964079, -29.8506269556485, >>> -29.8879227550404, -40.0438107454138 >>> GRINGPOINTLONGITUDE=-78.3244373487826, -69.2295648918179, >>> -57.4987926415755, -65.0468184138665 >>> GRINGPOINTSEQUENCENO=1, 2, 3, 4 >>> HDF-EOS file name=MOD16A2.A2000M01.h12v12.105.2010357170459.hdf >>> HDFEOSVersion=HDFEOS_V2.16 >>> HORIZONTALTILENUMBER=12 >>> INPUTPOINTER=MOD15A1.A2001008.h12v12.005.2006348070239.hdf, >>> MOD15A1.A2001007.h12v12.005.2006348001835.hdf, >>> MOD15A1.A2001006.h12v12.005.2006347202242.hdf, >>> MOD15A1.A2001005.h12v12.005.2006347152227.hdf, >>> MOD15A1.A2001004.h12v12.005.2006347093732.hdf, >>> MOD15A1.A2001003.h12v12.005.2006346135836.hdf, >>> MOD15A1.A2001002.h12v12.005.2006345172013.hdf, >>> MOD15A1.A2001001.h12v12.005.2006345093817.hdf, MOD15A2_ANC_RI4.hdf >>> LOCALGRANULEID=MOD15A2.A2001001.h12v12.005.2006351105708.hdf >>> LOCALVERSIONID=5.0.3 >>> long_name=MOD16A2 -- MODIS Gridded 1KM monthly Composite potential >>> Evapotranspiration (ET) >>> LONGNAME=MODIS/Terra Leaf Area Index/FPAR 8-Day L4 Global 1km SIN Grid >>> Maker=Qiaozhen Mu >>> MAXIMUMOBSERVATIONS=1 >>> MOD16 version=UM_collection5 >>> NADIRDATARESOLUTION=1km >>> NDAYS_COMPOSITED=8 >>> NORTHBOUNDINGCOORDINATE=-29.9999999973059 >>> NUMBEROFGRANULES=1 >>> PARAMETERNAME=MOD15A2 >>> PGEVERSION=5.0.5 >>> PROCESSINGCENTER=MODAPS >>> PROCESSINGENVIRONMENT=IRIX64 mtvs3 6.5 10070055 IP35 >>> PRODUCTIONDATETIME=2006-12-17T10:57:08.000Z >>> QAPERCENTCLOUDCOVER=5 >>> QAPERCENTEMPIRICALMODEL=2 >>> QAPERCENTGOODFPAR=98 >>> QAPERCENTGOODLAI=98 >>> QAPERCENTGOODQUALITY=98 >>> QAPERCENTINTERPOLATEDDATA=0 >>> QAPERCENTMAINMETHOD=98 >>> QAPERCENTMISSINGDATA=17 >>> QAPERCENTOTHERQUALITY=100 >>> QAPERCENTOUTOFBOUNDSDATA=17 >>> RANGEBEGINNINGDATE=2001-01-01 >>> RANGEBEGINNINGTIME=00:00:00 >>> RANGEENDINGDATE=2001-01-08 >>> RANGEENDINGTIME=23:59:59 >>> REPROCESSINGACTUAL=reprocessed >>> REPROCESSINGPLANNED=further update is anticipated >>> scale_factor=0.1 >>> scale_factor_err=0 >>> SCIENCEQUALITYFLAG=Not Investigated >>> SCIENCEQUALITYFLAGEXPLANATION=See >>> http://landweb.nascom/nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra >>> the product Science Quality status. >>> SHORTNAME=MOD15A2 >>> SOUTHBOUNDINGCOORDINATE=-39.9999999964079 >>> SPSOPARAMETERS=5367, 2680 >>> SYSTEMFILENAME=MOD15A1.A2001008.h12v12.005.2006348070239.hdf, >>> MOD15A1.A2001007.h12v12.005.2006348001835.hdf, >>> MOD15A1.A2001006.h12v12.005.2006347202242.hdf, >>> MOD15A1.A2001005.h12v12.005.2006347152227.hdf, >>> MOD15A1.A2001004.h12v12.005.2006347093732.hdf, >>> MOD15A1.A2001003.h12v12.005.2006346135836.hdf, >>> MOD15A1.A2001002.h12v12.005.2006345172013.hdf, >>> MOD15A1.A2001001.h12v12.005.2006345093817.hdf >>> TileID=51012012 >>> units=kg/m^2/mon >>> valid_range=-32767, 32700 >>> VERSIONID=5 >>> VERTICALTILENUMBER=12 >>> Webpage=http://www.ntsg.umt.edu >>> WESTBOUNDINGCOORDINATE=-78.3244373487825 >>> Corner Coordinates: >>> Upper Left (-6671703.118,-3335851.559) ( 69d16'55.32"W, 30d 0' 0.00"S) >>> Lower Left (-6671703.118,-4447802.079) ( 78d19'27.97"W, 40d 0' 0.00"S) >>> Upper Right (-5559752.598,-3335851.559) ( 57d44' 6.10"W, 30d 0' 0.00"S) >>> Lower Right (-5559752.598,-4447802.079) ( 65d16'13.31"W, 40d 0' 0.00"S) >>> Center (-6115727.858,-3891826.819) ( 67d 8'33.37"W, 35d 0' 0.00"S) >>> Band 1 Block=1200x833 Type=Int16, ColorInterp=Gray >>> Description = MOD16A2 -- MODIS Gridded 1KM monthly Composite >>> potential Evapotranspiration (ET) >>> NoData Value=32767 >>> Unit Type: kg/m^2/mon >>> Offset: 0, Scale:0.1 >>> >>> -- >>> ############################ >>> Chris Mutel >>> Ökologisches Systemdesign - Ecological Systems Design >>> Institut f.Umweltingenieurwissenschaften - Institute for Environmental >>> Engineering >>> ETH Zürich - HIF C 44 - Schafmattstr. 6 >>> 8093 Zürich >>> >>> Telefon: +41 44 633 71 45 - Fax: +41 44 633 10 61 >>> ############################ >>> _______________________________________________ >>> gdal-dev mailing list >>> [email protected] >>> http://lists.osgeo.org/mailman/listinfo/gdal-dev >> >> >> >> -- >> ---------------------------------------+-------------------------------------- >> I set the clouds in motion - turn up | Frank Warmerdam, [email protected] >> light and sound - activate the windows | http://pobox.com/~warmerdam >> and watch the world go round - Rush | Geospatial Software Developer > > > > -- > ############################ > Chris Mutel > Ökologisches Systemdesign - Ecological Systems Design > Institut f.Umweltingenieurwissenschaften - Institute for Environmental > Engineering > ETH Zürich - HIF C 44 - Schafmattstr. 6 > 8093 Zürich > > Telefon: +41 44 633 71 45 - Fax: +41 44 633 10 61 > ############################ > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev -- ---------------------------------------+-------------------------------------- I set the clouds in motion - turn up | Frank Warmerdam, [email protected] light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush | Geospatial Software Developer _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
