albert-decatur · August 29, 2015 14:19 · albert-decatur · Apr 24, 2015
diff --git a/example get polar sterographic coords from equirectangular b/example get polar sterographic coords from equirectangular
 select st_astext(st_transform(st_geomfromtext('POINT( 170.84215 -80.9338 )',4326),3031));
diff --git a/example get UL from fetched landsat-api json by scene b/example get UL from fetched landsat-api json by scene
 cat * |jq '.results[]|[.upperLeftCornerLongitude,.upperLeftCornerLatitude]|@csv' | sed "s:\"::g;s:,: :g;s:^:'POINT( :g;s:$: )',4326:g" | sed 's:^:ST_GeomFromText( :g;s:$: ):g;s:^:ST_Transform( :g;s:$:,3031 ):g;s:^:ST_AsText( :g;s:$: ):g;s:^:SELECT :g' | sed 's:^:COPY (:g;s:$: ) TO STDOUT;:g' | psql scratch
diff --git a/example WKT from UL landsat-api corner b/example WKT from UL landsat-api corner
 scene=LC80231222014362LGN00; curl -s "https://api.developmentseed.org/landsat?search=${scene}"  |jq '.results[]|[.upperLeftCornerLongitude,.upperLeftCornerLatitude]|@csv' | sed 's:"::g;s:,: :g;s:^:POINT( :g;s:$: ):g'
diff --git a/get antarctic landsat-api geojson from corners b/get antarctic landsat-api geojson from corners
 outdir=out/; a="'";db=scratch; antarctic_scenes=$( echo "copy ( select path,row from wrs2_asc_desc where row >= 120 and row <= 122 ) to stdout with csv header;" | psql $db | sed '1d' ); echo "$antarctic_scenes" | parallel --gnu 'landsat search -p {} -l 1 2>/dev/null | grep sceneID | grep -oE "LC8.*$" | sed "s:\",\s*$::g" | xargs -I :: curl -s "https://api.developmentseed.org/landsat?search=::" | jq '$a'.results[]|[.lowerLeftCornerLongitude,.lowerLeftCornerLatitude,"",.lowerRightCornerLongitude,.lowerRightCornerLatitude,"",.upperRightCornerLongitude,.upperRightCornerLatitude,"",.upperLeftCornerLongitude,.upperLeftCornerLatitude,"",.lowerLeftCornerLongitude,.lowerLeftCornerLatitude]|@csv'$a' | sed '$a's:\\::g;s:"::g;s:\(,\{2\}\):|:g;s:,: :g;s:|:,:g;s:^:POLYGON((:g;s:$:)):g'$a' | wellknown > '$outdir'/$( echo {} | sed "s:,::g" ).geojson'
diff --git a/get one browse url for each pathrow b/get one browse url for each pathrow
 # get wkt points one per line
 db=scratch; antarctic_scenes=$( echo "cop( select path,row from wrs2_asc_desc where row >= 120 and row <= 122 ) to stdout with csv header;" | psql $db | sed '1d' ); echo "$antarctic_scenes" | parallel --gnu 'landsat search -p {} -l 1 2>/dev/null | grep sceneID | grep -oE "LC8.*$" | sed "s:\",\s*$::g" | xargs -I :: curl -s "https://api.developmentseed.org/landsat?search=::" | grep -i browse | grep http | grep -oE "http.*\",$" | sed "s:\",$::g"'

 # get wkt points as single multipoint
 cat * |jq '.results[]|[.upperLeftCornerLongitude,.upperLeftCornerLatitude]|@csv' | sed "s:\"::g;s:,: :g;s:^:'POINT( :g;s:$: )',4326:g" | sed 's:^:ST_GeomFromText( :g;s:$: ):g;s:^:ST_Transform( :g;s:$:,3031 ):g;s:^:ST_AsText( :g;s:$: ):g;s:^:SELECT :g' | sed 's:^:COPY (:g;s:$: ) TO STDOUT;:g' | psql scratch | sed 's:^POINT::g' | tr '\n' ',' | sed 's:,$::g' | sed 's:^:MULTIPOINT (:g;s:$: ):g' | clipboard
	# get wkt points one per line
	db=scratch; antarctic_scenes=$( echo "cop( select path,row from wrs2_asc_desc where row >= 120 and row <= 122 ) to stdout with csv header;" \| psql $db \| sed '1d' ); echo "$antarctic_scenes" \| parallel --gnu 'landsat search -p {} -l 1 2>/dev/null \| grep sceneID \| grep -oE "LC8.$" \| sed "s:\",\s$::g" \| xargs -I :: curl -s "https://api.developmentseed.org/landsat?search=::" \| grep -i browse \| grep http \| grep -oE "http.*\",$" \| sed "s:\",$::g"'

	# get wkt points as single multipoint
	cat * \|jq '.results[]\|[.upperLeftCornerLongitude,.upperLeftCornerLatitude]\|@csv' \| sed "s:\"::g;s:,: :g;s:^:'POINT( :g;s:$: )',4326:g" \| sed 's:^:ST_GeomFromText( :g;s:$: ):g;s:^:ST_Transform( :g;s:$:,3031 ):g;s:^:ST_AsText( :g;s:$: ):g;s:^:SELECT :g' \| sed 's:^:COPY (:g;s:$: ) TO STDOUT;:g' \| psql scratch \| sed 's:^POINT::g' \| tr '\n' ',' \| sed 's:,$::g' \| sed 's:^:MULTIPOINT (:g;s:$: ):g' \| clipboard