MetOffice / MetOffice/aws-earth-examples
Add examples for working with Lambert Azimuthal projection (UKV, MOGREPS-UK)
Nobody has claimed this yet.
- Dominant language
- No language data
- Stars
- 45
- Forks
- 15
- PR merge metrics
- No merged PRs in 30d
Description
Code along the lines of the following can be used to turn the lambert azimuthal x / y co-ordinates into lat/lon. It'd be good to get a full example up to show one way of working with regional UK model data.
```
import pyproj
p = pyproj.Proj("+proj=laea +lon_0=-2.5 +lat_0=54.9 +x_0=0.0 +y_0=0.0 +ellps=GRS80")
# lon, lat = p(x, y, inverse=True)
lon, lat = p(-36000, -158000, inverse=True)
print (lat, lon)
```
53.47926025120268 -3.0422171310489534
_Originally posted by @boristyukin in https://github.com/MetOffice/aws-earth-examples/issues/19#issuecomment-643481826_
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
No target file or test is named. Start by locating the repository's existing dataset examples and review how they document regional model data; use the issue's pyproj Lambert Azimuthal snippet as the starting point. Done means a complete example showing how to work with UKV or MOGREPS-UK coordinates and convert x/y values to latitude/longitude.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- documentation
- Issue type
- Documentation
- Difficulty
- 2/5
- Estimated time
- Half a day
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 35/100