@geoscene/core 4.23.8 → 4.23.9

package/layers/support/imageryRendererUtils.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{clone as e}from"../../core/lang.js";import n from"./RasterFunction.js";import{schema as t}from"./rasterFunctionUtils.js";import{getRFxArgColorRampValue as r,convertColorRampToColormap as o,getColorRampName as a,rgbaConvertTo32Bit as i}from"../../renderers/support/colorRampUtils.js";import{stretchTypeJSONDict as s}from"../../renderers/support/stretchRendererUtils.js";import u from"../../renderers/visualVariables/SizeVariable.js";const c={u1:[0,1],u2:[0,3],u4:[0,15],u8:[0,255],s8:[-128,127],u16:[0,65535],s16:[-32768,32767]},l={simple_scalar:"Simple Scalar",wind_barb:"Wind Barb",single_arrow:"Single Arrow",beaufort_kn:"Beaufort Wind (Knots)",beaufort_m:"Beaufort Wind (MetersPerSecond)",ocean_current_m:"Ocean Current (MetersPerSecond)",ocean_current_kn:"Ocean Current (Knots)"},m=new Set(["raster-stretch","unique-value","class-breaks","raster-shaded-relief","vector-field","raster-colormap"]);function f(e){return m.has(e.type)}function p(n,t){if(!n||!t)return e(n||t);const r=e(n);if(t.rasterFunctionDefinition){const e=t.rasterFunctionDefinition;(e.thumbnail||e.thumbnailEx)&&(e.thumbnail=e.thumbnailEx=null),d(r.rasterFunctionDefinition.arguments,t)}else if("none"!==t.functionName.toLowerCase()){b(r.functionArguments).Raster=t}return r}function d(e,n){for(const t in e)"raster"===t.toLowerCase()&&("RasterFunctionVariable"===e[t].type?(e[t]=n.rasterFunctionDefinition,e[t].type="RasterFunctionTemplate"):"RasterFunctionTemplate"===e[t].type&&d(e[t].arguments,n))}function g(n){const r=e(t[n.functionName+"Function"]),o=n.functionArguments;for(const e in o)"raster"===e.toLowerCase()?(r.arguments[e]=g(o[e]),r.arguments[e].type="RasterFunctionTemplate"):"colormap"===e.toLowerCase()?(r.arguments[e].value=N(o[e]),r.arguments.ColorSchemeType.value=0):r.arguments[e].value=o[e];return r}function h(e,n){switch(n=n||{},e.type){case"raster-stretch":return w(e,n);case"class-breaks":return T(e,n);case"unique-value":return F(e,n);case"raster-colormap":return A(e,n);case"vector-field":return y(e,n);case"raster-shaded-relief":return v(e,n);case"flow":throw new Error("Unsupported rendering rule.")}}function b(e){const n=e.Raster;return n&&"esri.layers.support.RasterFunction"===n.declaredClass?b(n.functionArguments):e}const R={none:0,standardDeviation:3,histogramEqualization:4,minMax:5,percentClip:6,sigmoid:9};function y(e,t){const r=new n;r.functionName="VectorFieldRenderer";const{dataType:o,bandProperties:a}=t,i="vector-uv"===o;let s,c;a&&2===a.length&&(s=a.map((e=>e.BandName.toLowerCase())).indexOf("magnitude"),c=a.map((e=>e.BandName.toLowerCase())).indexOf("direction")),-1!==s&&null!==s||(s=0,c=1);const m="arithmetic"===e.rotationType?1:2,f="flow-from"===e.flowRepresentation?0:1,p=e.visualVariables?e.visualVariables.find((e=>"Magnitude"===e.field)):new u,d={magnitudeBandID:s,directionBandID:c,isUVComponents:i,referenceSystem:m,massFlowAngleRepresentation:f,symbolTileSize:50,symbolTileSizeUnits:100,calculationMethod:"Vector Average",symbologyName:l[e.style.toLowerCase().replace("-","_")],minimumMagnitude:p.minDataValue,maximumMagnitude:p.maxDataValue,minimumSymbolSize:p.minSize,maximumSymbolSize:p.maxSize};return r.functionArguments=d,t.convertToRFT?new n({rasterFunctionDefinition:g(r)}):r}function v(e,t){const a=t.convertToRFT;if("elevation"!==t.dataType)return new n;const i=new n;i.functionName="Hillshade";const s="traditional"===e.hillshadeType?0:1,u="none"===e.scalingType?1:3,c={HillshadeType:s,SlopeType:u,ZFactor:e.zFactor};return 0===s&&(c.Azimuth=e.azimuth,c.Altitude=e.altitude),3===u&&(c.PSPower=e.pixelSizePower,c.PSZFactor=e.pixelSizeFactor),i.functionArguments=c,i.variableName="Raster",e.colorRamp&&(i.functionName="ShadedRelief",a?c.ColorRamp=r(e.colorRamp):c.Colormap=o(e.colorRamp,256)),a?new n({rasterFunctionDefinition:g(i)}):i}function w(e,t){const i=t.convertToRFT,u=new n;u.functionName="Stretch";const c=R[s.toJSON(e.stretchType)],l="u8",m={StretchType:c,Statistics:x(e.statistics),DRA:e.dynamicRangeAdjustment,UseGamma:e.useGamma,Gamma:e.gamma,ComputeGamma:e.computeGamma};if(null!=e.outputMin&&(m.Min=e.outputMin),null!=e.outputMax&&(m.Max=e.outputMax),c===R.standardDeviation?(m.NumberOfStandardDeviations=e.numberOfStandardDeviations,u.outputPixelType=l):c===R.percentClip?(m.MinPercent=e.minPercent,m.MaxPercent=e.maxPercent,u.outputPixelType=l):c===R.minMax?u.outputPixelType=l:c===R.sigmoid&&(m.SigmoidStrengthLevel=e.sigmoidStrengthLevel),u.functionArguments=m,u.variableName="Raster",e.colorRamp){const s=e.colorRamp,c=new n;if(i)c.functionArguments={ColorRamp:r(s)};else{const n=a(s);n?c.functionArguments={colorRamp:n}:!t.convertColorRampToColormap||"algorithmic"!==s.type&&"multipart"!==s.type?c.functionArguments={colorRamp:e.colorRamp.toJSON()}:c.functionArguments={Colormap:o(s,256)}}return c.variableName="Raster",c.functionName="Colormap",c.functionArguments.Raster=u,i?new n({rasterFunctionDefinition:g(c)}):c}return i?new n({rasterFunctionDefinition:g(u)}):u}function T(e,t){const r=[],o=[],a=[],i=[],s=1e-6,{pixelType:u,rasterAttributeTable:c}=t,l=c&&c.features,m=C(c);if(l&&Array.isArray(l)&&e.classBreakInfos){e.classBreakInfos.forEach(((n,t)=>{const r=n.symbol.color;let o;r.a&&l.forEach((a=>{o=a.attributes[e.field],(o>=n.minValue&&o<n.maxValue||t===e.classBreakInfos.length-1&&o>=n.minValue)&&i.push([a.attributes[m],r.r,r.g,r.b])}))}));const r=u?S(i,u):i,o=new n;return o.functionName="Colormap",o.functionArguments={},o.functionArguments.Colormap=r,o.variableName="Raster",t.convertToRFT?new n({rasterFunctionDefinition:g(o)}):o}e.classBreakInfos.forEach(((e,n)=>{const t=e.symbol&&e.symbol.color;t.a?(0===n?r.push(e.minValue,e.maxValue+s):r.push(e.minValue+s,e.maxValue+s),o.push(n),i.push([n,t.r,t.g,t.b])):a.push(e.minValue,e.maxValue)}));const f=u?S(i,u):i,p=new n;p.functionName="Remap",p.functionArguments={InputRanges:r,OutputValues:o,NoDataRanges:a},p.variableName="Raster";const d=new n;return d.functionName="Colormap",d.functionArguments={Colormap:f,Raster:p},t.convertToRFT?new n({rasterFunctionDefinition:g(d)}):d}function S(e,n){const t=c[String(n).toLowerCase()];return t&&e.push([Math.floor(t[0]-1),0,0,0],[Math.ceil(t[1]+1),0,0,0]),e}function C(e){if(!e)return;const{fields:n}=e,t=n&&n.find((e=>e&&e.name&&"value"===e.name.toLowerCase()));return t&&t.name}function F(e,t){var r,o;const a=[],{pixelType:i,rasterAttributeTable:s}=t,u=s&&s.features,c=C(s),l=null==(r=e.defaultSymbol)||null==(o=r.color)?void 0:o.toRgb(),m=e.uniqueValueInfos;if(m)if(u){const n=new Map;m.forEach((e=>{var t;const r=e.value,o=null==(t=e.symbol.color)?void 0:t.toRgb();null!=r&&o&&n.set(String(r),o)}));const t=e.field;u.forEach((({attributes:e})=>{const r=String(e[t]),o=e[c];if(n.has(r)){const e=n.get(r);a.push([o,...e])}else l&&a.push([o,...l])}))}else for(let n=0;n<m.length;n++){const e=m[n],t=e.symbol.color,r=+e.value;if(null!=t&&t.a){if(isNaN(r))return null;a.push([r,t.r,t.g,t.b])}}const f=i&&a.length>0?S(a,i):a,p=new n;return p.functionName="Colormap",p.functionArguments={},p.functionArguments.Colormap=f,p.variableName="Raster",t.convertToRFT?new n({rasterFunctionDefinition:g(p)}):p}function A(e,t){const r=e.extractColormap();if(!r||0===r.length)return;const{pixelType:o}=t,a=o?S(r,o):r,i=new n;return i.functionName="Colormap",i.functionArguments={},i.functionArguments.Colormap=a,t.convertToRFT?new n({rasterFunctionDefinition:g(i)}):i}function x(e){const n=[];return e.forEach((e=>{const t=e;if(Array.isArray(t))n.push(t);else{if(null==t.min||null==t.max)return;const e=[t.min,t.max,t.avg||0,t.stddev||0];n.push(e)}})),n}function N(e){const n=[],t=[];return e.forEach((e=>{n.push(e[0]),t.push(i([...e.slice(1),255]))})),{type:"RasterColormap",values:n,colors:t}}export{p as combineRenderingRules,h as convertRendererToRenderingRule,g as convertRenderingRuleToRFT,f as isSupportedRendererType,c as pixelTypeRanges};
+import{clone as e}from"../../core/lang.js";import n from"./RasterFunction.js";import{schema as t}from"./rasterFunctionUtils.js";import{getRFxArgColorRampValue as r,convertColorRampToColormap as o,getColorRampName as a,rgbaConvertTo32Bit as i}from"../../renderers/support/colorRampUtils.js";import{stretchTypeJSONDict as s}from"../../renderers/support/stretchRendererUtils.js";import u from"../../renderers/visualVariables/SizeVariable.js";const c={u1:[0,1],u2:[0,3],u4:[0,15],u8:[0,255],s8:[-128,127],u16:[0,65535],s16:[-32768,32767]},l={simple_scalar:"Simple Scalar",wind_barb:"Wind Barb",single_arrow:"Single Arrow",beaufort_kn:"Beaufort Wind (Knots)",beaufort_m:"Beaufort Wind (MetersPerSecond)",ocean_current_m:"Ocean Current (MetersPerSecond)",ocean_current_kn:"Ocean Current (Knots)"},m=new Set(["raster-stretch","unique-value","class-breaks","raster-shaded-relief","vector-field","raster-colormap"]);function f(e){return m.has(e.type)}function p(n,t){if(!n||!t)return e(n||t);const r=e(n);if(t.rasterFunctionDefinition){const e=t.rasterFunctionDefinition;(e.thumbnail||e.thumbnailEx)&&(e.thumbnail=e.thumbnailEx=null),d(r.rasterFunctionDefinition.arguments,t)}else if("none"!==t.functionName.toLowerCase()){b(r.functionArguments).Raster=t}return r}function d(e,n){for(const t in e)"raster"===t.toLowerCase()&&("RasterFunctionVariable"===e[t].type?(e[t]=n.rasterFunctionDefinition,e[t].type="RasterFunctionTemplate"):"RasterFunctionTemplate"===e[t].type&&d(e[t].arguments,n))}function g(n){const r=e(t[n.functionName+"Function"]),o=n.functionArguments;for(const e in o)"raster"===e.toLowerCase()?(r.arguments[e]=g(o[e]),r.arguments[e].type="RasterFunctionTemplate"):"colormap"===e.toLowerCase()?(r.arguments[e].value=N(o[e]),r.arguments.ColorSchemeType.value=0):r.arguments[e].value=o[e];return r}function h(e,n){switch(n=n||{},e.type){case"raster-stretch":return w(e,n);case"class-breaks":return T(e,n);case"unique-value":return F(e,n);case"raster-colormap":return A(e,n);case"vector-field":return y(e,n);case"raster-shaded-relief":return v(e,n);case"flow":throw new Error("Unsupported rendering rule.")}}function b(e){const n=e.Raster;return n&&"geoscene.layers.support.RasterFunction"===n.declaredClass?b(n.functionArguments):e}const R={none:0,standardDeviation:3,histogramEqualization:4,minMax:5,percentClip:6,sigmoid:9};function y(e,t){const r=new n;r.functionName="VectorFieldRenderer";const{dataType:o,bandProperties:a}=t,i="vector-uv"===o;let s,c;a&&2===a.length&&(s=a.map((e=>e.BandName.toLowerCase())).indexOf("magnitude"),c=a.map((e=>e.BandName.toLowerCase())).indexOf("direction")),-1!==s&&null!==s||(s=0,c=1);const m="arithmetic"===e.rotationType?1:2,f="flow-from"===e.flowRepresentation?0:1,p=e.visualVariables?e.visualVariables.find((e=>"Magnitude"===e.field)):new u,d={magnitudeBandID:s,directionBandID:c,isUVComponents:i,referenceSystem:m,massFlowAngleRepresentation:f,symbolTileSize:50,symbolTileSizeUnits:100,calculationMethod:"Vector Average",symbologyName:l[e.style.toLowerCase().replace("-","_")],minimumMagnitude:p.minDataValue,maximumMagnitude:p.maxDataValue,minimumSymbolSize:p.minSize,maximumSymbolSize:p.maxSize};return r.functionArguments=d,t.convertToRFT?new n({rasterFunctionDefinition:g(r)}):r}function v(e,t){const a=t.convertToRFT;if("elevation"!==t.dataType)return new n;const i=new n;i.functionName="Hillshade";const s="traditional"===e.hillshadeType?0:1,u="none"===e.scalingType?1:3,c={HillshadeType:s,SlopeType:u,ZFactor:e.zFactor};return 0===s&&(c.Azimuth=e.azimuth,c.Altitude=e.altitude),3===u&&(c.PSPower=e.pixelSizePower,c.PSZFactor=e.pixelSizeFactor),i.functionArguments=c,i.variableName="Raster",e.colorRamp&&(i.functionName="ShadedRelief",a?c.ColorRamp=r(e.colorRamp):c.Colormap=o(e.colorRamp,256)),a?new n({rasterFunctionDefinition:g(i)}):i}function w(e,t){const i=t.convertToRFT,u=new n;u.functionName="Stretch";const c=R[s.toJSON(e.stretchType)],l="u8",m={StretchType:c,Statistics:x(e.statistics),DRA:e.dynamicRangeAdjustment,UseGamma:e.useGamma,Gamma:e.gamma,ComputeGamma:e.computeGamma};if(null!=e.outputMin&&(m.Min=e.outputMin),null!=e.outputMax&&(m.Max=e.outputMax),c===R.standardDeviation?(m.NumberOfStandardDeviations=e.numberOfStandardDeviations,u.outputPixelType=l):c===R.percentClip?(m.MinPercent=e.minPercent,m.MaxPercent=e.maxPercent,u.outputPixelType=l):c===R.minMax?u.outputPixelType=l:c===R.sigmoid&&(m.SigmoidStrengthLevel=e.sigmoidStrengthLevel),u.functionArguments=m,u.variableName="Raster",e.colorRamp){const s=e.colorRamp,c=new n;if(i)c.functionArguments={ColorRamp:r(s)};else{const n=a(s);n?c.functionArguments={colorRamp:n}:!t.convertColorRampToColormap||"algorithmic"!==s.type&&"multipart"!==s.type?c.functionArguments={colorRamp:e.colorRamp.toJSON()}:c.functionArguments={Colormap:o(s,256)}}return c.variableName="Raster",c.functionName="Colormap",c.functionArguments.Raster=u,i?new n({rasterFunctionDefinition:g(c)}):c}return i?new n({rasterFunctionDefinition:g(u)}):u}function T(e,t){const r=[],o=[],a=[],i=[],s=1e-6,{pixelType:u,rasterAttributeTable:c}=t,l=c&&c.features,m=C(c);if(l&&Array.isArray(l)&&e.classBreakInfos){e.classBreakInfos.forEach(((n,t)=>{const r=n.symbol.color;let o;r.a&&l.forEach((a=>{o=a.attributes[e.field],(o>=n.minValue&&o<n.maxValue||t===e.classBreakInfos.length-1&&o>=n.minValue)&&i.push([a.attributes[m],r.r,r.g,r.b])}))}));const r=u?S(i,u):i,o=new n;return o.functionName="Colormap",o.functionArguments={},o.functionArguments.Colormap=r,o.variableName="Raster",t.convertToRFT?new n({rasterFunctionDefinition:g(o)}):o}e.classBreakInfos.forEach(((e,n)=>{const t=e.symbol&&e.symbol.color;t.a?(0===n?r.push(e.minValue,e.maxValue+s):r.push(e.minValue+s,e.maxValue+s),o.push(n),i.push([n,t.r,t.g,t.b])):a.push(e.minValue,e.maxValue)}));const f=u?S(i,u):i,p=new n;p.functionName="Remap",p.functionArguments={InputRanges:r,OutputValues:o,NoDataRanges:a},p.variableName="Raster";const d=new n;return d.functionName="Colormap",d.functionArguments={Colormap:f,Raster:p},t.convertToRFT?new n({rasterFunctionDefinition:g(d)}):d}function S(e,n){const t=c[String(n).toLowerCase()];return t&&e.push([Math.floor(t[0]-1),0,0,0],[Math.ceil(t[1]+1),0,0,0]),e}function C(e){if(!e)return;const{fields:n}=e,t=n&&n.find((e=>e&&e.name&&"value"===e.name.toLowerCase()));return t&&t.name}function F(e,t){var r,o;const a=[],{pixelType:i,rasterAttributeTable:s}=t,u=s&&s.features,c=C(s),l=null==(r=e.defaultSymbol)||null==(o=r.color)?void 0:o.toRgb(),m=e.uniqueValueInfos;if(m)if(u){const n=new Map;m.forEach((e=>{var t;const r=e.value,o=null==(t=e.symbol.color)?void 0:t.toRgb();null!=r&&o&&n.set(String(r),o)}));const t=e.field;u.forEach((({attributes:e})=>{const r=String(e[t]),o=e[c];if(n.has(r)){const e=n.get(r);a.push([o,...e])}else l&&a.push([o,...l])}))}else for(let n=0;n<m.length;n++){const e=m[n],t=e.symbol.color,r=+e.value;if(null!=t&&t.a){if(isNaN(r))return null;a.push([r,t.r,t.g,t.b])}}const f=i&&a.length>0?S(a,i):a,p=new n;return p.functionName="Colormap",p.functionArguments={},p.functionArguments.Colormap=f,p.variableName="Raster",t.convertToRFT?new n({rasterFunctionDefinition:g(p)}):p}function A(e,t){const r=e.extractColormap();if(!r||0===r.length)return;const{pixelType:o}=t,a=o?S(r,o):r,i=new n;return i.functionName="Colormap",i.functionArguments={},i.functionArguments.Colormap=a,t.convertToRFT?new n({rasterFunctionDefinition:g(i)}):i}function x(e){const n=[];return e.forEach((e=>{const t=e;if(Array.isArray(t))n.push(t);else{if(null==t.min||null==t.max)return;const e=[t.min,t.max,t.avg||0,t.stddev||0];n.push(e)}})),n}function N(e){const n=[],t=[];return e.forEach((e=>{n.push(e[0]),t.push(i([...e.slice(1),255]))})),{type:"RasterColormap",values:n,colors:t}}export{p as combineRenderingRules,h as convertRendererToRenderingRule,g as convertRenderingRuleToRFT,f as isSupportedRendererType,c as pixelTypeRanges};

package/layers/support/labelFormatUtils.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import e from"../../core/Error.js";import t from"../../core/Logger.js";import{formatDate as r,convertDateFormatToIntlOptions as a}from"../../intl/date.js";import{formatNumber as n}from"../../intl/number.js";import{isNumericField as l}from"./fieldUtils.js";import{getLabelExpression as o,getSingleFieldArcadeExpression as i}from"./labelUtils.js";import{createLabelExpression as s}from"../../support/arcadeOnDemand.js";const u=t.getLogger("esri.layers.support.labelFormatUtils"),p={type:"simple",evaluate:()=>null},c={getAttribute:(e,t)=>e.field(t)};async function f(t,r,a){if(!t||!t.symbol)return p;const n=t.where,l=o(t),f=n?await import("../../core/sql/WhereClause.js"):null;let d;if("arcade"===l.type){const t=await s(l.expression,a,r);d={type:"arcade",evaluate(r){try{const e=t.evaluate({$feature:"attributes"in r?t.repurposeFeature(r):r});if(null!=e)return e.toString()}catch(a){u.error(new e("bad-arcade-expression","Encountered an error when evaluating label expression for feature",{feature:r,expression:l}))}return null},needsHydrationToEvaluate:()=>null==i(l.expression)}}else d={type:"simple",evaluate:e=>l.expression.replace(/{[^}]*}/g,(t=>{const a=t.slice(1,-1),n=r.get(a);if(!n)return t;let l=null;if("attributes"in e){e&&e.attributes&&(l=e.attributes[n.name])}else l=e.field(n.name);return null==l?"":m(l,n)}))};if(n){let e;try{e=f.WhereClause.create(n,r)}catch(y){return p}const t=d.evaluate;d.evaluate=r=>{const a="attributes"in r?void 0:c;return e.testFeature(r,a)?t(r):null}}return d}function m(e,t){if(null==e)return"";const o=t.domain;if(o)if("codedValue"===o.type||"coded-value"===o.type){const t=e;for(const e of o.codedValues)if(e.code===t)return e.name}else if("range"===o.type){const t=+e,r="range"in o?o.range[0]:o.minValue,a="range"in o?o.range[1]:o.maxValue;if(r<=t&&t<=a)return o.name}let i=e;return"date"===t.type||"esriFieldTypeDate"===t.type?i=r(i,a("short-date")):l(t)&&(i=n(+i)),i||""}export{f as createLabelFunction,m as formatField};
+import e from"../../core/Error.js";import t from"../../core/Logger.js";import{formatDate as r,convertDateFormatToIntlOptions as a}from"../../intl/date.js";import{formatNumber as n}from"../../intl/number.js";import{isNumericField as l}from"./fieldUtils.js";import{getLabelExpression as o,getSingleFieldArcadeExpression as i}from"./labelUtils.js";import{createLabelExpression as s}from"../../support/arcadeOnDemand.js";const u=t.getLogger("geoscene.layers.support.labelFormatUtils"),p={type:"simple",evaluate:()=>null},c={getAttribute:(e,t)=>e.field(t)};async function f(t,r,a){if(!t||!t.symbol)return p;const n=t.where,l=o(t),f=n?await import("../../core/sql/WhereClause.js"):null;let d;if("arcade"===l.type){const t=await s(l.expression,a,r);d={type:"arcade",evaluate(r){try{const e=t.evaluate({$feature:"attributes"in r?t.repurposeFeature(r):r});if(null!=e)return e.toString()}catch(a){u.error(new e("bad-arcade-expression","Encountered an error when evaluating label expression for feature",{feature:r,expression:l}))}return null},needsHydrationToEvaluate:()=>null==i(l.expression)}}else d={type:"simple",evaluate:e=>l.expression.replace(/{[^}]*}/g,(t=>{const a=t.slice(1,-1),n=r.get(a);if(!n)return t;let l=null;if("attributes"in e){e&&e.attributes&&(l=e.attributes[n.name])}else l=e.field(n.name);return null==l?"":m(l,n)}))};if(n){let e;try{e=f.WhereClause.create(n,r)}catch(y){return p}const t=d.evaluate;d.evaluate=r=>{const a="attributes"in r?void 0:c;return e.testFeature(r,a)?t(r):null}}return d}function m(e,t){if(null==e)return"";const o=t.domain;if(o)if("codedValue"===o.type||"coded-value"===o.type){const t=e;for(const e of o.codedValues)if(e.code===t)return e.name}else if("range"===o.type){const t=+e,r="range"in o?o.range[0]:o.minValue,a="range"in o?o.range[1]:o.maxValue;if(r<=t&&t<=a)return o.name}let i=e;return"date"===t.type||"esriFieldTypeDate"===t.type?i=r(i,a("short-date")):l(t)&&(i=n(+i)),i||""}export{f as createLabelFunction,m as formatField};

package/layers/support/layerUtils.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-function e(e){return e&&"base-tile"===e.type||"tile"===e.type||"elevation"===e.type||"imagery-tile"===e.type||"base-elevation"===e.type||"open-street-map"===e.type||"wcs"===e.type||"web-tile"===e.type||"wmts"===e.type||"vector-tile"===e.type}function t(e){return e&&"imagery-tile"===e.type}function n(e){return e&&"wmts"===e.type}function r(e){return e&&"voxel"===e.type}function i(e){return e.parent&&"esri.Basemap"===e.parent.declaredClass&&e.parent.baseLayers.indexOf(e)>-1}function a(e){const t=/[:;]/,n=[];let r=!1;if(e&&(e.forEach((function(e,i){n.push([i,e]),!r&&t.test(e)&&(r=!0)})),n.length>0)){let e;if(r){const t={};n.forEach((e=>{t[e[0]]=e[1]})),e=JSON.stringify(t)}else{const t=[];n.forEach((e=>{t.push(e[0]+":"+e[1])})),e=t.join(";")}return e}return null}function o(e){if(!e)return;const t=[];return e.forEach((function(e,n){t.push('"'+n+'":'+JSON.stringify(e))})),t.length?"{"+t.join(",")+"}":void 0}function l(e,t){const n=[];if(e>0&&t)for(let r=0;r<t.length;r++)if(!(t[r].parentLayerId>=0&&-1===n.indexOf(t[r].parentLayerId)&&t.some((function(e){return e.id===t[r].parentLayerId})))&&t[r].id>=0){let i=!0;const a=t[r].maxScale,o=t[r].minScale;(a>0||o>0)&&(a>0&&o>0?i=a<=e&&e<=o:a>0?i=a<=e:o>0&&(i=e<=o)),i&&n.push(t[r].id)}return n}function u(e){return!0===e.labelsVisible&&null!=e.labelingInfo&&e.labelingInfo.length>0}function c(e){if(e.activeLayer){const t=e.activeLayer.tileMatrixSet;if(t)return t;const n=e.activeLayer.tileMatrixSets;if(n)return n}return null}export{u as areLabelsVisible,l as getLayersForScale,c as getTileMaxtrixSetFromActiveLayer,i as isBaseLayer,t as isReprojectableTiledLayer,e as isTiledLayer,r as isVoxelLayer,n as isWMTSLayer,a as serializeLayerDefinitions,o as serializeTimeOptions};
+function e(e){return e&&"base-tile"===e.type||"tile"===e.type||"elevation"===e.type||"imagery-tile"===e.type||"base-elevation"===e.type||"open-street-map"===e.type||"wcs"===e.type||"web-tile"===e.type||"wmts"===e.type||"vector-tile"===e.type}function t(e){return e&&"imagery-tile"===e.type}function n(e){return e&&"wmts"===e.type}function r(e){return e&&"voxel"===e.type}function i(e){return e.parent&&"geoscene.Basemap"===e.parent.declaredClass&&e.parent.baseLayers.indexOf(e)>-1}function a(e){const t=/[:;]/,n=[];let r=!1;if(e&&(e.forEach((function(e,i){n.push([i,e]),!r&&t.test(e)&&(r=!0)})),n.length>0)){let e;if(r){const t={};n.forEach((e=>{t[e[0]]=e[1]})),e=JSON.stringify(t)}else{const t=[];n.forEach((e=>{t.push(e[0]+":"+e[1])})),e=t.join(";")}return e}return null}function o(e){if(!e)return;const t=[];return e.forEach((function(e,n){t.push('"'+n+'":'+JSON.stringify(e))})),t.length?"{"+t.join(",")+"}":void 0}function l(e,t){const n=[];if(e>0&&t)for(let r=0;r<t.length;r++)if(!(t[r].parentLayerId>=0&&-1===n.indexOf(t[r].parentLayerId)&&t.some((function(e){return e.id===t[r].parentLayerId})))&&t[r].id>=0){let i=!0;const a=t[r].maxScale,o=t[r].minScale;(a>0||o>0)&&(a>0&&o>0?i=a<=e&&e<=o:a>0?i=a<=e:o>0&&(i=e<=o)),i&&n.push(t[r].id)}return n}function u(e){return!0===e.labelsVisible&&null!=e.labelingInfo&&e.labelingInfo.length>0}function c(e){if(e.activeLayer){const t=e.activeLayer.tileMatrixSet;if(t)return t;const n=e.activeLayer.tileMatrixSets;if(n)return n}return null}export{u as areLabelsVisible,l as getLayersForScale,c as getTileMaxtrixSetFromActiveLayer,i as isBaseLayer,t as isReprojectableTiledLayer,e as isTiledLayer,r as isVoxelLayer,n as isWMTSLayer,a as serializeLayerDefinitions,o as serializeTimeOptions};

package/layers/support/rasterDatasets/BaseRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../request.js";import r from"../../../core/Error.js";import{JSONSupport as i}from"../../../core/JSONSupport.js";import o from"../../../core/Logger.js";import{isSome as n,isNone as s,unwrap as a}from"../../../core/maybe.js";import{EsriPromiseMixin as l}from"../../../core/Promise.js";import{onAbort as c}from"../../../core/promiseUtils.js";import{property as f}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as u}from"../../../core/accessorSupport/decorators/subclass.js";import{sanitizeUrl as h}from"../arcgisLayerUrl.js";import{url as m}from"../commonProperties.js";import d from"../RasterStorageInfo.js";import p from"../TileInfo.js";import{getRasterId as x,getBlock as g,putBlock as y,deleteBlock as I,decreaseRefCount as w}from"./RawBlockCache.js";import{decode as R}from"../rasterFormats/RasterCodec.js";import{mosaic as k,approximateTransform as v,getLocalArithmeticNorthRotations as b}from"../rasterFunctions/pixelUtils.js";import{load as T,getWorldWidth as S,projectPoint as B,snapPyramid as M,snapExtent as C,getWorldWrapCount as P,projectExtent as W,getProjectionOffsetGrid as _,getRasterDatasetAlignmentInfo as H,projectResolution as j}from"../rasterFunctions/rasterProjectionHelper.js";import{uvComponentToVector as D,convertToLocalDirections as F}from"../rasterFunctions/vectorFieldUtils.js";import E from"../../../geometry/Extent.js";import L from"../../../geometry/Point.js";const A=8;let z=class extends(l(i)){constructor(){super(...arguments),this.rasterJobHandler=null,this.datasetName=null,this.datasetFormat=null,this.rasterInfo=null,this.ioConfig={sampling:"closest"}}async init(){const e=T();this.addResolvingPromise(e),await this.when()}normalizeCtorArgs(e){return e&&e.ioConfig&&(e={...e,ioConfig:{resolution:null,bandIds:null,sampling:"closest",tileInfo:p.create(),...e.ioConfig}}),e}get _isGlobalWrappableSource(){const{rasterInfo:e}=this,t=S(e.spatialReference);return n(t)&&e.extent.width>=t/2}set url(e){this._set("url",h(e,o.getLogger(this.declaredClass)))}async open(e){throw new r("BaseRaster:open-not-implemented","open() is not implemented")}async fetchTile(e,t,r,i={}){const o=i.tileInfo||this.rasterInfo.storageInfo.tileInfo,n=this.getTileExtentFromTileInfo(e,t,r,o);return this.fetchPixels(n,o.size[0],o.size[1],i)}async identify(e,t={}){t=this._getRequestOptionsWithSliceId(t);const{spatialReference:r,extent:i}=this.rasterInfo,{datumTransformation:o}=t;let a=B(e,r,o);if(!i.intersects(a))return{location:a,value:null};if(n(this.rasterInfo.transform)){const e=this.rasterInfo.transform.inverseTransform(a);if(!this.rasterInfo.nativeExtent.intersects(e))return{location:e,value:null};a=e}let l=0;if(t.srcResolution){l=M(t.srcResolution,this.rasterInfo,this.ioConfig.sampling).pyramidLevel}else if(l=await this.computeBestPyramidLevelForLocation(e,t),null==l)return{location:a,value:null};const c=this.identifyPixelLocation(a,l,null);if(null===c)return{location:a,value:null};const{row:f,col:u,rowOffset:h,colOffset:m}=c,d=x(this.url,t.sliceId),p=`${l}/${f}/${u}`;let I=g(d,null,p);s(I)&&(I=this.fetchRawTile(l,f,u,t),y(d,null,p,I));const w=await I;if(s(w)||!w.pixels||0===w.pixels.length)return{location:a,value:null};const R=h*this.rasterInfo.storageInfo.blockHeight+m,k=!w.mask||w.mask[R]?w.pixels.map((e=>e[R])):null,v=this.rasterInfo.dataType;if(("vector-magdir"===v||"vector-uv"===v)&&(null==k?void 0:k.length)>1){return{location:a,value:k,magdirValue:"vector-magdir"===v?[k[0],k[1]]:D([k[0],k[1]]),pyramidLevel:l}}return{location:a,value:k,pyramidLevel:l}}async fetchPixels(e,t,r,i={}){if(e=C(e),(i=this._getRequestOptionsWithSliceId(i)).requestRawData)return this._fetchPixels(e,t,r,i);const o=S(e.spatialReference),n=P(e);if(s(o)||0===n||1===n&&this._isGlobalWrappableSource)return this._fetchPixels(e,t,r,i);if(n>=3)return{extent:e,pixelBlock:null};const a=[],{xmin:l,xmax:c}=e,f=Math.round(o/(c-l)*t),u=f-Math.round((o/2-l)/(c-l)*t);let h=0;const m=[];for(let s=0;s<=n;s++){const d=new E({xmin:0===s?l:-o/2,xmax:s===n?c-o*s:o/2,ymin:e.ymin,ymax:e.ymax,spatialReference:e.spatialReference}),p=0===s?f-u:s===n?t-h:f;h+=p,m.push(p);const x=i.disableWrapAround&&s>0?null:this._fetchPixels(d,p,r,i);a.push(x)}const d=(await Promise.all(a)).map((e=>null==e?void 0:e.pixelBlock));let p=null;const x={width:t,height:r};if(this.rasterJobHandler){p=(await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:d,srcMosaicSize:x,destDimension:null,coefs:null,sampleSpacing:null,interpolation:"nearest",alignmentInfo:null,blockWidths:m},i)).pixelBlock}else p=k(d,x,{blockWidths:m});return{extent:e,srcExtent:W(e,this.rasterInfo.spatialReference,i.datumTransformation),pixelBlock:p}}async fetchRawPixels(e,t,r,i={}){t={x:Math.floor(t.x),y:Math.floor(t.y)};const o=await this._fetchRawTiles(e,t,r,i),{nativeExtent:s,nativePixelSize:a,storageInfo:l}=this.rasterInfo,c=2**e,f=a.x*c,u=a.y*c,h=new E({xmin:s.xmin+f*t.x,xmax:s.xmin+f*(t.x+r.width-1),ymin:s.ymax-u*(t.y+r.height-1),ymax:s.ymax-u*t.y,spatialReference:s.spatialReference});if(!o)return{extent:h,srcExtent:h,pixelBlock:null};const{pixelBlocks:m,mosaicSize:d}=o;if(1===m.length&&n(m[0])&&m[0].width===r.width&&m[0].height===r.height)return{extent:h,srcExtent:h,pixelBlock:o.pixelBlocks[0]};const p=e>0?l.pyramidBlockWidth:l.blockWidth,x=e>0?l.pyramidBlockHeight:l.blockHeight,g={x:t.x%p,y:t.y%x};let y;if(this.rasterJobHandler){y=(await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:m,srcMosaicSize:d,destDimension:r,clipOffset:g,clipSize:r,coefs:null,sampleSpacing:null,interpolation:i.interpolation,alignmentInfo:null,blockWidths:null},i)).pixelBlock}else y=k(m,d,{clipOffset:g,clipSize:r});return{extent:h,srcExtent:h,pixelBlock:y}}fetchRawTile(e,t,i,o){throw new r("BaseRaster:read-not-implemented","fetchRawTile() is not implemented")}computeExtent(e){return W(this.rasterInfo.extent,e)}decodePixelBlock(e,t){return!this.rasterJobHandler||t.useCanvas?R(e,t):this.rasterJobHandler.decode({data:e,options:t})}async request(e,r,i){var o,n;const{customFetchParameters:s}=this.ioConfig,{range:a,query:l,headers:c}=r;i=null!=(o=null!=(n=i)?n:r.retryCount)?o:this.ioConfig.retryCount;const f=a?{Range:`bytes=${a.from}-${a.to}`}:null;try{return await t(e,{...r,query:{...l,...s},headers:{...c,...f}})}catch(u){if(i>0)return i--,this.request(e,r,i);throw u}}getSliceIndex(e){const{multidimensionalInfo:t}=this.rasterInfo;if(!n(t))return null;if(!n(e)||0===e.length)return null;let r=0;const i=e[0].variableName;for(let o=0;o<t.variables.length;o++){const n=t.variables[o],s=n.dimensions;if(n.name!==i){r+=s.map((e=>this._getDimensionValuesCount(e))).reduce(((e,t)=>e+t));break}const a=s.map((e=>this._getDimensionValuesCount(e))),l=s.length;for(let t=0;t<l;t++){const i=e.find((e=>e.dimensionName===s[t].name));if(null==i)return null;const o=Array.isArray(i.values[0])?i.values[0][0]:i.values[0],n=this._getIndexFromDimensions(o,s[t]);if(-1===n)return null;a.shift(),r+=t===l-1?n:n*a.reduce(((e,t)=>e+t))}}return r}getTileExtentFromTileInfo(e,t,r,i){const o=i.lodAt(e);return this.getTileExtent({x:o.resolution,y:o.resolution},t,r,i.origin,i.spatialReference,i.size)}updateTileInfo(){const{storageInfo:e,spatialReference:t,extent:r,pixelSize:i}=this.rasterInfo;if(!e.tileInfo){const o=[],n=e.maximumPyramidLevel||0;let s=Math.max(i.x,i.y),a=1/.0254*96*s;for(let e=0;e<=n;e++)o.push({level:n-e,resolution:s,scale:a}),s*=2,a*=2;const l=new L({x:r.xmin,y:r.ymax,spatialReference:t});e.tileInfo=new p({origin:l,size:[e.blockWidth,e.blockHeight],spatialReference:t,lods:o}),e.isVirtualTileInfo=!0}}createRemoteDatasetStorageInfo(e,t=512,r=512,i){const{width:o,height:n,nativeExtent:s,pixelSize:a,spatialReference:l}=e,c=new L({x:s.xmin,y:s.ymax,spatialReference:l});null==i&&(i=Math.max(0,Math.round(Math.log(Math.max(o,n))/Math.LN2-8)));const f=this.computeBlockBoundary(s,512,512,{x:s.xmin,y:s.ymax},[a],i);e.storageInfo=new d({blockWidth:t,blockHeight:r,pyramidBlockWidth:t,pyramidBlockHeight:r,origin:c,firstPyramidLevel:1,maximumPyramidLevel:i,blockBoundary:f})}async computeBestPyramidLevelForLocation(e,t={}){return 0}computeBlockBoundary(e,t,r,i,o,n=0,s=2){if(1===o.length&&n>0){o=[...o];let{x:e,y:t}=o[0];for(let r=0;r<n;r++)e*=s,t*=s,o.push({x:e,y:t})}const a=[],{x:l,y:c}=i;for(let f=0;f<o.length;f++){const{x:i,y:n}=o[f];a.push({minCol:Math.floor((e.xmin-l+.1*i)/t/i),maxCol:Math.floor((e.xmax-l-.1*i)/t/i),minRow:Math.floor((c-e.ymax+.1*n)/r/n),maxRow:Math.floor((c-e.ymin-.1*n)/r/n)})}return a}getPyramidPixelSize(e){const{nativePixelSize:t}=this.rasterInfo,{pyramidResolutions:r,pyramidScalingFactor:i}=this.rasterInfo.storageInfo;if(0===e)return t;if(n(r)&&r.length)return r[e-1];const o=i**e;return{x:t.x*o,y:t.y*o}}identifyPixelLocation(e,t,r){const{spatialReference:i,nativeExtent:o}=this.rasterInfo,{blockWidth:n,blockHeight:s,maximumPyramidLevel:a,origin:l}=this.rasterInfo.storageInfo,c=B(e,i,r);if(!o.intersects(c))return null;if(t<0||t>a)return null;const f=this.getPyramidPixelSize(t),{x:u,y:h}=f,m=(l.y-c.y)/h/s,d=(c.x-l.x)/u/n,p=Math.min(s-1,Math.floor((m-Math.floor(m))*s)),x=Math.min(n-1,Math.floor((d-Math.floor(d))*n));return{pyramidLevel:t,row:Math.floor(m),col:Math.floor(d),rowOffset:p,colOffset:x,srcLocation:c}}getTileExtent(e,t,r,i,o,n){const[s,a]=n,l=i.x+r*s*e.x,c=l+s*e.x,f=i.y-t*a*e.y,u=f-a*e.y;return new E({xmin:l,xmax:c,ymin:u,ymax:f,spatialReference:o})}getBlockWidthHeight(e){return{blockWidth:e>0?this.rasterInfo.storageInfo.pyramidBlockWidth:this.rasterInfo.storageInfo.blockWidth,blockHeight:e>0?this.rasterInfo.storageInfo.pyramidBlockHeight:this.rasterInfo.storageInfo.blockHeight}}isBlockOutside(e,t,r){const i=this.rasterInfo.storageInfo.blockBoundary[e];return!i||i.maxRow<t||i.maxCol<r||i.minRow>t||i.minCol>r}async _fetchPixels(e,t,r,i={}){let o=P(e);if(o>=2)return{extent:e,pixelBlock:null};const s=this._getSourceDataInfo(e,t,r,i),{pyramidLevel:l,pyramidResolution:c,srcResolution:f,srcExtent:u,srcWidth:h,srcHeight:m}=s;if(0===h||0===m)return{extent:e,srcExtent:u,pixelBlock:null};const d=a(this.rasterInfo.transform),p="gcs-shift"===(null==d?void 0:d.type),x=n(S(e.spatialReference));!p&&x||(o=P(s.srcExtent,p));const g=this.rasterInfo.storageInfo,y={x:Math.floor((u.xmin-g.origin.x)/c.x+.1),y:Math.floor((g.origin.y-u.ymax)/c.y+.1)},I=await this._fetchRawTiles(l,y,{width:h,height:m,wrapCount:o},i);if(!I)return{extent:e,srcExtent:u,pixelBlock:null};const w=l>0?g.pyramidBlockWidth:g.blockWidth,R=l>0?g.pyramidBlockHeight:g.blockHeight,B=w===h&&R===m&&y.x%w==0&&y.y%R==0,M=new L({x:(e.xmax-e.xmin)/t,y:(e.ymax-e.ymin)/r,spatialReference:e.spatialReference}),C=!e.spatialReference.equals(this.rasterInfo.spatialReference),{datumTransformation:W}=i;if(!C&&B&&1===I.pixelBlocks.length&&w===t&&R===r&&f.x===M.x&&f.y===M.y)return{extent:e,srcExtent:u,pixelBlock:I.pixelBlocks[0]};const H=x&&n(S(u.spatialReference)),j=i.requestProjectedLocalDirections&&this.rasterInfo.dataType.startsWith("vector");j&&!this.rasterJobHandler&&await T();const D=this.rasterJobHandler?await this.rasterJobHandler.getProjectionOffsetGrid({projectedExtent:e,srcBufferExtent:I.extent,pixelSize:M.toJSON(),datumTransformation:W,rasterTransform:d,hasWrapAround:o>0||H,isAdaptive:!1!==this.ioConfig.optimizeProjectionAccuracy,includeGCSGrid:j},i):_({projectedExtent:e,srcBufferExtent:I.extent,pixelSize:M,datumTransformation:W,rasterTransform:d,hasWrapAround:o>0||H,isAdaptive:!1,includeGCSGrid:j});let E;const A=!i.requestRawData,z={rows:D.spacing[0],cols:D.spacing[1]},O=a(this._getRasterTileAlignmentInfo(l,I.extent.xmin)),{pixelBlocks:U,mosaicSize:G,isPartiallyFilled:J}=I;let q=null;if(this.rasterJobHandler){const e=await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:U,srcMosaicSize:G,destDimension:A?{width:t,height:r}:null,coefs:A?D.coefficients:null,sampleSpacing:A?z:null,projectDirections:j,gcsGrid:j?D.gcsGrid:null,isUV:"vector-uv"===this.rasterInfo.dataType,interpolation:i.interpolation,alignmentInfo:O,blockWidths:null},i);({pixelBlock:E,localNorthDirections:q}=e)}else{const e=k(U,G,{alignmentInfo:O});E=A?v(e,{width:t,height:r},D.coefficients,z,i.interpolation):e,j&&D.gcsGrid&&(q=b({width:t,height:r},D.gcsGrid),E=F(E,this.rasterInfo.dataType,q))}return i.requestRawData||j?{srcExtent:u,pixelBlock:E,transformGrid:D,localNorthDirections:q,extent:e,isPartiallyFilled:J}:{srcExtent:u,extent:e,pixelBlock:E}}async _fetchRawTiles(e,t,r,i){const{origin:o,blockBoundary:s}=this.rasterInfo.storageInfo,{blockWidth:a,blockHeight:l}=this.getBlockWidthHeight(e);let{x:c,y:f}=t,{width:u,height:h,wrapCount:m}=r;const d=this._getRasterTileAlignmentInfo(e,0);i.buffer&&(c-=i.buffer.cols,f-=i.buffer.rows,u+=2*i.buffer.cols,h+=2*i.buffer.rows);let p=0,x=0,g=0;if(m&&n(d)){({worldColumnCountFromOrigin:x,originColumnOffset:g,rightPadding:p}=d);x*d.blockWidth-p>=c+u&&(p=0)}const y=Math.floor(c/a),I=Math.floor(f/l),w=Math.floor((c+u+p-1)/a),R=Math.floor((f+h+p-1)/l),k=s[e];if(!k)return null;const{minRow:v,minCol:b,maxCol:T,maxRow:S}=k;if(0===m&&(R<v||w<b||I>S||y>T))return null;const B=new Array;let M=!1;const C=null==this.ioConfig.allowPartialFill?i.allowPartialFill:this.ioConfig.allowPartialFill;for(let F=I;F<=R;F++)for(let t=y;t<=w;t++){let r=t;if(!i.disableWrapAround&&m&&n(d)&&x<=t&&(r=t-x-g),F>=v&&r>=b&&S>=F&&T>=r){const t=this._fetchRawTile(e,F,r,i);C?B.push(new Promise((e=>{t.then((t=>e(t))).catch((()=>{M=!0,e(null)}))}))):B.push(t)}else B.push(null)}if(0===B.length)return null;const P=await Promise.all(B),W={height:(R-I+1)*l,width:(w-y+1)*a},{spatialReference:_}=this.rasterInfo,H=this.getPyramidPixelSize(e),{x:j,y:D}=H;return{extent:new E({xmin:o.x+y*a*j,xmax:o.x+(w+1)*a*j,ymin:o.y-(R+1)*l*D,ymax:o.y-I*l*D,spatialReference:_}),pixelBlocks:P,mosaicSize:W,isPartiallyFilled:M}}_fetchRawTile(e,t,r,i){const o=this.rasterInfo.storageInfo.blockBoundary[e];if(!o)return Promise.resolve(null);const{minRow:n,minCol:a,maxCol:l,maxRow:f}=o;if(t<n||r<a||t>f||r>l)return Promise.resolve(null);const u=x(this.url,i.sliceId),h=`${e}/${t}/${r}`;let m=g(u,i.registryId,h);if(s(m)){const o=new AbortController;m=this.fetchRawTile(e,t,r,{...i,signal:o.signal}),y(u,i.registryId,h,m,o),m.catch((()=>I(u,i.registryId,h)))}return i.signal&&c(i,(()=>{w(u,i.registryId,h)})),m}_getIndexFromDimensions(e,t){const{extent:r,interval:i,unit:o,values:n}=t;if(null!=n&&n.length)return Array.isArray(n[0])?n.findIndex((t=>t[0]<=e&&t[1]>=e)):n.indexOf(e);if(e>r[1])return-1;const s=r[0];let a=-1;if("ISO8601"===o){var l;switch((null==(l=t.intervalUnit)?void 0:l.toLowerCase())||"seconds"){case"seconds":a=Math.round((e-s)/1e3/i);break;case"minutes":a=Math.round((e-s)/6e4/i);break;case"hours":a=Math.round((e-s)/36e5/i);break;case"days":a=Math.round((e-s)/864e5/i);break;case"years":a=Math.round((new Date(e).getUTCFullYear()-new Date(s).getUTCFullYear())/i);break;case"decades":a=Math.round((new Date(e).getUTCFullYear()-new Date(s).getUTCFullYear())/10/i)}return a}return Math.round((e-s)/i)}_getDimensionValuesCount(e){const{extent:t,interval:r,unit:i,values:o}=e;let n=(null==o?void 0:o.length)||0;if(n)return n;const s=t[0];if(0===n&&"ISO8601"===i){var a;switch((null==(a=e.intervalUnit)?void 0:a.toLowerCase())||"seconds"){case"seconds":n=Math.round((t[1]-t[0])/1e3/r);break;case"minutes":n=Math.round((t[1]-t[0])/6e4/r);break;case"hours":n=Math.round((t[1]-t[0])/36e5/r);break;case"days":n=Math.round((t[1]-t[0])/864e5/r);break;case"years":n=Math.round((new Date(t[1]).getUTCFullYear()-new Date(s).getUTCFullYear())/r);break;case"decades":n=Math.round((new Date(t[1]).getUTCFullYear()-new Date(s).getUTCFullYear())/10/r)}return n}return Math.round((t[1]-t[0])/r)}_getRasterTileAlignmentInfo(e,t){return null==this._rasterTileAlighmentInfo&&(this._rasterTileAlighmentInfo=H(this.rasterInfo)),n(this._rasterTileAlighmentInfo.pyramidsInfo)?{startX:t,halfWorldWidth:this._rasterTileAlighmentInfo.halfWorldWidth,hasGCSSShiftTransform:this._rasterTileAlighmentInfo.hasGCSSShiftTransform,...this._rasterTileAlighmentInfo.pyramidsInfo[e]}:null}_getSourceDataInfo(e,t,r,i={}){const o={datumTransformation:i.datumTransformation,pyramidLevel:0,pyramidResolution:null,srcExtent:null,srcHeight:0,srcResolution:null,srcWidth:0};i.srcResolution&&(o.srcResolution=i.srcResolution,this._updateSourceDataInfo(e,o));const n=this.rasterInfo.storageInfo.maximumPyramidLevel||0,{srcWidth:s,srcHeight:a,pyramidLevel:l}=o,c=s/t,f=a/r,u=l<n&&c*f>=16;if(u||l===n&&(c>A||f>A)||(0===s||0===a)){const s=new L({x:(e.xmax-e.xmin)/t,y:(e.ymax-e.ymin)/r,spatialReference:e.spatialReference});let a=j(s,this.rasterInfo.spatialReference,e,o.datumTransformation);const h=!a||i.srcResolution&&a.x+a.y<i.srcResolution.x+i.srcResolution.y;if(u&&i.srcResolution&&h){const e=Math.round(Math.log(Math.max(c,f))/Math.LN2)-1;if(n-l+3>=e){const t=2**e;a={x:i.srcResolution.x*t,y:i.srcResolution.y*t}}}a&&(o.srcResolution=a,this._updateSourceDataInfo(e,o))}return(o.srcWidth/t>A||o.srcHeight/r>A)&&(o.srcWidth=0,o.srcHeight=0),o}_updateSourceDataInfo(e,t){t.srcWidth=0,t.srcHeight=0;const r=this.rasterInfo.spatialReference,{srcResolution:i,datumTransformation:o}=t,{pyramidLevel:n,pyramidResolution:s,excessiveReading:l}=M(i,this.rasterInfo,this.ioConfig.sampling);if(l)return;let c=t.srcExtent||W(e,r,o);if(null==c)return;const f=a(this.rasterInfo.transform);f&&(c=f.inverseTransform(c)),t.srcExtent=c;const u=Math.ceil((c.xmax-c.xmin)/s.x-.1),h=Math.ceil((c.ymax-c.ymin)/s.y-.1);t.pyramidLevel=n,t.pyramidResolution=s,t.srcWidth=u,t.srcHeight=h}_getRequestOptionsWithSliceId(e){return n(this.rasterInfo.multidimensionalInfo)&&null==e.sliceId&&(e={...e,sliceId:this.getSliceIndex(e.multidimensionalDefinition)}),e}};e([f()],z.prototype,"_rasterTileAlighmentInfo",void 0),e([f({readOnly:!0})],z.prototype,"_isGlobalWrappableSource",null),e([f(m)],z.prototype,"url",null),e([f({type:String,json:{write:!0}})],z.prototype,"datasetName",void 0),e([f({type:String,json:{write:!0}})],z.prototype,"datasetFormat",void 0),e([f()],z.prototype,"rasterInfo",void 0),e([f()],z.prototype,"ioConfig",void 0),e([f()],z.prototype,"sourceJSON",void 0),z=e([u("esri.layers.support.rasterDatasets.BaseRaster")],z);const O=z;export{O as default};
+import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../request.js";import r from"../../../core/Error.js";import{JSONSupport as i}from"../../../core/JSONSupport.js";import o from"../../../core/Logger.js";import{isSome as n,isNone as s,unwrap as a}from"../../../core/maybe.js";import{EsriPromiseMixin as l}from"../../../core/Promise.js";import{onAbort as c}from"../../../core/promiseUtils.js";import{property as f}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as u}from"../../../core/accessorSupport/decorators/subclass.js";import{sanitizeUrl as h}from"../arcgisLayerUrl.js";import{url as m}from"../commonProperties.js";import d from"../RasterStorageInfo.js";import p from"../TileInfo.js";import{getRasterId as x,getBlock as g,putBlock as y,deleteBlock as I,decreaseRefCount as w}from"./RawBlockCache.js";import{decode as R}from"../rasterFormats/RasterCodec.js";import{mosaic as k,approximateTransform as v,getLocalArithmeticNorthRotations as b}from"../rasterFunctions/pixelUtils.js";import{load as T,getWorldWidth as S,projectPoint as B,snapPyramid as M,snapExtent as C,getWorldWrapCount as P,projectExtent as W,getProjectionOffsetGrid as _,getRasterDatasetAlignmentInfo as H,projectResolution as j}from"../rasterFunctions/rasterProjectionHelper.js";import{uvComponentToVector as D,convertToLocalDirections as F}from"../rasterFunctions/vectorFieldUtils.js";import E from"../../../geometry/Extent.js";import L from"../../../geometry/Point.js";const A=8;let z=class extends(l(i)){constructor(){super(...arguments),this.rasterJobHandler=null,this.datasetName=null,this.datasetFormat=null,this.rasterInfo=null,this.ioConfig={sampling:"closest"}}async init(){const e=T();this.addResolvingPromise(e),await this.when()}normalizeCtorArgs(e){return e&&e.ioConfig&&(e={...e,ioConfig:{resolution:null,bandIds:null,sampling:"closest",tileInfo:p.create(),...e.ioConfig}}),e}get _isGlobalWrappableSource(){const{rasterInfo:e}=this,t=S(e.spatialReference);return n(t)&&e.extent.width>=t/2}set url(e){this._set("url",h(e,o.getLogger(this.declaredClass)))}async open(e){throw new r("BaseRaster:open-not-implemented","open() is not implemented")}async fetchTile(e,t,r,i={}){const o=i.tileInfo||this.rasterInfo.storageInfo.tileInfo,n=this.getTileExtentFromTileInfo(e,t,r,o);return this.fetchPixels(n,o.size[0],o.size[1],i)}async identify(e,t={}){t=this._getRequestOptionsWithSliceId(t);const{spatialReference:r,extent:i}=this.rasterInfo,{datumTransformation:o}=t;let a=B(e,r,o);if(!i.intersects(a))return{location:a,value:null};if(n(this.rasterInfo.transform)){const e=this.rasterInfo.transform.inverseTransform(a);if(!this.rasterInfo.nativeExtent.intersects(e))return{location:e,value:null};a=e}let l=0;if(t.srcResolution){l=M(t.srcResolution,this.rasterInfo,this.ioConfig.sampling).pyramidLevel}else if(l=await this.computeBestPyramidLevelForLocation(e,t),null==l)return{location:a,value:null};const c=this.identifyPixelLocation(a,l,null);if(null===c)return{location:a,value:null};const{row:f,col:u,rowOffset:h,colOffset:m}=c,d=x(this.url,t.sliceId),p=`${l}/${f}/${u}`;let I=g(d,null,p);s(I)&&(I=this.fetchRawTile(l,f,u,t),y(d,null,p,I));const w=await I;if(s(w)||!w.pixels||0===w.pixels.length)return{location:a,value:null};const R=h*this.rasterInfo.storageInfo.blockHeight+m,k=!w.mask||w.mask[R]?w.pixels.map((e=>e[R])):null,v=this.rasterInfo.dataType;if(("vector-magdir"===v||"vector-uv"===v)&&(null==k?void 0:k.length)>1){return{location:a,value:k,magdirValue:"vector-magdir"===v?[k[0],k[1]]:D([k[0],k[1]]),pyramidLevel:l}}return{location:a,value:k,pyramidLevel:l}}async fetchPixels(e,t,r,i={}){if(e=C(e),(i=this._getRequestOptionsWithSliceId(i)).requestRawData)return this._fetchPixels(e,t,r,i);const o=S(e.spatialReference),n=P(e);if(s(o)||0===n||1===n&&this._isGlobalWrappableSource)return this._fetchPixels(e,t,r,i);if(n>=3)return{extent:e,pixelBlock:null};const a=[],{xmin:l,xmax:c}=e,f=Math.round(o/(c-l)*t),u=f-Math.round((o/2-l)/(c-l)*t);let h=0;const m=[];for(let s=0;s<=n;s++){const d=new E({xmin:0===s?l:-o/2,xmax:s===n?c-o*s:o/2,ymin:e.ymin,ymax:e.ymax,spatialReference:e.spatialReference}),p=0===s?f-u:s===n?t-h:f;h+=p,m.push(p);const x=i.disableWrapAround&&s>0?null:this._fetchPixels(d,p,r,i);a.push(x)}const d=(await Promise.all(a)).map((e=>null==e?void 0:e.pixelBlock));let p=null;const x={width:t,height:r};if(this.rasterJobHandler){p=(await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:d,srcMosaicSize:x,destDimension:null,coefs:null,sampleSpacing:null,interpolation:"nearest",alignmentInfo:null,blockWidths:m},i)).pixelBlock}else p=k(d,x,{blockWidths:m});return{extent:e,srcExtent:W(e,this.rasterInfo.spatialReference,i.datumTransformation),pixelBlock:p}}async fetchRawPixels(e,t,r,i={}){t={x:Math.floor(t.x),y:Math.floor(t.y)};const o=await this._fetchRawTiles(e,t,r,i),{nativeExtent:s,nativePixelSize:a,storageInfo:l}=this.rasterInfo,c=2**e,f=a.x*c,u=a.y*c,h=new E({xmin:s.xmin+f*t.x,xmax:s.xmin+f*(t.x+r.width-1),ymin:s.ymax-u*(t.y+r.height-1),ymax:s.ymax-u*t.y,spatialReference:s.spatialReference});if(!o)return{extent:h,srcExtent:h,pixelBlock:null};const{pixelBlocks:m,mosaicSize:d}=o;if(1===m.length&&n(m[0])&&m[0].width===r.width&&m[0].height===r.height)return{extent:h,srcExtent:h,pixelBlock:o.pixelBlocks[0]};const p=e>0?l.pyramidBlockWidth:l.blockWidth,x=e>0?l.pyramidBlockHeight:l.blockHeight,g={x:t.x%p,y:t.y%x};let y;if(this.rasterJobHandler){y=(await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:m,srcMosaicSize:d,destDimension:r,clipOffset:g,clipSize:r,coefs:null,sampleSpacing:null,interpolation:i.interpolation,alignmentInfo:null,blockWidths:null},i)).pixelBlock}else y=k(m,d,{clipOffset:g,clipSize:r});return{extent:h,srcExtent:h,pixelBlock:y}}fetchRawTile(e,t,i,o){throw new r("BaseRaster:read-not-implemented","fetchRawTile() is not implemented")}computeExtent(e){return W(this.rasterInfo.extent,e)}decodePixelBlock(e,t){return!this.rasterJobHandler||t.useCanvas?R(e,t):this.rasterJobHandler.decode({data:e,options:t})}async request(e,r,i){var o,n;const{customFetchParameters:s}=this.ioConfig,{range:a,query:l,headers:c}=r;i=null!=(o=null!=(n=i)?n:r.retryCount)?o:this.ioConfig.retryCount;const f=a?{Range:`bytes=${a.from}-${a.to}`}:null;try{return await t(e,{...r,query:{...l,...s},headers:{...c,...f}})}catch(u){if(i>0)return i--,this.request(e,r,i);throw u}}getSliceIndex(e){const{multidimensionalInfo:t}=this.rasterInfo;if(!n(t))return null;if(!n(e)||0===e.length)return null;let r=0;const i=e[0].variableName;for(let o=0;o<t.variables.length;o++){const n=t.variables[o],s=n.dimensions;if(n.name!==i){r+=s.map((e=>this._getDimensionValuesCount(e))).reduce(((e,t)=>e+t));break}const a=s.map((e=>this._getDimensionValuesCount(e))),l=s.length;for(let t=0;t<l;t++){const i=e.find((e=>e.dimensionName===s[t].name));if(null==i)return null;const o=Array.isArray(i.values[0])?i.values[0][0]:i.values[0],n=this._getIndexFromDimensions(o,s[t]);if(-1===n)return null;a.shift(),r+=t===l-1?n:n*a.reduce(((e,t)=>e+t))}}return r}getTileExtentFromTileInfo(e,t,r,i){const o=i.lodAt(e);return this.getTileExtent({x:o.resolution,y:o.resolution},t,r,i.origin,i.spatialReference,i.size)}updateTileInfo(){const{storageInfo:e,spatialReference:t,extent:r,pixelSize:i}=this.rasterInfo;if(!e.tileInfo){const o=[],n=e.maximumPyramidLevel||0;let s=Math.max(i.x,i.y),a=1/.0254*96*s;for(let e=0;e<=n;e++)o.push({level:n-e,resolution:s,scale:a}),s*=2,a*=2;const l=new L({x:r.xmin,y:r.ymax,spatialReference:t});e.tileInfo=new p({origin:l,size:[e.blockWidth,e.blockHeight],spatialReference:t,lods:o}),e.isVirtualTileInfo=!0}}createRemoteDatasetStorageInfo(e,t=512,r=512,i){const{width:o,height:n,nativeExtent:s,pixelSize:a,spatialReference:l}=e,c=new L({x:s.xmin,y:s.ymax,spatialReference:l});null==i&&(i=Math.max(0,Math.round(Math.log(Math.max(o,n))/Math.LN2-8)));const f=this.computeBlockBoundary(s,512,512,{x:s.xmin,y:s.ymax},[a],i);e.storageInfo=new d({blockWidth:t,blockHeight:r,pyramidBlockWidth:t,pyramidBlockHeight:r,origin:c,firstPyramidLevel:1,maximumPyramidLevel:i,blockBoundary:f})}async computeBestPyramidLevelForLocation(e,t={}){return 0}computeBlockBoundary(e,t,r,i,o,n=0,s=2){if(1===o.length&&n>0){o=[...o];let{x:e,y:t}=o[0];for(let r=0;r<n;r++)e*=s,t*=s,o.push({x:e,y:t})}const a=[],{x:l,y:c}=i;for(let f=0;f<o.length;f++){const{x:i,y:n}=o[f];a.push({minCol:Math.floor((e.xmin-l+.1*i)/t/i),maxCol:Math.floor((e.xmax-l-.1*i)/t/i),minRow:Math.floor((c-e.ymax+.1*n)/r/n),maxRow:Math.floor((c-e.ymin-.1*n)/r/n)})}return a}getPyramidPixelSize(e){const{nativePixelSize:t}=this.rasterInfo,{pyramidResolutions:r,pyramidScalingFactor:i}=this.rasterInfo.storageInfo;if(0===e)return t;if(n(r)&&r.length)return r[e-1];const o=i**e;return{x:t.x*o,y:t.y*o}}identifyPixelLocation(e,t,r){const{spatialReference:i,nativeExtent:o}=this.rasterInfo,{blockWidth:n,blockHeight:s,maximumPyramidLevel:a,origin:l}=this.rasterInfo.storageInfo,c=B(e,i,r);if(!o.intersects(c))return null;if(t<0||t>a)return null;const f=this.getPyramidPixelSize(t),{x:u,y:h}=f,m=(l.y-c.y)/h/s,d=(c.x-l.x)/u/n,p=Math.min(s-1,Math.floor((m-Math.floor(m))*s)),x=Math.min(n-1,Math.floor((d-Math.floor(d))*n));return{pyramidLevel:t,row:Math.floor(m),col:Math.floor(d),rowOffset:p,colOffset:x,srcLocation:c}}getTileExtent(e,t,r,i,o,n){const[s,a]=n,l=i.x+r*s*e.x,c=l+s*e.x,f=i.y-t*a*e.y,u=f-a*e.y;return new E({xmin:l,xmax:c,ymin:u,ymax:f,spatialReference:o})}getBlockWidthHeight(e){return{blockWidth:e>0?this.rasterInfo.storageInfo.pyramidBlockWidth:this.rasterInfo.storageInfo.blockWidth,blockHeight:e>0?this.rasterInfo.storageInfo.pyramidBlockHeight:this.rasterInfo.storageInfo.blockHeight}}isBlockOutside(e,t,r){const i=this.rasterInfo.storageInfo.blockBoundary[e];return!i||i.maxRow<t||i.maxCol<r||i.minRow>t||i.minCol>r}async _fetchPixels(e,t,r,i={}){let o=P(e);if(o>=2)return{extent:e,pixelBlock:null};const s=this._getSourceDataInfo(e,t,r,i),{pyramidLevel:l,pyramidResolution:c,srcResolution:f,srcExtent:u,srcWidth:h,srcHeight:m}=s;if(0===h||0===m)return{extent:e,srcExtent:u,pixelBlock:null};const d=a(this.rasterInfo.transform),p="gcs-shift"===(null==d?void 0:d.type),x=n(S(e.spatialReference));!p&&x||(o=P(s.srcExtent,p));const g=this.rasterInfo.storageInfo,y={x:Math.floor((u.xmin-g.origin.x)/c.x+.1),y:Math.floor((g.origin.y-u.ymax)/c.y+.1)},I=await this._fetchRawTiles(l,y,{width:h,height:m,wrapCount:o},i);if(!I)return{extent:e,srcExtent:u,pixelBlock:null};const w=l>0?g.pyramidBlockWidth:g.blockWidth,R=l>0?g.pyramidBlockHeight:g.blockHeight,B=w===h&&R===m&&y.x%w==0&&y.y%R==0,M=new L({x:(e.xmax-e.xmin)/t,y:(e.ymax-e.ymin)/r,spatialReference:e.spatialReference}),C=!e.spatialReference.equals(this.rasterInfo.spatialReference),{datumTransformation:W}=i;if(!C&&B&&1===I.pixelBlocks.length&&w===t&&R===r&&f.x===M.x&&f.y===M.y)return{extent:e,srcExtent:u,pixelBlock:I.pixelBlocks[0]};const H=x&&n(S(u.spatialReference)),j=i.requestProjectedLocalDirections&&this.rasterInfo.dataType.startsWith("vector");j&&!this.rasterJobHandler&&await T();const D=this.rasterJobHandler?await this.rasterJobHandler.getProjectionOffsetGrid({projectedExtent:e,srcBufferExtent:I.extent,pixelSize:M.toJSON(),datumTransformation:W,rasterTransform:d,hasWrapAround:o>0||H,isAdaptive:!1!==this.ioConfig.optimizeProjectionAccuracy,includeGCSGrid:j},i):_({projectedExtent:e,srcBufferExtent:I.extent,pixelSize:M,datumTransformation:W,rasterTransform:d,hasWrapAround:o>0||H,isAdaptive:!1,includeGCSGrid:j});let E;const A=!i.requestRawData,z={rows:D.spacing[0],cols:D.spacing[1]},O=a(this._getRasterTileAlignmentInfo(l,I.extent.xmin)),{pixelBlocks:U,mosaicSize:G,isPartiallyFilled:J}=I;let q=null;if(this.rasterJobHandler){const e=await this.rasterJobHandler.mosaicAndTransform({srcPixelBlocks:U,srcMosaicSize:G,destDimension:A?{width:t,height:r}:null,coefs:A?D.coefficients:null,sampleSpacing:A?z:null,projectDirections:j,gcsGrid:j?D.gcsGrid:null,isUV:"vector-uv"===this.rasterInfo.dataType,interpolation:i.interpolation,alignmentInfo:O,blockWidths:null},i);({pixelBlock:E,localNorthDirections:q}=e)}else{const e=k(U,G,{alignmentInfo:O});E=A?v(e,{width:t,height:r},D.coefficients,z,i.interpolation):e,j&&D.gcsGrid&&(q=b({width:t,height:r},D.gcsGrid),E=F(E,this.rasterInfo.dataType,q))}return i.requestRawData||j?{srcExtent:u,pixelBlock:E,transformGrid:D,localNorthDirections:q,extent:e,isPartiallyFilled:J}:{srcExtent:u,extent:e,pixelBlock:E}}async _fetchRawTiles(e,t,r,i){const{origin:o,blockBoundary:s}=this.rasterInfo.storageInfo,{blockWidth:a,blockHeight:l}=this.getBlockWidthHeight(e);let{x:c,y:f}=t,{width:u,height:h,wrapCount:m}=r;const d=this._getRasterTileAlignmentInfo(e,0);i.buffer&&(c-=i.buffer.cols,f-=i.buffer.rows,u+=2*i.buffer.cols,h+=2*i.buffer.rows);let p=0,x=0,g=0;if(m&&n(d)){({worldColumnCountFromOrigin:x,originColumnOffset:g,rightPadding:p}=d);x*d.blockWidth-p>=c+u&&(p=0)}const y=Math.floor(c/a),I=Math.floor(f/l),w=Math.floor((c+u+p-1)/a),R=Math.floor((f+h+p-1)/l),k=s[e];if(!k)return null;const{minRow:v,minCol:b,maxCol:T,maxRow:S}=k;if(0===m&&(R<v||w<b||I>S||y>T))return null;const B=new Array;let M=!1;const C=null==this.ioConfig.allowPartialFill?i.allowPartialFill:this.ioConfig.allowPartialFill;for(let F=I;F<=R;F++)for(let t=y;t<=w;t++){let r=t;if(!i.disableWrapAround&&m&&n(d)&&x<=t&&(r=t-x-g),F>=v&&r>=b&&S>=F&&T>=r){const t=this._fetchRawTile(e,F,r,i);C?B.push(new Promise((e=>{t.then((t=>e(t))).catch((()=>{M=!0,e(null)}))}))):B.push(t)}else B.push(null)}if(0===B.length)return null;const P=await Promise.all(B),W={height:(R-I+1)*l,width:(w-y+1)*a},{spatialReference:_}=this.rasterInfo,H=this.getPyramidPixelSize(e),{x:j,y:D}=H;return{extent:new E({xmin:o.x+y*a*j,xmax:o.x+(w+1)*a*j,ymin:o.y-(R+1)*l*D,ymax:o.y-I*l*D,spatialReference:_}),pixelBlocks:P,mosaicSize:W,isPartiallyFilled:M}}_fetchRawTile(e,t,r,i){const o=this.rasterInfo.storageInfo.blockBoundary[e];if(!o)return Promise.resolve(null);const{minRow:n,minCol:a,maxCol:l,maxRow:f}=o;if(t<n||r<a||t>f||r>l)return Promise.resolve(null);const u=x(this.url,i.sliceId),h=`${e}/${t}/${r}`;let m=g(u,i.registryId,h);if(s(m)){const o=new AbortController;m=this.fetchRawTile(e,t,r,{...i,signal:o.signal}),y(u,i.registryId,h,m,o),m.catch((()=>I(u,i.registryId,h)))}return i.signal&&c(i,(()=>{w(u,i.registryId,h)})),m}_getIndexFromDimensions(e,t){const{extent:r,interval:i,unit:o,values:n}=t;if(null!=n&&n.length)return Array.isArray(n[0])?n.findIndex((t=>t[0]<=e&&t[1]>=e)):n.indexOf(e);if(e>r[1])return-1;const s=r[0];let a=-1;if("ISO8601"===o){var l;switch((null==(l=t.intervalUnit)?void 0:l.toLowerCase())||"seconds"){case"seconds":a=Math.round((e-s)/1e3/i);break;case"minutes":a=Math.round((e-s)/6e4/i);break;case"hours":a=Math.round((e-s)/36e5/i);break;case"days":a=Math.round((e-s)/864e5/i);break;case"years":a=Math.round((new Date(e).getUTCFullYear()-new Date(s).getUTCFullYear())/i);break;case"decades":a=Math.round((new Date(e).getUTCFullYear()-new Date(s).getUTCFullYear())/10/i)}return a}return Math.round((e-s)/i)}_getDimensionValuesCount(e){const{extent:t,interval:r,unit:i,values:o}=e;let n=(null==o?void 0:o.length)||0;if(n)return n;const s=t[0];if(0===n&&"ISO8601"===i){var a;switch((null==(a=e.intervalUnit)?void 0:a.toLowerCase())||"seconds"){case"seconds":n=Math.round((t[1]-t[0])/1e3/r);break;case"minutes":n=Math.round((t[1]-t[0])/6e4/r);break;case"hours":n=Math.round((t[1]-t[0])/36e5/r);break;case"days":n=Math.round((t[1]-t[0])/864e5/r);break;case"years":n=Math.round((new Date(t[1]).getUTCFullYear()-new Date(s).getUTCFullYear())/r);break;case"decades":n=Math.round((new Date(t[1]).getUTCFullYear()-new Date(s).getUTCFullYear())/10/r)}return n}return Math.round((t[1]-t[0])/r)}_getRasterTileAlignmentInfo(e,t){return null==this._rasterTileAlighmentInfo&&(this._rasterTileAlighmentInfo=H(this.rasterInfo)),n(this._rasterTileAlighmentInfo.pyramidsInfo)?{startX:t,halfWorldWidth:this._rasterTileAlighmentInfo.halfWorldWidth,hasGCSSShiftTransform:this._rasterTileAlighmentInfo.hasGCSSShiftTransform,...this._rasterTileAlighmentInfo.pyramidsInfo[e]}:null}_getSourceDataInfo(e,t,r,i={}){const o={datumTransformation:i.datumTransformation,pyramidLevel:0,pyramidResolution:null,srcExtent:null,srcHeight:0,srcResolution:null,srcWidth:0};i.srcResolution&&(o.srcResolution=i.srcResolution,this._updateSourceDataInfo(e,o));const n=this.rasterInfo.storageInfo.maximumPyramidLevel||0,{srcWidth:s,srcHeight:a,pyramidLevel:l}=o,c=s/t,f=a/r,u=l<n&&c*f>=16;if(u||l===n&&(c>A||f>A)||(0===s||0===a)){const s=new L({x:(e.xmax-e.xmin)/t,y:(e.ymax-e.ymin)/r,spatialReference:e.spatialReference});let a=j(s,this.rasterInfo.spatialReference,e,o.datumTransformation);const h=!a||i.srcResolution&&a.x+a.y<i.srcResolution.x+i.srcResolution.y;if(u&&i.srcResolution&&h){const e=Math.round(Math.log(Math.max(c,f))/Math.LN2)-1;if(n-l+3>=e){const t=2**e;a={x:i.srcResolution.x*t,y:i.srcResolution.y*t}}}a&&(o.srcResolution=a,this._updateSourceDataInfo(e,o))}return(o.srcWidth/t>A||o.srcHeight/r>A)&&(o.srcWidth=0,o.srcHeight=0),o}_updateSourceDataInfo(e,t){t.srcWidth=0,t.srcHeight=0;const r=this.rasterInfo.spatialReference,{srcResolution:i,datumTransformation:o}=t,{pyramidLevel:n,pyramidResolution:s,excessiveReading:l}=M(i,this.rasterInfo,this.ioConfig.sampling);if(l)return;let c=t.srcExtent||W(e,r,o);if(null==c)return;const f=a(this.rasterInfo.transform);f&&(c=f.inverseTransform(c)),t.srcExtent=c;const u=Math.ceil((c.xmax-c.xmin)/s.x-.1),h=Math.ceil((c.ymax-c.ymin)/s.y-.1);t.pyramidLevel=n,t.pyramidResolution=s,t.srcWidth=u,t.srcHeight=h}_getRequestOptionsWithSliceId(e){return n(this.rasterInfo.multidimensionalInfo)&&null==e.sliceId&&(e={...e,sliceId:this.getSliceIndex(e.multidimensionalDefinition)}),e}};e([f()],z.prototype,"_rasterTileAlighmentInfo",void 0),e([f({readOnly:!0})],z.prototype,"_isGlobalWrappableSource",null),e([f(m)],z.prototype,"url",null),e([f({type:String,json:{write:!0}})],z.prototype,"datasetName",void 0),e([f({type:String,json:{write:!0}})],z.prototype,"datasetFormat",void 0),e([f()],z.prototype,"rasterInfo",void 0),e([f()],z.prototype,"ioConfig",void 0),e([f()],z.prototype,"sourceJSON",void 0),z=e([u("geoscene.layers.support.rasterDatasets.BaseRaster")],z);const O=z;export{O as default};

package/layers/support/rasterDatasets/CloudRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as r}from"../../../core/maybe.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import s from"../RasterInfo.js";import a from"../RasterStorageInfo.js";import n from"../TileInfo.js";import l from"./BaseRaster.js";import f from"./DBFParser.js";import{isTransformSupported as m,readTransform as c}from"../rasterTransforms/utils.js";import p from"../../../rest/support/FeatureSet.js";import d from"../../../geometry/SpatialReference.js";import u from"../../../geometry/Extent.js";import h from"../../../geometry/Point.js";const g=new Map;g.set("int16","esriFieldTypeSmallInteger"),g.set("int32","esriFieldTypeInteger"),g.set("int64","esriFieldTypeInteger"),g.set("float32","esriFieldTypeSingle"),g.set("float64","esriFieldTypeDouble"),g.set("text","esriFieldTypeString");const y=8;let x=class extends l{constructor(){super(...arguments),this.storageInfo=null,this.datasetFormat="CRF"}async open(e){await this.init();const{data:r}=await this.request(this.url+"/conf.json",{signal:null==e?void 0:e.signal});if(!this._validateHeader(r))throw new t("cloudraster:open","Invalid or unsupported conf.json.");this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const{storageInfo:o,rasterInfo:i}=this._parseHeader(r);if("thematic"===i.dataType){const e=await this._fetchAuxiliaryInformation();i.attributeTable=e}this._set("storageInfo",o),this._set("rasterInfo",i),this.ioConfig.retryCount=this.ioConfig.retryCount||0}async fetchRawTile(e,t,r,o={}){const i=this.rasterInfo.storageInfo.maximumPyramidLevel-e;if(i<0)return null;const s=this._buildCacheFilePath(i,t,r,o.multidimensionalDefinition),a=this._getIndexRecordFromBundle(t,r),n=await this.request(s,{range:{from:0,to:this.storageInfo.headerSize-1},responseType:"array-buffer",signal:o.signal});if(!n)return null;const l=new Uint8Array(n.data),f=this._getTileEndAndContentType(l,a);if(0===f.recordSize)return null;const m=await this.request(s,{range:{from:f.position,to:f.position+f.recordSize},responseType:"array-buffer",signal:o.signal});return m?this.decodePixelBlock(m.data,{width:this.rasterInfo.storageInfo.tileInfo.size[0],height:this.rasterInfo.storageInfo.tileInfo.size[1],planes:null,pixelType:null}):null}_validateHeader(e){const t=["origin","extent","geodataXform","LODInfos","blockWidth","blockHeight","bandCount","pixelType","pixelSizeX","pixelSizeY","format","packetSize"];return e&&"RasterInfo"===e.type&&!t.some((t=>!e[t]))}_parseHeader(e){var t,r;const o=["u1","u2","u4","u8","s8","u16","s16","u32","s32","f32","f64"][e.pixelType],{bandCount:i,histograms:l,colormap:f,blockWidth:m,blockHeight:c,firstPyramidLevel:p,maximumPyramidLevel:g}=e,x=e.statistics&&e.statistics.map((e=>({min:e.min,max:e.max,avg:e.mean,stddev:e.standardDeviation,median:e.median,mode:e.mode}))),I=e.extent.spatialReference,S=null==(t=e.geodataXform)?void 0:t.spatialReference,w=new d(null!=I&&I.wkid||null!=I&&I.wkt?I:S);let v=new u({xmin:e.extent.xmin,ymin:e.extent.ymin,xmax:e.extent.xmax,ymax:e.extent.ymax,spatialReference:w});const b=new h({x:e.pixelSizeX,y:e.pixelSizeY,spatialReference:w}),T=Math.round((v.xmax-v.xmin)/b.x),j=Math.round((v.ymax-v.ymin)/b.y),z=this._parseTransform(e.geodataXform),_=z?v:null;z&&(v=z.forwardTransform(v),b.x=(v.xmax-v.xmin)/T,b.y=(v.ymax-v.ymin)/j);const k=null!=(r=e.properties)?r:{},C=e.format.toLowerCase().replace("cache/",""),R=new h(e.origin.x,e.origin.y,w);let F,P,H,D;if(f&&f.colors)for(F=[],P=0;P<f.colors.length;P++)H=f.colors[P],D=f.values?f.values[P]:P,F.push([D,255&H,H<<16>>>24,H<<8>>>24,H>>>24]);const L=e.LODInfos,M=[];for(P=0;P<L.levels.length;P++)M.push({level:L.levels[P],resolution:L.resolutions[P],scale:96/.0254*L.resolutions[P]});const B=new n({dpi:96,lods:M,format:C,origin:R,size:[m,c],spatialReference:w}),O={recordSize:y,packetSize:e.packetSize,headerSize:e.packetSize*e.packetSize*y+64},$=[{maxCol:Math.ceil(T/m)-1,maxRow:Math.ceil(j/c)-1,minCol:0,minRow:0}];let q=2;if(g>0)for(P=0;P<g;P++)$.push({maxCol:Math.ceil(T/q/m)-1,maxRow:Math.ceil(j/q/c)-1,minCol:0,minRow:0}),q*=2;const A=e.mdInfo;return{storageInfo:O,rasterInfo:new s({width:T,height:j,pixelType:o,bandCount:i,extent:v,nativeExtent:_,transform:z,spatialReference:w,pixelSize:b,keyProperties:k,statistics:x,histograms:l,multidimensionalInfo:A,colormap:F,storageInfo:new a({blockWidth:m,blockHeight:c,pyramidBlockWidth:m,pyramidBlockHeight:c,origin:R,tileInfo:B,firstPyramidLevel:p,maximumPyramidLevel:g,blockBoundary:$})})}}_parseTransform(e){var r,o;if(!m(e))throw new t("cloudraster:open","the data contains unsupported geodata transform types");const i=c(e);if("identity"===i.type)return null;if("polynomial"!==i.type||null==(r=i.forwardCoefficients)||!r.length||null==(o=i.inverseCoefficients)||!o.length)throw new t("cloudraster:open","the data contains unsupported geodata transforms - both forward and inverse coefficients are required currently");return i}async _fetchAuxiliaryInformation(e){const t=this.request(this.url+"/conf.vat.json",{signal:e}).then((e=>e.data)).catch((()=>null)),r=this.request(this.url+"/conf.vat.dbf",{responseType:"array-buffer",signal:e}).then((e=>e.data)).catch((()=>null)),o=await Promise.all([t,r]);let i;if(o[0]){let e=o[0].fields;const t=o[0].values;if(e&&t){e=e.map((e=>({type:"OID"===e.name?"esriFieldTypeOID":g.get(e.type),name:e.name,alias:e.alias||e.name})));const r=t.map((e=>({attributes:e})));e&&t&&(i={fields:e,features:r})}}if(!i&&o[1]){i=f.parse(o[1]).recordSet}return p.fromJSON(i)}_buildCacheFilePath(e,t,o,i){const s=this.storageInfo.packetSize,a=Math.floor(t/s)*s,n=Math.floor(o/s)*s,l="R"+this._toHexString4(a)+"C"+this._toHexString4(n);let f="L";f+=e>=10?e.toString():"0"+e.toString();const{multidimensionalInfo:m}=this.rasterInfo,c=null==i?void 0:i[0];if(!r(m)||!c)return`${this.url}/_alllayers/${f}/${l}.bundle`;let p=m.variables.find((e=>e.name===c.variableName)).dimensions[0].values.indexOf(c.values[0]).toString(16);const d=4-p.length;for(let r=0;r<d;r++)p="0"+p;return p="S"+p,`${this.url}/_alllayers/${c.variableName}/${p}/${f}/${l}.bundle`}_getIndexRecordFromBundle(e,t){const r=this.storageInfo.packetSize,o=r*(e%r)+t%r;if(o<0)throw"Invalid level / row / col";return 20+o*this.storageInfo.recordSize+44}_getTileEndAndContentType(e,t){const r=e.subarray(t,t+8);let o,i=0;for(o=0;o<5;o++)i|=(255&r[o])<<8*o;const s=0xffffffffff&i;for(i=0,o=5;o<8;o++)i|=(255&r[o])<<8*(o-5);return{position:s,recordSize:0xffffffffff&i}}_toHexString4(e){let t=e.toString(16);if(4!==t.length){let e=4-t.length;for(;e-- >0;)t="0"+t}return t}};e([o({readOnly:!0})],x.prototype,"storageInfo",void 0),e([o({type:String,json:{write:!0}})],x.prototype,"datasetFormat",void 0),x=e([i("esri.layers.support.rasterDatasets.CloudRaster")],x);const I=x;export{I as default};
+import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as r}from"../../../core/maybe.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import s from"../RasterInfo.js";import a from"../RasterStorageInfo.js";import n from"../TileInfo.js";import l from"./BaseRaster.js";import f from"./DBFParser.js";import{isTransformSupported as m,readTransform as c}from"../rasterTransforms/utils.js";import p from"../../../rest/support/FeatureSet.js";import d from"../../../geometry/SpatialReference.js";import u from"../../../geometry/Extent.js";import h from"../../../geometry/Point.js";const g=new Map;g.set("int16","esriFieldTypeSmallInteger"),g.set("int32","esriFieldTypeInteger"),g.set("int64","esriFieldTypeInteger"),g.set("float32","esriFieldTypeSingle"),g.set("float64","esriFieldTypeDouble"),g.set("text","esriFieldTypeString");const y=8;let x=class extends l{constructor(){super(...arguments),this.storageInfo=null,this.datasetFormat="CRF"}async open(e){await this.init();const{data:r}=await this.request(this.url+"/conf.json",{signal:null==e?void 0:e.signal});if(!this._validateHeader(r))throw new t("cloudraster:open","Invalid or unsupported conf.json.");this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const{storageInfo:o,rasterInfo:i}=this._parseHeader(r);if("thematic"===i.dataType){const e=await this._fetchAuxiliaryInformation();i.attributeTable=e}this._set("storageInfo",o),this._set("rasterInfo",i),this.ioConfig.retryCount=this.ioConfig.retryCount||0}async fetchRawTile(e,t,r,o={}){const i=this.rasterInfo.storageInfo.maximumPyramidLevel-e;if(i<0)return null;const s=this._buildCacheFilePath(i,t,r,o.multidimensionalDefinition),a=this._getIndexRecordFromBundle(t,r),n=await this.request(s,{range:{from:0,to:this.storageInfo.headerSize-1},responseType:"array-buffer",signal:o.signal});if(!n)return null;const l=new Uint8Array(n.data),f=this._getTileEndAndContentType(l,a);if(0===f.recordSize)return null;const m=await this.request(s,{range:{from:f.position,to:f.position+f.recordSize},responseType:"array-buffer",signal:o.signal});return m?this.decodePixelBlock(m.data,{width:this.rasterInfo.storageInfo.tileInfo.size[0],height:this.rasterInfo.storageInfo.tileInfo.size[1],planes:null,pixelType:null}):null}_validateHeader(e){const t=["origin","extent","geodataXform","LODInfos","blockWidth","blockHeight","bandCount","pixelType","pixelSizeX","pixelSizeY","format","packetSize"];return e&&"RasterInfo"===e.type&&!t.some((t=>!e[t]))}_parseHeader(e){var t,r;const o=["u1","u2","u4","u8","s8","u16","s16","u32","s32","f32","f64"][e.pixelType],{bandCount:i,histograms:l,colormap:f,blockWidth:m,blockHeight:c,firstPyramidLevel:p,maximumPyramidLevel:g}=e,x=e.statistics&&e.statistics.map((e=>({min:e.min,max:e.max,avg:e.mean,stddev:e.standardDeviation,median:e.median,mode:e.mode}))),I=e.extent.spatialReference,S=null==(t=e.geodataXform)?void 0:t.spatialReference,w=new d(null!=I&&I.wkid||null!=I&&I.wkt?I:S);let v=new u({xmin:e.extent.xmin,ymin:e.extent.ymin,xmax:e.extent.xmax,ymax:e.extent.ymax,spatialReference:w});const b=new h({x:e.pixelSizeX,y:e.pixelSizeY,spatialReference:w}),T=Math.round((v.xmax-v.xmin)/b.x),j=Math.round((v.ymax-v.ymin)/b.y),z=this._parseTransform(e.geodataXform),_=z?v:null;z&&(v=z.forwardTransform(v),b.x=(v.xmax-v.xmin)/T,b.y=(v.ymax-v.ymin)/j);const k=null!=(r=e.properties)?r:{},C=e.format.toLowerCase().replace("cache/",""),R=new h(e.origin.x,e.origin.y,w);let F,P,H,D;if(f&&f.colors)for(F=[],P=0;P<f.colors.length;P++)H=f.colors[P],D=f.values?f.values[P]:P,F.push([D,255&H,H<<16>>>24,H<<8>>>24,H>>>24]);const L=e.LODInfos,M=[];for(P=0;P<L.levels.length;P++)M.push({level:L.levels[P],resolution:L.resolutions[P],scale:96/.0254*L.resolutions[P]});const B=new n({dpi:96,lods:M,format:C,origin:R,size:[m,c],spatialReference:w}),O={recordSize:y,packetSize:e.packetSize,headerSize:e.packetSize*e.packetSize*y+64},$=[{maxCol:Math.ceil(T/m)-1,maxRow:Math.ceil(j/c)-1,minCol:0,minRow:0}];let q=2;if(g>0)for(P=0;P<g;P++)$.push({maxCol:Math.ceil(T/q/m)-1,maxRow:Math.ceil(j/q/c)-1,minCol:0,minRow:0}),q*=2;const A=e.mdInfo;return{storageInfo:O,rasterInfo:new s({width:T,height:j,pixelType:o,bandCount:i,extent:v,nativeExtent:_,transform:z,spatialReference:w,pixelSize:b,keyProperties:k,statistics:x,histograms:l,multidimensionalInfo:A,colormap:F,storageInfo:new a({blockWidth:m,blockHeight:c,pyramidBlockWidth:m,pyramidBlockHeight:c,origin:R,tileInfo:B,firstPyramidLevel:p,maximumPyramidLevel:g,blockBoundary:$})})}}_parseTransform(e){var r,o;if(!m(e))throw new t("cloudraster:open","the data contains unsupported geodata transform types");const i=c(e);if("identity"===i.type)return null;if("polynomial"!==i.type||null==(r=i.forwardCoefficients)||!r.length||null==(o=i.inverseCoefficients)||!o.length)throw new t("cloudraster:open","the data contains unsupported geodata transforms - both forward and inverse coefficients are required currently");return i}async _fetchAuxiliaryInformation(e){const t=this.request(this.url+"/conf.vat.json",{signal:e}).then((e=>e.data)).catch((()=>null)),r=this.request(this.url+"/conf.vat.dbf",{responseType:"array-buffer",signal:e}).then((e=>e.data)).catch((()=>null)),o=await Promise.all([t,r]);let i;if(o[0]){let e=o[0].fields;const t=o[0].values;if(e&&t){e=e.map((e=>({type:"OID"===e.name?"esriFieldTypeOID":g.get(e.type),name:e.name,alias:e.alias||e.name})));const r=t.map((e=>({attributes:e})));e&&t&&(i={fields:e,features:r})}}if(!i&&o[1]){i=f.parse(o[1]).recordSet}return p.fromJSON(i)}_buildCacheFilePath(e,t,o,i){const s=this.storageInfo.packetSize,a=Math.floor(t/s)*s,n=Math.floor(o/s)*s,l="R"+this._toHexString4(a)+"C"+this._toHexString4(n);let f="L";f+=e>=10?e.toString():"0"+e.toString();const{multidimensionalInfo:m}=this.rasterInfo,c=null==i?void 0:i[0];if(!r(m)||!c)return`${this.url}/_alllayers/${f}/${l}.bundle`;let p=m.variables.find((e=>e.name===c.variableName)).dimensions[0].values.indexOf(c.values[0]).toString(16);const d=4-p.length;for(let r=0;r<d;r++)p="0"+p;return p="S"+p,`${this.url}/_alllayers/${c.variableName}/${p}/${f}/${l}.bundle`}_getIndexRecordFromBundle(e,t){const r=this.storageInfo.packetSize,o=r*(e%r)+t%r;if(o<0)throw"Invalid level / row / col";return 20+o*this.storageInfo.recordSize+44}_getTileEndAndContentType(e,t){const r=e.subarray(t,t+8);let o,i=0;for(o=0;o<5;o++)i|=(255&r[o])<<8*o;const s=0xffffffffff&i;for(i=0,o=5;o<8;o++)i|=(255&r[o])<<8*(o-5);return{position:s,recordSize:0xffffffffff&i}}_toHexString4(e){let t=e.toString(16);if(4!==t.length){let e=4-t.length;for(;e-- >0;)t="0"+t}return t}};e([o({readOnly:!0})],x.prototype,"storageInfo",void 0),e([o({type:String,json:{write:!0}})],x.prototype,"datasetFormat",void 0),x=e([i("geoscene.layers.support.rasterDatasets.CloudRaster")],x);const I=x;export{I as default};

package/layers/support/rasterDatasets/ImageAuxRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{unwrap as r}from"../../../core/maybe.js";import{eachAlways as s,createAbortError as a}from"../../../core/promiseUtils.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import n from"./BaseRaster.js";import l from"./InMemoryRaster.js";import{parsePAMInfo as m}from"./pamParser.js";import{getFormat as p}from"../rasterFormats/RasterCodec.js";import{estimateStatisticsFromHistograms as c}from"../rasterFunctions/pixelUtils.js";import f from"../rasterTransforms/PolynomialTransform.js";import u from"../../../geometry/SpatialReference.js";import h from"../../../geometry/Extent.js";let d=class extends n{async open(t){await this.init();const e=await this._fetchData(t);let{spatialReference:r,statistics:s,histograms:a,transform:o}=await this._fetchAuxiliaryData(t);const i=!r;i&&(r=new u({wkid:3857})),null!=a&&a.length&&null==s&&(s=c(a));const{width:n,height:m}=e;let p=new h({xmin:-.5,ymin:.5-m,xmax:n-.5,ymax:.5,spatialReference:r});const f=o?o.forwardTransform(p):p;let d=!0;if(o){const t=o.forwardCoefficients;d=t&&0===t[1]&&0===t[2],d&&(o=null,p=f)}const y=new l({data:{extent:f,nativeExtent:p,transform:o,pixelBlock:e,statistics:s,histograms:a,keyProperties:{DateType:"Processed"},isPseudoSpatialReference:i}});await y.open(),this._set("rasterInfo",y.rasterInfo),this._inMemoryRaster=y}fetchRawTile(t,e,r,s={}){return this._inMemoryRaster.fetchRawTile(t,e,r,s)}async _fetchData(t){const{data:r}=await this.request(this.url,{responseType:"array-buffer",signal:null==t?void 0:t.signal}),s=p(r).toUpperCase();if("JPG"!==s&&"PNG"!==s&&"GIF"!==s&&"BMP"!==s)throw new e("image-aux-raster:open","the data is not a supported format");this._set("datasetFormat",s);return await this.decodePixelBlock(r,{format:"jpg",width:1,height:1,useCanvas:!0})}async _fetchAuxiliaryData(t){var e,o;const i=r(null==t?void 0:t.signal),n=null!=(e=this.ioConfig.skipExtensions)?e:[],l=n.indexOf("aux.xml")>-1?null:this.request(this.url+".aux.xml",{responseType:"xml",signal:i}),p=this.datasetFormat,c="JPG"===p?"jgw":"PNG"===p?"pgw":"BMP"===p?"bpw":null,u=n.indexOf(c)>-1?null:this.request(this.url.slice(0,this.url.lastIndexOf("."))+"."+c,{responseType:"text",signal:i}),h=await s([l,u]);if(null!=i&&i.aborted)throw a();const d=m(null==(o=h[0].value)?void 0:o.data);if(!d.transform){const t=h[1].value?h[1].value.data.split("\n").slice(0,6).map((t=>Number(t))):null;d.transform=6===(null==t?void 0:t.length)?new f({forwardCoefficients:[t[4],t[5],t[0],-t[1],t[2],-t[3]]}):null}return d}};t([o({type:String,json:{write:!0}})],d.prototype,"datasetFormat",void 0),d=t([i("esri.layers.support.rasterDatasets.ImageAuxRaster")],d);const y=d;export{y as default};
+import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{unwrap as r}from"../../../core/maybe.js";import{eachAlways as s,createAbortError as a}from"../../../core/promiseUtils.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import n from"./BaseRaster.js";import l from"./InMemoryRaster.js";import{parsePAMInfo as m}from"./pamParser.js";import{getFormat as p}from"../rasterFormats/RasterCodec.js";import{estimateStatisticsFromHistograms as c}from"../rasterFunctions/pixelUtils.js";import f from"../rasterTransforms/PolynomialTransform.js";import u from"../../../geometry/SpatialReference.js";import h from"../../../geometry/Extent.js";let d=class extends n{async open(t){await this.init();const e=await this._fetchData(t);let{spatialReference:r,statistics:s,histograms:a,transform:o}=await this._fetchAuxiliaryData(t);const i=!r;i&&(r=new u({wkid:3857})),null!=a&&a.length&&null==s&&(s=c(a));const{width:n,height:m}=e;let p=new h({xmin:-.5,ymin:.5-m,xmax:n-.5,ymax:.5,spatialReference:r});const f=o?o.forwardTransform(p):p;let d=!0;if(o){const t=o.forwardCoefficients;d=t&&0===t[1]&&0===t[2],d&&(o=null,p=f)}const y=new l({data:{extent:f,nativeExtent:p,transform:o,pixelBlock:e,statistics:s,histograms:a,keyProperties:{DateType:"Processed"},isPseudoSpatialReference:i}});await y.open(),this._set("rasterInfo",y.rasterInfo),this._inMemoryRaster=y}fetchRawTile(t,e,r,s={}){return this._inMemoryRaster.fetchRawTile(t,e,r,s)}async _fetchData(t){const{data:r}=await this.request(this.url,{responseType:"array-buffer",signal:null==t?void 0:t.signal}),s=p(r).toUpperCase();if("JPG"!==s&&"PNG"!==s&&"GIF"!==s&&"BMP"!==s)throw new e("image-aux-raster:open","the data is not a supported format");this._set("datasetFormat",s);return await this.decodePixelBlock(r,{format:"jpg",width:1,height:1,useCanvas:!0})}async _fetchAuxiliaryData(t){var e,o;const i=r(null==t?void 0:t.signal),n=null!=(e=this.ioConfig.skipExtensions)?e:[],l=n.indexOf("aux.xml")>-1?null:this.request(this.url+".aux.xml",{responseType:"xml",signal:i}),p=this.datasetFormat,c="JPG"===p?"jgw":"PNG"===p?"pgw":"BMP"===p?"bpw":null,u=n.indexOf(c)>-1?null:this.request(this.url.slice(0,this.url.lastIndexOf("."))+"."+c,{responseType:"text",signal:i}),h=await s([l,u]);if(null!=i&&i.aborted)throw a();const d=m(null==(o=h[0].value)?void 0:o.data);if(!d.transform){const t=h[1].value?h[1].value.data.split("\n").slice(0,6).map((t=>Number(t))):null;d.transform=6===(null==t?void 0:t.length)?new f({forwardCoefficients:[t[4],t[5],t[0],-t[1],t[2],-t[3]]}):null}return d}};t([o({type:String,json:{write:!0}})],d.prototype,"datasetFormat",void 0),d=t([i("geoscene.layers.support.rasterDatasets.ImageAuxRaster")],d);const y=d;export{y as default};

package/layers/support/rasterDatasets/ImageServerRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as i,isNone as s,unwrap as a}from"../../../core/maybe.js";import{urlToObject as l}from"../../../core/urlUtils.js";import{property as r}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as n}from"../../../core/accessorSupport/decorators/subclass.js";import o from"../RasterInfo.js";import m from"../RasterStorageInfo.js";import{readServiceTileInfo as c}from"../serviceTileInfoProperty.js";import h from"../TileInfo.js";import{TilemapCache as u}from"../TilemapCache.js";import f from"./BaseRaster.js";import{setValidBoundary as p}from"../rasterFunctions/pixelUtils.js";import d from"../rasterTransforms/GCSShiftTransform.js";import y from"../../../rest/support/FeatureSet.js";import x from"../../../geometry/SpatialReference.js";import g from"../../../geometry/Extent.js";import v from"../../../geometry/Point.js";let S=class extends f{constructor(){super(...arguments),this._levelOffset=0,this._slices=null,this._tilemapCache=null,this.datasetFormat="RasterTileServer"}async open(e){await this.init();const s=e&&e.signal,a=this.sourceJSON?{data:this.sourceJSON}:await this.request(this.url,{query:{f:"json"},signal:s});a.ssl&&(this.url=this.url.replace(/^http:/i,"https:"));const r=a.data;if(this.sourceJSON=r,!r)throw new t("imageserverraster:open","cannot initialize tiled image service, missing service info");if(!r.tileInfo)throw new t("imageserverraster:open","use ImageryLayer to open non-tiled image services");this._fixScaleInServiceInfo();const n=["jpg","jpeg","png","png8","png24","png32","mixed"];this.tileType=r.cacheType,null==this.tileType&&(n.indexOf(r.tileInfo.format.toLowerCase())>-1?this.tileType="Map":"lerc"===r.tileInfo.format.toLowerCase()?this.tileType="Elevation":this.tileType="Raster"),this.datasetName=r.name.slice(r.name.indexOf("/")+1);const o=await this._fetchRasterInfo({signal:s});if(!i(o))throw new t("image-server-raster:open","cannot initialize image service");{const e="Map"===this.tileType?c(r.tileInfo,r):h.fromJSON(r.tileInfo),{extent:t,pixelSize:i}=o,s=.5/o.width*i.x;let a,l;const n=e.lodAt(Math.max.apply(null,e.lods.map((e=>e.level))));"Map"!==this.tileType&&0!==r.maxScale&&("Raster"===this.tileType?(a=e.lods.find((e=>e.resolution===i.x)),a||(a=e.lods[e.lods.length-1])):(a=e.lods.find((e=>Math.abs(e.scale-r.maxScale)<s)),a||(a=e.lods.filter((e=>e.scale>r.maxScale)).sort(((e,t)=>e.scale>t.scale?1:-1))[0])),i.x=i.y=a.resolution,o.width=Math.ceil((t.xmax-t.xmin)/i.x-.1),o.height=Math.ceil((t.ymax-t.ymin)/i.y-.1)),a||(a=n);const u=e.lodAt(Math.min.apply(null,e.lods.map((e=>e.level))));"Map"===this.tileType?this._levelOffset=e.lods[0].level:0!==r.minScale&&"Elevation"===this.tileType&&(l=e.lods.find((e=>Math.abs(e.scale-r.minScale)<s)),this._levelOffset=l.level-u.level),l||(l=u);const f=Math.max(i.x,i.y);(Math.abs(i.x-i.y)>s||!e.lods.some((e=>Math.abs(e.resolution-f)<s)))&&(i.x=i.y=a.resolution,o.width=Math.ceil((t.xmax-t.xmin)/i.x-.1),o.height=Math.ceil((t.ymax-t.ymin)/i.y-.1));const p=a.level-l.level,[d,y]=e.size,x=[];e.lods.forEach((e=>{e.level>=l.level&&e.level<=a.level&&x.push({x:e.resolution,y:e.resolution})})),x.sort(((e,t)=>e.x-t.x));const g=this.computeBlockBoundary(t,d,y,e.origin,x,p),v=x.length>1?x.slice(1):null;o.storageInfo=new m({blockWidth:e.size[0],blockHeight:e.size[1],pyramidBlockWidth:e.size[0],pyramidBlockHeight:e.size[1],pyramidResolutions:v,compression:e.format,origin:e.origin,firstPyramidLevel:1,maximumPyramidLevel:p,tileInfo:e,blockBoundary:g}),this._fixGCSShift(o),this._set("rasterInfo",o)}if(r.capabilities.toLowerCase().indexOf("tilemap")>-1){const e={tileInfo:o.storageInfo.tileInfo,parsedUrl:l(this.url),url:this.url,tileServers:[],type:"tile"};this._tilemapCache=new u({layer:e})}}async fetchRawTile(e,t,i,s={}){if(this._slices&&null==s.sliceId)return null;const{storageInfo:a,extent:l}=this.rasterInfo,r=a.maximumPyramidLevel-e+this._levelOffset,n=`${this.url}/tile/${r}/${t}/${i}`,o=this._slices?{sliceId:s.sliceId||0}:null,{data:m}=await this.request(n,{query:o,responseType:"array-buffer",signal:s.signal});if(!m)return null;const c=await this.decodePixelBlock(m,{width:a.tileInfo.size[0],height:a.tileInfo.size[1],planes:null,pixelType:null,isPoint:"Elevation"===this.tileType}),h=a.blockBoundary[e];if("jpg"!==a.compression||i>h.minCol&&i<h.maxCol&&t>h.minRow&&t<h.maxRow)return c;const{origin:u,blockWidth:f,blockHeight:d}=a,{x:y,y:x}=this.getPyramidPixelSize(e),g=Math.round((l.xmin-u.x)/y)%f,v=Math.round((l.xmax-u.x)/y)%f||f,S=Math.round((u.y-l.ymax)/x)%d,w=Math.round((u.y-l.ymin)/x)%d||d,I=i===h.minCol?g:0,b=t===h.minRow?S:0,j=i===h.maxCol?v:f,T=t===h.maxRow?w:d;return p(c,{x:I,y:b},{width:j-I,height:T-b}),c}getSliceIndex(e){if(!this._slices||s(e)||0===e.length)return null;const t=e;for(let i=0;i<this._slices.length;i++){const e=this._slices[i].multidimensionalDefinition;if(e.length===t.length&&!e.some((e=>{const i=t.find((t=>e.variableName===t.variableName&&t.dimensionName===e.dimensionName));if(!i)return!0;return(Array.isArray(e.values[0])?`${e.values[0][0]}-${e.values[0][1]}`:e.values[0])!==(Array.isArray(i.values[0])?`${i.values[0][0]}-${i.values[0][1]}`:i.values[0])})))return i}return null}async fetchVariableStatisticsHistograms(e,t){const i=this.request(this.url+"/statistics",{query:{variable:e,f:"json"},signal:t}).then((e=>{var t;return null==(t=e.data)?void 0:t.statistics})),s=this.request(this.url+"/histograms",{query:{variable:e,f:"json"},signal:t}).then((e=>{var t;return null==(t=e.data)?void 0:t.histograms})),a=await Promise.all([i,s]);return a[0]&&a[0].forEach((e=>{e.avg=e.mean,e.stddev=e.standardDeviation})),{statistics:a[0]||null,histograms:a[1]||null}}async computeBestPyramidLevelForLocation(e,t={}){if(!this._tilemapCache)return 0;let i=this.identifyPixelLocation(e,0,a(t.datumTransformation));if(null===i)return null;let s=0;const{maximumPyramidLevel:l}=this.rasterInfo.storageInfo;let r=l-s+this._levelOffset;const n=i.srcLocation;for(;r>=0;){try{if("available"===await this._tilemapCache.fetchAvailability(r,i.row,i.col,t))break}catch{}if(r--,s++,i=this.identifyPixelLocation(n,s,a(t.datumTransformation)),null===i)return null}return-1===r||null==i?null:s}async _fetchRasterInfo(e){const t=this.sourceJSON,i=Math.ceil((t.extent.xmax-t.extent.xmin)/t.pixelSizeX-.1),s=Math.ceil((t.extent.ymax-t.extent.ymin)/t.pixelSizeY-.1),a=x.fromJSON(t.spatialReference||t.extent.spatialReference);if("Map"===this.tileType)return new o({width:i,height:s,bandCount:3,extent:g.fromJSON(t.extent),spatialReference:a,pixelSize:new v({x:t.pixelSizeX,y:t.pixelSizeY,spatialReference:a}),pixelType:"u8",statistics:null,keyProperties:{DataType:"processed"}});const{slice:l,signal:r}=e,n=!!t.hasRasterAttributeTable&&this.request(this.url+"/rasterAttributeTable",{query:{slice:l,f:"json"},signal:r}).then((e=>y.fromJSON(e.data))).catch((()=>null)),m=!!t.hasColormap&&this.request(this.url+"/colormap",{query:{slice:l,f:"json"},signal:r}).then((e=>{var t;return null==(t=e.data)?void 0:t.colormap})),c=!!t.hasHistograms&&this.request(this.url+"/histograms",{query:{slice:l,f:"json"},signal:r}).then((e=>{var t;return null==(t=e.data)?void 0:t.histograms})),h=this.request(this.url+"/keyProperties",{query:{f:"json"},signal:r}).then((e=>e.data)).catch((()=>{})),u=!!t.hasMultidimensions&&this._fetchMultidimensionalInfo(),f=!!t.hasMultidimensions&&this.request(this.url+"/slices",{query:{f:"json"},signal:r}).then((e=>e.data&&e.data.slices)).catch((()=>{}));return Promise.all([n,m,c,h,u,f]).then((e=>{let l=null;if(t.minValues&&t.minValues.length===t.bandCount){l=[];for(let e=0;e<t.minValues.length;e++)l.push({min:t.minValues[e],max:t.maxValues[e],avg:t.meanValues[e],stddev:t.stdvValues[e]})}return this._slices=e[5]||null,new o({width:i,height:s,bandCount:t.bandCount,extent:g.fromJSON(t.extent),spatialReference:a,pixelSize:new v({x:t.pixelSizeX,y:t.pixelSizeY,spatialReference:a}),pixelType:t.pixelType.toLowerCase(),statistics:l,attributeTable:e[0]||null,colormap:e[1]||null,histograms:e[2]||null,keyProperties:e[3]||{},multidimensionalInfo:e[4]||null})}))}async _fetchMultidimensionalInfo(e){var t;const i=await this.request(this.url+"/multidimensionalInfo",{query:{f:"json"},signal:e}).then((e=>{var t;return null==(t=e.data)?void 0:t.multidimensionalInfo}));return null!=(t=i.variables)&&t.length&&i.variables.forEach((e=>{var t;null!=(t=e.statistics)&&t.length&&e.statistics.forEach((e=>{e.avg=e.mean,e.stddev=e.standardDeviation}))})),i}_fixScaleInServiceInfo(){const{sourceJSON:e}=this;e.minScale&&e.minScale<0&&(e.minScale=0),e.maxScale&&e.maxScale<0&&(e.maxScale=0)}_fixGCSShift(e){const{extent:t,spatialReference:i}=e;0===t.xmin&&360===t.xmax&&i.wkid&&i.isGeographic&&(e.nativeExtent=e.extent,e.transform=new d,e.extent=e.transform.forwardTransform(t))}};e([r({type:String,json:{write:!0}})],S.prototype,"datasetFormat",void 0),e([r()],S.prototype,"tileType",void 0),S=e([n("esri.layers.support.rasterDatasets.ImageServerRaster")],S);const w=S;export{w as default};
+import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as i,isNone as s,unwrap as a}from"../../../core/maybe.js";import{urlToObject as l}from"../../../core/urlUtils.js";import{property as r}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as n}from"../../../core/accessorSupport/decorators/subclass.js";import o from"../RasterInfo.js";import m from"../RasterStorageInfo.js";import{readServiceTileInfo as c}from"../serviceTileInfoProperty.js";import h from"../TileInfo.js";import{TilemapCache as u}from"../TilemapCache.js";import f from"./BaseRaster.js";import{setValidBoundary as p}from"../rasterFunctions/pixelUtils.js";import d from"../rasterTransforms/GCSShiftTransform.js";import y from"../../../rest/support/FeatureSet.js";import x from"../../../geometry/SpatialReference.js";import g from"../../../geometry/Extent.js";import v from"../../../geometry/Point.js";let S=class extends f{constructor(){super(...arguments),this._levelOffset=0,this._slices=null,this._tilemapCache=null,this.datasetFormat="RasterTileServer"}async open(e){await this.init();const s=e&&e.signal,a=this.sourceJSON?{data:this.sourceJSON}:await this.request(this.url,{query:{f:"json"},signal:s});a.ssl&&(this.url=this.url.replace(/^http:/i,"https:"));const r=a.data;if(this.sourceJSON=r,!r)throw new t("imageserverraster:open","cannot initialize tiled image service, missing service info");if(!r.tileInfo)throw new t("imageserverraster:open","use ImageryLayer to open non-tiled image services");this._fixScaleInServiceInfo();const n=["jpg","jpeg","png","png8","png24","png32","mixed"];this.tileType=r.cacheType,null==this.tileType&&(n.indexOf(r.tileInfo.format.toLowerCase())>-1?this.tileType="Map":"lerc"===r.tileInfo.format.toLowerCase()?this.tileType="Elevation":this.tileType="Raster"),this.datasetName=r.name.slice(r.name.indexOf("/")+1);const o=await this._fetchRasterInfo({signal:s});if(!i(o))throw new t("image-server-raster:open","cannot initialize image service");{const e="Map"===this.tileType?c(r.tileInfo,r):h.fromJSON(r.tileInfo),{extent:t,pixelSize:i}=o,s=.5/o.width*i.x;let a,l;const n=e.lodAt(Math.max.apply(null,e.lods.map((e=>e.level))));"Map"!==this.tileType&&0!==r.maxScale&&("Raster"===this.tileType?(a=e.lods.find((e=>e.resolution===i.x)),a||(a=e.lods[e.lods.length-1])):(a=e.lods.find((e=>Math.abs(e.scale-r.maxScale)<s)),a||(a=e.lods.filter((e=>e.scale>r.maxScale)).sort(((e,t)=>e.scale>t.scale?1:-1))[0])),i.x=i.y=a.resolution,o.width=Math.ceil((t.xmax-t.xmin)/i.x-.1),o.height=Math.ceil((t.ymax-t.ymin)/i.y-.1)),a||(a=n);const u=e.lodAt(Math.min.apply(null,e.lods.map((e=>e.level))));"Map"===this.tileType?this._levelOffset=e.lods[0].level:0!==r.minScale&&"Elevation"===this.tileType&&(l=e.lods.find((e=>Math.abs(e.scale-r.minScale)<s)),this._levelOffset=l.level-u.level),l||(l=u);const f=Math.max(i.x,i.y);(Math.abs(i.x-i.y)>s||!e.lods.some((e=>Math.abs(e.resolution-f)<s)))&&(i.x=i.y=a.resolution,o.width=Math.ceil((t.xmax-t.xmin)/i.x-.1),o.height=Math.ceil((t.ymax-t.ymin)/i.y-.1));const p=a.level-l.level,[d,y]=e.size,x=[];e.lods.forEach((e=>{e.level>=l.level&&e.level<=a.level&&x.push({x:e.resolution,y:e.resolution})})),x.sort(((e,t)=>e.x-t.x));const g=this.computeBlockBoundary(t,d,y,e.origin,x,p),v=x.length>1?x.slice(1):null;o.storageInfo=new m({blockWidth:e.size[0],blockHeight:e.size[1],pyramidBlockWidth:e.size[0],pyramidBlockHeight:e.size[1],pyramidResolutions:v,compression:e.format,origin:e.origin,firstPyramidLevel:1,maximumPyramidLevel:p,tileInfo:e,blockBoundary:g}),this._fixGCSShift(o),this._set("rasterInfo",o)}if(r.capabilities.toLowerCase().indexOf("tilemap")>-1){const e={tileInfo:o.storageInfo.tileInfo,parsedUrl:l(this.url),url:this.url,tileServers:[],type:"tile"};this._tilemapCache=new u({layer:e})}}async fetchRawTile(e,t,i,s={}){if(this._slices&&null==s.sliceId)return null;const{storageInfo:a,extent:l}=this.rasterInfo,r=a.maximumPyramidLevel-e+this._levelOffset,n=`${this.url}/tile/${r}/${t}/${i}`,o=this._slices?{sliceId:s.sliceId||0}:null,{data:m}=await this.request(n,{query:o,responseType:"array-buffer",signal:s.signal});if(!m)return null;const c=await this.decodePixelBlock(m,{width:a.tileInfo.size[0],height:a.tileInfo.size[1],planes:null,pixelType:null,isPoint:"Elevation"===this.tileType}),h=a.blockBoundary[e];if("jpg"!==a.compression||i>h.minCol&&i<h.maxCol&&t>h.minRow&&t<h.maxRow)return c;const{origin:u,blockWidth:f,blockHeight:d}=a,{x:y,y:x}=this.getPyramidPixelSize(e),g=Math.round((l.xmin-u.x)/y)%f,v=Math.round((l.xmax-u.x)/y)%f||f,S=Math.round((u.y-l.ymax)/x)%d,w=Math.round((u.y-l.ymin)/x)%d||d,I=i===h.minCol?g:0,b=t===h.minRow?S:0,j=i===h.maxCol?v:f,T=t===h.maxRow?w:d;return p(c,{x:I,y:b},{width:j-I,height:T-b}),c}getSliceIndex(e){if(!this._slices||s(e)||0===e.length)return null;const t=e;for(let i=0;i<this._slices.length;i++){const e=this._slices[i].multidimensionalDefinition;if(e.length===t.length&&!e.some((e=>{const i=t.find((t=>e.variableName===t.variableName&&t.dimensionName===e.dimensionName));if(!i)return!0;return(Array.isArray(e.values[0])?`${e.values[0][0]}-${e.values[0][1]}`:e.values[0])!==(Array.isArray(i.values[0])?`${i.values[0][0]}-${i.values[0][1]}`:i.values[0])})))return i}return null}async fetchVariableStatisticsHistograms(e,t){const i=this.request(this.url+"/statistics",{query:{variable:e,f:"json"},signal:t}).then((e=>{var t;return null==(t=e.data)?void 0:t.statistics})),s=this.request(this.url+"/histograms",{query:{variable:e,f:"json"},signal:t}).then((e=>{var t;return null==(t=e.data)?void 0:t.histograms})),a=await Promise.all([i,s]);return a[0]&&a[0].forEach((e=>{e.avg=e.mean,e.stddev=e.standardDeviation})),{statistics:a[0]||null,histograms:a[1]||null}}async computeBestPyramidLevelForLocation(e,t={}){if(!this._tilemapCache)return 0;let i=this.identifyPixelLocation(e,0,a(t.datumTransformation));if(null===i)return null;let s=0;const{maximumPyramidLevel:l}=this.rasterInfo.storageInfo;let r=l-s+this._levelOffset;const n=i.srcLocation;for(;r>=0;){try{if("available"===await this._tilemapCache.fetchAvailability(r,i.row,i.col,t))break}catch{}if(r--,s++,i=this.identifyPixelLocation(n,s,a(t.datumTransformation)),null===i)return null}return-1===r||null==i?null:s}async _fetchRasterInfo(e){const t=this.sourceJSON,i=Math.ceil((t.extent.xmax-t.extent.xmin)/t.pixelSizeX-.1),s=Math.ceil((t.extent.ymax-t.extent.ymin)/t.pixelSizeY-.1),a=x.fromJSON(t.spatialReference||t.extent.spatialReference);if("Map"===this.tileType)return new o({width:i,height:s,bandCount:3,extent:g.fromJSON(t.extent),spatialReference:a,pixelSize:new v({x:t.pixelSizeX,y:t.pixelSizeY,spatialReference:a}),pixelType:"u8",statistics:null,keyProperties:{DataType:"processed"}});const{slice:l,signal:r}=e,n=!!t.hasRasterAttributeTable&&this.request(this.url+"/rasterAttributeTable",{query:{slice:l,f:"json"},signal:r}).then((e=>y.fromJSON(e.data))).catch((()=>null)),m=!!t.hasColormap&&this.request(this.url+"/colormap",{query:{slice:l,f:"json"},signal:r}).then((e=>{var t;return null==(t=e.data)?void 0:t.colormap})),c=!!t.hasHistograms&&this.request(this.url+"/histograms",{query:{slice:l,f:"json"},signal:r}).then((e=>{var t;return null==(t=e.data)?void 0:t.histograms})),h=this.request(this.url+"/keyProperties",{query:{f:"json"},signal:r}).then((e=>e.data)).catch((()=>{})),u=!!t.hasMultidimensions&&this._fetchMultidimensionalInfo(),f=!!t.hasMultidimensions&&this.request(this.url+"/slices",{query:{f:"json"},signal:r}).then((e=>e.data&&e.data.slices)).catch((()=>{}));return Promise.all([n,m,c,h,u,f]).then((e=>{let l=null;if(t.minValues&&t.minValues.length===t.bandCount){l=[];for(let e=0;e<t.minValues.length;e++)l.push({min:t.minValues[e],max:t.maxValues[e],avg:t.meanValues[e],stddev:t.stdvValues[e]})}return this._slices=e[5]||null,new o({width:i,height:s,bandCount:t.bandCount,extent:g.fromJSON(t.extent),spatialReference:a,pixelSize:new v({x:t.pixelSizeX,y:t.pixelSizeY,spatialReference:a}),pixelType:t.pixelType.toLowerCase(),statistics:l,attributeTable:e[0]||null,colormap:e[1]||null,histograms:e[2]||null,keyProperties:e[3]||{},multidimensionalInfo:e[4]||null})}))}async _fetchMultidimensionalInfo(e){var t;const i=await this.request(this.url+"/multidimensionalInfo",{query:{f:"json"},signal:e}).then((e=>{var t;return null==(t=e.data)?void 0:t.multidimensionalInfo}));return null!=(t=i.variables)&&t.length&&i.variables.forEach((e=>{var t;null!=(t=e.statistics)&&t.length&&e.statistics.forEach((e=>{e.avg=e.mean,e.stddev=e.standardDeviation}))})),i}_fixScaleInServiceInfo(){const{sourceJSON:e}=this;e.minScale&&e.minScale<0&&(e.minScale=0),e.maxScale&&e.maxScale<0&&(e.maxScale=0)}_fixGCSShift(e){const{extent:t,spatialReference:i}=e;0===t.xmin&&360===t.xmax&&i.wkid&&i.isGeographic&&(e.nativeExtent=e.extent,e.transform=new d,e.extent=e.transform.forwardTransform(t))}};e([r({type:String,json:{write:!0}})],S.prototype,"datasetFormat",void 0),e([r()],S.prototype,"tileType",void 0),S=e([n("geoscene.layers.support.rasterDatasets.ImageServerRaster")],S);const w=S;export{w as default};

package/layers/support/rasterDatasets/InMemoryRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{isSome as s}from"../../../core/maybe.js";import{eachAlways as r}from"../../../core/promiseUtils.js";import{property as i}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as o}from"../../../core/accessorSupport/decorators/subclass.js";import a from"../RasterInfo.js";import m from"./BaseRaster.js";import{split as n,estimateStatisticsHistograms as l}from"../rasterFunctions/pixelUtils.js";import p from"../../../geometry/Extent.js";import c from"../../../geometry/SpatialReference.js";let h=class extends m{constructor(){super(...arguments),this.datasetFormat="MEMORY"}async open(t){var e;await this.init();const{pixelBlock:s,statistics:r,histograms:i,name:o,keyProperties:m,nativeExtent:n,transform:l}=this.data,{width:h,height:d,pixelType:f}=s,u=this.data.extent||new p({xmin:-.5,ymin:.5,xmax:h-.5,ymax:d-.5,spatialReference:new c({wkid:3857})}),y=null!=(e=this.data.isPseudoSpatialReference)?e:!this.data.extent,x={x:u.width/h,y:u.height/d},v=new a({width:h,height:d,pixelType:f,extent:u,nativeExtent:n,transform:l,pixelSize:x,spatialReference:u.spatialReference,bandCount:3,keyProperties:m||{},statistics:r,isPseudoSpatialReference:y,histograms:i});this.createRemoteDatasetStorageInfo(v,512,512),this._set("rasterInfo",v),this.updateTileInfo(),await this._buildInMemoryRaster(s,{width:512,height:512},t),this.datasetName=o,this.url="/InMemory/"+o}fetchRawTile(t,e,s,r={}){const i=this._pixelBlockTiles.get(`${t}/${e}/${s}`);return Promise.resolve(i)}async _buildInMemoryRaster(t,i,o){const a=this.rasterInfo.storageInfo.maximumPyramidLevel,m=this.rasterJobHandler?this.rasterJobHandler.split({pixelBlock:t,tileSize:i,maximumPyramidLevel:a},o):Promise.resolve(n(t,i,a)),p=s(this.rasterInfo.statistics),c=s(this.rasterInfo.histograms),h=p&&c?Promise.resolve({statistics:null,histograms:null}):this.rasterJobHandler?this.rasterJobHandler.estimateStatisticsHistograms({pixelBlock:t},o):Promise.resolve(l(t)),d=await r([m,h]);if(!d[0].value&&d[1].value)throw new e("inmemory-raster:open","failed to build in memory raster");var f,u;(this._pixelBlockTiles=d[0].value,p)||(this.rasterInfo.statistics=null==(f=d[1].value)?void 0:f.statistics);c&&(this.rasterInfo.histograms=null==(u=d[1].value)?void 0:u.histograms)}};t([i({type:String,json:{write:!0}})],h.prototype,"datasetFormat",void 0),t([i()],h.prototype,"data",void 0),h=t([o("esri.layers.support.rasterDatasets.InMemoryRaster")],h);const d=h;export{d as default};
+import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{isSome as s}from"../../../core/maybe.js";import{eachAlways as r}from"../../../core/promiseUtils.js";import{property as i}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as o}from"../../../core/accessorSupport/decorators/subclass.js";import a from"../RasterInfo.js";import m from"./BaseRaster.js";import{split as n,estimateStatisticsHistograms as l}from"../rasterFunctions/pixelUtils.js";import p from"../../../geometry/Extent.js";import c from"../../../geometry/SpatialReference.js";let h=class extends m{constructor(){super(...arguments),this.datasetFormat="MEMORY"}async open(t){var e;await this.init();const{pixelBlock:s,statistics:r,histograms:i,name:o,keyProperties:m,nativeExtent:n,transform:l}=this.data,{width:h,height:d,pixelType:f}=s,u=this.data.extent||new p({xmin:-.5,ymin:.5,xmax:h-.5,ymax:d-.5,spatialReference:new c({wkid:3857})}),y=null!=(e=this.data.isPseudoSpatialReference)?e:!this.data.extent,x={x:u.width/h,y:u.height/d},v=new a({width:h,height:d,pixelType:f,extent:u,nativeExtent:n,transform:l,pixelSize:x,spatialReference:u.spatialReference,bandCount:3,keyProperties:m||{},statistics:r,isPseudoSpatialReference:y,histograms:i});this.createRemoteDatasetStorageInfo(v,512,512),this._set("rasterInfo",v),this.updateTileInfo(),await this._buildInMemoryRaster(s,{width:512,height:512},t),this.datasetName=o,this.url="/InMemory/"+o}fetchRawTile(t,e,s,r={}){const i=this._pixelBlockTiles.get(`${t}/${e}/${s}`);return Promise.resolve(i)}async _buildInMemoryRaster(t,i,o){const a=this.rasterInfo.storageInfo.maximumPyramidLevel,m=this.rasterJobHandler?this.rasterJobHandler.split({pixelBlock:t,tileSize:i,maximumPyramidLevel:a},o):Promise.resolve(n(t,i,a)),p=s(this.rasterInfo.statistics),c=s(this.rasterInfo.histograms),h=p&&c?Promise.resolve({statistics:null,histograms:null}):this.rasterJobHandler?this.rasterJobHandler.estimateStatisticsHistograms({pixelBlock:t},o):Promise.resolve(l(t)),d=await r([m,h]);if(!d[0].value&&d[1].value)throw new e("inmemory-raster:open","failed to build in memory raster");var f,u;(this._pixelBlockTiles=d[0].value,p)||(this.rasterInfo.statistics=null==(f=d[1].value)?void 0:f.statistics);c&&(this.rasterInfo.histograms=null==(u=d[1].value)?void 0:u.histograms)}};t([i({type:String,json:{write:!0}})],h.prototype,"datasetFormat",void 0),t([i()],h.prototype,"data",void 0),h=t([o("geoscene.layers.support.rasterDatasets.InMemoryRaster")],h);const d=h;export{d as default};

package/layers/support/rasterDatasets/MRFRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{unwrap as r,isSome as s}from"../../../core/maybe.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import a from"../PixelBlock.js";import n from"../RasterInfo.js";import l from"../RasterStorageInfo.js";import f from"./BaseRaster.js";import{parseSpatialReference as c,parsePAMInfo as p}from"./pamParser.js";import{getElement as m,getElementValue as h}from"./xmlUtilities.js";import{isPlatformLittleEndian as u}from"../rasterFormats/utils.js";import{estimateStatisticsFromHistograms as g}from"../rasterFunctions/pixelUtils.js";import d from"../../../geometry/SpatialReference.js";import y from"../../../geometry/Extent.js";import x from"../../../geometry/Point.js";const w=new Map;w.set("Int8","s8"),w.set("UInt8","u8"),w.set("Int16","s16"),w.set("UInt16","u16"),w.set("Int32","s32"),w.set("UInt32","u32"),w.set("Float32","f32"),w.set("Float64","f32"),w.set("Double64","f32");const I=new Map;I.set("none",{blobExtension:".til",isOneSegment:!0,decoderFormat:"bip"}),I.set("lerc",{blobExtension:".lrc",isOneSegment:!1,decoderFormat:"lerc"}),I.set("deflate",{blobExtension:".pzp",isOneSegment:!0,decoderFormat:"deflate"}),I.set("jpeg",{blobExtension:".pjg",isOneSegment:!0,decoderFormat:"jpg"});let b=class extends f{constructor(){super(...arguments),this._files=null,this._storageIndex=null,this.datasetFormat="MRF"}async open(t){var e;await this.init(),this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const o=t?r(t.signal):null,i=await this.request(this.url,{responseType:"xml",signal:o}),{rasterInfo:a,files:n}=this._parseHeader(i.data);if(-1===(null==(e=this.ioConfig.skipExtensions)?void 0:e.indexOf("aux.xml"))){const e=await this._fetchAuxiliaryData(t);var l;if(null!=e)a.statistics=null!=(l=e.statistics)?l:a.statistics,a.histograms=e.histograms,e.histograms&&!s(a.statistics)&&(a.statistics=g(e.histograms))}this._set("rasterInfo",a),this._files=n;const f=await this.request(n.index,{responseType:"array-buffer",signal:o});this._storageIndex=this._parseIndex(f.data);const{blockWidth:c,blockHeight:p}=this.rasterInfo.storageInfo,m=this.rasterInfo.storageInfo.pyramidScalingFactor,{width:h,height:u}=this.rasterInfo,d=[],y=this._getBandSegmentCount();let x=0,w=-1;for(;x<this._storageIndex.length;){w++;const t=Math.ceil(h/c/m**w)-1,e=Math.ceil(u/p/m**w)-1;x+=(t+1)*(e+1)*y*4,d.push({maxRow:e,maxCol:t,minCol:0,minRow:0})}this.rasterInfo.storageInfo.blockBoundary=d,w>0&&(this.rasterInfo.storageInfo.firstPyramidLevel=1,this.rasterInfo.storageInfo.maximumPyramidLevel=w),this.updateTileInfo()}async fetchRawTile(t,e,r,o={}){const{blockWidth:i,blockHeight:n,blockBoundary:l}=this.rasterInfo.storageInfo,f=l[t];if(!f||f.maxRow<e||f.maxCol<r||f.minRow>e||f.minCol>r)return null;const{bandCount:c,pixelType:p}=this.rasterInfo,{ranges:m,actualTileWidth:h,actualTileHeight:u}=this._getTileLocation(t,e,r);if(!m||0===m.length)return null;if(0===m[0].from&&0===m[0].to){const t=new Uint8Array(i*n);return new a({width:i,height:n,pixels:null,mask:t,validPixelCount:0})}const{bandIds:g}=this.ioConfig,d=this._getBandSegmentCount(),y=[];let x=0;for(x=0;x<d;x++)(!g||g.indexOf[x]>-1)&&y.push(this.request(this._files.data,{range:{from:m[x].from,to:m[x].to},responseType:"array-buffer",signal:o.signal}));const w=await Promise.all(y),b=w.map((t=>t.data.byteLength)).reduce(((t,e)=>t+e)),A=new Uint8Array(b);let F=0;for(x=0;x<d;x++)A.set(new Uint8Array(w[x].data),F),F+=w[x].data.byteLength;const _=I.get(this.rasterInfo.storageInfo.compression).decoderFormat,R=await this.decodePixelBlock(A.buffer,{width:i,height:n,format:_,planes:(null==g?void 0:g.length)||c,pixelType:p});if(s(this.rasterInfo.noDataValue)&&"lerc"!==_&&!R.mask){const t=this.rasterInfo.noDataValue[0];if(null!=t){const e=R.width*R.height,r=new Uint8Array(e);if(Math.abs(t)>1e24)for(x=0;x<e;x++)Math.abs((R.pixels[0][x]-t)/t)>1e-6&&(r[x]=1);else for(x=0;x<e;x++)R.pixels[0][x]!==t&&(r[x]=1);R.mask=r}}let j=0,k=0;if(h!==i||u!==n){let t=R.mask;if(t)for(x=0;x<n;x++)if(k=x*i,x<u)for(j=h;j<i;j++)t[k+j]=0;else for(j=0;j<i;j++)t[k+j]=0;else for(t=new Uint8Array(i*n),R.mask=t,x=0;x<u;x++)for(k=x*i,j=0;j<h;j++)t[k+j]=1}return R}_parseIndex(t){if(t.byteLength%16>0)throw"invalid array buffer must be multiples of 16";let e,r,s,o,i,a;if(u){for(r=new Uint8Array(t),o=new ArrayBuffer(t.byteLength),s=new Uint8Array(o),i=0;i<t.byteLength/4;i++)for(a=0;a<4;a++)s[4*i+a]=r[4*i+3-a];e=new Uint32Array(o)}else e=new Uint32Array(t);return e}_getBandSegmentCount(){return I.get(this.rasterInfo.storageInfo.compression).isOneSegment?1:this.rasterInfo.bandCount}_getTileLocation(t,e,r){const{blockWidth:s,blockHeight:o,pyramidScalingFactor:i}=this.rasterInfo.storageInfo,{width:a,height:n}=this.rasterInfo,l=this._getBandSegmentCount();let f,c,p,m=0,h=0;for(p=0;p<t;p++)h=i**p,f=Math.ceil(a/s/h),c=Math.ceil(n/o/h),m+=f*c;h=i**t,f=Math.ceil(a/s/h),c=Math.ceil(n/o/h),m+=e*f+r,m*=4*l;const u=this._storageIndex.subarray(m,m+4*l);let g=0,d=0;const y=[];for(let x=0;x<l;x++)g=u[4*x+0]*2**32+u[4*x+1],d=g+u[4*x+2]*2**32+u[4*x+3],y.push({from:g,to:d});return{ranges:y,actualTileWidth:r<f-1?s:Math.ceil(a/h)-s*(f-1),actualTileHeight:e<c-1?o:Math.ceil(n/h)-o*(c-1)}}_parseHeader(t){const r=m(t,"MRF_META/Raster");if(!r)throw new e("mrf:open","not a valid MRF format");const s=m(r,"Size"),o=parseInt(s.getAttribute("x"),10),i=parseInt(s.getAttribute("y"),10),a=parseInt(s.getAttribute("c"),10),f=(h(r,"Compression")||"none").toLowerCase();if(!I.has(f))throw new e("mrf:open","currently does not support compression "+f);const p=h(r,"DataType")||"UInt8",u=w.get(p);if(null==u)throw new e("mrf:open","currently does not support pixel type "+p);const g=m(r,"PageSize"),b=parseInt(g.getAttribute("x"),10),A=parseInt(g.getAttribute("y"),10),F=m(r,"DataValues");let _,R;F&&(R=F.getAttribute("NoData"),null!=R&&(_=R.trim().split(" ").map((t=>parseFloat(t)))));if(m(t,"MRF_META/CachedSource"))throw new e("mrf:open","currently does not support MRF referencing other data files");const j=m(t,"MRF_META/GeoTags"),k=m(j,"BoundingBox");let M,S=!1;if(null!=k){const t=parseFloat(k.getAttribute("minx")),e=parseFloat(k.getAttribute("miny")),r=parseFloat(k.getAttribute("maxx")),s=parseFloat(k.getAttribute("maxy")),o=h(j,"Projection")||"";let i;if("LOCAL_CS[]"!==o)if(o.toLowerCase().startsWith("epsg:")){const t=Number(o.slice(5));isNaN(t)||0===t||(i=new d({wkid:t}))}else i=c(o);else S=!0,i=new d({wkid:3857});M=new y(t,e,r,s),M.spatialReference=i}else S=!0,M=new y({xmin:-.5,ymin:.5-i,xmax:o-.5,ymax:.5,spatialReference:new d({wkid:3857})});const T=m(t,"MRF_META/Rsets"),C=parseInt(T&&T.getAttribute("scale")||"2",10),U=M.spatialReference,B=new l({origin:new x({x:M.xmin,y:M.ymax,spatialReference:U}),blockWidth:b,blockHeight:A,pyramidBlockWidth:b,pyramidBlockHeight:A,compression:f,pyramidScalingFactor:C}),v=new x({x:M.width/o,y:M.height/i,spatialReference:U}),L=new n({width:o,height:i,extent:M,isPseudoSpatialReference:S,spatialReference:U,bandCount:a,pixelType:u,pixelSize:v,noDataValue:_,storageInfo:B}),E=h(t,"datafile"),P=h(t,"IndexFile");return{rasterInfo:L,files:{mrf:this.url,index:P||this.url.replace(".mrf",".idx"),data:E||this.url.replace(".mrf",I.get(f).blobExtension)}}}async _fetchAuxiliaryData(t){try{const{data:e}=await this.request(this.url+".aux.xml",{responseType:"xml",signal:null==t?void 0:t.signal});return p(e)}catch{return null}}};t([o()],b.prototype,"_files",void 0),t([o()],b.prototype,"_storageIndex",void 0),t([o({type:String,json:{write:!0}})],b.prototype,"datasetFormat",void 0),b=t([i("esri.layers.support.rasterIO.MRFRaster")],b);const A=b;export{A as default};
+import{_ as t}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import e from"../../../core/Error.js";import{unwrap as r,isSome as s}from"../../../core/maybe.js";import{property as o}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as i}from"../../../core/accessorSupport/decorators/subclass.js";import a from"../PixelBlock.js";import n from"../RasterInfo.js";import l from"../RasterStorageInfo.js";import f from"./BaseRaster.js";import{parseSpatialReference as c,parsePAMInfo as p}from"./pamParser.js";import{getElement as m,getElementValue as h}from"./xmlUtilities.js";import{isPlatformLittleEndian as u}from"../rasterFormats/utils.js";import{estimateStatisticsFromHistograms as g}from"../rasterFunctions/pixelUtils.js";import d from"../../../geometry/SpatialReference.js";import y from"../../../geometry/Extent.js";import x from"../../../geometry/Point.js";const w=new Map;w.set("Int8","s8"),w.set("UInt8","u8"),w.set("Int16","s16"),w.set("UInt16","u16"),w.set("Int32","s32"),w.set("UInt32","u32"),w.set("Float32","f32"),w.set("Float64","f32"),w.set("Double64","f32");const I=new Map;I.set("none",{blobExtension:".til",isOneSegment:!0,decoderFormat:"bip"}),I.set("lerc",{blobExtension:".lrc",isOneSegment:!1,decoderFormat:"lerc"}),I.set("deflate",{blobExtension:".pzp",isOneSegment:!0,decoderFormat:"deflate"}),I.set("jpeg",{blobExtension:".pjg",isOneSegment:!0,decoderFormat:"jpg"});let b=class extends f{constructor(){super(...arguments),this._files=null,this._storageIndex=null,this.datasetFormat="MRF"}async open(t){var e;await this.init(),this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const o=t?r(t.signal):null,i=await this.request(this.url,{responseType:"xml",signal:o}),{rasterInfo:a,files:n}=this._parseHeader(i.data);if(-1===(null==(e=this.ioConfig.skipExtensions)?void 0:e.indexOf("aux.xml"))){const e=await this._fetchAuxiliaryData(t);var l;if(null!=e)a.statistics=null!=(l=e.statistics)?l:a.statistics,a.histograms=e.histograms,e.histograms&&!s(a.statistics)&&(a.statistics=g(e.histograms))}this._set("rasterInfo",a),this._files=n;const f=await this.request(n.index,{responseType:"array-buffer",signal:o});this._storageIndex=this._parseIndex(f.data);const{blockWidth:c,blockHeight:p}=this.rasterInfo.storageInfo,m=this.rasterInfo.storageInfo.pyramidScalingFactor,{width:h,height:u}=this.rasterInfo,d=[],y=this._getBandSegmentCount();let x=0,w=-1;for(;x<this._storageIndex.length;){w++;const t=Math.ceil(h/c/m**w)-1,e=Math.ceil(u/p/m**w)-1;x+=(t+1)*(e+1)*y*4,d.push({maxRow:e,maxCol:t,minCol:0,minRow:0})}this.rasterInfo.storageInfo.blockBoundary=d,w>0&&(this.rasterInfo.storageInfo.firstPyramidLevel=1,this.rasterInfo.storageInfo.maximumPyramidLevel=w),this.updateTileInfo()}async fetchRawTile(t,e,r,o={}){const{blockWidth:i,blockHeight:n,blockBoundary:l}=this.rasterInfo.storageInfo,f=l[t];if(!f||f.maxRow<e||f.maxCol<r||f.minRow>e||f.minCol>r)return null;const{bandCount:c,pixelType:p}=this.rasterInfo,{ranges:m,actualTileWidth:h,actualTileHeight:u}=this._getTileLocation(t,e,r);if(!m||0===m.length)return null;if(0===m[0].from&&0===m[0].to){const t=new Uint8Array(i*n);return new a({width:i,height:n,pixels:null,mask:t,validPixelCount:0})}const{bandIds:g}=this.ioConfig,d=this._getBandSegmentCount(),y=[];let x=0;for(x=0;x<d;x++)(!g||g.indexOf[x]>-1)&&y.push(this.request(this._files.data,{range:{from:m[x].from,to:m[x].to},responseType:"array-buffer",signal:o.signal}));const w=await Promise.all(y),b=w.map((t=>t.data.byteLength)).reduce(((t,e)=>t+e)),A=new Uint8Array(b);let F=0;for(x=0;x<d;x++)A.set(new Uint8Array(w[x].data),F),F+=w[x].data.byteLength;const _=I.get(this.rasterInfo.storageInfo.compression).decoderFormat,R=await this.decodePixelBlock(A.buffer,{width:i,height:n,format:_,planes:(null==g?void 0:g.length)||c,pixelType:p});if(s(this.rasterInfo.noDataValue)&&"lerc"!==_&&!R.mask){const t=this.rasterInfo.noDataValue[0];if(null!=t){const e=R.width*R.height,r=new Uint8Array(e);if(Math.abs(t)>1e24)for(x=0;x<e;x++)Math.abs((R.pixels[0][x]-t)/t)>1e-6&&(r[x]=1);else for(x=0;x<e;x++)R.pixels[0][x]!==t&&(r[x]=1);R.mask=r}}let j=0,k=0;if(h!==i||u!==n){let t=R.mask;if(t)for(x=0;x<n;x++)if(k=x*i,x<u)for(j=h;j<i;j++)t[k+j]=0;else for(j=0;j<i;j++)t[k+j]=0;else for(t=new Uint8Array(i*n),R.mask=t,x=0;x<u;x++)for(k=x*i,j=0;j<h;j++)t[k+j]=1}return R}_parseIndex(t){if(t.byteLength%16>0)throw"invalid array buffer must be multiples of 16";let e,r,s,o,i,a;if(u){for(r=new Uint8Array(t),o=new ArrayBuffer(t.byteLength),s=new Uint8Array(o),i=0;i<t.byteLength/4;i++)for(a=0;a<4;a++)s[4*i+a]=r[4*i+3-a];e=new Uint32Array(o)}else e=new Uint32Array(t);return e}_getBandSegmentCount(){return I.get(this.rasterInfo.storageInfo.compression).isOneSegment?1:this.rasterInfo.bandCount}_getTileLocation(t,e,r){const{blockWidth:s,blockHeight:o,pyramidScalingFactor:i}=this.rasterInfo.storageInfo,{width:a,height:n}=this.rasterInfo,l=this._getBandSegmentCount();let f,c,p,m=0,h=0;for(p=0;p<t;p++)h=i**p,f=Math.ceil(a/s/h),c=Math.ceil(n/o/h),m+=f*c;h=i**t,f=Math.ceil(a/s/h),c=Math.ceil(n/o/h),m+=e*f+r,m*=4*l;const u=this._storageIndex.subarray(m,m+4*l);let g=0,d=0;const y=[];for(let x=0;x<l;x++)g=u[4*x+0]*2**32+u[4*x+1],d=g+u[4*x+2]*2**32+u[4*x+3],y.push({from:g,to:d});return{ranges:y,actualTileWidth:r<f-1?s:Math.ceil(a/h)-s*(f-1),actualTileHeight:e<c-1?o:Math.ceil(n/h)-o*(c-1)}}_parseHeader(t){const r=m(t,"MRF_META/Raster");if(!r)throw new e("mrf:open","not a valid MRF format");const s=m(r,"Size"),o=parseInt(s.getAttribute("x"),10),i=parseInt(s.getAttribute("y"),10),a=parseInt(s.getAttribute("c"),10),f=(h(r,"Compression")||"none").toLowerCase();if(!I.has(f))throw new e("mrf:open","currently does not support compression "+f);const p=h(r,"DataType")||"UInt8",u=w.get(p);if(null==u)throw new e("mrf:open","currently does not support pixel type "+p);const g=m(r,"PageSize"),b=parseInt(g.getAttribute("x"),10),A=parseInt(g.getAttribute("y"),10),F=m(r,"DataValues");let _,R;F&&(R=F.getAttribute("NoData"),null!=R&&(_=R.trim().split(" ").map((t=>parseFloat(t)))));if(m(t,"MRF_META/CachedSource"))throw new e("mrf:open","currently does not support MRF referencing other data files");const j=m(t,"MRF_META/GeoTags"),k=m(j,"BoundingBox");let M,S=!1;if(null!=k){const t=parseFloat(k.getAttribute("minx")),e=parseFloat(k.getAttribute("miny")),r=parseFloat(k.getAttribute("maxx")),s=parseFloat(k.getAttribute("maxy")),o=h(j,"Projection")||"";let i;if("LOCAL_CS[]"!==o)if(o.toLowerCase().startsWith("epsg:")){const t=Number(o.slice(5));isNaN(t)||0===t||(i=new d({wkid:t}))}else i=c(o);else S=!0,i=new d({wkid:3857});M=new y(t,e,r,s),M.spatialReference=i}else S=!0,M=new y({xmin:-.5,ymin:.5-i,xmax:o-.5,ymax:.5,spatialReference:new d({wkid:3857})});const T=m(t,"MRF_META/Rsets"),C=parseInt(T&&T.getAttribute("scale")||"2",10),U=M.spatialReference,B=new l({origin:new x({x:M.xmin,y:M.ymax,spatialReference:U}),blockWidth:b,blockHeight:A,pyramidBlockWidth:b,pyramidBlockHeight:A,compression:f,pyramidScalingFactor:C}),v=new x({x:M.width/o,y:M.height/i,spatialReference:U}),L=new n({width:o,height:i,extent:M,isPseudoSpatialReference:S,spatialReference:U,bandCount:a,pixelType:u,pixelSize:v,noDataValue:_,storageInfo:B}),E=h(t,"datafile"),P=h(t,"IndexFile");return{rasterInfo:L,files:{mrf:this.url,index:P||this.url.replace(".mrf",".idx"),data:E||this.url.replace(".mrf",I.get(f).blobExtension)}}}async _fetchAuxiliaryData(t){try{const{data:e}=await this.request(this.url+".aux.xml",{responseType:"xml",signal:null==t?void 0:t.signal});return p(e)}catch{return null}}};t([o()],b.prototype,"_files",void 0),t([o()],b.prototype,"_storageIndex",void 0),t([o({type:String,json:{write:!0}})],b.prototype,"datasetFormat",void 0),b=t([i("geoscene.layers.support.rasterIO.MRFRaster")],b);const A=b;export{A as default};

package/layers/support/rasterDatasets/TIFFRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{unwrap as r,isSome as i}from"../../../core/maybe.js";import{property as s}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as a}from"../../../core/accessorSupport/decorators/subclass.js";import n from"../RasterInfo.js";import o from"../RasterStorageInfo.js";import l from"./BaseRaster.js";import{parseSpatialReference as f,parsePAMInfo as u}from"./pamParser.js";import{parseSignature as m,getImageInfo as p,parseIFD as c,parseFieldValues as h,isBSQConfig as d}from"../rasterFormats/TiffDecoder.js";import y from"../rasterFormats/TiffTags.js";import{estimateStatisticsFromHistograms as g}from"../rasterFunctions/pixelUtils.js";import x from"../rasterTransforms/PolynomialTransform.js";import T from"../../../geometry/SpatialReference.js";import I from"../../../geometry/Extent.js";import w from"../../../geometry/Point.js";const _=function(e,t){const r=e.get(t);return r&&r.values},b=function(e,t){const r=e.get(t);return r&&r.values[0]};let E=class extends l{constructor(){super(...arguments),this._files=null,this._headerInfo=null,this._bufferSize=1048576,this.datasetFormat="TIFF"}async open(e){var s,a,l;await this.init();const u=e?r(e.signal):null,{data:c}=await this.request(this.url,{range:{from:0,to:this._bufferSize},responseType:"array-buffer",signal:u});if(!c)throw new t("tiffraster:open","failed to open url "+this.url);this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const{littleEndian:h,firstIFD:d,isBigTiff:y}=m(c),_=[];await this._readIFDs(_,c,h,d,0,y?8:4,u);const b=p(_),{width:E,height:S,tileWidth:F,tileHeight:v,planes:R,pixelType:k,compression:j,firstPyramidLevel:B,maximumPyramidLevel:D,pyramidBlockWidth:L,pyramidBlockHeight:z,tileBoundary:O,affine:P,metadata:H}=b,A=(null==(s=b.extent.spatialReference)?void 0:s.wkt)||(null==(a=b.extent.spatialReference)?void 0:a.wkid);let G=f(A),q=!1;null==G&&(q=!0,G=new T({wkid:3857}));const C=new I({...b.extent,spatialReference:G}),W=new w(C?{x:C.xmin,y:C.ymax,spatialReference:G}:{x:0,y:0}),U=new o({blockWidth:F,blockHeight:v,pyramidBlockWidth:L,pyramidBlockHeight:z,compression:j,origin:W,firstPyramidLevel:B,maximumPyramidLevel:D,blockBoundary:O}),Y=new w({x:(C.xmax-C.xmin)/E,y:(C.ymax-C.ymin)/S,spatialReference:G}),M=H?{BandProperties:H.bandProperties,DataType:H.dataType}:{},N=new n({width:E,height:S,bandCount:R,pixelType:k,compression:j,pixelSize:Y,storageInfo:U,spatialReference:G,isPseudoSpatialReference:q,keyProperties:M,extent:C,statistics:H?H.statistics:null});if(null!=P&&P.length&&(N.nativeExtent=new I({xmin:-.5,ymin:.5-S,xmax:E-.5,ymax:.5,spatialReference:G}),N.transform=new x({polynomialOrder:1,forwardCoefficients:[P[2]+P[0]/2,P[5]-P[3]/2,P[0],P[3],-P[1],-P[4]]}),N.extent=N.transform.forwardTransform(N.nativeExtent),N.pixelSize=new w({x:(C.xmax-C.xmin)/E,y:(C.ymax-C.ymin)/S,spatialReference:G}),U.origin.x=-.5,U.origin.y=.5),null==(l=this.ioConfig.skipExtensions)||!l.includes("aux.xml")){const t=await this._fetchAuxiliaryData(e);if(null!=t){var K;if(N.statistics=null!=(K=t.statistics)?K:N.statistics,N.histograms=t.histograms,t.histograms&&!i(N.statistics)&&(N.statistics=g(t.histograms)),t.transform&&!P){N.transform=t.transform,N.nativeExtent=N.extent;const e=N.transform.forwardTransform(N.nativeExtent);N.pixelSize=new w({x:(e.xmax-e.xmin)/E,y:(e.ymax-e.ymin)/S,spatialReference:G}),N.extent=e}N.spatialReference||(N.spatialReference=t.spatialReference)}}if(this._set("rasterInfo",N),this._headerInfo={littleEndian:h,isBigTiff:y,ifds:_,...b},!this._headerInfo.isSupported)throw new t("tiffraster:open","this tiff is not supported: "+this._headerInfo.message);this.updateTileInfo()}async fetchRawTile(e,t,r,i={}){var s;if(null==(s=this._headerInfo)||!s.isSupported||this.isBlockOutside(e,t,r))return null;const a=this._getTileLocation(e,t,r);if(!a)return null;const{ranges:n,actualTileWidth:o,actualTileHeight:l,ifd:f}=a,u=n.map((e=>this.request(this.url,{range:e,responseType:"array-buffer",signal:i.signal}))),m=await Promise.all(u),p=m.map((e=>e.data.byteLength)).reduce(((e,t)=>e+t)),c=1===m.length?m[0].data:new ArrayBuffer(p),h=[0],d=[0];if(m.length>1){const e=new Uint8Array(c);for(let t=0,r=0;t<m.length;t++){const i=m[t].data;e.set(new Uint8Array(i),r),h[t]=r,r+=i.byteLength,d[t]=i.byteLength}}const{blockWidth:y,blockHeight:g}=this.getBlockWidthHeight(e),x=await this.decodePixelBlock(c,{format:"tiff",customOptions:{headerInfo:this._headerInfo,ifd:f,offsets:h,sizes:d},width:y,height:g,planes:null,pixelType:null});let T,I,w;if(o!==y||l!==g){let e=x.mask;if(e)for(T=0;T<g;T++)if(w=T*y,T<l)for(I=o;I<y;I++)e[w+I]=0;else for(I=0;I<y;I++)e[w+I]=0;else for(e=new Uint8Array(y*g),x.mask=e,T=0;T<l;T++)for(w=T*y,I=0;I<o;I++)e[w+I]=1}return x}async _readIFDs(e,t,r,i,s,a=4,n){if(!i)return null;if(i>=t.byteLength||i<0){t=(await this.request(this.url,{range:{from:i+s,to:i+s+this._bufferSize},responseType:"array-buffer",signal:n})).data,s=i+s,i=0}const o=await this._readIFD(t,r,i,s,y.TIFF_TAGS,a,n);if(e.push(o.ifd),!o.nextIFD)return null;await this._readIFDs(e,t,r,o.nextIFD-s,s,a,n)}async _readIFD(e,t,r,i,s=y.TIFF_TAGS,a=4,n){if(!e)return null;const o=c(e,t,r,i,s,a);if(o.success){const r=[];if(o.ifd.forEach((e=>{e.values||r.push(e)})),r.length>0){const s=r.map((e=>e.offlineOffsetSize)),a=Math.min.apply(null,s.map((e=>e[0])));if(Math.min.apply(null,s.map((e=>e[0]+e[1])))-a<=this._bufferSize){const{data:s}=await this.request(this.url,{range:{from:a,to:a+this._bufferSize},responseType:"array-buffer",signal:n});e=s,i=a,r.forEach((r=>h(e,t,r,i)))}}if(o.ifd.has("GEOKEYDIRECTORY")){const r=o.ifd.get("GEOKEYDIRECTORY"),s=r.values;if(s&&s.length>4){const a=s[0]+"."+s[1]+"."+s[2],o=await this._readIFD(e,t,r.valueOffset+6-i,i,y.GEO_KEYS,2,n);r.data=o.ifd,r.data&&r.data.set("GEOTIFFVersion",{id:0,type:2,valueCount:1,valueOffset:null,values:[a]})}}return o}if(o.requiredBufferSize&&o.requiredBufferSize!==e.byteLength){const r=await this.request(this.url,{range:{from:i,to:i+o.requiredBufferSize+4},responseType:"array-buffer",signal:n});return(e=r.data).byteLength<o.requiredBufferSize?null:this._readIFD(e,t,0,i,y.TIFF_TAGS,4,n)}}_getTileLocation(e,t,r){const{firstPyramidLevel:i,blockBoundary:s}=this.rasterInfo.storageInfo,a=0===e?0:e-(i-1),n=this._headerInfo.ifds[a];if(!n)return null;const o=d(n,this._headerInfo),l=_(n,"TILEOFFSETS");if(void 0===l)return null;const f=_(n,"TILEBYTECOUNTS"),{minRow:u,minCol:m,maxRow:p,maxCol:c}=s[a];if(t>p||r>c||t<u||r<m)return null;const h=b(n,"IMAGEWIDTH"),y=b(n,"IMAGELENGTH"),g=b(n,"TILEWIDTH"),x=b(n,"TILELENGTH"),T=o?this.rasterInfo.bandCount:1,I=T*t*(c+1)+r,w=[{from:l[I],to:l[I+T-1]+f[I+T-1]-1}];if(o){let e=!0;for(let t=0;t<T;t++)if(l[I+t]+f[I+t]!==l[I+t+1]){e=!1;break}if(!e)for(let t=0;t<T;t++)w[t]={from:l[I+t],to:l[I+t]+f[I+t]-1}}const E=l[I],S=f[I];if(null==E||null==S)return null;return{ranges:w,ifd:n,actualTileWidth:r===c?h%g:g,actualTileHeight:t===p?y%x:x}}async _fetchAuxiliaryData(e){try{const{data:t}=await this.request(this.url+".aux.xml",{responseType:"xml",signal:null==e?void 0:e.signal});return u(t)}catch{return null}}};e([s()],E.prototype,"_files",void 0),e([s()],E.prototype,"_headerInfo",void 0),e([s()],E.prototype,"_bufferSize",void 0),e([s({type:String,json:{write:!0}})],E.prototype,"datasetFormat",void 0),E=e([a("esri.layers.support.rasterDatasets.TIFFRaster")],E);const S=E;export{S as default};
+import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{unwrap as r,isSome as i}from"../../../core/maybe.js";import{property as s}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as a}from"../../../core/accessorSupport/decorators/subclass.js";import n from"../RasterInfo.js";import o from"../RasterStorageInfo.js";import l from"./BaseRaster.js";import{parseSpatialReference as f,parsePAMInfo as u}from"./pamParser.js";import{parseSignature as m,getImageInfo as p,parseIFD as c,parseFieldValues as h,isBSQConfig as d}from"../rasterFormats/TiffDecoder.js";import y from"../rasterFormats/TiffTags.js";import{estimateStatisticsFromHistograms as g}from"../rasterFunctions/pixelUtils.js";import x from"../rasterTransforms/PolynomialTransform.js";import T from"../../../geometry/SpatialReference.js";import I from"../../../geometry/Extent.js";import w from"../../../geometry/Point.js";const _=function(e,t){const r=e.get(t);return r&&r.values},b=function(e,t){const r=e.get(t);return r&&r.values[0]};let E=class extends l{constructor(){super(...arguments),this._files=null,this._headerInfo=null,this._bufferSize=1048576,this.datasetFormat="TIFF"}async open(e){var s,a,l;await this.init();const u=e?r(e.signal):null,{data:c}=await this.request(this.url,{range:{from:0,to:this._bufferSize},responseType:"array-buffer",signal:u});if(!c)throw new t("tiffraster:open","failed to open url "+this.url);this.datasetName=this.url.slice(this.url.lastIndexOf("/")+1);const{littleEndian:h,firstIFD:d,isBigTiff:y}=m(c),_=[];await this._readIFDs(_,c,h,d,0,y?8:4,u);const b=p(_),{width:E,height:S,tileWidth:F,tileHeight:v,planes:R,pixelType:k,compression:j,firstPyramidLevel:B,maximumPyramidLevel:D,pyramidBlockWidth:L,pyramidBlockHeight:z,tileBoundary:O,affine:P,metadata:H}=b,A=(null==(s=b.extent.spatialReference)?void 0:s.wkt)||(null==(a=b.extent.spatialReference)?void 0:a.wkid);let G=f(A),q=!1;null==G&&(q=!0,G=new T({wkid:3857}));const C=new I({...b.extent,spatialReference:G}),W=new w(C?{x:C.xmin,y:C.ymax,spatialReference:G}:{x:0,y:0}),U=new o({blockWidth:F,blockHeight:v,pyramidBlockWidth:L,pyramidBlockHeight:z,compression:j,origin:W,firstPyramidLevel:B,maximumPyramidLevel:D,blockBoundary:O}),Y=new w({x:(C.xmax-C.xmin)/E,y:(C.ymax-C.ymin)/S,spatialReference:G}),M=H?{BandProperties:H.bandProperties,DataType:H.dataType}:{},N=new n({width:E,height:S,bandCount:R,pixelType:k,compression:j,pixelSize:Y,storageInfo:U,spatialReference:G,isPseudoSpatialReference:q,keyProperties:M,extent:C,statistics:H?H.statistics:null});if(null!=P&&P.length&&(N.nativeExtent=new I({xmin:-.5,ymin:.5-S,xmax:E-.5,ymax:.5,spatialReference:G}),N.transform=new x({polynomialOrder:1,forwardCoefficients:[P[2]+P[0]/2,P[5]-P[3]/2,P[0],P[3],-P[1],-P[4]]}),N.extent=N.transform.forwardTransform(N.nativeExtent),N.pixelSize=new w({x:(C.xmax-C.xmin)/E,y:(C.ymax-C.ymin)/S,spatialReference:G}),U.origin.x=-.5,U.origin.y=.5),null==(l=this.ioConfig.skipExtensions)||!l.includes("aux.xml")){const t=await this._fetchAuxiliaryData(e);if(null!=t){var K;if(N.statistics=null!=(K=t.statistics)?K:N.statistics,N.histograms=t.histograms,t.histograms&&!i(N.statistics)&&(N.statistics=g(t.histograms)),t.transform&&!P){N.transform=t.transform,N.nativeExtent=N.extent;const e=N.transform.forwardTransform(N.nativeExtent);N.pixelSize=new w({x:(e.xmax-e.xmin)/E,y:(e.ymax-e.ymin)/S,spatialReference:G}),N.extent=e}N.spatialReference||(N.spatialReference=t.spatialReference)}}if(this._set("rasterInfo",N),this._headerInfo={littleEndian:h,isBigTiff:y,ifds:_,...b},!this._headerInfo.isSupported)throw new t("tiffraster:open","this tiff is not supported: "+this._headerInfo.message);this.updateTileInfo()}async fetchRawTile(e,t,r,i={}){var s;if(null==(s=this._headerInfo)||!s.isSupported||this.isBlockOutside(e,t,r))return null;const a=this._getTileLocation(e,t,r);if(!a)return null;const{ranges:n,actualTileWidth:o,actualTileHeight:l,ifd:f}=a,u=n.map((e=>this.request(this.url,{range:e,responseType:"array-buffer",signal:i.signal}))),m=await Promise.all(u),p=m.map((e=>e.data.byteLength)).reduce(((e,t)=>e+t)),c=1===m.length?m[0].data:new ArrayBuffer(p),h=[0],d=[0];if(m.length>1){const e=new Uint8Array(c);for(let t=0,r=0;t<m.length;t++){const i=m[t].data;e.set(new Uint8Array(i),r),h[t]=r,r+=i.byteLength,d[t]=i.byteLength}}const{blockWidth:y,blockHeight:g}=this.getBlockWidthHeight(e),x=await this.decodePixelBlock(c,{format:"tiff",customOptions:{headerInfo:this._headerInfo,ifd:f,offsets:h,sizes:d},width:y,height:g,planes:null,pixelType:null});let T,I,w;if(o!==y||l!==g){let e=x.mask;if(e)for(T=0;T<g;T++)if(w=T*y,T<l)for(I=o;I<y;I++)e[w+I]=0;else for(I=0;I<y;I++)e[w+I]=0;else for(e=new Uint8Array(y*g),x.mask=e,T=0;T<l;T++)for(w=T*y,I=0;I<o;I++)e[w+I]=1}return x}async _readIFDs(e,t,r,i,s,a=4,n){if(!i)return null;if(i>=t.byteLength||i<0){t=(await this.request(this.url,{range:{from:i+s,to:i+s+this._bufferSize},responseType:"array-buffer",signal:n})).data,s=i+s,i=0}const o=await this._readIFD(t,r,i,s,y.TIFF_TAGS,a,n);if(e.push(o.ifd),!o.nextIFD)return null;await this._readIFDs(e,t,r,o.nextIFD-s,s,a,n)}async _readIFD(e,t,r,i,s=y.TIFF_TAGS,a=4,n){if(!e)return null;const o=c(e,t,r,i,s,a);if(o.success){const r=[];if(o.ifd.forEach((e=>{e.values||r.push(e)})),r.length>0){const s=r.map((e=>e.offlineOffsetSize)),a=Math.min.apply(null,s.map((e=>e[0])));if(Math.min.apply(null,s.map((e=>e[0]+e[1])))-a<=this._bufferSize){const{data:s}=await this.request(this.url,{range:{from:a,to:a+this._bufferSize},responseType:"array-buffer",signal:n});e=s,i=a,r.forEach((r=>h(e,t,r,i)))}}if(o.ifd.has("GEOKEYDIRECTORY")){const r=o.ifd.get("GEOKEYDIRECTORY"),s=r.values;if(s&&s.length>4){const a=s[0]+"."+s[1]+"."+s[2],o=await this._readIFD(e,t,r.valueOffset+6-i,i,y.GEO_KEYS,2,n);r.data=o.ifd,r.data&&r.data.set("GEOTIFFVersion",{id:0,type:2,valueCount:1,valueOffset:null,values:[a]})}}return o}if(o.requiredBufferSize&&o.requiredBufferSize!==e.byteLength){const r=await this.request(this.url,{range:{from:i,to:i+o.requiredBufferSize+4},responseType:"array-buffer",signal:n});return(e=r.data).byteLength<o.requiredBufferSize?null:this._readIFD(e,t,0,i,y.TIFF_TAGS,4,n)}}_getTileLocation(e,t,r){const{firstPyramidLevel:i,blockBoundary:s}=this.rasterInfo.storageInfo,a=0===e?0:e-(i-1),n=this._headerInfo.ifds[a];if(!n)return null;const o=d(n,this._headerInfo),l=_(n,"TILEOFFSETS");if(void 0===l)return null;const f=_(n,"TILEBYTECOUNTS"),{minRow:u,minCol:m,maxRow:p,maxCol:c}=s[a];if(t>p||r>c||t<u||r<m)return null;const h=b(n,"IMAGEWIDTH"),y=b(n,"IMAGELENGTH"),g=b(n,"TILEWIDTH"),x=b(n,"TILELENGTH"),T=o?this.rasterInfo.bandCount:1,I=T*t*(c+1)+r,w=[{from:l[I],to:l[I+T-1]+f[I+T-1]-1}];if(o){let e=!0;for(let t=0;t<T;t++)if(l[I+t]+f[I+t]!==l[I+t+1]){e=!1;break}if(!e)for(let t=0;t<T;t++)w[t]={from:l[I+t],to:l[I+t]+f[I+t]-1}}const E=l[I],S=f[I];if(null==E||null==S)return null;return{ranges:w,ifd:n,actualTileWidth:r===c?h%g:g,actualTileHeight:t===p?y%x:x}}async _fetchAuxiliaryData(e){try{const{data:t}=await this.request(this.url+".aux.xml",{responseType:"xml",signal:null==e?void 0:e.signal});return u(t)}catch{return null}}};e([s()],E.prototype,"_files",void 0),e([s()],E.prototype,"_headerInfo",void 0),e([s()],E.prototype,"_bufferSize",void 0),e([s({type:String,json:{write:!0}})],E.prototype,"datasetFormat",void 0),E=e([a("geoscene.layers.support.rasterDatasets.TIFFRaster")],E);const S=E;export{S as default};

package/layers/support/rasterDatasets/WCSRaster.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as i,isNone as s,unwrap as n}from"../../../core/maybe.js";import{isAbortError as o}from"../../../core/promiseUtils.js";import{property as r}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as a}from"../../../core/accessorSupport/decorators/subclass.js";import{isAxesOrderReversedForWkid as l}from"../../ogc/crsUtils.js";import c from"../DimensionalDefinition.js";import d from"./BaseRaster.js";import{parse as p}from"./multipartParser.js";import{parseCapabilities as h}from"./wcsCapabilitiesParser.js";import{standardizeInterpolations as f,parseCoverages as u}from"./wcsCoverageParser.js";import{clip as g}from"../rasterFunctions/pixelUtils.js";import m from"../../../geometry/Extent.js";const v=["nearest neighbor","bilinear","bicubic"],x=["nearest","linear","cubic"];let w=class extends d{constructor(){super(...arguments),this.datasetFormat="WCSServer"}async open(e){await this.init();const i=null==e?void 0:e.signal,s=await this._getCapabilities(i);if(this.capabilities=s,!this.version){var n;let e=null==(n=s.capabilitiesVersion)?void 0:n.slice(0,3);"2.0"===e||"1.1"===e||"1.0"===e?this.version=s.capabilitiesVersion:(e=s.supportedVersions.find((e=>"2.0.1"===e))||s.supportedVersions.find((e=>"2.0"===e.slice(0,3)))||s.supportedVersions.find((e=>"1.1"===e.slice(0,3)))||s.supportedVersions.find((e=>"1.0"===e.slice(0,3)))||"1.0.0",this.version=e)}null==this.coverageId&&(this.coverageId=s.coverages[0].id);const o=s.coverages.find((e=>e.id===this.coverageId));if(null==o)throw new t("wcsraster-open",`the coverageId ${this.coverageId} does not exist in capabilities`);const r=await this._describeCoverage(i);this.coverageInfo=r[0],"2.0"===this.version.slice(0,3)&&(this.coverageInfo.lonLatEnvelope=o.lonLatEnvelope,this.coverageInfo.supportedInterpolations=f(s.supportedInterpolations)),this.datasetName=this.coverageInfo.title;const{rasterInfo:a}=this.coverageInfo;if(this.createRemoteDatasetStorageInfo(a,512,512),this._set("rasterInfo",a),null==a.spatialReference)throw new t("wcsraster-open",`coverage without spatial referernce is not supported: ${this.coverageId}`);const{pixelType:l,bandCount:c}=await this._getPixelTypeAndBandCount(i);a.pixelType=l,1===a.bandCount&&c>1&&(a.bandCount=c),this.updateTileInfo()}async fetchRawTile(e,n,o,r={}){if(this.isBlockOutside(e,n,o))return null;const{nativePixelSize:a,spatialReference:l}=this.rasterInfo,c=2**e,d=a.x*c,p=a.y*c,{blockWidth:h,blockHeight:f}=this.getBlockWidthHeight(e),{origin:u}=this.rasterInfo.storageInfo.tileInfo,m=this.getTileExtent({x:d,y:p},n,o,u,l,[h,f]),v=this.rasterInfo.extent,x=m.xmax>v.xmax||m.ymin<v.ymin;let w=m,y=h,I=f;if(x&&(w=m.clone().intersection(v),i(w)&&(y=Math.floor((w.xmax-w.xmin)/d),I=Math.floor((w.ymax-w.ymin)/p),w.xmax=w.xmin+d*y,w.ymin=w.ymax-p*I)),s(w)||y<=1||I<=1)return null;const b=await this._getCoverage(w,y,I,c,r);if(!b)return null;const C=await this.decodePixelBlock(b,{width:y,height:I,planes:null,pixelType:null});if(C&&(C.width!==y||C.height!==I))throw new t("wcsraster-fetch",`the reponse has unexpected dimension width: ${C.width}, height: {pixelBlock.height}`);return x?g(C,{x:0,y:0},{width:f,height:f}):C}async _getCapabilities(e){const i={service:"WCS",request:"GetCapabilities"};this.version&&(i.version=this.version,i.acceptVersions=this.version);try{const{data:t}=await this.request(this.url,{query:i,responseType:"xml",signal:e});return h(t)}catch(s){if(!o(s))throw new t("wcslayer:open","wcs capabilities is not valid or supported");throw s}}async _describeCoverage(e){const i={service:"WCS",request:"DescribeCoverage",version:this.version},s=this.version.slice(0,3);"1.0"===s?i.coverage=this.coverageId:"1.1"===s?i.identifiers=this.coverageId:"2.0"===s&&(i.coverageId=this.coverageId);try{const{data:t}=await this.request(this.url,{query:i,responseType:"xml",signal:e});return u(t,this.version)}catch(n){if(!o(n))throw new t("wcslayer:open","wcs coverage description is not valid or supported");throw n}}async _getPixelTypeAndBandCount(e){const{pixelSize:s,extent:o,multidimensionalInfo:r}=this.rasterInfo,a=o.center,l=new m({xmin:a.x-s.x,xmax:a.x+s.x,ymin:a.y-s.y,ymax:a.y+s.y,spatialReference:o.spatialReference});let d;if(i(r)){const e=r.variables[0];d=[],e.dimensions.forEach((t=>{d.push(new c({variableName:e.name,dimensionName:t.name,values:t.hasRegularIntervals?t.extent[0]:t.values[0],isSlice:!0}))}))}const p=await this._getCoverage(l,2,2,1,{multidimensionalDefinition:d,signal:n(e)});if(!p)throw new t("wcsraster-open","unable to determine pixel type");const h=await this.decodePixelBlock(p,{width:2,height:2,planes:null,pixelType:null});return{pixelType:h.pixelType,bandCount:h.getPlaneCount()}}async _getCoverage(e,i,s,n,o){const{coverageDescription:r}=this.coverageInfo,{version:a}=r,l="2.0"===r.version?this._getCoverage201Parameters(e,i,s,n,o,r):"1.1"===r.version?this._getCoverage110Parameters(e,i,s,o,r):this._getCoverage100Parameters(e,i,s,o),c="2.0"===r.version?await this.request(this._constructWCS201Url(l),{signal:o.signal,responseType:"array-buffer"}):await this.request(this.url,{query:l,signal:o.signal,responseType:"array-buffer"});if("1.0"===a)return c.data;const d=p(c);if(d.isMultipart&&d.data){const e=d.data.find((e=>{var t;return(null==(t=e.contentType)?void 0:t.toLowerCase().includes("image"))&&null!=e.contentData}));return null==e?void 0:e.contentData}if(this.ioConfig.allowAnyMediaType)return c.data;throw new t("wcsraster:getcoverage","response is not a valid coverage, multipart response is expected")}_getInterpolationIndex(e){return-1===this.coverageInfo.supportedInterpolations.indexOf(e)||"nearest"===e?0:"bilinear"===e?1:"cubic"===e?2:0}_getCoverage100Parameters(e,t,i,s){const n=`${e.xmin},${e.ymin},${e.xmax},${e.ymax}`,o=e.spatialReference.wkid,r=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"GEOTIFF",{bandIds:a,interpolation:l}=s,c=this._getInterpolationIndex(l),d=a?a.map((e=>this.coverageInfo.bandNames[e])):null,p=v[c];return{service:"WCS",request:"GetCoverage",version:this.version,coverage:this.coverageId,format:r,crs:`EPSG:${o}`,bbox:n,width:t,height:i,interpolation:p,band:null==d?void 0:d.join(",")}}_getCoverage110Parameters(e,t,s,n,o){var r;const{multidimensionalDefinition:a,bandIds:c,interpolation:d}=n,p=e.spatialReference.wkid,h=`urn:ogc:def:crs:EPSG::${p}`,f=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"image/tiff",u=this._getInterpolationIndex(d),g=x[u],m=null==d||0===(null==(r=this.coverageInfo.supportedInterpolations)?void 0:r.indexOf(d)),v=o.domain.spatialDomain,w=v.origin.x<=v.envelope.xmin&&v.origin.y<=v.envelope.ymin,y=e.width/t,I=e.height/s*(w?1:-1),b=w?[e.xmin,e.ymin]:[e.xmin,e.ymax],C=v.useEPSGAxis&&l(p),$=C?`${b[1]},${b[0]}`:`${b[0]},${b[1]}`,S=C?`${I},0,0,${y}`:`${y},0,0,${I}`,j=y/2,T=e.xmin+j,_=e.xmax-j,O=Math.abs(I)/2,P=e.ymin+O,L=e.ymax-O,E=C?`${P},${T},${L},${_},${h}`:`${T},${P},${_},${L},${h}`,R=o.range.find((e=>e.axis.some((e=>e.identifier.toLowerCase().indexOf("band")>-1))));let q,k=R&&g&&c?m?`${R.identifier}[${R.axis[0].identifier}[${c.join(",")}]]`:`${R.identifier}:${g}[${R.axis[0].identifier}[${c.join(",")}]]`:null;if(i(a)&&a.length)for(let i=0;i<a.length;i++){let e=a[i].values;const t=a[i].dimensionName.toLowerCase(),s=a[i].variableName.toLowerCase();if(e.length>0)if(e[0]instanceof Array&&(e=e[0]),"stdtime"===t)q=e.map((e=>this._convertToISOTime(e))).join(",");else{const i=o.range.find((e=>e.identifier.toLowerCase()===s));if(i){const s=i.axis.find((e=>e.identifier.toLowerCase()===t));s&&(k=m?i.identifier+"["+s.identifier+"["+e.join(",")+"]]":i.identifier+":"+g+"["+s.identifier+"["+e.join(",")+"]]")}}}return{service:"WCS",request:"GetCoverage",version:this.version,identifier:this.coverageId,format:f,crs:`EPSG:${p}`,boundingbox:E,gridCS:"urn:ogc:def:cs:OGC:0.0:Grid2dSquareCS",gridType:"urn:ogc:def:method:WCS:1.1:2dGridIn2dCrs",gridOrigin:$,gridOffsets:S,gridBaseCRS:h,timeSequence:q,rangeSubset:k}}_convertToISOTime(e,t=!1){return(t?new Date(24*(e-25569)*60*60*1e3):new Date(e)).toISOString()}_getCoverage201Parameters(e,t,s,n,o,r){const{multidimensionalDefinition:a,interpolation:l}=o,c=this._getInterpolationIndex(l);let d=null;const{supportedInterpolations:p}=this.capabilities;if(null!=p&&p.length)switch(c){case 0:d=p.find((e=>e.includes("nearest")));break;case 1:d=p.find((e=>e.includes("linear")));break;case 2:d=p.find((e=>e.includes("cubic")||e.includes("quadratic")))}const h=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"image/tiff",{bandNames:f}=this.coverageInfo,{boundedBy:u,domainSet:g,rangeType:m}=r,v=u.isEastFirst?0:1,x=1-v,{axisLabels:w}=u,y=w[v],I=w[x],b=`http://www.opengis.net/def/crs/EPSG/0/${e.spatialReference.wkid}`,C=b,$=[];$.push(`${y}(${e.xmin},${e.xmax})`),$.push(`${I}(${e.ymin},${e.ymax})`);const S=[];if(w.length>2)for(let i=2;i<w.length;i++){const e=g.origin[i];if(w[i].toLowerCase().indexOf("time")>-1){let t=e.toString();u.uomLabels[i].toLowerCase().indexOf("ole")>-1&&(S.push(w[i]),t=this._convertToISOTime(e,!0)),$.push(w[i]+",http://www.opengis.net("+t+")")}else $.push(w[i]+",http://www.opengis.net("+e+")")}let j=null;if(i(a)&&a.length){const e=[];m.forEach((t=>t.forEach((t=>e.push(t.name)))));const t=[];for(let i=0;i<a.length;i++){const s=w.find((e=>e===a[i].dimensionName)),n=e.find((e=>e===a[i].variableName));if(-1===t.indexOf(n)&&t.push(n),s){let e=a[i].values;if(e.length>0){Array.isArray(e[0])&&(e=e[0]);let t="";t=s.toLowerCase().indexOf("time")>-1?e.map((e=>this._convertToISOTime(e))).join(","):e.join(",");const i=$.findIndex((e=>0===e.indexOf(s+",http://www.opengis.net")));-1===i&&$.push(s+",http://www.opengis.net("+t+")"),-1!==i&&-1===$[i].indexOf("("+t+")")&&$.splice(i,1,s+",http://www.opengis.net("+t+")")}}}t.length&&(j=t.join(","))}else if(null!=f&&f.length){j=(o.bandIds?o.bandIds.map((e=>f[e])):f).join(",")}const T=$.join("&subset="),_=!(r.domainSet.axisLabels.join("")===r.boundedBy.axisLabels.join(""))&&!1!==this.ioConfig.allowScaleFactor,O=_?null:`${y}(${t}),${I}(${s})`,P=_?1/n:null;return{service:"WCS",request:"GetCoverage",version:this.version,coverageId:this.coverageId,rangesubset:j,interpolation:d,scaleSize:O,scaleFactor:P,subset:T,format:h,mediaType:this.ioConfig.allowAnyMediaType?null:"multipart/related",outputcrs:b,subsettingcrs:C}}_constructWCS201Url(e){const t={...this.ioConfig.customFetchParameters,...e},i=[];Object.keys(t).forEach((e=>{const s=t[e];null!=s&&("subset"===e?s.split("&subset=").forEach((e=>{e&&i.push(`subset=${encodeURIComponent(e)}`)})):i.push(`${e}=${encodeURIComponent(s)}`))}));return`${encodeURI(this.url)}?${i.join("&")}`}};e([r({type:String,json:{write:!0}})],w.prototype,"datasetFormat",void 0),e([r()],w.prototype,"tileType",void 0),e([r({type:String,json:{write:!0}})],w.prototype,"version",void 0),e([r({type:String,json:{write:!0}})],w.prototype,"coverageId",void 0),w=e([a("esri.layers.support.rasterDatasets.ImageServerRaster")],w);const y=w;export{y as default};
+import{_ as e}from"../../../chunks/tslib.es6.js";import"../../../geometry.js";import t from"../../../core/Error.js";import{isSome as i,isNone as s,unwrap as n}from"../../../core/maybe.js";import{isAbortError as o}from"../../../core/promiseUtils.js";import{property as r}from"../../../core/accessorSupport/decorators/property.js";import"../../../core/arrayUtils.js";import"../../../core/has.js";import"../../../core/accessorSupport/ensureType.js";import{subclass as a}from"../../../core/accessorSupport/decorators/subclass.js";import{isAxesOrderReversedForWkid as l}from"../../ogc/crsUtils.js";import c from"../DimensionalDefinition.js";import d from"./BaseRaster.js";import{parse as p}from"./multipartParser.js";import{parseCapabilities as h}from"./wcsCapabilitiesParser.js";import{standardizeInterpolations as f,parseCoverages as u}from"./wcsCoverageParser.js";import{clip as g}from"../rasterFunctions/pixelUtils.js";import m from"../../../geometry/Extent.js";const v=["nearest neighbor","bilinear","bicubic"],x=["nearest","linear","cubic"];let w=class extends d{constructor(){super(...arguments),this.datasetFormat="WCSServer"}async open(e){await this.init();const i=null==e?void 0:e.signal,s=await this._getCapabilities(i);if(this.capabilities=s,!this.version){var n;let e=null==(n=s.capabilitiesVersion)?void 0:n.slice(0,3);"2.0"===e||"1.1"===e||"1.0"===e?this.version=s.capabilitiesVersion:(e=s.supportedVersions.find((e=>"2.0.1"===e))||s.supportedVersions.find((e=>"2.0"===e.slice(0,3)))||s.supportedVersions.find((e=>"1.1"===e.slice(0,3)))||s.supportedVersions.find((e=>"1.0"===e.slice(0,3)))||"1.0.0",this.version=e)}null==this.coverageId&&(this.coverageId=s.coverages[0].id);const o=s.coverages.find((e=>e.id===this.coverageId));if(null==o)throw new t("wcsraster-open",`the coverageId ${this.coverageId} does not exist in capabilities`);const r=await this._describeCoverage(i);this.coverageInfo=r[0],"2.0"===this.version.slice(0,3)&&(this.coverageInfo.lonLatEnvelope=o.lonLatEnvelope,this.coverageInfo.supportedInterpolations=f(s.supportedInterpolations)),this.datasetName=this.coverageInfo.title;const{rasterInfo:a}=this.coverageInfo;if(this.createRemoteDatasetStorageInfo(a,512,512),this._set("rasterInfo",a),null==a.spatialReference)throw new t("wcsraster-open",`coverage without spatial referernce is not supported: ${this.coverageId}`);const{pixelType:l,bandCount:c}=await this._getPixelTypeAndBandCount(i);a.pixelType=l,1===a.bandCount&&c>1&&(a.bandCount=c),this.updateTileInfo()}async fetchRawTile(e,n,o,r={}){if(this.isBlockOutside(e,n,o))return null;const{nativePixelSize:a,spatialReference:l}=this.rasterInfo,c=2**e,d=a.x*c,p=a.y*c,{blockWidth:h,blockHeight:f}=this.getBlockWidthHeight(e),{origin:u}=this.rasterInfo.storageInfo.tileInfo,m=this.getTileExtent({x:d,y:p},n,o,u,l,[h,f]),v=this.rasterInfo.extent,x=m.xmax>v.xmax||m.ymin<v.ymin;let w=m,y=h,I=f;if(x&&(w=m.clone().intersection(v),i(w)&&(y=Math.floor((w.xmax-w.xmin)/d),I=Math.floor((w.ymax-w.ymin)/p),w.xmax=w.xmin+d*y,w.ymin=w.ymax-p*I)),s(w)||y<=1||I<=1)return null;const b=await this._getCoverage(w,y,I,c,r);if(!b)return null;const C=await this.decodePixelBlock(b,{width:y,height:I,planes:null,pixelType:null});if(C&&(C.width!==y||C.height!==I))throw new t("wcsraster-fetch",`the reponse has unexpected dimension width: ${C.width}, height: {pixelBlock.height}`);return x?g(C,{x:0,y:0},{width:f,height:f}):C}async _getCapabilities(e){const i={service:"WCS",request:"GetCapabilities"};this.version&&(i.version=this.version,i.acceptVersions=this.version);try{const{data:t}=await this.request(this.url,{query:i,responseType:"xml",signal:e});return h(t)}catch(s){if(!o(s))throw new t("wcslayer:open","wcs capabilities is not valid or supported");throw s}}async _describeCoverage(e){const i={service:"WCS",request:"DescribeCoverage",version:this.version},s=this.version.slice(0,3);"1.0"===s?i.coverage=this.coverageId:"1.1"===s?i.identifiers=this.coverageId:"2.0"===s&&(i.coverageId=this.coverageId);try{const{data:t}=await this.request(this.url,{query:i,responseType:"xml",signal:e});return u(t,this.version)}catch(n){if(!o(n))throw new t("wcslayer:open","wcs coverage description is not valid or supported");throw n}}async _getPixelTypeAndBandCount(e){const{pixelSize:s,extent:o,multidimensionalInfo:r}=this.rasterInfo,a=o.center,l=new m({xmin:a.x-s.x,xmax:a.x+s.x,ymin:a.y-s.y,ymax:a.y+s.y,spatialReference:o.spatialReference});let d;if(i(r)){const e=r.variables[0];d=[],e.dimensions.forEach((t=>{d.push(new c({variableName:e.name,dimensionName:t.name,values:t.hasRegularIntervals?t.extent[0]:t.values[0],isSlice:!0}))}))}const p=await this._getCoverage(l,2,2,1,{multidimensionalDefinition:d,signal:n(e)});if(!p)throw new t("wcsraster-open","unable to determine pixel type");const h=await this.decodePixelBlock(p,{width:2,height:2,planes:null,pixelType:null});return{pixelType:h.pixelType,bandCount:h.getPlaneCount()}}async _getCoverage(e,i,s,n,o){const{coverageDescription:r}=this.coverageInfo,{version:a}=r,l="2.0"===r.version?this._getCoverage201Parameters(e,i,s,n,o,r):"1.1"===r.version?this._getCoverage110Parameters(e,i,s,o,r):this._getCoverage100Parameters(e,i,s,o),c="2.0"===r.version?await this.request(this._constructWCS201Url(l),{signal:o.signal,responseType:"array-buffer"}):await this.request(this.url,{query:l,signal:o.signal,responseType:"array-buffer"});if("1.0"===a)return c.data;const d=p(c);if(d.isMultipart&&d.data){const e=d.data.find((e=>{var t;return(null==(t=e.contentType)?void 0:t.toLowerCase().includes("image"))&&null!=e.contentData}));return null==e?void 0:e.contentData}if(this.ioConfig.allowAnyMediaType)return c.data;throw new t("wcsraster:getcoverage","response is not a valid coverage, multipart response is expected")}_getInterpolationIndex(e){return-1===this.coverageInfo.supportedInterpolations.indexOf(e)||"nearest"===e?0:"bilinear"===e?1:"cubic"===e?2:0}_getCoverage100Parameters(e,t,i,s){const n=`${e.xmin},${e.ymin},${e.xmax},${e.ymax}`,o=e.spatialReference.wkid,r=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"GEOTIFF",{bandIds:a,interpolation:l}=s,c=this._getInterpolationIndex(l),d=a?a.map((e=>this.coverageInfo.bandNames[e])):null,p=v[c];return{service:"WCS",request:"GetCoverage",version:this.version,coverage:this.coverageId,format:r,crs:`EPSG:${o}`,bbox:n,width:t,height:i,interpolation:p,band:null==d?void 0:d.join(",")}}_getCoverage110Parameters(e,t,s,n,o){var r;const{multidimensionalDefinition:a,bandIds:c,interpolation:d}=n,p=e.spatialReference.wkid,h=`urn:ogc:def:crs:EPSG::${p}`,f=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"image/tiff",u=this._getInterpolationIndex(d),g=x[u],m=null==d||0===(null==(r=this.coverageInfo.supportedInterpolations)?void 0:r.indexOf(d)),v=o.domain.spatialDomain,w=v.origin.x<=v.envelope.xmin&&v.origin.y<=v.envelope.ymin,y=e.width/t,I=e.height/s*(w?1:-1),b=w?[e.xmin,e.ymin]:[e.xmin,e.ymax],C=v.useEPSGAxis&&l(p),$=C?`${b[1]},${b[0]}`:`${b[0]},${b[1]}`,S=C?`${I},0,0,${y}`:`${y},0,0,${I}`,j=y/2,T=e.xmin+j,_=e.xmax-j,O=Math.abs(I)/2,P=e.ymin+O,L=e.ymax-O,E=C?`${P},${T},${L},${_},${h}`:`${T},${P},${_},${L},${h}`,R=o.range.find((e=>e.axis.some((e=>e.identifier.toLowerCase().indexOf("band")>-1))));let q,k=R&&g&&c?m?`${R.identifier}[${R.axis[0].identifier}[${c.join(",")}]]`:`${R.identifier}:${g}[${R.axis[0].identifier}[${c.join(",")}]]`:null;if(i(a)&&a.length)for(let i=0;i<a.length;i++){let e=a[i].values;const t=a[i].dimensionName.toLowerCase(),s=a[i].variableName.toLowerCase();if(e.length>0)if(e[0]instanceof Array&&(e=e[0]),"stdtime"===t)q=e.map((e=>this._convertToISOTime(e))).join(",");else{const i=o.range.find((e=>e.identifier.toLowerCase()===s));if(i){const s=i.axis.find((e=>e.identifier.toLowerCase()===t));s&&(k=m?i.identifier+"["+s.identifier+"["+e.join(",")+"]]":i.identifier+":"+g+"["+s.identifier+"["+e.join(",")+"]]")}}}return{service:"WCS",request:"GetCoverage",version:this.version,identifier:this.coverageId,format:f,crs:`EPSG:${p}`,boundingbox:E,gridCS:"urn:ogc:def:cs:OGC:0.0:Grid2dSquareCS",gridType:"urn:ogc:def:method:WCS:1.1:2dGridIn2dCrs",gridOrigin:$,gridOffsets:S,gridBaseCRS:h,timeSequence:q,rangeSubset:k}}_convertToISOTime(e,t=!1){return(t?new Date(24*(e-25569)*60*60*1e3):new Date(e)).toISOString()}_getCoverage201Parameters(e,t,s,n,o,r){const{multidimensionalDefinition:a,interpolation:l}=o,c=this._getInterpolationIndex(l);let d=null;const{supportedInterpolations:p}=this.capabilities;if(null!=p&&p.length)switch(c){case 0:d=p.find((e=>e.includes("nearest")));break;case 1:d=p.find((e=>e.includes("linear")));break;case 2:d=p.find((e=>e.includes("cubic")||e.includes("quadratic")))}const h=(this.coverageInfo.supportedFormats||[]).find((e=>e.toLowerCase().indexOf("tiff")>-1))||"image/tiff",{bandNames:f}=this.coverageInfo,{boundedBy:u,domainSet:g,rangeType:m}=r,v=u.isEastFirst?0:1,x=1-v,{axisLabels:w}=u,y=w[v],I=w[x],b=`http://www.opengis.net/def/crs/EPSG/0/${e.spatialReference.wkid}`,C=b,$=[];$.push(`${y}(${e.xmin},${e.xmax})`),$.push(`${I}(${e.ymin},${e.ymax})`);const S=[];if(w.length>2)for(let i=2;i<w.length;i++){const e=g.origin[i];if(w[i].toLowerCase().indexOf("time")>-1){let t=e.toString();u.uomLabels[i].toLowerCase().indexOf("ole")>-1&&(S.push(w[i]),t=this._convertToISOTime(e,!0)),$.push(w[i]+",http://www.opengis.net("+t+")")}else $.push(w[i]+",http://www.opengis.net("+e+")")}let j=null;if(i(a)&&a.length){const e=[];m.forEach((t=>t.forEach((t=>e.push(t.name)))));const t=[];for(let i=0;i<a.length;i++){const s=w.find((e=>e===a[i].dimensionName)),n=e.find((e=>e===a[i].variableName));if(-1===t.indexOf(n)&&t.push(n),s){let e=a[i].values;if(e.length>0){Array.isArray(e[0])&&(e=e[0]);let t="";t=s.toLowerCase().indexOf("time")>-1?e.map((e=>this._convertToISOTime(e))).join(","):e.join(",");const i=$.findIndex((e=>0===e.indexOf(s+",http://www.opengis.net")));-1===i&&$.push(s+",http://www.opengis.net("+t+")"),-1!==i&&-1===$[i].indexOf("("+t+")")&&$.splice(i,1,s+",http://www.opengis.net("+t+")")}}}t.length&&(j=t.join(","))}else if(null!=f&&f.length){j=(o.bandIds?o.bandIds.map((e=>f[e])):f).join(",")}const T=$.join("&subset="),_=!(r.domainSet.axisLabels.join("")===r.boundedBy.axisLabels.join(""))&&!1!==this.ioConfig.allowScaleFactor,O=_?null:`${y}(${t}),${I}(${s})`,P=_?1/n:null;return{service:"WCS",request:"GetCoverage",version:this.version,coverageId:this.coverageId,rangesubset:j,interpolation:d,scaleSize:O,scaleFactor:P,subset:T,format:h,mediaType:this.ioConfig.allowAnyMediaType?null:"multipart/related",outputcrs:b,subsettingcrs:C}}_constructWCS201Url(e){const t={...this.ioConfig.customFetchParameters,...e},i=[];Object.keys(t).forEach((e=>{const s=t[e];null!=s&&("subset"===e?s.split("&subset=").forEach((e=>{e&&i.push(`subset=${encodeURIComponent(e)}`)})):i.push(`${e}=${encodeURIComponent(s)}`))}));return`${encodeURI(this.url)}?${i.join("&")}`}};e([r({type:String,json:{write:!0}})],w.prototype,"datasetFormat",void 0),e([r()],w.prototype,"tileType",void 0),e([r({type:String,json:{write:!0}})],w.prototype,"version",void 0),e([r({type:String,json:{write:!0}})],w.prototype,"coverageId",void 0),w=e([a("geoscene.layers.support.rasterDatasets.ImageServerRaster")],w);const y=w;export{y as default};

package/layers/support/rasterFormats/Lerc.js

@@ -2,4 +2,4 @@
 All material copyright GeoScene, All Rights Reserved, unless otherwise specified.
 See https://js.geoscene.cn/4.23/geoscene/copyright.txt for details.
 */
-import{getAssetUrl as e}from"../../../assets.js";const t=[{pixelType:"S8",size:1,ctor:Int8Array,range:[-128,128]},{pixelType:"U8",size:1,ctor:Uint8Array,range:[0,255]},{pixelType:"S16",size:2,ctor:Int16Array,range:[-32768,32767]},{pixelType:"U16",size:2,ctor:Uint16Array,range:[0,65536]},{pixelType:"S32",size:4,ctor:Int32Array,range:[-2147483648,2147483647]},{pixelType:"U32",size:4,ctor:Uint32Array,range:[0,4294967296]},{pixelType:"F32",size:4,ctor:Float32Array,range:[-34027999387901484e22,34027999387901484e22]},{pixelType:"F64",size:8,ctor:Float64Array,range:[-17976931348623157e292,17976931348623157e292]}];let n=null,r=!1;function a(){return n||(n=import("../../../chunks/lerc-wasm.js").then((e=>e.l)).then((({default:t})=>t({locateFile:t=>e(`esri/layers/support/rasterFormats/${t}`)}))).then((e=>{o(e),r=!0})),n)}function l(){return r}const s={getBlobInfo:null,decode:null};function i(e){return 16+(e>>3<<3)}function o(e){const{_malloc:n,_free:r,_lerc_getBlobInfo:a,_lerc_getDataRanges:l,_lerc_decode_4D:o,asm:u}=e;let c;const f=Object.values(u).find((t=>t&&"buffer"in t&&t.buffer===e.HEAPU8.buffer)),y=e=>{const t=e.map((e=>i(e))),r=t.reduce(((e,t)=>e+t)),a=n(r);c=new Uint8Array(f.buffer);let l=t[0];t[0]=a;for(let n=1;n<t.length;n++){const e=t[n];t[n]=t[n-1]+l,l=e}return t};s.getBlobInfo=e=>{const t=12,n=3,s=new Uint8Array(4*t),i=new Uint8Array(8*n),[o,u,h]=y([e.length,s.length,i.length]);c.set(e,o),c.set(s,u),c.set(i,h);let p=a(o,e.length,u,h,t,n);if(p)throw r(o),`lerc-getBlobInfo: error code is ${p}`;c=new Uint8Array(f.buffer),s.set(c.slice(u,u+s.length)),i.set(c.slice(h,h+i.length));const g=new Uint32Array(s.buffer),d=new Float64Array(i.buffer),[b,m,A,w,U,x,T,C,V,,D]=g,z={version:b,dimCount:A,width:w,height:U,validPixelCount:T,bandCount:x,blobSize:C,maskCount:V,dataType:m,minValue:d[0],maxValue:d[1],maxZerror:d[2],statistics:[],bandCountWithNoData:D};if(D)return z;if(1===A&&1===x)return r(o),z.statistics.push({minValue:d[0],maxValue:d[1]}),z;const I=A*x*8,F=new Uint8Array(I),_=new Uint8Array(I);let k,v=o,B=0,S=!1;if(c.byteLength<o+2*I?(r(o),S=!0,[v,B,k]=y([e.length,I,I]),c.set(e,v)):[B,k]=y([I,I]),c.set(F,B),c.set(_,k),p=l(o,e.length,A,x,B,k),p)throw r(v),S||r(B),`lerc-getDataRanges: error code is ${p}`;c=new Uint8Array(f.buffer),F.set(c.slice(B,B+I)),_.set(c.slice(k,k+I));const O=new Float64Array(F.buffer),$=new Float64Array(_.buffer),j=z.statistics;for(let r=0;r<x;r++)if(A>1){const e=O.slice(r*A,(r+1)*A),t=$.slice(r*A,(r+1)*A),n=Math.min.apply(null,e),a=Math.max.apply(null,t);j.push({minValue:n,maxValue:a,dimStats:{minValues:e,maxValues:t}})}else j.push({minValue:O[r],maxValue:$[r]});return r(v),S||r(B),z},s.decode=(e,n)=>{const{maskCount:a,dimCount:l,bandCount:s,width:i,height:u,dataType:h,bandCountWithNoData:p}=n,g=t[h],d=i*u,b=new Uint8Array(d*s),m=d*l*s*g.size,A=new Uint8Array(m),w=new Uint8Array(s),U=8*s,x=new Uint8Array(U),[T,C,V,D,z]=y([e.length,b.length,A.length,w.length,x.length]);c.set(e,T),c.set(b,C),c.set(A,V),c.set(w,D),c.set(x,z);const I=o(T,e.length,a,C,l,i,u,s,h,V,D,z);if(I)throw r(T),`lerc-decode: error code is ${I}`;c=new Uint8Array(f.buffer),A.set(c.slice(V,V+m)),b.set(c.slice(C,C+d));let F=null;if(p){w.set(c.slice(D,D+s)),x.set(c.slice(z,z+U)),F=[];for(let e=0;e<w.length;e++)F.push(w[e]?x[e]:null)}return r(T),{data:A,maskData:b,noDataValues:F}}}function u(e,t,n,r,a){if(n<2)return e;const l=new r(t*n);if(a)for(let s=0,i=0;s<t;s++)for(let r=0,a=s;r<n;r++,a+=t)l[a]=e[i++];else for(let s=0,i=0;s<t;s++)for(let r=0,a=s;r<n;r++,a+=t)l[i++]=e[a];return l}function c(e,n={}){var r;const a=null!=(r=n.inputOffset)?r:0,l=e instanceof Uint8Array?e.subarray(a):new Uint8Array(e,a),i=s.getBlobInfo(l),{data:o,maskData:c}=s.decode(l,i),{width:f,height:y,bandCount:h,dimCount:p,dataType:g,maskCount:d,statistics:b}=i,m=t[g],A=new m.ctor(o.buffer),w=[],U=[],x=f*y,T=x*p;for(let t=0;t<h;t++){const e=A.subarray(t*T,(t+1)*T);if(n.returnPixelInterleavedDims)w.push(e);else{const t=u(e,x,p,m.ctor,!0);w.push(t)}U.push(c.subarray(t*T,(t+1)*T))}const C=0===d?null:1===d?U[0]:new Uint8Array(x);if(d>1){C.set(U[0]);for(let e=1;e<U.length;e++){const t=U[e];for(let e=0;e<x;e++)C[e]=C[e]&t[e]}}const{noDataValue:V}=n,D=null!=V&&m.range[0]<=V&&m.range[1]>=V;if(d>0&&D){const{noDataValue:e}=n;for(let t=0;t<h;t++){const n=w[t],r=U[t]||C;for(let t=0;t<x;t++)0===r[t]&&(n[t]=e)}}const z=d===h&&h>1?U:null;return{width:f,height:y,bandCount:h,pixelType:n.pixelType&&0===i.version?n.pixelType:m.pixelType,dimCount:p,statistics:b,pixels:w,mask:C,bandMasks:z}}function f(e,t={}){var n,r;const a=e instanceof Uint8Array?e.subarray(null!=(n=t.inputOffset)?n:0):new Uint8Array(e,null!=(r=t.inputOffset)?r:0);return s.getBlobInfo(a)}export{c as decode,f as getBlobInfo,l as isLoaded,a as load};
+import{getAssetUrl as e}from"../../../assets.js";const t=[{pixelType:"S8",size:1,ctor:Int8Array,range:[-128,128]},{pixelType:"U8",size:1,ctor:Uint8Array,range:[0,255]},{pixelType:"S16",size:2,ctor:Int16Array,range:[-32768,32767]},{pixelType:"U16",size:2,ctor:Uint16Array,range:[0,65536]},{pixelType:"S32",size:4,ctor:Int32Array,range:[-2147483648,2147483647]},{pixelType:"U32",size:4,ctor:Uint32Array,range:[0,4294967296]},{pixelType:"F32",size:4,ctor:Float32Array,range:[-34027999387901484e22,34027999387901484e22]},{pixelType:"F64",size:8,ctor:Float64Array,range:[-17976931348623157e292,17976931348623157e292]}];let n=null,r=!1;function a(){return n||(n=import("../../../chunks/lerc-wasm.js").then((e=>e.l)).then((({default:t})=>t({locateFile:t=>e(`geoscene/layers/support/rasterFormats/${t}`)}))).then((e=>{o(e),r=!0})),n)}function l(){return r}const s={getBlobInfo:null,decode:null};function i(e){return 16+(e>>3<<3)}function o(e){const{_malloc:n,_free:r,_lerc_getBlobInfo:a,_lerc_getDataRanges:l,_lerc_decode_4D:o,asm:u}=e;let c;const f=Object.values(u).find((t=>t&&"buffer"in t&&t.buffer===e.HEAPU8.buffer)),y=e=>{const t=e.map((e=>i(e))),r=t.reduce(((e,t)=>e+t)),a=n(r);c=new Uint8Array(f.buffer);let l=t[0];t[0]=a;for(let n=1;n<t.length;n++){const e=t[n];t[n]=t[n-1]+l,l=e}return t};s.getBlobInfo=e=>{const t=12,n=3,s=new Uint8Array(4*t),i=new Uint8Array(8*n),[o,u,h]=y([e.length,s.length,i.length]);c.set(e,o),c.set(s,u),c.set(i,h);let p=a(o,e.length,u,h,t,n);if(p)throw r(o),`lerc-getBlobInfo: error code is ${p}`;c=new Uint8Array(f.buffer),s.set(c.slice(u,u+s.length)),i.set(c.slice(h,h+i.length));const g=new Uint32Array(s.buffer),d=new Float64Array(i.buffer),[b,m,A,w,U,x,T,C,V,,D]=g,z={version:b,dimCount:A,width:w,height:U,validPixelCount:T,bandCount:x,blobSize:C,maskCount:V,dataType:m,minValue:d[0],maxValue:d[1],maxZerror:d[2],statistics:[],bandCountWithNoData:D};if(D)return z;if(1===A&&1===x)return r(o),z.statistics.push({minValue:d[0],maxValue:d[1]}),z;const I=A*x*8,F=new Uint8Array(I),_=new Uint8Array(I);let k,v=o,B=0,S=!1;if(c.byteLength<o+2*I?(r(o),S=!0,[v,B,k]=y([e.length,I,I]),c.set(e,v)):[B,k]=y([I,I]),c.set(F,B),c.set(_,k),p=l(o,e.length,A,x,B,k),p)throw r(v),S||r(B),`lerc-getDataRanges: error code is ${p}`;c=new Uint8Array(f.buffer),F.set(c.slice(B,B+I)),_.set(c.slice(k,k+I));const O=new Float64Array(F.buffer),$=new Float64Array(_.buffer),j=z.statistics;for(let r=0;r<x;r++)if(A>1){const e=O.slice(r*A,(r+1)*A),t=$.slice(r*A,(r+1)*A),n=Math.min.apply(null,e),a=Math.max.apply(null,t);j.push({minValue:n,maxValue:a,dimStats:{minValues:e,maxValues:t}})}else j.push({minValue:O[r],maxValue:$[r]});return r(v),S||r(B),z},s.decode=(e,n)=>{const{maskCount:a,dimCount:l,bandCount:s,width:i,height:u,dataType:h,bandCountWithNoData:p}=n,g=t[h],d=i*u,b=new Uint8Array(d*s),m=d*l*s*g.size,A=new Uint8Array(m),w=new Uint8Array(s),U=8*s,x=new Uint8Array(U),[T,C,V,D,z]=y([e.length,b.length,A.length,w.length,x.length]);c.set(e,T),c.set(b,C),c.set(A,V),c.set(w,D),c.set(x,z);const I=o(T,e.length,a,C,l,i,u,s,h,V,D,z);if(I)throw r(T),`lerc-decode: error code is ${I}`;c=new Uint8Array(f.buffer),A.set(c.slice(V,V+m)),b.set(c.slice(C,C+d));let F=null;if(p){w.set(c.slice(D,D+s)),x.set(c.slice(z,z+U)),F=[];for(let e=0;e<w.length;e++)F.push(w[e]?x[e]:null)}return r(T),{data:A,maskData:b,noDataValues:F}}}function u(e,t,n,r,a){if(n<2)return e;const l=new r(t*n);if(a)for(let s=0,i=0;s<t;s++)for(let r=0,a=s;r<n;r++,a+=t)l[a]=e[i++];else for(let s=0,i=0;s<t;s++)for(let r=0,a=s;r<n;r++,a+=t)l[i++]=e[a];return l}function c(e,n={}){var r;const a=null!=(r=n.inputOffset)?r:0,l=e instanceof Uint8Array?e.subarray(a):new Uint8Array(e,a),i=s.getBlobInfo(l),{data:o,maskData:c}=s.decode(l,i),{width:f,height:y,bandCount:h,dimCount:p,dataType:g,maskCount:d,statistics:b}=i,m=t[g],A=new m.ctor(o.buffer),w=[],U=[],x=f*y,T=x*p;for(let t=0;t<h;t++){const e=A.subarray(t*T,(t+1)*T);if(n.returnPixelInterleavedDims)w.push(e);else{const t=u(e,x,p,m.ctor,!0);w.push(t)}U.push(c.subarray(t*T,(t+1)*T))}const C=0===d?null:1===d?U[0]:new Uint8Array(x);if(d>1){C.set(U[0]);for(let e=1;e<U.length;e++){const t=U[e];for(let e=0;e<x;e++)C[e]=C[e]&t[e]}}const{noDataValue:V}=n,D=null!=V&&m.range[0]<=V&&m.range[1]>=V;if(d>0&&D){const{noDataValue:e}=n;for(let t=0;t<h;t++){const n=w[t],r=U[t]||C;for(let t=0;t<x;t++)0===r[t]&&(n[t]=e)}}const z=d===h&&h>1?U:null;return{width:f,height:y,bandCount:h,pixelType:n.pixelType&&0===i.version?n.pixelType:m.pixelType,dimCount:p,statistics:b,pixels:w,mask:C,bandMasks:z}}function f(e,t={}){var n,r;const a=e instanceof Uint8Array?e.subarray(null!=(n=t.inputOffset)?n:0):new Uint8Array(e,null!=(r=t.inputOffset)?r:0);return s.getBlobInfo(a)}export{c as decode,f as getBlobInfo,l as isLoaded,a as load};