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

#South
# collection_name="JAXA.JASMES_GCOM-W.AMSR2_ic0.v201_south_daily"
#North
collection_name="JAXA.JASMES_GCOM-W.AMSR2_ic0.v201_north_daily"
band="IC0"
dlim = ["2021-01-01T00:00:00","2021-01-01T00:00:00"]
ppu  = 4
bbox = [-4388608, -4388608, 4388608, 4388608]


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=(7,7))
plt.imshow(array[:,:,0], cmap="rainbow")
plt.colorbar()
plt.show()

-------------Process Start---------------
 - Collection : JAXA.JASMES_GCOM-W.AMSR2_ic0.v201_north_daily
 - Date : 2021-01/01/, 
 - Resolution : 4.0 pixels per 32768 m 
 - Bounds : [-4388608, -4388608, 4388608, 4388608]
 - Band : IC0
 - Loading images No.0 : 2021-01/01/
   ------10------20------30------40------50------60------70------80------90-----100
   ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
 - ROI mask : masked
 - Show images : 

showed
shape=> (1072, 1072, 1)