使用遥感图像来判断是否适合赏枫

1
//VERSION=3
2
function setup() {
3
  return {
4
    input: ["B04", "B08", "SCL", "dataMask"],
57 collapsed lines
5
    output: [
6
      { id: "default", bands: 4 },
7
      { id: "index", bands: 1, sampleType: "FLOAT32" },
8
      { id: "eobrowserStats", bands: 2, sampleType: "FLOAT32" },
9
      { id: "dataMask", bands: 1 }
10
    ]
11
  };
12
}
13

14
// 参数：根据需要微调这三个阈值
15
const NDVI_MIN = 0.25;  // 低于此值基本认为非有效植被（透明）
16
const FALL_LO = 0.30;   // 落叶/变色窗口下界（偏黄）
17
const FALL_HI = 0.55;   // 落叶/变色窗口上界（偏红）
18

19
function isVegetation(scl) {
20
  return scl === 4; // SCL=4: Vegetation
21
}
22

23
function lerp(a,b,t){ return a + (b-a)*t; }
24
function clamp01(x){ return Math.max(0, Math.min(1, x)); }
25

26
function evaluatePixel(s) {
27
  const ndvi = (s.B08 - s.B04) / (s.B08 + s.B04);
28
  const veg = isVegetation(s.SCL) && s.dataMask === 1;
29

30
  // 非植被或无效像素：透明
31
  if (!veg || ndvi < NDVI_MIN) {
32
    return {
33
      default: [0,0,0,0],
34
      index: [NaN],
35
      eobrowserStats: [0,0],
36
      dataMask: [0]
37
    };
38
  }
39

40
  let rgba;
41

42
  if (ndvi >= FALL_LO && ndvi <= FALL_HI) {
43
    // 在“变色窗口”内：黄(#FFC800) → 红(#FF0000) 渐变
44
    const t = clamp01((ndvi - FALL_LO) / (FALL_HI - FALL_LO)); // 0..1
45
    const r = lerp(1.00, 1.00, t);
46
    const g = lerp(0.78, 0.00, t);
47
    const b = lerp(0.00, 0.00, t);
48
    rgba = [r, g, b, 1];   // 高亮标红
49
  } else if (ndvi > FALL_HI) {
50
    // 常绿/高NDVI：降饱和的深绿，避免喧宾夺主
51
    const g = clamp01((ndvi - FALL_HI) / (0.85 - FALL_HI)); // 0..1
52
    rgba = [0.05, 0.25 + 0.5*g, 0.08, 0.9];
53
  } else {
54
    // 介于 NDVI_MIN 和 FALL_LO：稀疏植被，弱黄棕
55
    rgba = [0.75, 0.65, 0.35, 0.9];
56
  }
57

58
  return {
59
    default: rgba,
60
    index: [ndvi],
61
    eobrowserStats: [ndvi, 0],
62
    dataMask: [1]
63
  };
64
}

1
//VERSION=3
2
function setup() {
50 collapsed lines
3
 return {
4
   input: ["B04","B08","SCL","dataMask"],
5
   output: [{ id:"default", bands:4 }],
6
   mosaicking: "ORBIT"  // 允许多时相
7
 };
8
}
9

10
// ---- 参数可调 ----
11
const NDVI_MIN = 0.25;         // 植被下限（低于则透明）
12
const DROP_LO  = 0.05;         // 认为“开始变色”的NDVI下降阈值
13
const DROP_HI  = 0.25;         // 强变色上限（越接近越红）
14
const CLOUDS   = new Set([7,8,9,10]); // SCL 云类（排除）
15

16
function ndvi(s){ return (s.B08 - s.B04) / (s.B08 + s.B04); }
17
function clamp01(x){ return Math.max(0, Math.min(1, x)); }
18

19
function goodVeg(s){
20
 return s.dataMask === 1 && s.SCL === 4; // 仅植被
21
}
22
function notCloud(s){ return !CLOUDS.has(s.SCL); }
23

24
function evaluatePixel(samples, scenes) {
25
 if (samples.length < 2) return [0,0,0,0];
26

27
 // 约定：samples[0] 较新，samples[samples.length-1] 较早
28
 const sNew = samples[0];
29
 const sOld = samples[samples.length-1];
30

31
 // 基本质量控制：两期都得是植被且非云
32
 if (!(goodVeg(sNew) && goodVeg(sOld) && notCloud(sNew) && notCloud(sOld))) {
33
   return [0,0,0,0];
34
 }
35

36
 const nNew = ndvi(sNew);
37
 const nOld = ndvi(sOld);
38
 if (nNew < NDVI_MIN && nOld < NDVI_MIN) return [0,0,0,0];
39

40
 const drop = nOld - nNew;  // NDVI 下降为正值 = 疑似变色/落叶
41
 if (drop >= DROP_LO) {
42
   // 黄(#FFC800) -> 红(#FF0000) 渐变
43
   const t = clamp01((drop - DROP_LO) / (DROP_HI - DROP_LO));
44
   const r = 1.0;
45
   const g = 0.78*(1.0 - t);  // 0.78 -> 0
46
   const b = 0.0;
47
   return [r, g, b, 1.0];
48
 } else {
49
   // 变化不大：用低饱和绿显示“仍然茂盛的林地”
50
   // 这样红色区域会非常醒目
51
   const g = clamp01((nNew - NDVI_MIN) / (0.85 - NDVI_MIN));
52
   return [0.05, 0.25 + 0.5*g, 0.08, 0.6];
53
 }
54
}