code: def get_geostationary_mask(area, chunks=None):
      """Compute a mask of the earth's shape as seen by a geostationary satellite.
  
      Args:
          area (pyresample.geometry.AreaDefinition) : Corresponding area
                                                      definition
          chunks (int or tuple): Chunk size for the 2D array that is generated.
  
      Returns:
          Boolean mask, True inside the earth's shape, False outside.
  
      """
      # Compute projection coordinates at the earth's limb
      h = _get_geostationary_height(area)
      xmax, ymax = get_geostationary_angle_extent(area)
      xmax *= h
      ymax *= h
  
      # Compute projection coordinates at the centre of each pixel
      x, y = area.get_proj_coords(chunks=chunks or CHUNK_SIZE)
  
      # Compute mask of the earth's elliptical shape
      return ((x / xmax) ** 2 + (y / ymax) ** 2) <= 1
  
imports: from __future__ import annotations
  import bz2
  import logging
  import os
  import shutil
  import tempfile
  import warnings
  from contextlib import closing, contextmanager
  from io import BytesIO
  from shutil import which
  from subprocess import PIPE, Popen
  import numpy
  import pyproj
  import xarray
  from pyresample.geometry import AreaDefinition
  from satpy import config
  from satpy.readers import FSFile
  from satpy.utils import get_legacy_chunk_size
label_desc: Use image node_id 6 for calling a global function with example usage: get_geostationary_mask(area, chunks) and returns: unknown