Does someone have information on how can I convert a bathymetric 250m data in Geotiff with utm26 projection, to .grd in latitude and longitude coordinates? I was trying the following two commands, but it does not work.
No. There is no such a thing as GMT grids (well there are those old native binary grids but GDAL knows nothing about them). Just give the output name with a .nc extension to select netCDF. The alternative is use grdproject.
The distinction between convert (we convert formats) and (re)project must be understood.
Some GDAL builds have a “GMT” raster driver that is a recent netCDF. Not all builds have it though. I used to use gdal_translate -of GMT all the time to convert other formats to a GMT netCDF.
gdalinfo --formats | grep GMT
returns nothing for me. There used to be a GMT netCDF driver but that was the old GMT3 netCDF grids that were made up with a vector instead of a matrix. That diver was dropped in the driver cleaning that GDAL performed 1 or 2 years ago.
I used to use the GDAL from MacPorts that did include a “GMT” raster driver. Now MacPorts no longer builds GDAL correctly because the maintainer has to fix something for newer compilers. The gdal that can be installed from conda-forge does not include that driver. I must have been using the older GDAL.
Thanks. I got using a webpage. I am trying now to plot a spatial-temporal graphic of seismicity below the map. I am using the psbasemap, however, the plot its not going out of the map. For example, I used -Y-18.5. It does shows the window. However, if I decrease to -Y-6, the new psbasemap is shown inside the map window. Do you know why? Here is the command of the spatial-temporal plot:
The spatial-temporal plot should be plotting in a separate figure below the map.
It shows that message:
psbasemap [INFORMATION]: Constructing the basemap
psbasemap [INFORMATION]: Linear projection implies x-axis distance exaggeration relative to the y-axis by a factor of 2.02012e-09
psbasemap [INFORMATION]: Map scale is 0.00027 km per cm or 1:27.