геопривязка данных Modis

[jj_kan]

Уважаемые еще вопрос такого плана.

Имеется тематическая карта на базе MODIS в формате tif
ее геопривязка:
DATUM - D_Sphere_ARC_INFO
Sp. Referens - GCS_Sphere_ARC_INFO

Попытки изменить геопривязку в ArcGis и ENVI ни к чему не привели.

Имеется пакет MRT. В нем файлы видны с расширением hdr.
В каком пакете это лучше сделать и
Каким образом можно изменить имеющуюся привязку на следующую:
Projected Coordinate System: Pulkovo_1942_GK_Zone
Geographic Coordinate System: GCS_Pulkovo_1942
Datum: D_Pulkovo_1942

2. Если можно ссылки где база тематических карт по Modis (Terra)или Aster или еще других спутн.систем.

Заранее благодарен.

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

DATUM и все? обычно это синусоидальная проекция
см. например http://gis-lab.info/qa/mod14a1-import.html

что такое "база тематических карт по MODIS"?

[Александр Мурый]

Вопрос №1.
MRT здесь ни к чему --- менять проекцию в HDF не надо, вырезать часть HDF не надо, перегонять в GeoTIFF тоже не надо. Значит, пользуйтесь перепроецированием растров в том софте, что у вас есть (ArcGis / ENVI / etc.)
Главное --- определите чётко, какие проекции на входе и на выходе. У исходного растра я увидел только датум.

Вопрос №2.
http://lmgtfy.com/?q=modis+terra+products и далее по этому шаблону

[jj_kan]

amuriy СПАСИБО конечно, но это тем карта, а не снимок MODIS.
вСЕ ПОПЫТКИ перепроецировать ее или в GCS_WGS_1984 или в Pulkovo 42 ни к чему не привели.

Не хотелось бы привязываться по точкам, это очень грубо.

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

[Александр Мурый]

amuriy СПАСИБО конечно, но это тем карта, а не снимок MODIS.
вСЕ ПОПЫТКИ перепроецировать ее или в GCS_WGS_1984 или в Pulkovo 42 ни к чему не привели.

Любая "тем. карта" в *.tif --- это по-любому растр, так? Или вы изобрели свой "фирменный" велосипед тип геоданных?
Поэтому вопросом на вопрос: что мешает спроецировать этот растр как надо? Опишите попытки подробнее.

Файлы *.tif не отображаются во вложениях на форуме. Как вариант --> заархивировать и выложить архив.

[jj_kan]

Уважаемый amuriy,
ничего я не пытаюсь изобретать.
Процедура стандартная.
Есть тематический растр модис, загружаем его в ArcMap c привязкой на лету, не получается.

Тогда пытаемся принудительно присвоить этой картинке пространственную привязку в ArcCatalog один в один одинаковую с привязкой которая есть в проекте у других полигонов или снимков, не получается. Лежит в сторонке.

И третий вариант, снова принудительно присваивая из базы имеющиеся проекции. Снова картинка в тороне.

В принципе все...

[Александр Мурый]

Можно попробовать задать проекцию в чём-то кроме ArcGIS или ENVI. Проще всего с помощью GDAL --> gdalwarp. Интересно "пощупать" сам растр, было бы недурно, если бы вы его прикрепили.

[gamm]

Есть тематический растр модис, загружаем его в ArcMap c привязкой на лету, не получается.

откуда взялся "тематический растр модис", и что это вообще такое? Вы не могли бы определить это общепринятыми терминами? иначе вас невозможно понять ... настораживает употребление слова картинка - [telepat]ее руками нарисовали в фотошопе (и это объясняет отсутствие привязки)[/telepat] ?

[jj_kan]

уважаемый gamm,
у меня есть тематическая карта в виде растра (сам не раз такие делал).
Она сделана на базе MOdis.

мне необходимо наложить ее на другую карту. Поэтому я пытаюсь выяснить как это сделать.
Т.к. перечисленными выше способами этого сделать не получилось. К сожалению.

[gamm]

уважаемый gamm,
у меня есть тематическая карта в виде растра (сам не раз такие делал).
Она сделана на базе MOdis.

к сожалению, вы так и не ответили на вопрос, а мои телепатические способности исчерпаны ... успехов.

[jj_kan]

Уважаемые,
мне все таки необходимо решить вопрос о привязке тематической карты.
Попробую еще раз.

Envi выдает следующие данные о привязке:
Proj: Arbitrary (Map Based)
Datum: D_Unknown

в ArcGis привязка вообще не определяется.
Spatial Reference: Undefined

Я понимаю, что это достаточно сложно. Карта на базе MODIS.
Но, можно ли, что то предпринять в этом случае?


Ув. gamm извините, но я возможно не понял, что конкретно Вы хотели выяснить.
Разве thematic map это не общепринятый термин?

Ув. amuriy спасибо за помощь, попытаемся перепроецировать данные в указанном Вами пакете.

Заранее благодарен.

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

пример дадите?
с вопросом определитесь?
у вас то сфероид, то не СК определена, определитесь.

приложенный выше tif - абсолютно нормален и должен распознаваться без проблем.

[Александр Мурый]

Только что догадался, что неотображаемый тиф можно просто скачать путём "сохранить как.." // facepalm.tif //

По теме.
Вывод

Код: Выделить всё

gdalinfo nelda2_h20v03_resize.tif

подтверждает слова Максима, что растр нормальный:

Код: Выделить всё

Driver: GTiff/GeoTIFF
Files: nelda2_h20v03_resize.tif
Size is 473, 288
Coordinate System is:
PROJCS["WGS 84 / UTM zone 39N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",51],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32639"]]
Origin = (2822501.119644165039062,6301979.610885620117188)
Pixel Size = (463.312713623046875,-463.312713623046875)
Metadata:
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 2822501.120, 6301979.611) ( 85d 7'20.92"E, 51d42'36.48"N)
Lower Left  ( 2822501.120, 6168545.549) ( 84d18'49.61"E, 50d42'40.09"N)
Upper Right ( 3041648.033, 6301979.611) ( 87d35'41.30"E, 50d49'40.05"N)
Lower Right ( 3041648.033, 6168545.549) ( 86d46'51.14"E, 49d51'24.29"N)
Center      ( 2932074.576, 6235262.580) ( 85d58'12.86"E, 50d46'56.58"N)
Band 1 Block=473x1 Type=Byte, ColorInterp=Red
Band 2 Block=473x1 Type=Byte, ColorInterp=Green
Band 3 Block=473x1 Type=Byte, ColorInterp=Blue

Для перепроецирования растра в соотв-ю зону Pulkovo42 / GK (видимо, это она) надо выполнить:

Код: Выделить всё

gdalwarp -t_srs "EPSG:28409" nelda2_h20v03_resize.tif nelda2_h20v03_resize_GKzone9.tif

[jj_kan]

Уважаемый amuriy извините все что Вы показали это пишется в командной строке?

Т.Е. выделяю файл, далее ввожу тот скрипт, что Вы прислали и потом для перевода в зону пулково последнюю строчку, правильно?
Еще вопрос для этих целей (или вместо Workstation) в ArcGis 9.3 был вставлен модуль Command Line, а в ArcGis 10 эту роль играет Piton?
Спрашиваю, потому, что никогда этим не занимался.

Для Максима выкладываю файл.
Эти файлы были мне присланы, поэтому произошла путаница с СК.

Спасибо за помощь.

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

1. Разберитесь окончательно в какой СК у вас данные, пока это не сделаете, все остальное бестолку. Научитесь пользоваться нормальными инструментами получения информации об СК.

2. Если те кто вам прислал данные неудосужились нормально прописать СК в геотифах, мысленно пошлите им луч ненависти и напишите Кранкиной или ее коллабораторам, Nelda - ее проект, попросите описание проекции.

3. Если писать не хотите - читайте мой самый первый ответ в этой теме.

PS: И последнее, пожалуйста, не называйте это "тематической картой MODIS". Называйте например "результатом классификации данных MODIS".