I can code this, given figures/a formula; but I just flat out don't have the G.I.S knowledge for those figures/formula.

I have to come up with some Python which parses the SRTM elevstion data at 30 arc second resolution, and divide it into (approximately) 1km squares.

I then have to examine each elevation in each square and find the highest point in that square.

I am vague aware that I can't just multiply by a constant to convert - but that's my limit.

How can I divide Europe (the political entity) into roughly km squares using Python?


You need to convert your lat-long bounding box into a bounding box in a projected coordinate system in metres. You then chop this up into 1,000 unit boxes, giving you kilometre squares.

There are pan-Europe data products that use squares like this, so you might just be able to download a raster with that coverage and get the coordinates of each cell - that's a win especially if you might be overlaying your data with anything else on that system. For example the EU CORINE land cover data uses EPSG:3035 or EPSG:3034 (I forget which) for projections to a metric grid system.

There also appears to be a proposal for a Euro grid system based on one of those CRS: https://en.wikipedia.org/wiki/European_grid

For converting between coordinate systems in Python various solutions exist: How to convert projected coordinates to lat/lon using Python?


Your Answer

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

Not the answer you're looking for? Browse other questions tagged or ask your own question.