Refactoring SpatialRefSysMixin.srs for efficiency and better error handling

[Arnaldo Govene]

The srs property in the SpatialRefSysMixin class has several issues that can be addressed to improve code efficiency, maintainability, and adherence to best practices. The following points highlight the concerns:

1. Recursive Calls:
The property uses return self.srs within the try blocks. This introduces recursion, which could lead to stack overflow in edge cases or unexpected behavior.

Current Code:

try:
    self._srs = gdal.SpatialReference(self.wkt)
    return self.srs  # Recursive call
except Exception as e:
    msg = e

2. Error Handling:
The msg variable is overwritten in the second try block without being utilized after the first block, making error reporting ambiguous. This results in the final error message only reflecting the second exception, losing context from the first one.

Current Code:

try:
    self._srs = gdal.SpatialReference(self.wkt)
    return self.srs
except Exception as e:
    msg = e  # First exception caught

try:
    self._srs = gdal.SpatialReference(self.proj4text)
    return self.srs
except Exception as e:
    msg = e  # Second exception overwrites the first

Proposal:

1. Improve Error Reporting:
Capture and differentiate between errors raised during the initialization from WKT and PROJ.4. Ensure that both error messages are included in the final exception for better clarity and debugging.

2. Raplace Manual Caching with @cached_propert and Remove Recursive Calls:
To further optimize the srs property, replace the manual caching mechanism with @cached_property. This simplifies the code while retaining the benefits of caching, which is particularly important given the computational cost of creating a gdal.SpatialReference object.

Refactored Code:

from django.utils.functional import cached_property

class SpatialRefSysMixin:

    @cached_property
    def srs(self):
        """
        Return a GDAL SpatialReference object.
        """
        try:
            return gdal.SpatialReference(self.wkt)
        except Exception as e:
            wkt_error = f"Error initializing from WKT: {str(e)}"

        try:
            return gdal.SpatialReference(self.proj4text)
        except Exception as e:
            proj4_error = f"Error initializing from PROJ.4: {str(e)}"

        raise Exception(f"Could not get OSR SpatialReference. WKT Error: {wkt_error} | PROJ.4 Error: {proj4_error}")

Rationale for Using cached_property:

• Caches the computationally expensive creation of gdal.SpatialReference, improving efficiency for repeated access.
• Removes manual caching logic, making the code cleaner and easier to maintain.
• Ideal for mostly static spatial reference data, avoiding unnecessary recomputation.

[Sarah Boyce]

Thank you, agree the error overwrite could be improved, the caching logic is also due a refresh