<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Thanks Even,<div>So doing that command, and then looking at;</div><div>"gdalinfo -hist out.vrt" gives me a list of wrong looking buckets;</div><div><br></div><div><div>Band 1 Block=128x128 Type=Int16, ColorInterp=Undefined</div><div>0...10...20...30...40...50...60...70...80...90...100 - done.</div><div> 256 buckets from 0 to 0:</div><div> 6483600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</div></div><div><br></div><div>(and the resulting tif from the gdal_warp of that also doesn't display properly in qgis).</div><div><br></div><div>For info, the original netcdf gdalinfo is;</div><div><br></div><div><div>Driver: netCDF/Network Common Data Format</div><div>Files: MACC-IFS_forecasts_of_total_column_CO_2018-08-05.nc</div><div> MACC-IFS_forecasts_of_total_column_CO_2018-08-05.nc.aux.xml</div><div>Size is 3600, 1801</div><div>Coordinate System is `'</div><div>Origin = (-0.049999999152054,90.049999999999997)</div><div>Pixel Size = (0.099999998304108,-0.100000000000000)</div><div>Metadata:</div><div> latitude#long_name=latitude</div><div> latitude#units=degrees_north</div><div> longitude#long_name=longitude</div><div> longitude#units=degrees_east</div><div> NC_GLOBAL#Conventions=CF-1.6</div><div> NC_GLOBAL#history=2018-10-16 12:30:21 GMT by grib_to_netcdf-2.8.0: grib_to_netcdf /data/data01/scratch/64/10/_mars-atls14-a82bacafb5c306db76464bc7e824bb75-stvecY.grib -o /data/data03/scratch/20/51/_grib2netcdf-atls13-98f536083ae965b31b0d04811be6f4c6-1XurHN.nc -utime</div><div> NETCDF_DIM_EXTRA={time}</div><div> NETCDF_DIM_time_DEF={1,4}</div><div> NETCDF_DIM_time_VALUES=1039560</div><div> tcco#add_offset=0.005068144666631755</div><div> tcco#long_name=Total column Carbon monoxide</div><div> tcco#missing_value=-32767</div><div> tcco#scale_factor=1.473959627507736e-07</div><div> tcco#units=kg m**-2</div><div> tcco#_FillValue=-32767</div><div> time#calendar=gregorian</div><div> time#long_name=time</div><div> time#units=hours since 1900-01-01 00:00:0.0</div><div>Corner Coordinates:</div><div>Upper Left ( -0.0500000, 90.0500000)</div><div>Lower Left ( -0.0500000, -90.0500000)</div><div>Upper Right ( 359.950, 90.050)</div><div>Lower Right ( 359.950, -90.050)</div><div>Center ( 179.9499969, -0.0000000)</div><div>Band 1 Block=3600x1 Type=Int16, ColorInterp=Undefined</div><div> Min=-32766.000 Max=32767.000</div><div> Minimum=-32766.000, Maximum=32767.000, Mean=-29973.406, StdDev=1279.310</div><div> NoData Value=-32767</div><div> Unit Type: kg m**-2</div><div> Offset: 0.00506814466663176, Scale:1.47395962750774e-07</div><div> Metadata:</div><div> add_offset=0.005068144666631755</div><div> long_name=Total column Carbon monoxide</div><div> missing_value=-32767</div><div> NETCDF_DIM_time=1039560</div><div> NETCDF_VARNAME=tcco</div><div> scale_factor=1.473959627507736e-07</div><div> STATISTICS_MAXIMUM=32767</div><div> STATISTICS_MEAN=-29973.406142881</div><div> STATISTICS_MINIMUM=-32766</div><div> STATISTICS_STDDEV=1279.3102503388</div><div> units=kg m**-2</div><div> _FillValue=-32767</div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, 25 Oct 2018 at 19:11, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Tobias,<br>
<br>
> <br>
> Is there a way to pass the scale / offset information from the original NC<br>
> file to the output of gdal_warp?<br>
<br>
You can 'unscale' priorly with "gdal_translate -unscale <a href="http://in.nc" rel="noreferrer" target="_blank">in.nc</a> out.vrt -of VRT" <br>
and gdalwarp the VRT.<br>
<br>
Even<br>
<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div>