> 影像相對配準方法及殘餘誤差見 [配準方法與檢查](COREGISTRATION_RESULTS.md)。

# 台北2024／2025／2026衛星樹覆分類變化（探索版）

更新：2026-09-24。衛星分析呈現同季樹冠面積估計；政府登錄資料另列背景參考。

本版回答同季衛星光譜分類呈現的增減；未經2024–2026高解析參考獨立驗證，數字是10m樹覆類別像元面積，**不是實際樹冠投影面積、砍樹面積或株數**。

## 主結果：市街內部＋交界帶

三年都使用7/1–9/23；下表固定在三年有效日期>=3、分類分數>=0.7，且候選翻轉通過多日期一致性檢查的同一共同可判範圍C。缺測或低信心不判為非樹。

| 期間 | C面積 ha | C占U＋I | 類別增加 ha | 類別減少 ha | 類別淨變化 ha | >=0.09ha成片淨變化 ha |
|---|---:|---:|---:|---:|---:|---:|
| 2024→2025 | 10,569.04 | 73.7% | 1.42 | 10.60 | -9.18 | -3.76 |
| 2025→2026 | 10,569.04 | 73.7% | 0.95 | 16.16 | -15.21 | -7.69 |
| 2024→2026 | 10,569.04 | 73.7% | 1.48 | 25.87 | -24.39 | -11.55 |

成片面積只加總符合斑塊門檻的候選；不必等於全部類別面積差。主圖預設成片候選，可切換全部像元候選；光譜支持不再作讀者頁驗證選項。

## 結論穩定性

各區以medoid／逐波段中位數、有效日期2／3／5天、分數0.6／0.7／0.8做18種設定；另以相同的嚴格共同支持範圍比較兩種合成法。不同設定的可判範圍會改變；範圍不是95%信賴區間，也不涵蓋共同模型的系統性誤差。

- U＋I 2024→2025：18種設定淨值範圍 -10.32 至 -5.50 ha；各設定方向一致，但仍是未驗證分類結果。
- U＋I 2025→2026：18種設定淨值範圍 -16.96 至 -9.80 ha；各設定方向一致，但仍是未驗證分類結果。

## 資料與處理

1. 原始影像：Copernicus Sentinel-2 L2A，經Microsoft Planetary Computer公開STAC及原數值子集取得。2024／2025／2026主窗口候選16／22／22個日期。沒有以全景雲量先剔除；台北SCL合格面積低於0.05%（約13.43ha）的日期略過波段下載並登錄，並非保證該景完全沒有無雲像元。來源ID、平台、相對軌道、處理版本及遮罩覆蓋均留存。
2. QA：SCL4／5／6，雲影及雲類別3／8／9／10擴張20m（10m格網四鄰接膨脹2次）；所有必要波段有效、反射率範圍檢查、極暗像元列未知。原DN對PB>=04套(DN−1000)/10000一次。本版未加入Cloud Score+，SCL可能殘留薄雲與建物陰影。
3. B2/B3/B4/B8為10m，B5/B8A/B11/B12原生20m後雙線性重投影至固定EPSG:3826。影像已完成逐景相對配準，並以保留區塊檢查年度殘餘位移，詳見配準方法文件；這不是獨立高解析樹冠精度驗證。
4. 每年選距逐波段中位數最近的完整8波段實際觀測作medoid；另存中位數版本及NDVI/NDRE/NDMI逐日期中位數、有效日期數與觀測日期分布。所選時間代表季節狀態，不是9/23當日。
5. 免登入分類器：2021 WorldCover弱標籤，搭配同年既有6–9月S2合成；輸入8波段及6光譜指數，固定RF模型套三年。只用光譜作預測，沒有用WRI舊樹冠、2026名冊或分區遮罩強制限制樹類。2021訓練合成與本版端點合成不完全相同，是可能的移轉誤差來源。
6. WorldCover樹覆類別定義含樹覆>=10%，本地光譜轉移模型沒有量測樹高，也沒有證明50%冠覆。模型分數不是冠覆比例。2km空間留出衡量與WorldCover參照的一致性，不能當當代獨立精度或面積誤差。
7. 候選翻轉回看每個原始日期：每年須至少2個日期、且占有效日期至少60%，以分數>=0.7支持該年樹／非樹狀態；不符合者列低信心。C為三年都可觀測、分類分數合格且通過此時間一致性檢查的固定範圍。此檢查沿用同一模型，不是獨立驗證。每格0.01ha；增加與減少分開，再求淨值。所有區的未知觀測、低信心面積及C覆蓋率列在CSV。原始格點增減恆等於兩端樹覆類別面積差。
8. 9像元（0.09ha）為成片顯示門檻，不是已驗證的最小可偵測事件；1／4／9／25像元另報。同源光譜支持不作獨立驗證選項。
9. 本版是三年同季端點篩查；沒有聲稱完成2019起CCDC/BFAST、事件日期歸因、獨立高解析驗證或誤差校正面積。颱風、修剪、含水量、物候、混合像元仍可能造成分類翻轉。

## 固定範圍及顯示

U市街內部8,742.62ha、I交界帶5,600.74ha、F連續林地背景11,351.66ha、O其他1,160.18ha；行政區與分區交叉。市街範圍沿用2021 WorldCover建成地斑塊>=1ha及100m鄰近範圍、填補封閉孔洞<=50ha、排除水面。I以連續樹覆>=10ha且其中>=5ha在市街範圍外者為種子，包含其市街重疊區及市街內100m鄰近帶。F為這些林地在市街範圍外的部分，不稱已證實天然林；分區不是產權或管理權。變化圖預設顯示U＋I候選，森林以淡色分區背景共同呈現；F變化可由分析範圍另選。朱紅表示候選減少、藍染表示候選增加；未知斜線可開啟，隱藏斜線不代表無變化。三年都市主題圖以若竹、山吹、淡苔分別表示都市、交界帶及森林背景的樹覆類別，另保留原始RGB。政府半年報僅在補充連結。衛星變化不能歸因於市府。

## 同源影像目視檢查

已檢視U＋I的前三大增加及前三大減少斑塊。較大減少候選可見由綠色轉裸露／鋪面的光譜差異；前三個增加候選位於同一片開放地，可能混入草地或農地返綠，不能當作已證實新增樹冠。此檢查僅用同一批10m影像，不是獨立驗證，也不因此改動既定門檻。

## 可重現檔案

- `district_zone_change.csv`：12區×分區與全市，兩段及累積分類增減、共同面積、日期與未知。
- `yearly_tree_class_area.csv`：同一C內三年樹覆類別面積。
- `classification_sensitivity.csv`、`patch_sensitivity.csv`、`composite_comparison_same_support.csv`：敏感度及同範圍合成比較。
- `s2_YYYY_jul01_sep23_medoid.tif`：31層合成、指數、觀測與日期、中位數光譜；`tree_proxy_YYYY.tif`：類別及模型分數；`change_YYYY_YYYY_10m.tif`：五類候選圖。
- `temporal_votes_YYYY.tif`、`temporal_check_metadata.json`：候選像元逐日有效觀測、樹／非樹高分票數及時間一致性規則。
- `proxy_model_metadata.json`：參照標籤、模型設定、訓練／測試樣本與限制。
- 原始STAC與DN在 `data/interim/change_2024_2026/satellite_v1/`，未覆寫既有timeseries。

來源：[Copernicus產品說明](https://sentiwiki.copernicus.eu/web/s2-products)、[SCL處理說明](https://sentiwiki.copernicus.eu/web/s2-processing)、[Microsoft Planetary Computer Sentinel-2](https://planetarycomputer.microsoft.com/dataset/sentinel-2-l2a)、[ESA WorldCover 2021使用手冊](https://esa-worldcover.s3.eu-central-1.amazonaws.com/v200/2021/docs/WorldCover_PUM_V2.0.pdf)。Contains modified Copernicus Sentinel data 2021, 2024–2026; ESA WorldCover 2021 v200 / CC BY 4.0。


發布說明：本網站提供公開表格與展示圖；文中提到的原始波段、GeoTIFF、模型和研究快取不包含在網站中。