Определение расстояния от Земли до Солнца в момент съемки

[Максим Дубинин]

Обсуждение темы:
Определение расстояния от Земли до Солнца в момент получения снимка

http://gis-lab.info/qa/earthsundist.html

[Svetla_GIS]

Ссылка на статью с сайта NASA "Перевод данных в показатели reflectance" изменила свое местоположение.
Старый адрес:
http://ltpwww.gsfc.nasa.gov/IAS/handboo ... ter11.html
Новый адрес:
http://landsathandbook.gsfc.nasa.gov/ha ... ter11.html

И измените, пожалуйста, название статьи ("расстояние" на "расстояния").

[Svetla_GIS]

Очень полезная информация. Спасибо.

[Максим Дубинин]

спасибо за внимательность, все исправлено

[Marina TS]

Добрый день. Подскажите, пожалуйста, в метаданных Landsat 8 есть параметр EARTH_SUN_DISTANCE. Можно ли его использовать напрямую для атмосферной коррекции? При расчете расстояния от Земли до Солнца по методике, описанной в статье, расстояние равно 1.01229973227962, в метаданных указано 1.0122853. Какое значение лучше использовать, рассчитанное по описанной методике или можно сразу брать из метаданных?

[bolotoved]

При ответе на этот вопрос я бы руководствовался следующими соображениями:

1. Я бы оценил насколько сильно влияет различие в 6-м знаке после запятой, подставив соотв. значения в расчетные формулы, может сказаться, что результаты будут отличаться незначительно.

2. В обоих случаях мы не можем установить точность определения расстояния - это просто данность. В этом случае я предпочел бы использовать значение из метаданных, т.к. именно оно используется по умолчанию автоматическими алгоритмами коррекции (например в GRASS) и на это значение проще сослаться в публикации.

[Marina TS]

Спасибо большое за ответ!:)

[Vedmed]

Уважаемые, а как найти координаты центра снимка, чтобы ввести его в Horizons? Координаты угловых точек снимка в файле метаданных есть, а центра не вижу, может не туда смотрю?

метаданныеShow

GROUP = L1_METADATA_FILE
GROUP = METADATA_FILE_INFO
ORIGIN = "Image courtesy of the U.S. Geological Survey"
REQUEST_ID = "0101404309991_00001"
LANDSAT_SCENE_ID = "LT51450252011249IKR01"
FILE_DATE = 2014-04-30T07:50:48Z
STATION_ID = "IKR"
PROCESSING_SOFTWARE_VERSION = "LPGS_12.4.1"
DATA_CATEGORY = "NOMINAL"
END_GROUP = METADATA_FILE_INFO
GROUP = PRODUCT_METADATA
DATA_TYPE = "L1T"
DATA_TYPE_L0RP = "TMR_L0RP"
ELEVATION_SOURCE = "GLS2000"
OUTPUT_FORMAT = "GEOTIFF"
EPHEMERIS_TYPE = "DEFINITIVE"
SPACECRAFT_ID = "LANDSAT_5"
SENSOR_ID = "TM"
SENSOR_MODE = "BUMPER"
WRS_PATH = 145
WRS_ROW = 025
DATE_ACQUIRED = 2011-09-06
SCENE_CENTER_TIME = 04:55:19.8670630Z
CORNER_UL_LAT_PRODUCT = 51.25264
CORNER_UL_LON_PRODUCT = 84.53809
CORNER_UR_LAT_PRODUCT = 51.27370
CORNER_UR_LON_PRODUCT = 88.06372
CORNER_LL_LAT_PRODUCT = 49.24973
CORNER_LL_LON_PRODUCT = 84.63924
CORNER_LR_LAT_PRODUCT = 49.26936
CORNER_LR_LON_PRODUCT = 88.01998
CORNER_UL_PROJECTION_X_PRODUCT = 328200.000
CORNER_UL_PROJECTION_Y_PRODUCT = 5680800.000
CORNER_UR_PROJECTION_X_PRODUCT = 574200.000
CORNER_UR_PROJECTION_Y_PRODUCT = 5680800.000
CORNER_LL_PROJECTION_X_PRODUCT = 328200.000
CORNER_LL_PROJECTION_Y_PRODUCT = 5457900.000
CORNER_LR_PROJECTION_X_PRODUCT = 574200.000
CORNER_LR_PROJECTION_Y_PRODUCT = 5457900.000
REFLECTIVE_LINES = 7431
REFLECTIVE_SAMPLES = 8201
THERMAL_LINES = 7431
THERMAL_SAMPLES = 8201
FILE_NAME_BAND_1 = "LT51450252011249IKR01_B1.TIF"
FILE_NAME_BAND_2 = "LT51450252011249IKR01_B2.TIF"
FILE_NAME_BAND_3 = "LT51450252011249IKR01_B3.TIF"
FILE_NAME_BAND_4 = "LT51450252011249IKR01_B4.TIF"
FILE_NAME_BAND_5 = "LT51450252011249IKR01_B5.TIF"
FILE_NAME_BAND_6 = "LT51450252011249IKR01_B6.TIF"
FILE_NAME_BAND_7 = "LT51450252011249IKR01_B7.TIF"
GROUND_CONTROL_POINT_FILE_NAME = "LT51450252011249IKR01_GCP.txt"
REPORT_VERIFY_FILE_NAME = "LT51450252011249IKR01_VER.txt"
BROWSE_VERIFY_FILE_NAME = "LT51450252011249IKR01_VER.jpg"
METADATA_FILE_NAME = "LT51450252011249IKR01_MTL.txt"
CPF_NAME = "L5CPF20110701_20110930.10"
END_GROUP = PRODUCT_METADATA
GROUP = IMAGE_ATTRIBUTES
CLOUD_COVER = 0.00
IMAGE_QUALITY = 9
SUN_AZIMUTH = 152.86738887
SUN_ELEVATION = 43.34897272
GROUND_CONTROL_POINTS_MODEL = 147
GEOMETRIC_RMSE_MODEL = 5.171
GEOMETRIC_RMSE_MODEL_Y = 4.333
GEOMETRIC_RMSE_MODEL_X = 2.823
GROUND_CONTROL_POINTS_VERIFY = 3056
GEOMETRIC_RMSE_VERIFY = 0.214
GEOMETRIC_RMSE_VERIFY_QUAD_UL = 0.169
GEOMETRIC_RMSE_VERIFY_QUAD_UR = 0.197
GEOMETRIC_RMSE_VERIFY_QUAD_LL = 0.242
GEOMETRIC_RMSE_VERIFY_QUAD_LR = 0.241
END_GROUP = IMAGE_ATTRIBUTES
GROUP = MIN_MAX_RADIANCE
RADIANCE_MAXIMUM_BAND_1 = 193.000
RADIANCE_MINIMUM_BAND_1 = -1.520
RADIANCE_MAXIMUM_BAND_2 = 365.000
RADIANCE_MINIMUM_BAND_2 = -2.840
RADIANCE_MAXIMUM_BAND_3 = 264.000
RADIANCE_MINIMUM_BAND_3 = -1.170
RADIANCE_MAXIMUM_BAND_4 = 221.000
RADIANCE_MINIMUM_BAND_4 = -1.510
RADIANCE_MAXIMUM_BAND_5 = 30.200
RADIANCE_MINIMUM_BAND_5 = -0.370
RADIANCE_MAXIMUM_BAND_6 = 15.303
RADIANCE_MINIMUM_BAND_6 = 1.238
RADIANCE_MAXIMUM_BAND_7 = 16.500
RADIANCE_MINIMUM_BAND_7 = -0.150
END_GROUP = MIN_MAX_RADIANCE
GROUP = MIN_MAX_PIXEL_VALUE
QUANTIZE_CAL_MAX_BAND_1 = 255
QUANTIZE_CAL_MIN_BAND_1 = 1
QUANTIZE_CAL_MAX_BAND_2 = 255
QUANTIZE_CAL_MIN_BAND_2 = 1
QUANTIZE_CAL_MAX_BAND_3 = 255
QUANTIZE_CAL_MIN_BAND_3 = 1
QUANTIZE_CAL_MAX_BAND_4 = 255
QUANTIZE_CAL_MIN_BAND_4 = 1
QUANTIZE_CAL_MAX_BAND_5 = 255
QUANTIZE_CAL_MIN_BAND_5 = 1
QUANTIZE_CAL_MAX_BAND_6 = 255
QUANTIZE_CAL_MIN_BAND_6 = 1
QUANTIZE_CAL_MAX_BAND_7 = 255
QUANTIZE_CAL_MIN_BAND_7 = 1
END_GROUP = MIN_MAX_PIXEL_VALUE
GROUP = PRODUCT_PARAMETERS
CORRECTION_GAIN_BAND_1 = "CPF"
CORRECTION_GAIN_BAND_2 = "CPF"
CORRECTION_GAIN_BAND_3 = "CPF"
CORRECTION_GAIN_BAND_4 = "CPF"
CORRECTION_GAIN_BAND_5 = "CPF"
CORRECTION_GAIN_BAND_6 = "INTERNAL_CALIBRATION"
CORRECTION_GAIN_BAND_7 = "CPF"
CORRECTION_BIAS_BAND_1 = "CPF"
CORRECTION_BIAS_BAND_2 = "CPF"
CORRECTION_BIAS_BAND_3 = "CPF"
CORRECTION_BIAS_BAND_4 = "CPF"
CORRECTION_BIAS_BAND_5 = "CPF"
CORRECTION_BIAS_BAND_6 = "CPF"
CORRECTION_BIAS_BAND_7 = "CPF"
END_GROUP = PRODUCT_PARAMETERS
GROUP = RADIOMETRIC_RESCALING
RADIANCE_MULT_BAND_1 = 0.766
RADIANCE_MULT_BAND_2 = 1.448
RADIANCE_MULT_BAND_3 = 1.044
RADIANCE_MULT_BAND_4 = 0.876
RADIANCE_MULT_BAND_5 = 0.120
RADIANCE_MULT_BAND_6 = 0.055
RADIANCE_MULT_BAND_7 = 0.066
RADIANCE_ADD_BAND_1 = -2.28583
RADIANCE_ADD_BAND_2 = -4.28819
RADIANCE_ADD_BAND_3 = -2.21398
RADIANCE_ADD_BAND_4 = -2.38602
RADIANCE_ADD_BAND_5 = -0.49035
RADIANCE_ADD_BAND_6 = 1.18243
RADIANCE_ADD_BAND_7 = -0.21555
END_GROUP = RADIOMETRIC_RESCALING
GROUP = PROJECTION_PARAMETERS
MAP_PROJECTION = "UTM"
DATUM = "WGS84"
ELLIPSOID = "WGS84"
UTM_ZONE = 45
GRID_CELL_SIZE_REFLECTIVE = 30.00
GRID_CELL_SIZE_THERMAL = 30.00
ORIENTATION = "NORTH_UP"
RESAMPLING_OPTION = "CUBIC_CONVOLUTION"
MAP_PROJECTION_L0RA = "NA"
END_GROUP = PROJECTION_PARAMETERS
END_GROUP = L1_METADATA_FILE
END

[Ariki]

Вычислить среднее из координат углов
Но если вы укажете координаты одного из углов, думаю, существенной разницы не будет.

[Vedmed]

То есть, погрешность некритична? Просто думал, что центр нужно считать с учетом сфероидности Земли... Ну, так даже проще, спасибо!

[Ariki]

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