Opened 9 years ago
Closed 5 years ago
#5935 closed defect (wontfix)
GTiff to netcdf conversion fails
Reported by: | aaschwanden | Owned by: | warmerdam |
---|---|---|---|
Priority: | normal | Milestone: | closed_because_of_github_migration |
Component: | default | Version: | 1.11.2 |
Severity: | normal | Keywords: | netcdf, conversion |
Cc: |
Description
Hi,
I've come across a very strange behavior of the gdal netcdf driver when I try to remap a 500MB GTiff to new domain as a netCDF file. See bash script below.
When trying to directly convert the GTiff to a netcdf file of resolution $GRID, the resulting file has missing values on the whole left half of the domain as if only the right half was processed. This happens for both resolutions, 450m and 600m.
Now my first workaround was to first remap the file onto the new domain as a GTiff, then convert the GTiff to netcdf. The intermediate GTiff file looks correct for both resolutions. The resulting netcdf file is correct for 600m resolution, but, again, the right half is shown for 450m.
I've tried different creation options for the netcdf file (nc,nc2,nc4) and two gdal versions: GDAL 1.10.1, released 2013/08/26 and GDAL 1.11.2, released 2015/02/10. Always the same behavior. See attached png file.
I'm quite puzzled. Any ideas would be greatly appreciated.
Andy
gdal is a fantastic tool that I use on a daily basis, keep up the good work.
#!/bin/bash IBCAO=IBCAO_V3_500m_RR wget -nc http://www.ngdc.noaa.gov/mgg/bathymetry/arctic/grids/version3_0/${IBCAO}.tif # Create a buffer that is a multiple of the grid resolution # and works for grid resolutions up to 36km. buffer_x=40650 buffer_y=22000 xmin=$((-638000 - $buffer_x)) ymin=$((-3349600 - $buffer_y)) xmax=$((864700 + $buffer_x)) ymax=$((-657600 + $buffer_y)) # choose netcdf type. Doesn't matter, behavior is the same for all nc/nc2/nc4/nc4c CO='-co "FORMAT=NC"' # First attempt, direct conversion from GTiff to netCDF # This works for neither resolution for GRID in 600 450; do rm -f g${GRID}m_ibcao.* gdalwarp $CO -overwrite -r average -t_srs EPSG:3413 -te $xmin $ymin $xmax $ymax -tr $GRID $GRID -of netCDF ${IBCAO}.tif g${GRID}m_ibcao.nc done # Second attempt, create an intermediate GTiff at desired resolution, # then convert to netCDF. # The intermediate GTiff file looks ok for both resolutions. # This works for 600m but not for 450m for GRID in 600 450; do rm -f g${GRID}m_ibcao-2s.* gdalwarp -overwrite -r average -t_srs EPSG:3413 -te $xmin $ymin $xmax $ymax -tr $GRID $GRID -of GTiff ${IBCAO}.tif g${GRID}m_ibcao-2s.tif gdalwarp $CO -overwrite -of netCDF -s_srs EPSG:3413 -t_srs EPSG:3413 g${GRID}m_ibcao-2s.tif g${GRID}m_ibcao-2s.nc done
Attachments (1)
Change History (2)
by , 9 years ago
Attachment: | g450m_ibcao.png added |
---|
comment:1 by , 5 years ago
Milestone: | → closed_because_of_github_migration |
---|---|
Resolution: | → wontfix |
Status: | new → closed |
This ticket has been automatically closed because Trac is no longer used for GDAL bug tracking, since the project has migrated to GitHub. If you believe this ticket is still valid, you may file it to https://github.com/OSGeo/gdal/issues if it is not already reported there.
converted file. Only right side was processed?