Ищем на Луне черные дырки

Слушай, давай снимем эту сухую «техно-прозу» и переложим её на нормальный рабочий язык — так, чтобы ты сразу видел, где тут физика, где подвох, и что с этим делать в нашем «институте кратеров».

---

## Что это вообще за файл

Это глобальная карта высот Луны (DEM) с разрешением **118 метров на пиксель**, сделанная по данным лазерного альтиметра LOLA с зонда LRO. По сути — это «цифровой слепок» всей Луны: для каждой точки знаешь, на сколько метров она выше или ниже условного шара радиусом 1737,4 км.

## Главное, что важно именно тебе (и по «Правилу три Пащенко»)

1. **Точность по высоте — примерно 1 метр.** А сам лазерный замер — ещё круче: вертикальное разрешение около 10 сантиметров. То есть перепады на дне кратера и тонкие «одеяла» выброса (ejecta blanket) ты реально увидишь. Это как твой лидар: не просто «есть яма», а «ямка с нюансом».
2. **Горизонтальная точность — лучше 20 метров.** Этого хватает, чтобы уверенно различать зоны выброса и не путать их с фоном.
3. **Есть пробелы между трассами.** Вдоль экватора между лазерными замерами бывают щели по 1–2 км, местами до 4 км. В этих дырах высоту дорисовали интерполяцией — то есть «по среднему». Для грубой оценки рельефа сойдёт, но если ты ищешь тонкую аномалию (тот самый маркер экзотического удара), эти места лучше не брать как главные доказательства.
4. **Формат и масштаб — вот где ловушка.** Файл хранит высоты как 16;битные целые числа, и там стоит множитель: **значение в файле надо делить на 2**, чтобы получить метры. То есть:
   $$
   \text{Высота (м)} = \text{пиксель} \times 0{,}5
   $$
   В QGIS это часто делают автоматом, а в ArcGIS — нет. Если не учесть, получится, что кратер глубиной 2 км, а на деле — 1 км. Это как перепутать масштаб на чертеже: вся физика сразу ломается.
5. **Размер — 8 ГБ.** Это не «картинка для соцсетей», это серьёзный массив. Обрабатывать его лучше не в «простых» просмотрщиках, а в ГИС (QGIS, ArcGIS) или через GDAL.

---

## Как быстро привести к виду «можно работать»

Если ты хочешь сразу получить нормальные метры, без ручных пересчётов, самый надёжный способ — через GDAL:

```bash
gdal_translate -ot Float32 -unscale Lunar_LRO_LOLA_Global_LDEM_118m_Mar2014.tif LOLA_LDEM_118m_Mar2014_unscale.tif
```

После этого у тебя будет GeoTIFF, где каждый пиксель — это реальная высота в метрах. Дальше можно резать на кусочки (clip) под конкретные кратеры, строить профили вдоль радиуса выброса, считать шероховатость и т.д.

---

## Где брать и как искать

- **USGS Astrogeology Science Center** — это главный «склад» таких карт.
- **PDS (NASA Planetary Data System)** — первоисточник сырых и производных данных.
- Для быстрого подбора и нарезки есть **Map a Planet** — там можно выбрать область, формат и скачать уже обрезанное.

---

## Как это использовать по твоей методике

Помнишь про «на пальцах» и про аэрозольную оптику? Вот как это ложится на реальные данные:

1. **Берёшь кратер.** Вырезаешь кусок DEM вокруг него (например, радиус 20–30 км).
2. **Строишь профили** от центра к краю по нескольким направлениям. Смотришь, как плавно или резко спадает толщина выброса. Для обычного удара характерен один профиль, для экзотического (PBH;сценарий) — другой: более крутой, с «провалом» или, наоборот, с узкой аномальной зоной.
3. **Считаешь шероховатость** на разных масштабах. Это как дисперсность аэрозоля: если поверхность «мелкая», почти гладкая — значит, выброс мелкодисперсный; если «рваная» — много крупных обломков.
4. **Накладываешь на снимки LROC** (они есть в тех же каталогах). Видишь: где на карте «гладко», там на фото может быть не россыпь валунов, а тонкая пыль. Это и есть маркер.

---

## Практический совет

Не пытайся сразу «съесть весь файл». Начни с 3–4 кратеров:
- один «эталонный» — классический, хорошо изученный;
- два «подозрительных» — с необычной формой или странным распределением выброса;
- один заведомо «шумный» (с интерполяцией в центре) — чтобы понять, где данные врут.

Так ты быстро поймёшь, где физика, а где артефакты обработки. И это как раз то самое «получить ответ на пальцах», прежде чем гонять тяжёлые модели.

---

Если скажешь, какой кратер хочешь взять первым (или хотя бы район: Море Дождей, Южный полюс и т.п.), я подскажу, как именно нарезать данные и какие профили строить, чтобы сразу увидеть «не тот» выброс. Или можем вообще пока не трогать файлы, а просто прикинуть, как из такого DEM вытащить тот самый аэрозольный параметр, который ты ищешь. Что ближе? ;
USGS: Science for a Changing WorldNASAPDS Cartography and Imaging Science Node
Astropedia AnnexAstropediaLunar and Planetary Cartographic CatalogReturn to search

Sample (jpg) 1024px wide
Download (8 GB)
open
Moon LRO LOLA DEM 118m
Primary AuthorsLOLA Science TeamOriginatorsGoddard Space Flight CenterPublisherUSGS Astrogeology Science CenterPublication Date2014-03-11Abstract
   
    Product Information:

This digital elevation model (DEM) is based on data from the Lunar Orbiter Laser Altimeter (LOLA; Smith et al., 2010), an instrument on the National Aeronautics and Space Agency (NASA) Lunar Reconnaissance Orbiter (LRO) spacecraft (Tooley et al., 2010). The created DEM represents more than 6.5 billion measurements gathered between July 2009 and July 2013, adjusted for consistency in the coordinate system described below, and then converted to lunar radii (Mazarico et al., 2012).

Elevations were computed by subtracting the lunar reference radius of 1737.4 km from the surface radius measurements (LRO Project and LGCWG, 2008; Archinal et al., 2011). Thus elevation values are the distance above or below the reference sphere. The average accuracy of each point after crossover correction is better than 20 meters in horizontal position and ~1 meter in radius (Mazarico et al., 2012). The measurements were converted into a DEM (Neumann et al., 2011) using Generic Mapping Tools software (Wessel & Smith, 2008), with a resolution of 256 pixels per degree. In projection, the pixels are 118 meters in size at the equator. Gaps between tracks of 1-2 km are common, and some gaps of up to 4 km occur near the equator. DEM points located in these gaps in LOLA data were filled by interpolation Smith et al., 2010).

Mission and Instrument Information:

The U.S. National Aeronautics and Space Administration (NASA) launched the Lunar Reconnaissance Orbiter (LRO) spacecraft to the Moon in June 2009 (Tooley et al., 2010) carrying a variety of instruments that continue to return high-resolution images of the lunar surface from its eccentric polar mapping orbit (Petro et al., 2019).

The LOLA has collected over more than 6.5 billion measurements of global surface height with a vertical precision of ~10 cm and an accuracy of ~1m (Mazarico et al., 2013). With such highly accurate global coverage, the resulting topographic map has become the reference geodetic framework for the lunar community and has led to the highest resolution and most accurate polar DEMs to date (Barker et al., 2016).

Note this file is stored using 16bit signed integer with a defined "scaler=0.5" and "offset=0" in the GeoTIFF label. This means all pixel values should be divided by 2. So:
Lunar elevation (m) = pixel_DN * 0.5 + 0.0

Some applications like QGIS will automatically scale the data for the user. Other applications, like ArcGIS Pro will not. You can apply the scale in the raster calculator or using GDAL tools. Using GDAL, you could run:
> gdal_translate -ot Float32 -unscale Lunar_LRO_LOLA_Global_LDEM_118m_Mar2014.tif LOLA_LDEM_118m_Mar2014_unscale.tif

References:

Archinal, B. A., A'Hearn, M. F., Bowell, E., Conrad, A., Consolmagno, G. J., Courtin, R., Fukushima, T., et al. (2011). Report of the IAU Working Group on cartographic coordinates and rotational elements: 2009. Celestial Mechanics and Dynamical Astronomy, 109 (2), 101-135. https://doi.org/10.1007/s10569-010-9320-4

Barker, M. K., Mazarico, E., Neumann, G. A., Zuber, M. T., Haruyama, J., & Smith, D. E. (2016). A new lunar digital elevation model from the Lunar Orbiter Laser Altimeter and SELENE Terrain Camera. Icarus, 273, 346-355. https://doi.org/10.1016/j.icarus.2015.07.039

Davies, M. E., & Colvin, T. R. (2000). Lunar coordinates in the regions of the Apollo landers. Journal of Geophysical Research, 105(E8), 20277-20280. https://doi.org/10.1029/1999JE001165

Folkner, W. M., Williams, J. G., & Boggs, D. H. (2008). The Planetary and Lunar Ephemeris DE 421. JPL Memorandum IOM 343R-08-003, 31 March.
Folkner, W. M., Williams, J. G., & Boggs, D. H. (2009). The Planetary and Lunar Ephemeris DE 421. IPN Progress Report 42-178, August 15. http://ipnpr.jpl.nasa.gov/progress_report/42-178/178C.pdf

Greeley, R., & Batson, R. M. (1990). Planetary mapping. (ISBN 0-521-30774-0, pp. 274-275). New York, NY: Cambridge University Press.

LRO Project and Lunar Geodesy and Cartography Working Group (2008). A standardized lunar coordinate system for the Lunar Reconnaissance Orbiter and lunar datasets: LRO Project and LGCWG White Paper, version 5 of October 1.
Mazarico, E., Rowlands, D. D., Neumann, G. A., Smith, D. E., Torrence, M. H., Lemoine, F. G., & Zuber, M. T. (2012). Orbit determination of the Lunar Reconnaissance Orbiter. Journal of Geodesy, 86-(3), 193-207. https://doi.org/10.1007/s00190-011-0509-4

Neumann, G. A. (2011). Lunar Reconnaissance Orbiter Lunar Orbiter Laser Altimeter reduced data record and derived products software interface specification, version 2.42, LRO-L-LOLA-4- GDR-V1.0, NASA Planetary Data System (PDS).
Petro, N. E., Keller, J. W., Cohen, B. A., & McClanahan, T. P. (2019). Ten years of the Lunar Reconnaissance Orbiter: Advancing lunar science and context for future lunar exploration. Paper presented at the 50th Lunar and Planetary Science Conference, Lunar and Planetary Institute, Houston, TX. https://www.hou.usra.edu/meetings/lpsc2019/pdf/2780.pdf

Smith, D. E., Zuber, M. T., Neumann, G. A., Lemoine, F. G., Torrence, M. H., McGarry, J. F., Rowlands, D. D., et al. (2010). Initial observations from the Lunar Orbiter Laser Altimeter (LOLA). Geophysical Research Letters, 37(18). https://doi.org/10.1029/2010GL043751

Smith, D. E., Zuber, M. T., Neumann, G. A., Mazarico, E., Head III, J., Torrence, M. H., & LOLA Science Team (2011). Results from the Lunar Orbiter Laser Altimeter (LOLA): global, high-resolution topographic mapping of the Moon.  Paper presented at the 42nd Lunar and Planetary Science Conference, Lunar and Planetary Institute, Houston, TX. https://www.lpi.usra.edu/meetings/lpsc2011/pdf/2350.pdf

Tooley, C. R., Houghton, M. B., Saylor Jr., S. S., Peddie, C., Everett, D. F., Baker, C. L., & Safdie, K. N. (2010). Lunar Reconnaissance Orbiter mission and spacecraft design. Space Science Review, 150, 23;;;62. https://doi.org/10.1007/s11214-009-9624-4

Wessel, P., & Smith, W. H. F. (1998). New, improved version of Generic Mapping Tools released. Eos, Transactions of the American Geophysical Union, 79(47), 579. https://doi.org/10.1002/2013EO450001

Williams, J. G., Boggs, D. H., & Folkner, W. M. (2008). DE421 Lunar Orbit, Physical Librations, and Surface Coordinates. JPL Interoffice Memorandum IOM 335-JW,DB,WF-20080314-001, 14 March.     Purpose
   
    To create a preliminary elevation model for the Moon.
   
Contact and Distribution
FormatDigital Elevation Model, Remote-sensing Data, Topographic MapAccess Constraintspublic domainAccess ScopePDSUse ConstraintsPlease cite authorsEditionMarch 11, 2014Supplemental Informationhttps://lola.gsfc.nasa.gov/, https://ode.rsl.wustl.edu/moon/, http://pds-geosciences.wustl.edu/missions/lro/lola_faq.htm, http://pds-geosciences.wustl.edu/missions/lro/default.htmNative Data Set EnvironmentISIS v3Astrogeology ThemeImage Processing, Remote Sensing, TopographyMission NamesLunar Reconnaissance OrbiterInstrument NamesLOLAOnline Package File Size8 GBOnline File Address2255 N. Gemini DriveContact CityFlagstaffContact StateAZContact Postal Code86001Contact Emailastroweb@usgs.gov
Data Status and Quality
Time Period of Content (start)2009-07-13Time Period of Content (stop)2013-07-18Currentness ReferenceGround conditionProgressIn WorkUpdate FrequencyAs neededLogical ConsistencyThe LOLA data were initially referenced to an internally consistent inertial coordinate system, derived from tracking of the LRO spacecraft. By adopting appropriate values for the orientation of the Moon as defined by the International Astronomical Union (IAU, Archinal and others, 2011), these inertial coordinates were converted into the planet-fixed coordinates (longitude and latitude) used on this map. The coordinate system defined for this product is the mean Earth/polar axis (ME) system. The ME system, sometimes called the mean Earth/rotation axis system, is the method most often used for cartographic products of the past (Davies and Colvin, 2000). Values for the orientation of the Moon were derived from the JPL DE/LE 421 planetary ephemeris (Williams and others, 2008; Folkner and others, 2008, 2009), rotated into the ME system. Longitude increases to the east and latitude is planetocentric as allowed in accordance with current international and NASA standards (Archinal and others, 2011; LRO Project and LGCWG, 2008). The intersection of the lunar equator and prime meridian occurs at what can be called the Moon's mean sub-Earth point. The concept of a lunar sub-Earth point derives from the fact that the Moon's rotation is tidally locked to the Earth. The actual sub-Earth point on the Moon varies slightly due to orbital eccentricity, inclination, and other factors, so a mean sub-Earth point is used to define the point on the lunar surface where longitude equals 0°. This point does not coincide with any prominent crater or other lunar surface feature (Archinal and others, 2011; LRO Project and LGCWG, 2008).Completeness ReportGaps between tracks of 1-2 km are common, and some gaps of up to 4 km occur near the equator. DEM points located in these gaps in LOLA data were filled by interpolation (Smith and others, 2011). The high and low elevations in the DEM are approximate because of the interpolation process used.Process DescriptionThe LOLA measurements were converted into a digital elevation model (DEM; Neumann and others, 2010) using Generic Mapping Tools software (Wessel and Smith, 1998), with a resolution of 256 pixels per degrees. In projection, the pixels are 118 meters in size at the equator. PRODUCT_CREATION_TIME = 2014-03-11T00:00:00Source TitleThe LOLA Data NodeSource Online  StatusPDS 3 ArchivedAttribute Accuracy ReportBest EffortHorizontal Positional Accuracy ReportBest EffortVertical Positional Accuracy ReportBest Effort
Geospatial Information
TargetMoonSystemEarthMinimum Latitude-90Maximum Latitude90Minimum Longitude-180Maximum Longitude180Direct Spatial Reference MethodRasterObject TypeGrid CellRaster Row Count (lines)46080Raster Column Count (samples)92160Bit Type (8, 16, 32)16Quad NameRadius A1737400Radius C1737400Bands1Pixel Resolution (meters/pixel)118.4505876Scale (pixels/degree)256Map Projection NameSimple CylindricalLatitude TypePlanetocentricLongitude DirectionPositive EastLongitude Domain-180 to 180

Project, clip, convert format (e.g. PDS, GeoTIFF, Jpeg) and download this product with Map a Planet

FGDC: XML Metadata
Ancillary Data
PDS3 Label
ISIS3 Label
thumb image
browse image
full image
Related Products

Moon LRO LOLA Color Shaded Relief 388m

Moon LRO South Pole DEM

Moon LRO LROC WAC Global Morphology Mosaic 100m

Moon Apollo 17 LROC NAC Landing Site Orthomosaic 50cm

Moon LRO LOLA - SELENE Kaguya TC DEM Merge 60N60S 59m
ISIS Logo
ISIS
Nomenclature Logo
Planetary Nomenclature
Map-a-Planet Logo
Map a Planet 2
PDS Logo
PDS Cartography and Imaging Science Node
AstroLink Logo
AstroLink
Photogrammetry Guest Faciltiy Logo
Photogrammetry Guest Facility
Pilot Logo
PILOT
MRCTR GIS Lab Logo
MRCTR GIS Lab
About
USGS Astrogeology Science Center