Help needed

[sgindeed]

import geemap
import geopandas as gpd
import ee
import pandas as pd
from shapely.geometry import mapping
import geemap.foliumap as geemap
import pandas as pd
import numpy as np

ee.Initialize()

region3_gdf = gpd.read_file("Untitled_layer.geojson")

aoi_total_bounds = region3_gdf.total_bounds
center_lat = (aoi_total_bounds[1] + aoi_total_bounds[3]) / 2
center_lon = (aoi_total_bounds[0] + aoi_total_bounds[2]) / 2

Map = geemap.Map(center=(center_lat, center_lon), zoom=10)
Map.add_gdf(region3_gdf, layer_name="Region AOI", style_dict={'color': 'blue', 'fillColor': 'none'})

geojson = region3_gdf.geo_interface
roi = ee.Geometry(geojson['features'][0]['geometry'])
modis_aod = ee.ImageCollection("MODIS/061/MCD19A2_GRANULES")

start_date = '2023-11-01'
end_date = '2023-11-11'
modisaod_filtered = modis_aod.filterDate(start_date, end_date).filterBounds(roi)

aod_mean_image = modisaod_filtered.mean().clip(roi)
aod_viz_params = {
'min': 0.0,
'max': 1.0,
'palette': ['blue', 'green', 'yellow', 'red']
}
Map.addLayer(aod_mean_image.select('Optical_Depth_047'), aod_viz_params, 'MODIS AOD (Clipped to Polygon)')

grid_generator1k = grids.SquareGridGenerator(1_000)
grid_gdf1k = grid_generator1k.generate_grid(region3_gdf)
Map.add_gdf(grid_gdf1k, layer_name="Grid", style_dict={'color': '#95ddf1', 'fillColor': 'black', 'opacity': 0.5})
Map.addLayerControl()
Map

i want to extract date, time, latitudes, longitudes and Optical_Density_047 values for every cell inside the grid but i am having trouble doing it. please help as it either provides a client side error or query aborts after 5000 elements.