Gibs python use cases

[grallewellyn]

Relevant use cases we identified in our ticket that I implemented in this PR:

• Converting an image to GeoTIFF. Can just elaborated on existing: https://nasa-gibs.github.io/gibs-api-docs/python-usage/#bulk-downloading-of-gibs-imagery
• Rasterio gives access to properties of a geospatial raster file. https://github.com/rasterio/rasterio?tab=readme-ov-file#api-overview
• Overlaying. Useful link geopandas: https://geopandas.org/en/stable/docs/user_guide/set_operations.html   And also: https://gis.stackexchange.com/questions/345917/overlay-two-raster-geotiffs-using-rasterio-mask  
• Masking with a shapefile in Rasterio: https://github.com/rasterio/rasterio/blob/main/docs/topics/masking-by-shapefile.rst
• Reprojection with rasterio: https://github.com/rasterio/rasterio/blob/main/docs/topics/reproject.rst
• Masks rasterio (including no data mask)- https://github.com/rasterio/rasterio/blob/main/docs/topics/masks.rst . Current example in the docs is only for MODIS Terra: https://nasa-gibs.github.io/gibs-api-docs/python-usage/#using-a-no-data-mask DID NOT IMPLEMENT rasterio.read_masks just returns 255 for all values for our data so this function doesn't work for us 
• Use Pillow to create an animation like we do in this notebook
• Include interactive map examples: cartopy
  https://scitools.org.uk/cartopy/docs/v0.15/examples/wmts.html ipyleaflet
  https://ipyleaflet.readthedocs.io/en/latest/layers/wms_layer.html (WAS ALREADY IMPLEMENTED) folium
  https://python-visualization.github.io/folium/latest/user_guide/raster_layers/wms_tile_layer.html

Other notes:

I verified that mapping wmts raster data with cartopy gives the same result as creating the jpg with gdal_translate from the .xml file.

Unfortunately, the animated gif example cannot be rendered in the docs page like the existing ipyleaflet example

I know you can’t see the changes from Git in this PR, but I like ReviewNB for looking at changes to notebooks

@jtroberts I wasn't sure who else to add as a reviewer for this PR

[grallewellyn]

We decided in the GITC standup meeting that Alex and Matt would be the best to review
No rush @grabermtw with AGU this week

[review-notebook-app]

Check out this pull request on 

See visual diffs & provide feedback on Jupyter Notebooks.

---

Powered by ReviewNB

[grallewellyn]

I made the following updates suggested by Joe:

• Should change reproject to instead use EPSG:9822
• Note for ipyleaflet and animated example that it won't render
• swap Visualize WMTS Raster Data By Cartopy to be above Visualize WMTS Raster Data By GDAL
• For overlaying example, use GHRSST_L4_MUR_Sea_Surface_Temperature instead of VIIRS
• Make raster cropping better resolution by increasing size=(1200, 600), for global extents tiff
• "NOTE The available options for the format parameter are: image/png, image/tiff, image/jpeg" should be in a comment instead of a markdown cell

@jtroberts I tried to fix the properties about a tiff file cell with rasterio, but when I download the request to

https://gibs.earthdata.nasa.gov/wms/epsg4326/best/wms.cgi?version=1.3.0&service=WMS&request=GetMap&
format=image/tiff&STYLE=default&bbox=-90,-180,90,180&CRS=EPSG:4326&
HEIGHT=512&WIDTH=512&TIME=2021-11-25&layers=MODIS_Terra_SurfaceReflectance_Bands143

rasterio.open gives me height, width as 512 x 512 which is different than the 163840 x 81920 you said it should be. This is probably because the url requests height and width as 512. Do we still want to include the rasterio properties example?

Also, this is what the new reprojected image looks like:

It is EPSG:3857 to EPSG:9921

[jtroberts]

rasterio.open gives me height, width as 512 x 512 which is different than the 163840 x 81920 you said it should be.

That makes sense given that you specified HEIGHT=512&WIDTH=512. I would change WIDTH=1024 to keep the image in the correct proportion. In that case, the expected output size should be 1024 x 512.

It is EPSG:3857 to EPSG:9921

How about we use EASE Grid instead since it is a more common equal-area global projection? So EPSG:3857 to EPSG:6933.

[grallewellyn]

That makes sense given that you specified HEIGHT=512&WIDTH=512. I would change WIDTH=1024 to keep the image in the correct proportion. In that case, the expected output size should be 1024 x 512.

Okay made that change

Okay, I changed to EPSG:6933 and this is the reprojected image

[jtroberts]

Changes look good to me! My last suggestion is to render that reprojected image in the notebook for reference but not a big deal.

[grallewellyn]

Added showing the reprojected image to the notebook

[grabermtw]

Looks good to me!