code: def _get_geostationary_reference_longitude(geos_area):
      params = geos_area.crs.coordinate_operation.params
      lon_0_params = [
          p
          for p in params
          if "longitude of natural origin" in p.name.lower()
      ]
      if not lon_0_params:
          return 0
      elif len(lon_0_params) != 1:
          raise ValueError(
              "Not sure how to get reference longitude "
              "information from AreaDefinition."
          )
      return lon_0_params[0].value
  
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 3 for calling a global function with example usage: _get_geostationary_reference_longitude(geos_area) and returns: lon_0_params, int