<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>Hi,</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><br>
</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>I am struggling with an issue with the NetCDF file with dealing with data from ESA Copernicus, Proba-V 1KM Level 3.</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><br>
</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>If I read the data using the HDF5 driver I lose the georeferenced but the values seems to more aceturate, since the pixels values over the ocean are supposed t be flagged as 254. See the Python printout:</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><br>
</span></div>
<blockquote style="padding: 0px; margin-top: 0px; margin-right: 0px; margin-bottom: 0px;" dir="ltr">
<blockquote style="padding: 0px; margin-top: 0px; margin-right: 0px; margin-bottom: 0px;" dir="ltr">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><span style="font-family: "Courier New",monospace; font-size: 11pt;">from osgeo import gdal</span><span><br>
</span></span></div>
</blockquote>
</blockquote>
<span>
<blockquote style="padding: 0px; margin-top: 0px; margin-right: 0px; margin-bottom: 0px;" dir="ltr">
<blockquote style="padding: 0px; margin-top: 0px; margin-right: 0px; margin-bottom: 0px;" dir="ltr">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">import numpy</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"></span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ds = gdal.Open('HDF5:"c_gls_NDVI_201907110000_GLOBE_PROBAV_V2.2.1.nc"://NDVI')</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">bd = ds.GetRasterBand(1)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ar = bd.ReadAsArray(25109,9083,10,10)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">print(ar)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ds = None</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"></span><span style="font-family: "Courier New",monospace; font-size: 11pt;">[[254 254 254 254 254 254 254 254 254 254]</span><span><br>
</span>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]</span></div>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;"> [254 254 254 254 254 254 254 254 254 254]]</span></div>
<div><br>
</div>
<span></span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ds = gdal.Open('NETCDF:"c_gls_NDVI_201907110000_GLOBE_PROBAV_V2.2.1.nc":NDVI')</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">bd = ds.GetRasterBand(1)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ar = bd.ReadAsArray(25109,9083,10,10)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">print(ar)</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">ds = None</span></div>
<span></span>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">[[255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]</span></div>
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> [255 255 255 255 255 255 255 255 255 255]]</span></div>
</blockquote>
</blockquote>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
By looking at the gdalinfo information it clear that netCDF gathers more information, like I said, georeferenced but both drivers are not capable to show the list of flagged values on the NDVI dataset:</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<blockquote style="padding: 0px; margin-top: 0px; margin-right: 0px; margin-bottom: 0px;" dir="ltr">
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">$ gdalinfo </span><span style="font-family: "Courier New",monospace; font-size: 11pt; display: inline !important; font-size-adjust: none; background-color: rgb(255, 255, 255);">HDF5:"c_gls_NDVI_201907110000_GLOBE_PROBAV_V2.2.1.nc"://NDVI</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">Driver: HDF5Image/HDF5 Dataset</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">. . .</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> NDVI_flag_meanings=Missing cloud snow sea background</span><span><br>
</span><span style="font-family: "Courier New",monospace; font-size: 11pt;"> NDVI_flag_values=√ⁿ²■ </span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">$ gdalinfo </span><span style="font-family: "Courier New",monospace; font-size: 11pt; display: inline !important; font-size-adjust: none; background-color: rgb(255, 255, 255);">NETCDF:"c_gls_NDVI_201907110000_GLOBE_PROBAV_V2.2.1.nc":NDVI</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">Warning 1: NetCDF driver detected file type=5, but libnetcdf detected type=3</span><span><br>
</span>
<div><span style="font-family: "Courier New",monospace; font-size: 11pt;">Warning 1: Unsupported netCDF datatype (7), treat as Float32.</span></div>
<span style="font-family: "Courier New",monospace; font-size: 11pt;">Driver: netCDF/Network Common Data Format</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;">. . .</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span style="font-family: "Courier New",monospace; font-size: 11pt;"> NDVI#flag_meanings=Missing cloud snow sea background</span><span><br>
</span><span style="font-family: "Courier New",monospace; font-size: 11pt;"> NDVI#flag_values={}</span></div>
</blockquote>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
There is this warning but that doesn't seems to be affecting the pixels values on the land.</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Does anybody has a clue on what is going on here?</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Thanks,</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Ivan</div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<br>
</div>
</span>
</body>
</html>