<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><br>
</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>Hi Folks,</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 Downloaded NetCDF files from <a id="LPlnk648854" href="http://land.copernicus.eu/global/products/ndvi">http://land.copernicus.eu/global/products/ndvi</a> and I am having problems with the GeoTransform that the driver is getting.</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>The image extents, as you can see in following gdalinfo report, goes beyond the -180 and +180 longitude degrees (half resolution, half pixel) but that's OK according to page 27 in the PDF documentation:
</span><span><span 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"><a title="https://earth.esa.int/documents/700255/1929094/PROBAV-Products_User_Manual_v1.3.pdf/fd5e30f4-5305-4d41-86aa-fab7b9568251" 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" href="https://earth.esa.int/documents/700255/1929094/PROBAV-Products_User_Manual_v1.3.pdf/fd5e30f4-5305-4d41-86aa-fab7b9568251">PROBAV-Products_User_Manual_v1.3.pdf</a> </span></span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span>
<div><br>
</div>
<div>But if I open that image on QGIS and overlay a vector (GADM from gadm.org) or any other global layer it is clear that -180x80 should be in the upper left corner of the image, not In the center. There is a half pixel missing on the left and bottom of the
image. This in in Fiji:</div>
<div><br>
</div>
<div><img tabindex="-1" style="max-width: 100%; -ms-user-select: none;" size="10134" contenttype="image/png" unselectable="on" data-outlook-trace="F:1|T:1" src="cid:fe6be096-c2db-42a9-a1c0-15358ccf71bb"><br>
</div>
<div><br>
</div>
<div>The gdalinfo report shows that the crs:GeoTransform is -180.0000000000 0.0029761905 0.0 80.0000000000 0.0 -0.0029761905 but it seems like the driver decided to take that to as center of the pixel based. Looking at the NetCDF source I can see that there
this possibility exists.</div>
</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;">
Of course, I am not waiting for a solution. I just used gdal_translate "- a_ullr" and fix that into the target format. It takes centuries to converted that image to geotiff. But I am concerned that others might get in trouble because of that.</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;">
I tried to contact the data provider but did not received an answer yet.</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;">
<div style="border-bottom-color: currentColor; border-left-color: currentColor; border-right-color: currentColor; border-top-color: currentColor; font-family: Calibri,Helvetica,sans-serif; font-size: 12pt; font-size-adjust: none; margin-bottom: 0px; margin-top: 0px">
Does anybody has experience with that dataset and/or NetCDF that can help figure out what is going on. It is a bug in the driver or an issue with the data.</div>
<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;">
-</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>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
Here is the gdalinfo report:</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;">
Don't worry about the warnings in regards to the pixel values. It is fine.</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><span>gdalinfo netcdf:c_gls_NDVI300_201910110000_GLOBE_PROBAV_V1.0.1.nc --debug on</span><br>
<div>GDAL_netCDF: driver detected file type=5, libnetcdf detected type=3</div>
<div>Warning 1: NetCDF driver detected file type=5, but libnetcdf detected type=3</div>
<div>GDAL_netCDF: setting file type to 3, was 5</div>
<div>GDAL_netCDF: dim_count = 2</div>
<div>GDAL_netCDF: var_count = 4</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute _FillValue</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute missing_value</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute flag_values</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute valid_range</div>
<div>GDAL_netCDF:</div>
<div>=====</div>
<div>SetProjectionFromVar( 3)</div>
<div><b>GDAL_netCDF: got grid_mapping crs</b></div>
<div>GDAL_netCDF: bIsGdalFile=1 bIsGdalCfFile=0 bBottomUp=0</div>
<div>GDAL_netCDF: got spheroid from CF: (6378137.000000 , 298.257224)</div>
<div>GDAL_netCDF: set bBottomUp = 0 from Y axis</div>
<div>GDAL_netCDF: setting WKT from CF</div>
<div>GDAL_netCDF: SetProjection, WKT = GEOGCS["unknown",DATUM["unknown",SPHEROID["Spheroid",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]]]</div>
<div>GDAL_netCDF: xdim: 120960 nSpacingBegin: 3 nSpacingMiddle: 3 nSpacingLast: 3</div>
<div>GDAL_netCDF: ydim: 47040 nSpacingBegin: -3 nSpacingMiddle: -3 nSpacingLast: -3</div>
<div>GDAL_netCDF: setting WKT from GDAL</div>
<div>GDAL_netCDF: SetProjection, WKT = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]</div>
<div><b>GDAL_netCDF: SetGeoTransform(-180.001488,0.002976,0.000000,80.001488,0.000000,-0.002976)</b></div>
<div>GDAL_netCDF: bGotGeogCS=1 bGotCfSRS=1 bGotCfGT=1 bGotGdalSRS=1 bGotGdalGT=0</div>
<div>Warning 1: Unsupported netCDF datatype (7), treat as Float32.</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute _FillValue</div>
<div>GDAL_netCDF: netcdf type=5 gdal type=6 signedByte=1</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute _FillValue</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute missing_value</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute flag_values</div>
<div>GDAL_netCDF: NCDFGetAttr unsupported type 7 for attribute valid_range</div>
<div>GDAL_netCDF: got add_offset=-0.07999999821186066, status=0</div>
<div>GDAL_netCDF: got scale_factor=0.004000000189989805, status=0</div>
<div>GDAL: GDALOpen(netcdf:c_gls_NDVI300_201910110000_GLOBE_PROBAV_V1.0.1.nc, this=00000219E2D291F0) succeeds as netCDF.</div>
<div>Driver: netCDF/Network Common Data Format</div>
<div>GDAL: GDALDefaultOverviews::OverviewScan()</div>
<div>Files: none associated</div>
<div>Size is 120960, 47040</div>
<div>Coordinate System is:</div>
<div>GEOGCS["WGS 84",</div>
<div> DATUM["WGS_1984",</div>
<div> SPHEROID["WGS 84",6378137,298.257223563,</div>
<div> AUTHORITY["EPSG","7030"]],</div>
<div> TOWGS84[0,0,0,0,0,0,0],</div>
<div> AUTHORITY["EPSG","6326"]],</div>
<div> PRIMEM["Greenwich",0,</div>
<div> AUTHORITY["EPSG","8901"]],</div>
<div> UNIT["degree",0.0174532925199433,</div>
<div> AUTHORITY["EPSG","9108"]],</div>
<div> AUTHORITY["EPSG","4326"]]</div>
<div><b>Origin = (-180.001488095238102,80.001488095238088)</b></div>
<div>Pixel Size = (0.002976190476204,-0.002976190476190)</div>
<div>Metadata:</div>
<div><b> crs#GeoTransform=-180.0000000000 0.0029761905 0.0 80.0000000000 0.0 -0.0029761905</b></div>
<div> crs#grid_mapping_name=latitude_longitude</div>
<div> crs#inverse_flattening=298.257223563</div>
<div> crs#longitude_of_prime_meridian=0</div>
<div> crs#long_name=coordinate reference system</div>
<div> crs#semi_major_axis=6378137</div>
<div> crs#spatial_ref=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]</div>
<div> crs#_CoordinateAxisTypes=GeoX GeoY</div>
<div> crs#_CoordinateTransformType=Projection</div>
<div> lat#axis=Y</div>
<div> lat#DIMENSION_LABELS=lat</div>
<div> lat#long_name=latitude</div>
<div> lat#standard_name=latitude</div>
<div> lat#units=degrees_north</div>
<div> lat#_CoordinateAxisType=Lat</div>
<div> lon#axis=X</div>
<div> lon#DIMENSION_LABELS=lon</div>
<div> lon#long_name=longitude</div>
<div> lon#standard_name=longitude</div>
<div> lon#units=degrees_east</div>
<div> lon#_CoordinateAxisType=Lon</div>
<div> NC_GLOBAL#archive_facility=VITO</div>
<div> NC_GLOBAL#Conventions=CF-1.6</div>
<div> NC_GLOBAL#copyright=Copernicus Service information 2019</div>
<div> NC_GLOBAL#history=Processing line NDVI: 2019-10-22</div>
<div> NC_GLOBAL#identifier=urn:cgls:global:ndvi300_v1_333m:NDVI300_201910110000_GLOBE_PROBAV_V1.0.1</div>
<div> NC_GLOBAL#institution=VITO NV</div>
<div> NC_GLOBAL#long_name=Normalized Difference Vegetation Index</div>
<div> NC_GLOBAL#orbit_type=LEO</div>
<div> NC_GLOBAL#parent_identifier=urn:cgls:global:ndvi300_v1_333m</div>
<div> NC_GLOBAL#platform=Proba-V</div>
<div> NC_GLOBAL#processing_level=L3</div>
<div> NC_GLOBAL#processing_mode=Near Real Time</div>
<div> NC_GLOBAL#product_version=V1.0.1</div>
<div> NC_GLOBAL#references=http://land.copernicus.eu/global/products/ndvi</div>
<div> NC_GLOBAL#sensor=VEGETATION</div>
<div> NC_GLOBAL#source=Derived from EO satellite imagery</div>
<div> NC_GLOBAL#time_coverage_end=2019-10-20T23:59:59Z</div>
<div> NC_GLOBAL#time_coverage_start=2019-10-11T00:00:00Z</div>
<div> NC_GLOBAL#title=10-daily Normalized Difference Vegetation Index 333M: GLOBE 2019-10-11T00:00:00Z</div>
<div> NDVI#add_offset=-0.079999998</div>
<div> NDVI#flag_meanings=Missing cloud snow sea background</div>
<div> NDVI#flag_values={}</div>
<div> NDVI#grid_mapping=crs</div>
<div> NDVI#long_name=Normalized Difference Vegetation Index 333M</div>
<div> NDVI#missing_value=</div>
<div> NDVI#scale_factor=0.0040000002</div>
<div> NDVI#standard_name=normalized_difference_vegetation_index</div>
<div> NDVI#units=</div>
<div> NDVI#valid_range={}</div>
<div> NDVI#_FillValue=</div>
<div>OGRCT: PROJ >= 4.8.0 features enabled</div>
<div>OGRCT: Using locale-safe proj version</div>
<div>OGRCT: Source: +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs</div>
<div>OGRCT: Target: +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs</div>
<div>Corner Coordinates:</div>
<div><b>Upper Left (-180.0014881, 80.0014881) (180d 0' 5.36"W, 80d 0' 5.36"N)</b></div>
<div><b>Lower Left (-180.0014881, -59.9985119) (180d 0' 5.36"W, 59d59'54.64"S)</b></div>
<div><b>Upper Right ( 179.9985119, 80.0014881) (179d59'54.64"E, 80d 0' 5.36"N)</b></div>
<div><b>Lower Right ( 179.9985119, -59.9985119) (179d59'54.64"E, 59d59'54.64"S)</b></div>
<div>Center ( -0.0014881, 10.0014881) ( 0d 0' 5.36"W, 10d 0' 5.36"N)</div>
<div>Band 1 Block=120960x1 Type=Float32, ColorInterp=Undefined</div>
<div> NoData Value=0</div>
<div> Offset: -0.0799999982118607, Scale:0.00400000018998981</div>
<div> Metadata:</div>
<div> add_offset=-0.079999998</div>
<div> flag_meanings=Missing cloud snow sea background</div>
<div> flag_values={}</div>
<div> grid_mapping=crs</div>
<div> long_name=Normalized Difference Vegetation Index 333M</div>
<div> missing_value=</div>
<div> NETCDF_VARNAME=NDVI</div>
<div> scale_factor=0.0040000002</div>
<div> standard_name=normalized_difference_vegetation_index</div>
<div> units=</div>
<div> valid_range={}</div>
<div> _FillValue=</div>
<div>GDAL: GDALClose(netcdf:c_gls_NDVI300_201910110000_GLOBE_PROBAV_V1.0.1.nc, this=00000219E2D291F0)</div>
<span></span><br>
</span></div>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Helvetica,sans-serif; font-size: 12pt;">
<span><br>
</span></div>
</body>
</html>