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

collection_name="JAXA.EORC_ALOS.PRISM_AW3D30.v3.2_global"
band="DSM"
ppu  = 4
bbox = [0, -90, 360, 90]


print("-------------Process Start---------------")
# Get an image
data = je.ImageCollection(collection=collection_name,ssl_verify=True)\
         .filter_date()\
         .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_ALOS.PRISM_AW3D30.v3.2_global
 - Date : 2021-02/, 
 - Resolution : 5.0 pixels per 1 degree 
 - Bounds : [0, -90, 360, 90]
 - Band : DSM
 - Loading images No.0 : 2021-02/
   ------10------20------30------40------50------60------70------80------90-----100
   ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
 - ROI mask : masked
 - Show images : 

showed
shape=> (900, 1800, 1)