Repository: obspy/obspy

Issue:
Cartopy local projection for south hemisphere. 
Hi, 
There is a misbehavior on the local projection option on the South hemisphere using cartopy:

Test code as follow:

```
import matplotlib.pyplot as plt
import obspy
from obspy.clients.fdsn import Client
auspass = Client("http://auspass.edu.au:8080")
station = auspass.get_stations(network='*',station='*',starttime=obspy.UTCDateTime('1950-01-01'),endtime=obspy.UTCDateTime())
station.plot(projection='local', color_per_network=True)
```

**The fix is the following:**
In obspy/imaging/maps.py (line 618 --  620)
```
        proj_kwargs['central_latitude'] = lat_0
        proj_kwargs['central_longitude'] = lon_0
        proj_kwargs['standard_parallels'] = [lat_0,lat_0] #Line to add to have the same behavior than Basemap
        proj = ccrs.AlbersEqualArea(**proj_kwargs)
```


Thanks

Sebastien


Cartopy local projection for south hemisphere. 
Hi, 
There is a misbehavior on the local projection option on the South hemisphere using cartopy:

Test code as follow:

```
import matplotlib.pyplot as plt
import obspy
from obspy.clients.fdsn import Client
auspass = Client("http://auspass.edu.au:8080")
station = auspass.get_stations(network='*',station='*',starttime=obspy.UTCDateTime('1950-01-01'),endtime=obspy.UTCDateTime())
station.plot(projection='local', color_per_network=True)
```

**The fix is the following:**
In obspy/imaging/maps.py (line 618 --  620)
```
        proj_kwargs['central_latitude'] = lat_0
        proj_kwargs['central_longitude'] = lon_0
        proj_kwargs['standard_parallels'] = [lat_0,lat_0] #Line to add to have the same behavior than Basemap
        proj = ccrs.AlbersEqualArea(**proj_kwargs)
```


Thanks

Sebastien




Discussion:
Thanks for the report! @QuLogic Can you have a look at this? I don't really have a lot of experience with cartopy. Thanks!
Side question.. is this australian FDSNWS a long term installment? If yes, you might want to let FDSN know about it (https://www.fdsn.org/webservices/datacenters/) and we could also add it to our list (http://docs.obspy.org/master/packages/obspy.clients.fdsn.html#basic-fdsn-client-usage)..
Hi

The FDSNWS server will be a long term server. We are just starting to import data into it. We will inform FDSN about it in a near future. 

Cheers, 

Sebastien 

> On 17 Aug 2018, at 19:25, Tobias Megies <notifications@github.com> wrote:
> 
> Side question.. is this australian FDSNWS a long term installment? If yes, you might want to let FDSN know about it (https://www.fdsn.org/webservices/datacenters/) and we could also add it to our list (http://docs.obspy.org/master/packages/obspy.clients.fdsn.html#basic-fdsn-client-usage)..
> 
> —
> You are receiving this because you authored the thread.
> Reply to this email directly, view it on GitHub, or mute the thread.

Fix looks good to me on a quick glance..

There's still behavior discrepancy in default of station labels and missing lat/lon ticklabels.. but no idea about that, never touched the cartopy versions of the map code..

## Basemap
![basemap](https://user-images.githubusercontent.com/1842780/44259614-4d03b600-a212-11e8-831c-9dacedfffd26.png)

## Cartopy Before
![cartopy](https://user-images.githubusercontent.com/1842780/44259634-5db42c00-a212-11e8-9440-51439c7d4641.png)

## Cartopy After
![cartopy](https://user-images.githubusercontent.com/1842780/44259610-49702f00-a212-11e8-80d2-2e4234a3705a.png)

Thanks for the report! @QuLogic Can you have a look at this? I don't really have a lot of experience with cartopy. Thanks!
Side question.. is this australian FDSNWS a long term installment? If yes, you might want to let FDSN know about it (https://www.fdsn.org/webservices/datacenters/) and we could also add it to our list (http://docs.obspy.org/master/packages/obspy.clients.fdsn.html#basic-fdsn-client-usage)..
Hi

The FDSNWS server will be a long term server. We are just starting to import data into it. We will inform FDSN about it in a near future. 

Cheers, 

Sebastien 

> On 17 Aug 2018, at 19:25, Tobias Megies <notifications@github.com> wrote:
> 
> Side question.. is this australian FDSNWS a long term installment? If yes, you might want to let FDSN know about it (https://www.fdsn.org/webservices/datacenters/) and we could also add it to our list (http://docs.obspy.org/master/packages/obspy.clients.fdsn.html#basic-fdsn-client-usage)..
> 
> —
> You are receiving this because you authored the thread.
> Reply to this email directly, view it on GitHub, or mute the thread.

Fix looks good to me on a quick glance..

There's still behavior discrepancy in default of station labels and missing lat/lon ticklabels.. but no idea about that, never touched the cartopy versions of the map code..

## Basemap
![basemap](https://user-images.githubusercontent.com/1842780/44259614-4d03b600-a212-11e8-831c-9dacedfffd26.png)

## Cartopy Before
![cartopy](https://user-images.githubusercontent.com/1842780/44259634-5db42c00-a212-11e8-9440-51439c7d4641.png)

## Cartopy After
![cartopy](https://user-images.githubusercontent.com/1842780/44259610-49702f00-a212-11e8-80d2-2e4234a3705a.png)


Resolving patch:
diff --git a/obspy/imaging/maps.py b/obspy/imaging/maps.py
--- a/obspy/imaging/maps.py
+++ b/obspy/imaging/maps.py
@@ -627,6 +627,7 @@ def plot_cartopy(lons, lats, size, color, labels=None, projection='global',
 
         proj_kwargs['central_latitude'] = lat_0
         proj_kwargs['central_longitude'] = lon_0
+        proj_kwargs['standard_parallels'] = [lat_0, lat_0]
         proj = ccrs.AlbersEqualArea(**proj_kwargs)
 
     # User-supplied projection.


Test patch: