import gdal

import numpy

import pickle

with open("can_be_seen_counts.pickle", "rb") as f:

data = pickle.load(f)

rows = 1201

cols = 1201

base_lat = 40.064531

base_long = -7.462507

with open("can_be_seen_counts.asc", "w") as f:

f.write(f"ncols {cols}\n")

f.write(f"nrows {rows}\n")

f.write(f"xllcorner {base_long}\n")

f.write(f"yllcorner {base_lat}\n")

f.write(f"cellsize {1/1200}\n")

f.write(f"NODATA_value {-1}")

for n, v in enumerate(data.values()):

if n % 1201 == 0:

f.write("\n")

else:

f.write(" ")

f.write(str(v))

import cProfile

from srtm.utilities import haversine

def do_it():

for _ in range(1, 1000000):

haversine(1.23456, 2.23456, 3.23456, 4.23456)

cProfile.run('do_it()', filename='haversine.cprof')

Metadata-Version: 2.1

Name: python-srtm

Version: 0.2.0

Summary:

Home-page: https://github.com/adamcharnock/python-srtm

License: MIT

Keywords: nasa,geospatial,altitude,elevation-profile

Author: Adam Charnock

Author-email: [email protected]

Requires-Python: >=3.8,<4.0

Classifier: License :: OSI Approved :: MIT License

Classifier: Programming Language :: Python :: 3

Classifier: Programming Language :: Python :: 3.8

Project-URL: Repository, https://github.com/adamcharnock/python-srtm

Description-Content-Type: text/markdown

# NASA SRTM altitude data parsing in Python

Provides an API onto SRTM `.hgt` or `.hgt.zip` files.

Requires Python 3.8, **may** work with Python 3.6 & 3.7.

## Installation

```

pip install python-srtm

export SRTM1_DIR=/path/to/srtm1/

export SRTM3_DIR=/path/to/srtm3/

```

## Use

You can access either SRTM1 or SRTM3 data. SRTM 1, for example:

```python

# SRTM1 - 30m resolution

>>> from srtm import Srtm1HeightMapCollection

>>> srtm1_data = Srtm1HeightMapCollection()

>>> srtm1_data.get_altitude(latitude=40.123, longitude=-7.456)

615

>>> Srtm1HeightMapCollection().get_elevation_profile(40.123, -7.456, 40.129, -7.460)

[615, 620, 618, 620, 616, 603, 593, 582, 575, 579, 580, 589, 589, 581, 565, 553, 545, 541, 534, 533, 529, 520, 514]

```

Or SRTM3:

```python

# SRTM3 - 90m resolution

>>> from srtm import Srtm3HeightMapCollection

>>> srtm3_data = Srtm3HeightMapCollection()

>>> srtm3_data.get_altitude(latitude=40.123, longitude=-7.456)

608

>>> Srtm3HeightMapCollection().get_elevation_profile(40.123, -7.456, 40.129, -7.460)

[626, 616, 585, 593, 577, 548, 528, 514]

```

## Profiling

```python

import cProfile

cProfile.run('function_to_profile()', filename='output.cprof')

```

```bash

brew install qcachegrind

pip install pyprof2calltree

pyprof2calltree -k -i /pythonprofiling/profiler/first_iteration.cprof

```

## Release process

For internal reference:

```

# Run the tests

pytest

# Update the setup.py

dephell convert

black setup.py

# Ensure poetry.lock is up to date

poetry lock

export VERSION="VERSION HERE"

# Version bump

poetry version $VERSION

# Commit

git add .

git commit -m "Releasing version $VERSION"

# Tagging and branching

git tag "v$VERSION"

git branch "v$VERSION"

git push origin \

refs/tags/"v$VERSION" \

refs/heads/"v$VERSION" \

master

poetry publish --build

```

Wheel-Version: 1.0

Generator: poetry 1.0.5

Root-Is-Purelib: true

Tag: py3-none-any