# Data discovered via CMR with granule UR SC:NSIDC-0725.004:243900064 and collection concept ID C2386646586-NSIDC_ECS

# Data processed using the following steps:

gdal_merge.py -a_nodata -2e+09 -separate -o {output_dir}/merged.tif {input_dir}/*vx*.tif {input_dir}/*vy*.tif

gdalwarp -tr 1500 1500 {input_dir}/merged.tif {output_dir}/downsampled.tif

gdal2xyz.py -skipnodata -csv -allbands {input_dir}/downsampled.tif {output_dir}/as_xyz.xyz && echo "x,y,vx,vy" > {output_dir}/data_with_header.csv && cat {output_dir}/as_xyz.xyz >> {output_dir}/data_with_header.csv

ogr2ogr -lco ENCODING=UTF-8 -t_srs EPSG:3413 -clipdst {assets_dir}/latitude_shape_40_degrees.geojson -makevalid -oo X_POSSIBLE_NAMES=x -oo Y_POSSIBLE_NAMES=y -s_srs EPSG:3413 -oo KEEP_GEOM_COLUMNS=NO -oo AUTODETECT_TYPE=YES {output_dir}/vectors.gpkg CSV:{input_dir}/data_with_header.csv