To follow-up on my latest email, I confirmed that GDAL is interpreting the
pixels in the augmented image file from Octave as uint16 rather than int16,
and this results in incorrect elevation values after cubic interpolation
and subsampling (e.g., a pixel elevation of -28 m is misinterpreted as
65,5
Thanks Joaquim and Even,
Per my email I augmented the original GLOBE image and created a properly
referenced image using gdal_translate as follows: I did an imread from
QGIS/GDAL into Octave, added the appropriate 24 rows and columns on each
side, and did an imwrite back to QGIS/GDAL to continue p
Alternatively, use GMT's grdsample with -fg option to force the knowledge
(if it's not already in the file) that the Earth is round.
http://gmt.soest.hawaii.edu/doc/5.4.2/grdsample.html
Harvey,
I subsampled a GLOBE image of elevation data by a subsampling factor of
12
to a resol
Harvey,
>
> I subsampled a GLOBE image of elevation data by a subsampling factor of 12
> to a resolution of 0.1 deg x 0.1 deg on integer multiples of 0.1 deg. The
> original GLOBE image is: Latitude from +90-1/240 deg to -90+1/240 deg in
> 1/120 deg steps, and Longitude from -180+1/240 deg to +18
Hi,
I subsampled a GLOBE image of elevation data by a subsampling factor of 12
to a resolution of 0.1 deg x 0.1 deg on integer multiples of 0.1 deg. The
original GLOBE image is: Latitude from +90-1/240 deg to -90+1/240 deg in
1/120 deg steps, and Longitude from -180+1/240 deg to +180-1/240 deg in
