<div dir="ltr"><div>Hi,</div><div><br></div><div>I would like to convert the next NetCDF file to a GeoTIFF WGS 84.</div><br>NetCDF file `gdalinfo` output:<br><br>Driver: netCDF/Network Common Data Format<br>Files: /home/me/<a href="http://netcdf_file.nc">netcdf_file.nc</a><br>Size is 512, 512<br>Metadata:<br>  NC_GLOBAL#altimetry=altiCERSAT<br>  NC_GLOBAL#altitude_resolution=n/a<br>  NC_GLOBAL#area=Ireland 2 min wave grid<br>  NC_GLOBAL#BETAMAX=1.52<br>  NC_GLOBAL#CICE0=0.25<br>  NC_GLOBAL#CICEN=0.75<br>  NC_GLOBAL#comment=<br>  NC_GLOBAL#contact=<a href="mailto:waves@mycompany.com">waves@mycompany.com</a><br>  NC_GLOBAL#Conventions=CF-1.5<br>  NC_GLOBAL#creation_date=2017-07-27T12:20:30Z<br>  NC_GLOBAL#data_type=SWASH straight grid<br>  NC_GLOBAL#easternmost_longitude=-1.0333328<br>  NC_GLOBAL#easting=longitude<br>  NC_GLOBAL#field_type=hourly<br>  NC_GLOBAL#FLAGTR=4<br>  NC_GLOBAL#forcing_ice=iceCFSR<br>  NC_GLOBAL#forcing_wind=wndCFSR<br>  NC_GLOBAL#forecast_range=analysis<br>  NC_GLOBAL#forecast_type=hindcast<br>  NC_GLOBAL#grid=irluk-2m<br>  NC_GLOBAL#grid_projection=n/a<br>  NC_GLOBAL#history=2017-07-27T12:20:30Z : Creation<br>  NC_GLOBAL#institution=My Company<br>  NC_GLOBAL#institution_references=<a href="http://www.mycompany.com/">http://www.mycompany.com/</a><br>  NC_GLOBAL#latitude_resolution=   0.0333333<br>  NC_GLOBAL#longitude_resolution=   0.0333333<br>  NC_GLOBAL#maximum_altitude=9000 m<br>  NC_GLOBAL#minimum_altitude=-12000 m<br>  NC_GLOBAL#netcdf_version=netCDF 4.3.0<br>  NC_GLOBAL#northernmost_latitude=59.5000000<br>  NC_GLOBAL#northing=latitude<br>  NC_GLOBAL#product_name=<a href="http://ww3.irluk-2m.199501_hs.nc">ww3.irluk-2m.199501_hs.nc</a><br>  NC_GLOBAL#product_version=1.0<br>  NC_GLOBAL#REFCOAST=0.05<br>  NC_GLOBAL#references=<a href="http://www.mycompany.com/waves/">http://www.mycompany.com/waves/</a><br>  NC_GLOBAL#REFICEBERG=0.3<br>  NC_GLOBAL#REFSUBGRID=0.05<br>  NC_GLOBAL#run_time=2017-07-27T12:20:30Z<br>  NC_GLOBAL#SDS4 namelist parameter WHITECAPWIDTH=0.30000001<br>  NC_GLOBAL#source=TWM WAVEWATCH III (R) v4.18<br>  NC_GLOBAL#southernmost_latitude=49.5000000<br>  NC_GLOBAL#start_date=1995-01-01T00:00:00Z<br>  NC_GLOBAL#stop_date=1995-01-31T23:00:00Z<br>  NC_GLOBAL#title=TWM SWASH HINDCAST<br>  NC_GLOBAL#WAVEWATCH_III_switches=F90 NOGRB NOPA LRB4 SHRD PR3 UQ FLX0 LN1 ST4 STAB0 NL1 BT4 DB1 MLIM TR0 BS0 IC0 REF1 XX0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 011 NC4 TRKNC IS0<br>  NC_GLOBAL#WAVEWATCH_III_version_number=4.18b<br>  NC_GLOBAL#westernmost_longitude=-12.0000000<br>Subdatasets:<br>  SUBDATASET_1_NAME=NETCDF:"/home/me/<a href="http://netcdf_file.nc">netcdf_file.nc</a>":MAPSTA<br>  SUBDATASET_1_DESC=[301x330] status map (16-bit integer)<br>  SUBDATASET_2_NAME=NETCDF:"/home/me/<a href="http://netcdf_file.nc">netcdf_file.nc</a>":hs<br>  SUBDATASET_2_DESC=[744x301x330] sea_surface_wave_significant_height (32-bit floating-point)<br>Corner Coordinates:<br>Upper Left  (    0.0,    0.0)<br>Lower Left  (    0.0,  512.0)<br>Upper Right (  512.0,    0.0)<br>Lower Right (  512.0,  512.0)<br>Center      (  256.0,  256.0)<br><br><br>I'm expecting as output a file with WGS 84 projection like this: <br><br>`gdalinfo` output:<br><br><br>Driver: GTiff/GeoTIFF<br>Files: tiff_file.tif<br>Size is 4000, 4000<br>Coordinate System is:<br>PROJCRS["WGS 84 / Pseudo-Mercator",<br>    BASEGEOGCRS["WGS 84",<br>        DATUM["World Geodetic System 1984",<br>            ELLIPSOID["WGS 84",6378137,298.257223563,<br>                LENGTHUNIT["metre",1]]],<br>        PRIMEM["Greenwich",0,<br>            ANGLEUNIT["degree",0.0174532925199433]],<br>        ID["EPSG",4326]],<br>    CONVERSION["Popular Visualisation Pseudo-Mercator",<br>        METHOD["Popular Visualisation Pseudo Mercator",<br>            ID["EPSG",1024]],<br>        PARAMETER["Latitude of natural origin",0,<br>            ANGLEUNIT["degree",0.0174532925199433],<br>            ID["EPSG",8801]],<br>        PARAMETER["Longitude of natural origin",0,<br>            ANGLEUNIT["degree",0.0174532925199433],<br>            ID["EPSG",8802]],<br>        PARAMETER["False easting",0,<br>            LENGTHUNIT["metre",1],<br>            ID["EPSG",8806]],<br>        PARAMETER["False northing",0,<br>            LENGTHUNIT["metre",1],<br>            ID["EPSG",8807]]],<br>    CS[Cartesian,2],<br>        AXIS["easting (X)",east,<br>            ORDER[1],<br>            LENGTHUNIT["metre",1]],<br>        AXIS["northing (Y)",north,<br>            ORDER[2],<br>            LENGTHUNIT["metre",1]],<br>    USAGE[<br>        SCOPE["unknown"],<br>        AREA["World - 85°S to 85°N"],<br>        BBOX[-85.06,-180,85.06,180]],<br>    ID["EPSG",3857]]<br>Data axis to CRS axis mapping: 1,2<br>Origin = (-1160000.000000000000000,7799999.999999998137355)<br>Pixel Size = (300.000000000000000,-300.000000000000000)<br>Metadata:<br>  AREA_OR_POINT=Area<br>Image Structure Metadata:<br>  COMPRESSION=LZW<br>  INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left  (-1160000.000, 7800000.000) ( 10d25'13.65"W, 57d11'40.60"N)<br>Lower Left  (-1160000.000, 6600000.000) ( 10d25'13.65"W, 50d52'46.07"N)<br>Upper Right (   40000.000, 7800000.000) (  0d21'33.57"E, 57d11'40.60"N)<br>Lower Right (   40000.000, 6600000.000) (  0d21'33.57"E, 50d52'46.07"N)<br>Center      ( -560000.000, 7200000.000) (  5d 1'50.04"W, 54d 9'26.21"N)<br>Band 1 Block=4000x1 Type=Byte, ColorInterp=Red<br>  Mask Flags: PER_DATASET ALPHA <br>Band 2 Block=4000x1 Type=Byte, ColorInterp=Green<br>  Mask Flags: PER_DATASET ALPHA <br>Band 3 Block=4000x1 Type=Byte, ColorInterp=Blue<br>  Mask Flags: PER_DATASET ALPHA <br>Band 4 Block=4000x1 Type=Byte, ColorInterp=Alpha<br><br><br><div>Any help will be appreciated.</div><div><br></div><div>Javier<br></div></div>