АСТАНА — ПЕРЕНОС МЕТОДИКИ MILAN JANOSOV

Источники
Полный текст статьи прочитан через доступ пользователя в Chrome:
https://milanjanosov.substack.com/p/mapping-summer-heat-from-space
Интерфейс и режимы сверены с открытой картой:
https://bejewelled-dango-fb95d2.netlify.app/hungary_lst_map.html

Что повторено из опубликованной статьи
• Сезон каждого года: 15 июня — 15 августа включительно.
• Граница из предоставленного GeoJSON; прямоугольное окно вокруг неё
  расширено на 0,08° во все стороны. Растр не обрезается по границе города.
• Landsat Collection 2 Level-2, Tier 1, через Microsoft Planetary Computer.
• Четыре канала: готовая температура поверхности, QA_PIXEL, QA_RADSAT, ST_QA.
• LST = DN × 0,00341802 + 149 − 273,15; нулевые DN исключены.
• Исключены биты QA_PIXEL 0–5, все ненулевые QA_RADSAT и ST_QA > 4 K.
• Маски применены на исходной сетке до билинейного приведения температуры.
• Единая сетка 30 м. Для Астаны используется UTM 42N / EPSG:32642 вместо
  венгерской EPSG:23700, как предписано в статье для другого города.
• Медиана по отдельным сценам. Сцены одной даты не объединяются отдельно:
  это соответствует циклу и временному стеку в опубликованном примере.
• Пиксель может опираться на одно пригодное наблюдение; это не прежний
  вариант с требованием минимум трёх дат. Годовой кадр создаётся при
  наличии минимум двух сцен с пригодными данными в расчётном окне.
• Остаточные пробелы заполняются средним известных соседей в окне 3 × 3;
  не более 60 проходов. Используется scipy.ndimage.uniform_filter, как
  в статье. Маска интерполированных пикселей сохранена отдельно.
• Для здания берётся температура пикселя под representative_point;
  для дороги — под серединой линии. Дороги рисуются буфером 6 м.
• Классы дорог из примера: primary, secondary, tertiary, residential.

Выбор сцен: статья и полная карта
В учебном примере автор сокращает набор до 10 наименее облачных сцен
и прямо указывает, что полная обработка использует все сцены.
Интерактивная карта Астаны использует все подходящие сцены.
Дополнительно для 2026 года сохранён вариант *_2026_tutorial.tif,
рассчитанный по 10 наименее облачным сценам — для воспроизведения
именно короткого учебного примера. Наборы нельзя смешивать.

Исторические режимы открытой карты
Шкала времени — 1984–2026; годы с недостаточным архивом не дорисовываются.
Recent — попиксельная медиана годовых карт 2020–2026.
Baseline — попиксельная медиана доступных годовых карт 1990–1999.
Δ — Recent минус Baseline. Временной интерполяции между годами нет.
Окраска изменена 04.10.2026 по запросу пользователя.
Для годовых карт и Recent: цвет = LST минус среднее арифметическое
всех конечных значений растра внутри предоставленной границы Астаны.
Все пиксели имеют одинаковую площадь в UTM; учитываются также вода,
незастроенные территории и заполненные пропуски. Это не среднее
по зданиям, не медиана и не средняя температура воздуха.
Для Δ цвет по-прежнему означает Recent минус Baseline, без вычитания
городского среднего изменения. Эти два смысла явно подписаны в легенде.
Палитра заменена 04.10.2026 на девять ступеней с пользовательского
скриншота «SUMMER HEAT (LST)». Цвета считаны после преобразования
встроенного ICC-профиля изображения в sRGB:
#313696 #4575b4 #75acd0 #aad9e9 #ffffbf #fee090 #fdae61 #f46d43 #a20025.
Границы ступеней отклонения: −18, −14, −10, −6, −2, +2, +6, +10, +14, +18 °C.
Центральная светло-жёлтая ступень — ±2 °C относительно среднего.
Крайние ступени включают значения вне ±18; знаки ≤ / ≥ в легенде
обозначают насыщение окраски. Фактические диапазоны показаны отдельно.
Числа 33–48 °C с образца не перенесены: пользователь запросил цвета,
а ранее выбранный смысл окраски — отклонение от среднего — сохранён.
Непрозрачность стартует с 85%. Граница белая, толщина 1 px.

Проверка диапазона Δ
Внутри Астаны: минимум −11,5006, максимум +22,4008 °C;
среднее +4,6983; 2-й и 98-й процентили −2,4367 и +12,1667 °C.
9,5703% пикселей выходили за старые пределы окраски ±9 °C.
Во всём окне карты: −14,0306…+22,4008 °C.
В точках выборки зданий внутри города: −4,8195…+21,5921 °C;
для дорог: −3,7245…+20,5765 °C. Интерфейс показывает диапазон
активного слоя; опорное среднее всегда вычислено по растру города.
Проверено попиксельное равенство delta = recent − baseline.
Контроль без пространственного заполнения пропусков даёт диапазон
−17,9625…+22,2776 °C (иной набор наблюдений и покрытия).
Следовательно, ±9 не описывали фактические экстремумы исходного расчёта.
Экстремумы не следует интерпретировать как достоверный климатический тренд:
число наблюдений, даты, покрытие и землепользование различаются.
Медианное число лет с непосредственно измеренной температурой в пикселе:
4 для периода 1990–1999, 7 для 2020–2026; местами их нет и данные заполнены.

Неопределённые в статье детали полной интерактивной версии
Статья не содержит полный генератор 1984–2026 и не определяет загородную
референтную область для показанного на сайте числа SUHI. SUHI
теперь рассчитан отдельно для Астаны по явному определению ниже;
это не восстановленная неизвестная референтная область автора.
Исторические годы здесь рассчитаны повторением опубликованного
годового алгоритма, а многолетние периоды — по определениям интерфейса.
Статья показывает здания и дороги в демонстрационном круге 2,5 км;
интерактивный сайт показывает их во всём окне. Здесь для режима карты
использовано полное окно, а не ограниченный учебный круг.
Сетка фиксирована для всех лет, чтобы обеспечить попиксельное сравнение.
Положение её начала может отличаться от автоматического начала сетки
первой сцены в учебном notebook; шаг и способ пересчёта сохранены.

Оформление и подложка
Сохранены расположение панелей, переключатели слоёв и периодов,
ползунки, клавиши стрелок и белый контур границы. Подписи переведены
на русский. Палитра заменена по запросу пользователя; ориентир:
https://www.ventusky.com/ru
Здания и дороги рисуются из векторных контуров на Canvas с учётом
плотности экрана, без растягивания уменьшенных PNG. Геометрия упрощена
с допуском 0,3 м в Web Mercator. Цвет берётся из исходного 30-м растра.
Пиксельный слой экспортируется с длинной стороной 2481 пиксель
(раньше 800), ближайшим соседом, без визуального сглаживания при зуме.
Повышение чёткости изображения не повышает детализацию теплового сенсора.
Подложка: Esri World Dark Gray Base. Подписи: отдельный встроенный слой OpenStreetMap / Geofabrik 2026
из gis_osm_places_free, gis_osm_places_a_free и gis_osm_roads_free.
Рисуются поверх тепловой карты, выбор названий зависит от масштаба;
пересечения подписей подавляются. Проверен публичный сервис подложки:
https://services.arcgisonline.com/arcgis/rest/services/Canvas/World_Dark_Gray_Base/MapServer
© Esri, HERE, Garmin, OpenStreetMap contributors, GIS User Community.
Подложка загружается через интернет; автономный HTML содержит все
тепловые данные, но картографические плитки не встроены. Это историческая
растровая подложка Esri; стандартный слой её подписей заменён
локальными названиями OSM из-за обнаруженного устаревшего «Нур-Султан».

SUHI: отдельное определение для Астаны
SUHI = медиана LST внутри административной границы города минус
медиана LST открытых пригородных территорий. Городская медиана использует
всю предоставленную границу, включая незастроенные территории и воду.
Пригородный фон: кольцо от 1 до 5 км вне границы, только классы OSM
farmland, grass, meadow, scrub. Исключены вода с буфером 90 м,
здания с буфером 100 м и residential / industrial / commercial / retail /
quarry / landfill / military с буфером 100 м. Все расстояния в EPSG:32642.
Точная область сохранена в suhi_reference.geojson. Маска одинакова для
всех лет и основана на OSM 2026, а не исторической реконструкции.
Используются те же температурные композиты после заполнения пропусков,
что и на карте. При покрытии менее 50% референтной маски SUHI не показан.
Положительное число: административная территория теплее фона;
отрицательное: холоднее. Знак не принуждается к положительному.
В режиме Recent медианы рассчитываются по многолетнему композиту.
В режиме Δ SUHI = SUHI Recent − SUHI Baseline; это разность двух индексов.
Scenes: число сцен с пригодными данными в расчётном окне для годового
кадра; для многолетнего режима — сумма по использованным годам.
Подход сравнения с открытыми пригородными территориями описан, например,
в отчёте World Bank, но параметры маски выше выбраны для этой карты:
https://documents1.worldbank.org/curated/en/099400002252239104/pdf/P17117609533dd04a0aa2503707c6e7527d.pdf

Геометрия зданий и дорог
OpenStreetMap / Geofabrik, локальный файл Kazakhstan от 16.09.2026,
пространственная выборка по тому же окну. Геометрия одинакова для всех
лет: старые кадры не являются реконструкцией существовавшей тогда застройки.
Тепловая детализация Landsat — примерно 100 м (для ранних сенсоров иная);
30-метровая сетка и контуры отдельных зданий не дают температуру крыши.
Температура относится к поверхности, а не к воздуху и не к суточному максимуму.

Файлы
index.html — рабочая версия с внешними локальными кадрами.
astana_lst_map.html — автономный HTML, когда выполнена окончательная сборка.
rasters/lst_YEAR.tif — годовая температура после заполнения пропусков.
rasters/raw_YEAR.tif — медиана только измеренных значений.
rasters/count_YEAR.tif — число пригодных сцен в пикселе.
rasters/filled_mask_YEAR.tif — 1: пространственно интерполированный пиксель.
rasters/YEAR.json — список сцен, параметры и доля интерполяции.
rasters/recent.tif, baseline.tif, delta.tif — многолетние расчёты.
buildings.gpkg, roads.gpkg — использованные геометрии.
catalog/ — метаданные источников; cache/ — обработанные сцены для повторного запуска.

Запуск
Python с зависимостями requirements.txt.
python pipeline.py --catalog-only --years all
python pipeline.py --years all
python pipeline.py --tutorial --years 2026
python render.py
python package_map.py

Источники данных и лицензии
USGS Landsat 8–9 C2 L2: https://doi.org/10.5066/P9OGBGM6
USGS Landsat 7 C2 L2: https://doi.org/10.5066/P9C7I13B
USGS Landsat 4–5 C2 L2: https://doi.org/10.5066/P9IAXOVV
Граница, здания, дороги: © OpenStreetMap contributors, ODbL.
Данные предоставлены USGS, доступ через Microsoft Planetary Computer.
Leaflet 1.9.4: Vladimir Agafonkin / CloudMade, BSD-2-Clause.
Методика и образец интерфейса: Milan Janosov, Mapping Summer Heat from Space.
HTML интерфейса Астаны реализован отдельно; венгерские температурные
изображения и персональные данные браузера в результат не включены.

Итог выполнения
Проверены все 43 года шкалы. Готовы 38 кадров с данными внутри Астаны.
Без кадра: 1984, 1988, 1990, 1995, 1998 (нет достаточных пригодных данных).
2026: 35 отобранных сцен, 34 с пригодными пикселями в окне.
Внутри границы города: 99,0391% пикселей имеют измеренные значения;
0,9609% заполнены методом автора; после заполнения пропусков не осталось.
Медиана поверхности по городу: 41,1134 °C.
Геометрии: 76 187 зданий и 8 048 дорожных объектов.
Архив ZIP содержит автономную карту, код, метаданные и основные GeoTIFF;
полный комплект годовых GeoTIFF остаётся в соседней папке rasters/.

Проверка SUHI и новой визуализации (04.10.2026)
Референтная область: 95,4387 км² на 30-метровой сетке.
Для 2026 года все её ячейки имеют непосредственно измеренную
температуру в годовом композите, без пространственного заполнения.
City median 41,1134 °C; rural median 45,1420 °C; SUHI −4,0286 °C.
Проверены 38 годовых значений, Recent и разность SUHI двух периодов.
Проверено соответствие всех 40 цветовых кадров и выборок зданий/дорог
исходным растрам. Изображение подложки и подписи проверены в браузере.


АРХИТЕКТУРА РАССКАЗА И ОТБОР ТЕРРИТОРИЙ (07.10.2026)
История использует ANALYSIS: медианы измеренных пикселей, без заполнения
пропусков. Критерий показа годовой медианы: покрытие города >= 50%.
Предварительный отбор: rule333_score 0–1 и измеренная температура ячейки
30 м под representative_point здания выше медианы города 2020–2026.
Порог 38.85 °C — правило проекта, не норматив опасной температуры.
Знаменатели районных долей показаны в таблице; здания имеют равный вес.
14 зданий вне районных полигонов исключены из районной таблицы.
Баллы KBH требуют проверки даты и методики. Код публичного сайта объясняет
их как число выполненных критериев, но архив не содержит компонентов.
Дополнительный green-330300-astana.geojson содержит 2115 объектов, все
mock: true. Он исключён из аналитики. Проверка записана в source_audit.json.
Руководство: https://trekontoret.no/wp-content/uploads/2025/02/The-330300-Handbook-1.pdf
Численная прогнозная модель и величина эффекта озеленения не рассчитаны.
В разделе будущего описаны необходимые сценарии и проверка результата.