2f:Tb22,s = ({ evalscript : ` //VERSION=3 function setup() { return ({ input: [{ bands:["B04", "B03", "B02"], }], output: [{ id: "default", bands: 3, sampleType: SampleType.${image_config.format.out_type} }] }); } const typeFactor = ('${image_config.format.out_type}'=='UINT16') ? 65535 : 255; function stretchType(arr, off, scale) { if (!arr) { return arr;} switch (arr.length) { case 4: arr[3] = typeFactor * (off + arr[3]*scale); case 3: arr[2] = typeFactor * (off + arr[2]*scale); case 2: arr[1] = typeFactor * (off + arr[1]*scale); default: arr[0] = typeFactor * (off + arr[0]*scale); } return arr; } //== PARAMETERS =========================== const c0r = 0.031; const tx = 0.1; const ty0 = 0.1; const ty1 = 1.0; const max = 0.4; const sat = 1.35; const debug = false; // Set to 'true' to highlight out-of-range values const atmRatios = [0.9, 1.6, 2.7]; // Rayleigh-derived consts for automated atmosphere offsets //== FUNCTIONS ============================ const sRGBenc = C => C < 0.0031308 ? (12.92 * C) : (1.055 * Math.pow(C, 0.41666) - 0.055); // atmospheric adjustment const atm2p = (a, c0, c1) => (a - c0) / c1; const atm1p = (a, c0) => atm2p(a, c0, (1 - c0)**2); const atm = (a, ii) => (typeof cManual !== 'undefined') ? (cManual[ii] instanceof Array) ? atm2p(a, cManual[ii][0], cManual[ii][1]) : atm1p(a, cManual[ii]) : atm1p(a, c0r * atmRatios[ii]); //contrast enhancement const adjFun = (a, cx, cy, cmax) => { const ar = a / cmax; const txr = cx / cmax; const bot = (2 * txr - 1) * ar - txr; return ar * (1 + (txr - cy) * (1 - ar) / bot); }; const adj2 = (a, cx, cy0, cy1, max) => { const blend = adjFun(a, 0.5 * max, (max - cx)/max, max); return (1 - blend) * adjFun(a, cx, cy0, max) + blend * adjFun(a, cx, cy1, max); } const adj = a => adj2(a, tx, ty0, ty1, max); const satEnh = rgbArr => { const avg = (rgbArr[0]+rgbArr[1]+rgbArr[2])/3.0; return [ avg * (1 - sat) + rgbArr[0] * sat, avg * (1 - sat) + rgbArr[1] * sat, avg * (1 - sat) + rgbArr[2] * sat ] }; const checkDebug = arr => { if (!debug) { return arr; } const minC = Math.min(arr[0], arr[1], arr[2]); if (minC < 0) { return [ arr[0] < 0 ? 1 : 0, arr[1] < 0 ? 1 : 0, arr[2] < 0 ? 1 : 0 ]; } const maxC = Math.max(arr[0], arr[1], arr[2]); if (maxC > 1) { return [ arr[0] > 1 ? 1 : 0, arr[1] > 1 ? 1 : 0, arr[2] > 1 ? 1 : 0 ]; } return arr; }; //== SCRIPT ============================ function evaluatePixel(sample) { const rgb = satEnh([ adj(atm(sample.B04, 0)), adj(atm(sample.B03, 1)), adj(atm(sample.B02, 2)) ]); const retLin = checkDebug(rgb); return stretchType([ sRGBenc(retLin[0]), sRGBenc(retLin[1]), sRGBenc(retLin[2]) ], 0, 1); } ` })30:[] 2d:{"title":"Sentinel Hub large 16-bit image in UTM","cells":[{"id":0,"value":"md`# Sentinel Hub large 16-bit image in UTM`","pinned":false,"hidden":false,"mode":"ojs"},{"id":64,"value":"viewof login_client_id = Text({type: \"login\", label: \"Client ID:\", placeholder: \"XXX-XXX\", width: 400})","pinned":false,"hidden":false,"mode":"ojs"},{"id":72,"value":"viewof login_client_secret = Text({label: \"Client secret:\", type: \"text\", width: 400})","pinned":false,"hidden":false,"mode":"ojs"},{"id":238,"value":"html`${(logged_in && \"LOGGED_IN\") || \"NOT LOGGED IN\"}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":170,"value":"$2f","pinned":true,"hidden":false,"mode":"ojs"},{"id":101,"value":"image_config = ({\n format: PNG, // PNG (8-bit) or TIF (16-bit)\n date: '2019-11-16',\n // center\n lat: 38.66,\n lon: -9.26,\n // real-world pixel size for the final image\n resolution: 20,\n // size of the final image in pixels \n img_w: 3400,\n img_h: Math.round(3400*594/420),\n img_split: 2500, //max size of the tile\n preview_scale: 8 //e.g. 8 for thumbnail preview. Use 1 for the final download\n});","pinned":true,"hidden":false,"mode":"ojs"},{"id":244,"value":"img = {\n const header = [html`\n \n \n ${tiles_arr[0].map((a,x) => html`${x}`)}\n `\n ];\n \n const rows = tiles_arr.map((row, y) => {\n const cells = row.map((tileURL, x) => {\n return html`\n

`\n });\n \n return html`${y}${cells}`;\n });\n \n return html`${header.concat(rows)}

`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":328,"value":"md`## How to combine tiles\n- Install GDAL, and\n- download all tiles into a directory (click image to start download).\n- In the directory, create a VRT by using the ${'`'}gdalbuildvrt${'`'} command.\n- Then use ${'`gdal_translate`'} to convert a VRT to a single TIF file:\n${'```'}bash\ngdalbuildvrt merged.vrt *.tif\ngdal_translate merged.vrt merged.tif\n${'```'}\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":357,"value":"html`

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":523,"value":"tiles_arr = {\n const res = image_params.resolution;\n const tilesX = Math.ceil(image_params.img_w / image_params.img_split);\n const tilesY = Math.ceil(image_params.img_h / image_params.img_split);\n const rows = [];\n let topYpx=0;\n for (let y = 0; y < tilesY; y++) {\n const cells = [];\n const botYpx = Math.round((y+1)/tilesY * image_params.img_h);\n let leftXpx = 0;\n for (let x = 0; x < tilesX; x++) {\n const rightXpx = Math.round((x+1)/tilesX * image_params.img_w);\n const min_e = image_params.bbox[0] + res * leftXpx;\n const max_n = image_params.bbox[3] - res * topYpx;\n const tile_h = botYpx - topYpx;\n const tileURL = await fetch_tile(min_e, max_n, res, rightXpx - leftXpx, tile_h);\n cells.push(tileURL);\n leftXpx = rightXpx;\n }\n topYpx = botYpx;\n rows.push(cells);\n }\n return rows;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":117,"value":"image_params = {\n const utm_zone = Math.ceil(image_config.lon/6) + 30;\n const utm = proj4(`+proj=utm +zone=${utm_zone} +datum=WGS84 +units=m +no_defs`);\n \n const preview = image_config.preview_scale;\n \n const en = utm.forward([image_config.lon, image_config.lat]);\n const res = image_config.resolution * preview;\n const half_w = Math.round(image_config.img_w / 2 / preview);\n const half_h = Math.round(image_config.img_h / 2 / preview);\n const center_e = Math.round(en[0]/res) * res;\n const center_n = Math.round(en[1]/res) * res;\n return ({\n ...image_config.format,\n ...image_config,\n resolution: res,\n utm_zone: utm_zone,\n img_w: 2 * half_w,\n img_h: 2 * half_h,\n img_split: image_config.img_split / preview,\n center_e: center_e,\n center_n: center_n,\n bbox: [center_e - half_w * res, center_n - half_h * res, center_e + half_w * res, center_n + half_h * res],\n evalscript: s.evalscript\n })\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":92,"value":"function fetch_tile(min_e, max_n, res, w, h) {\n return instance.post('https://services.sentinel-hub.com/api/v1/process',\n {\n \"input\": {\n \"bounds\": {\n \"bbox\": [min_e, max_n - res * h, min_e + res * w, max_n],\n \"properties\": {\n \"crs\": \"http://www.opengis.net/def/crs/EPSG/0/326\"+image_params.utm_zone\n }\n },\n \"data\": [{\n \"type\": \"S2L1C\",\n \"dataFilter\": {\n \"timeRange\": {\n \"from\": image_params.date + \"T00:00:00Z\",\n \"to\": image_params.date + \"T23:59:59Z\"\n }\n }\n }]\n },\n \"output\": {\n \"width\": w,\n \"height\": h,\n \"responses\": [\n {\n \"identifier\": \"default\",\n \"format\": {\n \"type\": image_params.mime_type\n }\n }\n ]\n },\n \"evalscript\": image_params.evalscript\n },\n {\n responseType: 'blob'\n }\n).then(resp => URL.createObjectURL(resp.data))\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":15,"value":"/*instance.post('https://cors-anywhere.herokuapp.com/services.sentinel-hub.com/api/v1/catalog/search',\n {\n \"bbox\": [13,45,14,46],\n \"datetime\": \"2019-12-10T00:00:00Z/2019-12-10T23:59:59Z\",\n \"collections\": [\"sentinel-1-grd\"],\n \"limit\": 5\n })*/","pinned":true,"hidden":false,"mode":"ojs"},{"id":285,"value":"TIF = ({\n out_type: 'UINT16',\n mime_type: 'image/tiff',\n image_suffix: 'tif'\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":287,"value":"PNG = ({\n out_type: 'UINT8',\n mime_type: 'image/png',\n image_suffix: 'png'\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":228,"value":"logged_in = (instance && true) || false;","pinned":false,"hidden":false,"mode":"ojs"},{"id":2,"value":"instance = {\n const instance = axios.create({\n baseURL: \"https://services.sentinel-hub.com\"\n });\n const config = {\n headers: {\n 'Content-Type': 'application/x-www-form-urlencoded;charset=utf-8'\n }\n };\n const body = qs.stringify({\n client_id: login_client_id,\n client_secret: login_client_secret,\n grant_type: \"client_credentials\"\n });\n \n // All requests using this instance will have an access token automatically added\n await instance.post(\"/oauth/token\", body, config).then(resp => {\n Object.assign(instance.defaults, {headers: {authorization: `Bearer ${resp.data.access_token}`}})\n });\n return instance;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":66,"value":"import {Text} from \"@observablehq/inputs\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":34,"value":"axios = require('https://bundle.run/axios@0.21.1')","pinned":false,"hidden":false,"mode":"ojs"},{"id":37,"value":"qs = require('qs@6.10.1/dist/qs.js')","pinned":false,"hidden":false,"mode":"ojs"},{"id":109,"value":"proj4 = require(\"proj4\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W30","stdlib":"1","titling":"auto"}