grdimage error

Discussion:

grdimage error

Felix Bostel

2014-08-12 19:06:40 UTC

Hello,

I want to visualize a slope (format: netCDF) using grdimage.
I get the following error message:
"error 104: execution failed: GMT Fatal Error: grdimage requesting
negative n_items (-64887781) - exceeding 32-bit counting?"
My grid file is a 32-bit file. I also tried to process a smaller part of
the grid using the -R option.

Has anyone an idea what could help to slove this problem?

Best regards,
Felix

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
To unsubscribe, send the message "signoff gmt-help" to ***@lists.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Paul Wessel

2014-08-12 20:58:34 UTC

Permalink

Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.

Cheers, Paul

Sent from my iPad

Post by Felix Bostel
Hello,
I want to visualize a slope (format: netCDF) using grdimage.
"error 104: execution failed: GMT Fatal Error: grdimage requesting negative n_items (-64887781) - exceeding 32-bit counting?"
My grid file is a 32-bit file. I also tried to process a smaller part of the grid using the -R option.
Has anyone an idea what could help to slove this problem?
Best regards,
Felix
Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
To unsubscribe, send the message "signoff gmt-help" to ***@lists.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Felix Bostel

2014-08-13 17:18:08 UTC

Permalink

command line: (grdimage)
-Ju32T/1:50000 -R8.5/9/46.5/47r -K -O
(I'm using gmt via a service which runs on a university server. grdinfo
doesn't run there). Is it possible to get the gird information in other
ways (e.g.: QGIS)? I tried to install gmt locally, but haven't
succeeded, yet.

Cheers, Felix

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Paul Wessel

2014-08-13 17:38:18 UTC

Permalink

Your -R is probably wrong. The trailing "r" means that your area is defined by the lower left (8.5/9) and upper right (46.5/47) corners you specified, but I suspect you instead mean to say

-R8.5/46.5/9/47r

The other selection would requires an enormous image grid and hence the memory message.

-p

Post by Felix Bostel
command line: (grdimage)
-Ju32T/1:50000 -R8.5/9/46.5/47r -K -O
(I'm using gmt via a service which runs on a university server. grdinfo doesn't run there). Is it possible to get the gird information in other ways (e.g.: QGIS)? I tried to install gmt locally, but haven't succeeded, yet.
Cheers, Felix

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Felix Bostel

2014-08-13 19:24:42 UTC

Permalink

Of course! You're totally right. But I still get the same error message...

Post by Paul Wessel
Your -R is probably wrong. The trailing "r" means that your area is defined by the lower left (8.5/9) and upper right (46.5/47) corners you specified, but I suspect you instead mean to say
-R8.5/46.5/9/47r
The other selection would requires an enormous image grid and hence the memory message.
-p

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Paul Wessel

2014-08-13 19:30:23 UTC

Permalink

Seems odd if you can run grdimage on the serve but not grdinfo? I doubt they have disabled some of the programs. Unless you have a web interface that only allows a few things. We also don't know what GMT version you are running so it is still pretty unconstrained. Without info on file area etc I cannot really help.

-p

Post by Felix Bostel
Of course! You're totally right. But I still get the same error message...

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Felix Bostel

2014-08-13 19:52:29 UTC

Permalink

Yes, it is. I just can run grdimage, grdclip, grd2xyz. grdclip works
with my grid file, grdimage and grd2xyz don't.

This is the file information derived from QGIS metadata.

Driver:

GDAL provider

netCDF
Network Common Data Format

Dataset Description

/home/felix/Dokumente/jABC/Test/data/slope_uri_utm.nc

NC_GLOBAL#Conventions=CF-1.0

NC_GLOBAL#Source_Software=Esri ArcGIS

slope_uri_utm.tif#coordinates=x y

slope_uri_utm.tif#esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]

slope_uri_utm.tif#grid_mapping=transverse_mercator

slope_uri_utm.tif#long_name=slope_uri_utm.tif

slope_uri_utm.tif#missing_value=-3.4028235e+38

slope_uri_utm.tif#units=Meter

transverse_mercator#false_easting=500000

transverse_mercator#false_northing=0

transverse_mercator#grid_mapping_name=transverse_mercator

transverse_mercator#latitude_of_projection_origin=0

transverse_mercator#longitude_of_central_meridian=9

transverse_mercator#scale_factor_at_central_meridian=0.9996

x#long_name=x coordinate of projection

x#standard_name=projection_x_coordinate

x#units=Meter

y#long_name=y coordinate of projection

y#standard_name=projection_y_coordinate

y#units=Meter

Band 1

coordinates=x y

esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]

grid_mapping=transverse_mercator

long_name=slope_uri_utm.tif

missing_value=-3.4028235e+38

NETCDF_VARNAME=slope_uri_utm.tif

units=Meter

Dimensions:

X: 1001 Y: 1442 Bands: 1

Origin:

422535,5.2057e+06

Pixel Size:

77.467,-77.467

No Data Value

-3.40282e+38

Data Type:

GDT_Float32 - Thirty two bit floating point

Pyramid overviews:

Layer Spatial Reference System:

+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs

Layer Extent (layer original source projection):

422534.9726139035774395,5093988.6136936303228140 :
500079.4581443184870295,5205696.0543880462646484

Band

Band 1

Band No

1

No Stats

No stats collected yet

Post by Paul Wessel
Seems odd if you can run grdimage on the serve but not grdinfo? I doubt they have disabled some of the programs. Unless you have a web interface that only allows a few things. We also don't know what GMT version you are running so it is still pretty unconstrained. Without info on file area etc I cannot really help.
-p

Post by Felix Bostel
Of course! You're totally right. But I still get the same error message...

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Paul Wessel

2014-08-13 19:56:40 UTC

Permalink

So the grid is already in UTM hence you cannot apply -Ju again. You'll need to see example 28 for how to mix UTM grid and geo annotations. This may perhaps not be possible through the web interface. Basically you would plot the grid using a linear scale -Jx1:50000 and overlay a psbasemap with the -Ju projection.

-p

Yes, it is. I just can run grdimage, grdclip, grd2xyz. grdclip works with my grid file, grdimage and grd2xyz don't.
This is the file information derived from QGIS metadata.
GDAL
provider
netCDF
Network Common Data Format
Dataset Description
/home/felix/Dokumente/jABC/Test/data/slope_uri_utm.nc
NC_GLOBAL#Conventions=CF-1.0
NC_GLOBAL#Source_Software=Esri
ArcGIS
slope_uri_utm.tif#coordinates=x
y
slope_uri_utm.tif#esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]
slope_uri_utm.tif#grid_mapping=transverse_mercator
slope_uri_utm.tif#long_name=slope_uri_utm.tif
slope_uri_utm.tif#missing_value=-3.4028235e+38
slope_uri_utm.tif#units=Meter
transverse_mercator#false_easting=500000
transverse_mercator#false_northing=0
transverse_mercator#grid_mapping_name=transverse_mercator
transverse_mercator#latitude_of_projection_origin=0
transverse_mercator#longitude_of_central_meridian=9
transverse_mercator#scale_factor_at_central_meridian=0.9996
x#long_name=x
coordinate of projection
x#standard_name=projection_x_coordinate
x#units=Meter
y#long_name=y
coordinate of projection
y#standard_name=projection_y_coordinate
y#units=Meter
Band 1
coordinates=x
y
esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]
grid_mapping=transverse_mercator
long_name=slope_uri_utm.tif
missing_value=-3.4028235e+38
NETCDF_VARNAME=slope_uri_utm.tif
units=Meter
1442 Bands: 1
422535,5.2057e+06
77.467,-77.467
No Data Value
-3.40282e+38
GDT_Float32
- Thirty two bit floating point
+proj=utm
+zone=32 +datum=WGS84 +units=m +no_defs
Layer Extent (layer original source
422534.9726139035774395,5093988.6136936303228140
: 500079.4581443184870295,5205696.0543880462646484
Band
Band 1
Band No
1
No Stats
No stats
collected yet

Post by Felix Bostel
Of course! You're totally right. But I still get the same error message...

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue on gmt.soest.hawaii.edu
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.

Felix Bostel

2014-08-13 21:36:46 UTC

Permalink

Thank you very much!
psbasemap with -Ju and grid -Jx1:50000 works perfectly :)

Cheers, Felix

Post by Paul Wessel
So the grid is already in UTM hence you cannot apply -Ju again.
You'll need to see example 28 for how to mix UTM grid and geo
annotations. This may perhaps not be possible through the web
interface. Basically you would plot the grid using a linear scale
-Jx1:50000 and overlay a psbasemap with the -Ju projection.
-p

Post by Felix Bostel
Yes, it is. I just can run grdimage, grdclip, grd2xyz. grdclip works
with my grid file, grdimage and grd2xyz don't.
This is the file information derived from QGIS metadata.
GDAL provider
netCDF
Network Common Data Format
Dataset Description
/home/felix/Dokumente/jABC/Test/data/slope_uri_utm.nc
NC_GLOBAL#Conventions=CF-1.0
NC_GLOBAL#Source_Software=Esri ArcGIS
slope_uri_utm.tif#coordinates=x y
slope_uri_utm.tif#esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]
slope_uri_utm.tif#grid_mapping=transverse_mercator
slope_uri_utm.tif#long_name=slope_uri_utm.tif
slope_uri_utm.tif#missing_value=-3.4028235e+38
slope_uri_utm.tif#units=Meter
transverse_mercator#false_easting=500000
transverse_mercator#false_northing=0
transverse_mercator#grid_mapping_name=transverse_mercator
transverse_mercator#latitude_of_projection_origin=0
transverse_mercator#longitude_of_central_meridian=9
transverse_mercator#scale_factor_at_central_meridian=0.9996
x#long_name=x coordinate of projection
x#standard_name=projection_x_coordinate
x#units=Meter
y#long_name=y coordinate of projection
y#standard_name=projection_y_coordinate
y#units=Meter
Band 1
coordinates=x y
esri_pe_string=PROJCS["WGS_1984_UTM_zone_32N",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",9.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0]]
grid_mapping=transverse_mercator
long_name=slope_uri_utm.tif
missing_value=-3.4028235e+38
NETCDF_VARNAME=slope_uri_utm.tif
units=Meter
X: 1001 Y: 1442 Bands: 1
422535,5.2057e+06
77.467,-77.467
No Data Value
-3.40282e+38
GDT_Float32 - Thirty two bit floating point
+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs
500079.4581443184870295,5205696.0543880462646484
Band
Band 1
Band No
1
No Stats
No stats collected yet

Post by Felix Bostel
Of course! You're totally right. But I still get the same error message...

Post by Paul Wessel
Insufficient info. Please post your grdimage command line as well as output of grdinfo of your grid.
Cheers, Paul
Sent from my iPad

Post by Felix Bostel
Hello,
I want to visualize a slope (format: netCDF) using grdimage.
"error 104: execution failed: GMT Fatal Error: grdimage requesting negative n_items (-64887781) - exceeding 32-bit counting?"
My grid file is a 32-bit file. I also tried to process a smaller part of the grid using the -R option.
Has anyone an idea what could help to slove this problem?
Best regards,
Felix
Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue ongmt.soest.hawaii.edu <http://gmt.soest.hawaii.edu>
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum ongmt.soest.hawaii.edu <http://gmt.soest.hawaii.edu> instead.

Mailing list for GMT discussions of all kinds. If you are not sure you have found a bug, discuss it here first.
To formally report bugs or request features, please register and add New Issue ongmt.soest.hawaii.edu <http://gmt.soest.hawaii.edu>
Note: gmt-help will become obsolete on Sept 1, 2014 - please use forum ongmt.soest.hawaii.edu <http://gmt.soest.hawaii.edu> instead.

Mailing list for GMT discussions of all kinds. If you are not sure
you have found a bug, discuss it here first. To formally report bugs
or request features, please register and add New Issue on
gmt.soest.hawaii.edu <http://gmt.soest.hawaii.edu> To unsubscribe,
obsolete on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu
<http://gmt.soest.hawaii.edu> instead.

Mailing list for GMT discussions of all kinds. If you are not sure you
have found a bug, discuss it here first. To formally report bugs or
request features, please register and add New Issue on
gmt.soest.hawaii.edu To unsubscribe, send the message "signoff
on Sept 1, 2014 - please use forum on gmt.soest.hawaii.edu instead.