Using ContourWriterAGG from pycoast

[B. Schitter]

Hi all,

when I use ContourWriterAGG from pycoast to add coastlines, borders, ... to a with satpy created image, I have to define i.e. the following things:

proj4_string = '+proj=stere +lon_0=8.00 +lat_0=50.00 +lat_ts=50.00 +ellps=WGS84'

area_extent = (-3363403.31,-2291879.85,2630596.69,2203620.1)

area_def = (proj4_string, area_extent)

In satpy I use an area for resample, is it possible to get this values from the area used in satpy so that I don't have to define this values separately?

Best regards, Bernhard

[David Hoese]

Hi Bernhard,

You should be able to pass the area definition object to pycoast directly. Anywhere you pass `(proj4_string, area_extent)` you should be able to pass an AreaDefinition. Note this is the object, not the name of the area that you may have stored in a YAML file of areas. If you want to use an area in a YAML file then there is a way to extract that. However, if you have an image generated from Satpy, is this a geotiff or a PNG? If a geotiff then there are pyresample utility functions for creating an AreaDefinition from a geotiff.

Dave

On 4/13/23 23:49, B. Schitter wrote:
> Hi all,
>
> when I use ContourWriterAGG from pycoast to add coastlines, borders, ... to a with satpy created image, I have to define i.e. the following things:
>

> proj4_string='+proj=stere +lon_0=8.00 +lat_0=50.00 +lat_ts=50.00 +ellps=WGS84'

> area_extent=(-3363403.31,-2291879.85,2630596.69,2203620.1)
> area_def=(proj4_string,area_extent)
>
> In satpy I use an area for resample, is it possible to get this values from the area used in satpy so that I don't have to define this values separately?
>
> Best regards, Bernhard
>

> --
> You received this message because you are subscribed to the Google Groups "pytroll" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to pytroll+u...@googlegroups.com <mailto:pytroll+u...@googlegroups.com>.
> To view this discussion on the web, visit https://groups.google.com/d/msgid/pytroll/cecdb3a7-dbee-4570-b3f7-dd3ee03f5081n%40googlegroups.com <https://groups.google.com/d/msgid/pytroll/cecdb3a7-dbee-4570-b3f7-dd3ee03f5081n%40googlegroups.com?utm_medium=email&utm_source=footer>.

[lobsiger...@gmail.com]

Hi,

when you create the image in satpy you can directly add all the (PyCoast and PyDecorate) overlay stuff in the save dataset call.

Ernst