# Load module
from jaxa.earth import je
import numpy as np
import matplotlib.pyplot as plt

#collection_name="JAXA.EORC_GSMaP_standard.Gauge.00Z-23Z.v6_monthly" #for monthly data
#collection_name="JAXA.EORC_GSMaP_standard.Gauge.00Z-23Z.v6_half-monthly" # for half-monthly data
collection_name="JAXA.EORC_GSMaP_standard.Gauge.00Z-23Z.v6_daily"
band="PRECIP"
dlim = ["2021-01-01T00:00:00","2021-01-01T00:00:00"]
ppu  = 4
bbox = [0, -90, 360, 90]


print("-------------Process Start---------------")
# Get an image
data = je.ImageCollection(collection=collection_name,ssl_verify=True)\
         .filter_date(dlim=dlim)\
         .filter_resolution(ppu=ppu)\
         .filter_bounds(bbox=bbox)\
         .select(band=band)\
         .get_images()

# Process and show an image
img = je.ImageProcess(data)\
        .show_images()

# Convert to numpy array
array=data.raster.img[0]

# Print array shape
print("shape=>",array.shape)

# Show image by matplotlib
plt.figure(figsize=(10,10))
plt.imshow(array[:,:,0], cmap="rainbow")
plt.show()

-------------Process Start---------------
 - Collection : JAXA.EORC_GSMaP_standard.Gauge.00Z-23Z.v6_daily
 - Date : 2021-01/01/, 
 - Resolution : 5.0 pixels per 1 degree 
 - Bounds : [0, -90, 360, 90]
 - Band : PRECIP
 - Loading images No.0 : 2021-01/01/
   ------10------20------30------40------50------60------70------80------90-----100
   ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
 - ROI mask : masked
 - Show images : 

showed
shape=> (900, 1800, 1)