gdalinfo on your data returned below:
Warning 1: Recode from UTF-8 to CP_ACP failed with the error: "Invalid argument".
Driver: netCDF/Network Common Data Format
Size is 5000, 3000
Coordinate System is:
EXTENSION["PROJ4","+proj=geos +lon_0=-75 +h=35786023 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs +sweep=x"]]
Origin = (-3627271.328111211769283,1583173.512170388596132)
Pixel Size = (1002.008657743770527,1002.008657743770527)
Geostationary_Satellite... I honestly don't know anything about this projection... So I tried
gdalwarp giving proj4 string explicitly to your nc (Rad layer) as below:
gdalwarp -of netCDF
-s_srs "+proj=geos +lon_0=-75 +h=35786023 +x_0=0 +y_0=0
+ellps=GRS80 +units=m +no_defs +sweep=x"
-t_srs EPSG:3857 -tr 1002.008657743770527 1002.008657743770527
[Below] Both images are underlain by Bing image (in QGIS...Left image is displayed on On-the-fly projection).
- Following this post we started working on a new scope to create True color image from three GOES-16 images (Band1, 2, and 3).
- Initial idea was to reproject the provided NetCDF files along with all variables (Rad for radiance and DQF as QC value). Then with three reprojected bands David's software should be able to convert them into color enhanced Jpeg image... Using GDAL tool I could reproject one variable (Rad) but not two variables together.
- Instead, I extracted
Rad variable from each band and created RGB composite by SAGA. There was a nice article Bah et al. (2018) how to create "Green" band by mixing red, blue and NIR bands. Following this paper I mixed 45% of Band2 (red), 45% of Band1 (blue) and 10% of Band3 (Nir) by raster calculator. Output composite image is as below: