mdsumner’s gists

mdsumner / HLS.md

Created May 16, 2024 21:53

https://cmr.earthdata.nasa.gov/stac/LPCLOUD/search?collections=HLSL30.v2.0,HLSS30.v2.0&bbox=0,0,10,10&datetime=2024-04-16T00:00:00Z/2024-05-14T23:59:59Z

https://git.earthdata.nasa.gov/projects/LPDUR/repos/lpdaac_cloud_data_access/browse/notebooks/Topic_3__Data_Discovery_STAC_CMR-STAC_API.ipynb

mdsumner / python.md

Created May 15, 2024 10:16

Python 3.10.12 (main, Nov 20 2023, 15:14:05) [GCC 11.4.0] on linux Type "help", "copyright", "credits" or "license" for more information.

from osgeo import gdal Traceback (most recent call last): File "", line 1, in ModuleNotFoundError: No module named 'osgeo'

mdsumner / url_extract.md

Last active May 10, 2024 15:33

    set.seed(1)

url = "https://naciscdn.org/naturalearth/10m/raster/SR_LR.zip/SR_LR.tif"
url = paste0("/vsizip/vsicurl/", url)
#url <- "/vsicurl/https://gebco2023.s3.valeria.science/gebco_2023_land_cog.tif"

n = 44000
pts = cbind(x = runif(n, -180, 180), y = runif(n, -90, 90))

  

mdsumner / no-comply.md

Last active May 2, 2024 00:21

ATLAS/ICESat-2 L3B Daily and Monthly Gridded Sea Ice Freeboard

https://nsidc.org/data/atl20/versions/4

https://n5eil01u.ecs.nsidc.org/ATLAS/ATL20.004/2018.10.14/ATL20-01_20181001010615_00370101_004_01.h5

there’s a problem with the grid_mapping, the ‘crs’ declaration is non-compliant and prevents downstream tools picking up the projection and georeferenced metadata.

ncdump -h ATL20-02_20181001010615_00370101_004_01.h5 | grep grid_mapping

mdsumner / SGAT_mask.md

Last active April 30, 2024 11:47

    ## whereever SGAT code does this: 


library(maptools)
data("wrld_simpl")
## Cells per degree
n <- 1
grid <- raster(nrows=n*diff(ylim),ncols=n*diff(xlim),
  xmn=xlim[1],xmx=xlim[2],ymn=ylim[1],ymx=ylim[2],

  

mdsumner / gist:621c3a37720b22edb92b1240c7b2837a

Created April 30, 2024 00:54

    gdalinfo "/vsicurl/https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5" -sd 26
Driver: HDF5Image/HDF5 Dataset
Files: /vsicurl/https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5
Size is 608, 896
Metadata:
  Conventions=CF-1.6
  HDFEOS_INFORMATION_HDFEOSVersion=HDFEOS_5.1.15
  history=This version of the Sea Ice processing code contains updates provided by the science team on September 16, 2019. For details on these updates, see the release notes provided in the DAP.
  institution=NASA's AMSR Science Investigator-led Processing System (SIPS)


  

mdsumner / npole_seaice_25km_h5.md

Last active April 29, 2024 19:34

user guide: https://nsidc.org/sites/default/files/au_si12-v001-userguide_1.pdf (lists the extent for the grid, it's the same as the 25km ones)

https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/

Can read this with your earthdata token set up

header <- "Authorization: Bearer <token>"
u <- "/vsicurl/https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5"
Sys.setenv("GDAL_HTTP_HEADERS" = header)

mdsumner / DM_simulations.md

Created April 26, 2024 17:49

See

> ## Number of chunks and runs per chunk
> chunks <- 250
> runs <- 40
> plan(multisession, workers= availableCores()-1)
Error: Cannot create 127 parallel PSOCK nodes. Each node needs one connection, but there are only 124 connections left out of the maximum 128 available on this R installation

mdsumner / geobox_dsn.md

Created April 19, 2024 00:32

GeoBox from a DSN

def extent(dimension, transform): 
  return [transform[0], transform[0] + dimension[0] * transform[1], 
  transform[3] + dimension[1] * transform[5], transform[3]]
  
  
def bbox(dimension, transform):

mdsumner / parquet.md

Created April 17, 2024 07:57

    library(duckdbfs)
library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 

  

Michael Sumner mdsumner

GeoBox from a DSN