Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SRTM driver: AREA_OR_POINT=Point #1505

Closed
egouden opened this issue May 3, 2019 · 5 comments
Closed

SRTM driver: AREA_OR_POINT=Point #1505

egouden opened this issue May 3, 2019 · 5 comments

Comments

@egouden
Copy link

egouden commented May 3, 2019

Dear GDAL developers,

I rely heavily on GDAL for developing hydrometeorological products.
Thank you for your amazing work.

After reading the documentation ( here is a summary), it appears that most SRTM height data are based on averaging original satellite radar measurement data and should therefore be interpreted as areas and not points.
Do I understand correctly that the AREA_OR_POINT metadata should be set accordingly to "Area"?
Is this metada used in any of the resampling algorithms?

Best regards,

Edouard Goudenhoofdt

@jratike80
Copy link
Collaborator

jratike80 commented May 3, 2019

Gdal-dev list is the place to ask questions and there you can reach much wider audience. Citation from the page https://github.com/OSGeo/gdal/issues where you created the ticket:

Questions should go to the gdal-dev mailing list (https://lists.osgeo.org/mailman/listinfo/gdal-dev)
or other support forums. GitHub issues are for bug reports and suggestions
for new features.

@rouault
Copy link
Member

rouault commented May 3, 2019

Is this metada used in any of the resampling algorithms?

No, this is informational only.
It appears that AREA_OR_POINT=Point was introduced per 332ecba as a result of https://trac.osgeo.org/gdal/ticket/1884

@egouden
Copy link
Author

egouden commented May 4, 2019

Thank you for your reply.
I have read this ticket before submitting my bug report.
The reporter maintains that the SRTM data should be interpreted as points, which is wrong for the NASA datasets following the documentation.
Apart from that the reader is correct for NASA SRTM 2.1 and 3.0 !

@jratike80
Copy link
Collaborator

Because it has been AREA_OR_POINT=Point for 9 years now and there must be other users who deal with SRTM data I think it should not be changed before discussing about the matter on the mailing list first. If it is only informative metadata for SRTM driver then it should not make any difference, but in some other cases pixel-is-point vs. pixel-is-area leads to half a pixel shifts in georeferencing.

@piyushrpt
Copy link
Contributor

piyushrpt commented May 5, 2019

My understanding is that the AREA_OR_POINT field is meant to convey information about the grid on which the data is laid out rather than the contents themselves. SRTM uses a pixel (center) convention, as opposed to most other DEMs (https://www.usgs.gov/land-resources/eros/topochange/science/srtm-ned-vertical-differencing?qt-science_center_objects=0#qt-science_center_objects) and GDAL's internal (area) model. Other tools also have to use a similar mechanism to handle grid definitions - e.g, grid line registration or pixel registration in GMT.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants