Breaking News: Grepper is joining You.com. Read the official announcement!

Now I have 12 raster for 12 month for each rainfall, temperature, tmax and tmin. for the rainfall i want sum all raster to get the annual and for temperature i want to take the average to find the annual average. please give me the correct code to do this

Green Team answered on June 17, 2023 Popularity 1/10 Helpfulness 1/10

answer Now I have 12 raster for 12 month for each rainfall, temperature, tmax and tmin. for the rainfall i want sum all raster to get the annual and for temperature i want to take the average to find the annual average. please give me the correct code to do this

Now I have 12 raster for 12 month for each rainfall, temperature, tmax and tmin. for the rainfall i want sum all raster to get the annual and for temperature i want to take the average to find the annual average. please give me the correct code to do this

Comment

import rasterio
from rasterio import merge
import os
import numpy as np

# Paths to the folders containing the raster files
folders = ["Pre", "Tave", "Tmx", "Tmn"]

# Directory where the result files will be saved
output_directory = "F:/Suitability area/All final data/"

for folder in folders:
    # List to store the data from each raster file
    data_list = []
    # Get a list of all the .tif files in the current folder
    files = [f for f in os.listdir(f"F:/Suitability area/All final data/{folder}") if f.endswith('.tif')]
    
    for file in files:
        # Open the raster file
        with rasterio.open(f"F:/Suitability area/All final data/{folder}/{file}") as src:
            # Read the raster data and add it to the list
            data_list.append(src.read(1))

    # Convert the list to a 3D numpy array
    data = np.array(data_list)

    if folder == "Pre":
        # For rainfall, we want to calculate the sum
        result = np.sum(data, axis=0)
    else:
        # For temperature, we want to calculate the average
        result = np.mean(data, axis=0)

    # Save the result to a new raster file
    with rasterio.open(f"{output_directory}/{folder}_annual.tif", 'w', driver='GTiff', 
                       height=result.shape[0], width=result.shape[1], count=1, dtype=result.dtype, 
                       crs=src.crs, transform=src.transform) as dst:
        dst.write(result, 1)

Popularity 1/10 Helpfulness 1/10 Language python

Source: chat.openai.com

Tags: average each find get python raster sum take temperature this

Link to this answer
Share Copy Link

Contributed on Jun 17 2023

Green Team

0 Answers Avg Quality 2/10

Now I have 12 raster for 12 month for each rainfall, temperature, tmax and tmin. for the rainfall i want sum all raster to get the annual and for temperature i want to take the average to find the annual average. please give me the correct code to do this

Contents

More Related Answers

Now I have 12 raster for 12 month for each rainfall, temperature, tmax and tmin. for the rainfall i want sum all raster to get the annual and for temperature i want to take the average to find the annual average. please give me the correct code to do this

Grepper

Documentation

Social

Legal

Contact

Oops, You will need to install Grepper and log-in to perform this action.