Gents - I have exported a .dat file (plain xyz format) from a seismic interpretation package thinking “let’s see if I can plot this feature in 3D with grdview”. Grdview didn’t like “.dat” so I thought let’s do xyz2grd first. It’s a relatively small file (153KB, about 4k rows).
It seems that the increment in the source data is 0.00044966 in some units. I guess that the unit is degree that is 113 km long at the equator. 113 km x 0.00044966 makes about 50 m.
With your -R the size of the output is 2300 x 1480 units. The intended unit is perhaps meter. However, the extents divided by 0.00044966 make a raster that has 5114975 by 3291375 pixels. That’s a big raster.
The problem is that 'e'. It makes xyz2grd think that your limits were given in degrees and that you want increments at 50 meters. Just drop the 'e'. Alternatively, if your my_grid.dat is really already gridded, rename it to my_grid.xyz and let GDAL (indirectly) figure out -R and -I. That is, you should be able to do
Yep, its was the “e”. I wanted to make sure I was getting meters which is why I appended it. xyz2grd worked just fine and then so did grdview.
I tried the other thing you recommended just for fun - the gmt convert function but it didn’t like the fact there was “a change in direction” in the file so it clearly doesn’t think my exported file is a regular grid.