Pick standard parallels automatically?

[QuLogic]

Currently, we have default values for standard parallels hard-coded in both LambertConformal and AlbersEqualArea.

With basemap though, it automatically picks standard parallels based on central latitude:

        if projection in ['lcc', 'eqdc', 'aea']:
            if projection == 'lcc' and k_0 is not None:
                projparams['k_0']=k_0
            # if lat_0 is given, but not lat_1,
            # set lat_1=lat_0
            if lat_1 is None and lat_0 is not None:
                lat_1 = lat_0
                projparams['lat_1'] = lat_1
            if lat_1 is None or lon_0 is None:
                raise ValueError('must specify lat_1 or lat_0 and lon_0 for %s basemap (lat_2 is optional)' % _projnames[projection])
            if lat_2 is None:
                projparams['lat_2'] = lat_1

Is this perhaps something we should be emulating?

[dopplershift]

That seems reasonable. It might also make people migrating from Basemap happy.

[pelson]

It will definitely make #1023 harder though, so I'm 👍, but not until we have got a good sense of how #1023 would be impacted.