to_netcdf() seems to work fine, but there is one small problem: when I use grdinfo, it doesn't provide the maximum and minimum values for the z coordinate. Do you know if there is a way to do this?$ grdinfo test.grd
grdinfo: Warning: Parameter TIME_LANGUAGE is deprecated. Use GMT_LANGUAGE instead.
Please see http://gmt.soest.hawaii.edu/doc/5.4.5/GMT_Docs.html#new-features-in-gmt-5 for more information.
test.grd: Title:
test.grd: Command:
test.grd: Remark:
test.grd: Gridline node registration used [Geographic grid]
test.grd: Grid file format: nf = GMT netCDF format (32-bit float), COARDS, CF-1.5
test.grd: x_min: 0 x_max: 359.875 x_inc: 0.125 name: Longitude n_columns: 2880
test.grd: y_min: -89.875 y_max: 90 y_inc: 0.125 name: Latitude n_rows: 1440
test.grd: z_min: 0 z_max: 0 name: Mars Topography
test.grd: scale_factor: 1 add_offset: 0
test.grd: format: classic-L0 option to get GMT to actually look at the data: http://www.generic-mapping-tools.org/gmt/latest/grdinfo.html
str values only in the dict to avoid errors. And xarray doesn't preserve that metadata for many computations on the data because they can't guarantee that the metadata is still correct.
@MarkWieczorek Okay, I think I've determined how you might modify the z_min and z_max values. MWE:
import xarray as xr
ds = xr.open_dataset('foo.grd')
ds.z.attrs['actual_range'] = [-9.0, 9.0]
ds.to_netcdf('bar.grd')The key here is setting the attributes of the z variable (or whatever it's called), not those of the overall Dataset.
gxr = xr.DataArray(grid.to_array(), dims=('Latitude', 'Longitude'), coords={'Latitude': grid.lats(), 'Longitude': grid.lons()}, name='Mars Topography', attrs={'actual_range': [grid.min(), grid.max()]})
gxr.to_netcdf('test.grd')
Hi all! I'm new to pygmt. Is there a way to add titles to plots? I've dug around the other Gitter page and Github, but apologies if I missed this somewhere. Also my axis labels don't seem to show up. For example:
figd.plot(x=lon,y=lat,J="W0/6i",R=region,color=colors,style='c0.1c',frame=['xg60+l"longitude"', 'yf5g30a60+l"latitude"'],cmap="viridis")I'm not seeing "latitude" or "longitude"
frame=['xg60+l"longitude"', 'yf5g30a60+l"latitude"', '+t"This is title"'] to add a title to your plot. GMT cannot add any labels to geographic plots, because the X axis is always longitude and Y axis is always latitude for geographic projections. It's a "feature" of GMT, not a bug. If you really want to add "latitude" or "longitude", you have to use pstext to add them manually. However, the text function hasn't been implemented in pygmt.
data = pygmt.datasets.load_japan_quakes() and I got an error message. I notice that gmt.soest.hawaii.edu is offline at the moment. Does all the data come from that server?`
Hello. I'm wondering if there's any plan to add some require version info on requirement.txt requirement-dev.txt to let dependencies more clarify,
because recently I'd like to trying package pygmt via PKGBUILD for makepkg tools in ArchLinux,
https://github.com/sean0921/pygmt-aur.git
as another python package built by makepkg, in build(), package() procedure, it seems ok,
but when I try to use pygmt.test() in check() procedure,
it encouters some error,
related message: https://gist.github.com/sean0921/66eb9e45220ad20b3e1b0d355f67d199
although it seems like that only Anaconda(Miniconda) Python is supported, I still curious that if it's possible to packaging it in the way of distribution package manager, let it more friendly and accessible for more users and developers. :)
Hi! Very nice and useful initiative. I'm implementing laser altimetry simulations in python and happen to need the GMT grdtrack function to extract elevation data from a topography basemap. Is this function already implemented in pygmt or do you plan to implement it at some point? Thanks!
Hi @steo85it, I'm using grdtrack myself and would be happy to implement it. Already have some code and will submit an issue and pull request which you're more than welcome to contribute to :)
