wgs84_to_ha.py
def area_of_pixel(pixel_size, center_lat):
    """Calculate m^2 area of a wgs84 square pixel.

    Adapted from: https://gis.stackexchange.com/a/127327/2397

    Args:
        pixel_size (float): length of side of pixel in degrees.
        center_lat (float): latitude of the center of the pixel. Note this
            value +/- half the `pixel-size` must not exceed 90/-90 degrees
            latitude or an invalid area will be calculated.

    Returns:
        Area of square pixel of side length `pixel_size` centered at
        `center_lat` in m^2.

    """
    a = 6378137  # meters
    b = 6356752.3142  # meters
    e = math.sqrt(1 - (b/a)**2)
    area_list = []
    for f in [center_lat+pixel_size/2, center_lat-pixel_size/2]:
        zm = 1 - e*math.sin(math.radians(f))
        zp = 1 + e*math.sin(math.radians(f))
        area_list.append(
            math.pi * b**2 * (
                math.log(zp/zm) / (2*e) +
                math.sin(math.radians(f)) / (zp*zm)))
    return abs(pixel_size / 360. * (area_list[0] - area_list[1]))
Calculate m^2 area of a wgs84 square pixel.

    Adapted from: https://gis.stackexchange.com/a/127327/2397

    Args:
        pixel_size (float): length of side of pixel in degrees.
        center_lat (float): latitude of the center of the pixel. Note this
            value +/- half the `pixel-size` must not exceed 90/-90 degrees
            latitude or an invalid area will be calculated.

    Returns:
        Area of square pixel of side length `pixel_size` centered at
        `center_lat` in m^2.

    
Calculate m^2 area of a wgs84 square pixel.
Calculate m^2 area of a wgs84 square pixel .
def area_of_pixel(pixel_size, center_lat):
    a = 6378137  
    b = 6356752.3142  
    e = math.sqrt(1 - (b/a)**2)
    area_list = []
    for f in [center_lat+pixel_size/2, center_lat-pixel_size/2]:
        zm = 1 - e*math.sin(math.radians(f))
        zp = 1 + e*math.sin(math.radians(f))
        area_list.append(
            math.pi * b**2 * (
                math.log(zp/zm) / (2*e) +
                math.sin(math.radians(f)) / (zp*zm)))
    return abs(pixel_size / 360. * (area_list[0] - area_list[1]))
def area_of_pixel ( pixel_size , center_lat ) : a = 6378137 b = 6356752.3142 e = math . sqrt ( 1 - ( b / a ) ** 2 ) area_list = [ ] for f in [ center_lat + pixel_size / 2 , center_lat - pixel_size / 2 ] : zm = 1 - e * math . sin ( math . radians ( f ) ) zp = 1 + e * math . sin ( math . radians ( f ) ) area_list . append ( math . pi * b ** 2 * ( math . log ( zp / zm ) / ( 2 * e ) + math . sin ( math . radians ( f ) ) / ( zp * zm ) ) ) return abs ( pixel_size / 360. * ( area_list [ 0 ] - area_list [ 1 ] ) )
Calculate m^2 area of a wgs84 square pixel.
Calculate m^2 area of a wgs84 square pixel .
"""Calculate m^2 area of a wgs84 square pixel.

    Adapted from: https://gis.stackexchange.com/a/127327/2397

    Args:
        pixel_size (float): length of side of pixel in degrees.
        center_lat (float): latitude of the center of the pixel. Note this
            value +/- half the `pixel-size` must not exceed 90/-90 degrees
            latitude or an invalid area will be calculated.

    Returns:
        Area of square pixel of side length `pixel_size` centered at
        `center_lat` in m^2.

    """ # meters
 # meters

Area of square pixel of side length `pixel_size` centered at
`center_lat` in m^2.
Area of square pixel of side length ` pixel_size ` centered at ` center_lat ` in m^2 .
pixel_size
length of side of pixel in degrees.
length of side of pixel in degrees .
center_lat
latitude of the center of the pixel. Note this
value +/- half the `pixel-size` must not exceed 90/-90 degrees
latitude or an invalid area will be calculated.
latitude of the center of the pixel . Note this value + / - half the ` pixel - size ` must not exceed 90 / - 90 degrees latitude or an invalid area will be calculated .