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

collection_name="JAXA.EORC_ALOS-2.PALSAR-2_FNF.v2.1.0_global_yearly"
band="FNF"
dlim = ["2020-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_ALOS-2.PALSAR-2_FNF.v2.1.0_global_yearly
 - Date : 2020/, 
 - Resolution : 5.0 pixels per 1 degree 
 - Bounds : [0, -90, 360, 90]
 - Band : FNF
 - Loading images No.0 : 2020/
   ------10------20------30------40------50------60------70------80------90-----100
   ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
 - ROI mask : masked
 - Show images : 

showed
shape=> (900, 1800, 1)