I'm trying to generate a simple profile for some gridded topography data in GMT. However, even though the end of the segment produced by project is clearly in the NaN region of the data set (I usually mask it but included NaN as grey here). When I extract values with grdtrack it exports values for the entire length of the segment. The beginning of the segment also starts in NaN region successfully exports NaN until valid values are reached (second picture). I've double checked that my projections and regions are consistent throughout. Anyone know what might be causing this behaviour

table end table begining profile

gmt begin lake png
    gmt grdimage $DEM $region -JM6i -Qc 
    gmt project -C-66.225/18.330 -A105 -L0/0.9 -G0.001 -Q | gmt grdtrack -G$DEM > $depth
    cut -f3,4 $depth > $depth1
    gmt plot $depth  -W0.1c,red
gmt end

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.