GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


Импорт продуктов MODIS уровней 2G, 3, 4 с помощью MRT

Обсуждение материалов сайта: вопросы, замечания, предложения

Импорт продуктов MODIS уровней 2G, 3, 4 с помощью MRT

Сообщение Максим Дубинин » 11 фев 2006, 22:32

Обсуждение темы "Импорт продуктов MODIS уровней 2G, 3, 4 с помощью MRT"

http://gis-lab.info/qa/modisimport.html
  • 0

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6

Re: Импорт продуктов MODIS уровней 2G, 3, 4 с помощью MRT

Сообщение manisha » 24 авг 2006, 18:47

С чем может быть связана ошибка Invalid field ERS?

resample.exe -p d:\tmp\3\modis09q1.prm -g d:\tmp\3\out.log
***************************************************

MODIS Reprojection Tool (v3.3a January 2006)
Start Time: Thu Aug 24 18:31:51 2006

------------------------------------------------------------------

Error: ReadParameterFile : Reading Input Parameter File
: Invalid field ERS
Fatal Error, Terminating...


.prm

INPUT_FILENAME = D:\tmp\3\MOD09Q1.A2006217.h20v02.004.2006228110437.hdf
SPECTRAL_SUBSET = (1 1 1)
OUTPUT_FILENAME = D:\tmp\3\MOD09Q1.A2006217.h20v02.004.2006228110437.out.hdf
RESAMPLING_TYPE = NEAREST_NEIGHBOR
OUTPUT_PROJECTION_TYPE = GEO
OUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )
DATUM = WGS84
OUTPUT_PIXEL_SIZE = 250


Спасибо
  • 0

manisha
Завсегдатай
 
Зарегистрирован: 24 авг 2006
 
Сообщения: 410
Репутация: 10
Статьи: 2
Проекты: 1

Сообщение Максим Дубинин » 25 авг 2006, 00:55

OUTPUT_PROJECTION_TYPE = GEO
...
OUTPUT_PIXEL_SIZE = 250


А как это может быть, географическая проекция и 250 (м?) пиксел?
  • 0

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6

Сообщение manisha » 25 авг 2006, 15:32

Спасибо, исправил, осталось от попыток выкинуть в альберса. MRT тоже видимо удивляется как такое может быть :)
Тем не менее если исправить например на OUTPUT_PIXEL_SIZE = 0.02, то
Код: Выделить всё
Error: ReadParameterFile : Reading Input Parameter File
     : Invalid field RS

Работает если из файла с параметрами убрать и DATUM и OUTPUT_PIXEL_SIZE
Код: Выделить всё
MODIS Reprojection Tool (v3.3a January 2006)
Start Time:  Fri Aug 25 13:58:40 2006

------------------------------------------------------------------


Input image and reprojection info
---------------------------------
input_filename:          D:\tmp\3\MOD09Q1.A2006217.h20v02.004.2006228110437.hdf
output_filename:         D:\tmp\3\MOD09Q1.A2006217.h20v02.004.2006228110437.out.hdf
input_filetype:          HDF-EOS
output_filetype:         HDF-EOS
input_projection_type:   SIN
input_datum:             No Datum (use projection parameters)
output_projection_type:  GEO
output_datum:            No Datum (use projection parameters)
resampling_type:         NN
input projection parameters:  6371007.18 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
output projection parameters: 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00

input image corners (lat/lon):
    UL:  70.00 58.48
    UR:  70.00 87.71
    LL:  60.00 40.00
    LR:  60.00 60.00

input image spatial subset corners (lat/lon):
    UL:  70.00 58.48
    UR:  70.00 87.71
    LL:  60.00 40.00
    LR:  60.00 60.00

     band          select   type lines smpls pixsiz      min    max   fill
  1) sur_refl_b01       1  INT16  4800  4800 231.6564   -100  16000 -28672
  2) sur_refl_b02       1  INT16  4800  4800 231.6564   -100  16000 -28672
  3) sur_refl_qc_250m   1 UINT16  4800  4800 231.6564      0  32767  65535

SINUSOIDAL PROJECTION PARAMETERS:
 Radius of Sphere              = 6371007.181000 meters
 Longitude of C. Meridian      = 0 0 0.000000
 False Easting                 = 0.000000 meters
 False Northing                = 0.000000 meters

GEOGRAPHIC COORDINATES: None

NNResample : processing band sur_refl_b01

NNResample : processing band sur_refl_b02

NNResample : processing band sur_refl_qc_250m

Output image info
-----------------
output image corners (lat/lon):
    UL:  69.999999993717 39.999999989886
    UR:  69.999999993717 87.718398568823
    LL:  59.997868721334 39.999999989886
    LR:  59.997868721334 87.718398568823

output image corners (X-Y projection units):
    UL:  39.999999989887 69.999999993717
    UR:  87.718398568823 69.999999993717
    LL:  39.999999989887 59.997868721334
    LR:  87.718398568823 59.997868721334

     band               type lines smpls pixsiz      min    max   fill
  1) sur_refl_b01      INT16  2029  9680 0.0049   -100  16000 -28672
  2) sur_refl_b02      INT16  2029  9680 0.0049   -100  16000 -28672
  3) sur_refl_qc_250m UINT16  2029  9680 0.0049      0  32767  65535

End Time:  Fri Aug 25 14:02:56 2006

Finished processing!


Если импортировать в европейский альберс (как в примере в статье, изменены только пути/имена файлов), то
Код: Выделить всё
Error: ReadParameterFile : Reading Input Parameter File
     : Invalid field ERS


Где-то еще ошибка видимо
В итоге втянул через ердас, но...
  • 0

manisha
Завсегдатай
 
Зарегистрирован: 24 авг 2006
 
Сообщения: 410
Репутация: 10
Статьи: 2
Проекты: 1

Сообщение Максим Дубинин » 25 авг 2006, 19:18

Хмм, странно, надо еще попробовать и доразобраться все-таки в чем проблема:

Содержание моих файлов (они сначала мозаицируются из двух tile'ов, обрабатываются только нужные каналы из множества). И датум и pixel-size на месте, нужно быть внимательным со spectral_subset:

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

Prm:

Код: Выделить всё
INPUT_FILENAME = Z:\test\A2000097\A2000097_mos_b1.hdf
SPECTRAL_SUBSET = (1 0 0 0 0 0 0 0 0 0 0)
SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG
OUTPUT_FILENAME = Z:\test\A2000097\A2000097_mos_b1_res.tif
RESAMPLING_TYPE = NEAREST_NEIGHBOR
OUTPUT_PROJECTION_TYPE = GEO
OUTPUT_PROJECTION_PARAMETERS = (
  0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )
DATUM = WGS84
OUTPUT_PIXEL_SIZE = 0.002


Log:
Код: Выделить всё
MODIS Reprojection Tool (v3.3a January 2006)
Start Time:  Fri Aug 25 10:16:29 2006

------------------------------------------------------------------

Warning: ReadParameterFile : No Error
     : Only 11 out of 1 bands were addressed (0 or 1) in the SPECTRAL_SUBSET field. The remaining bands will not be processed (assumed to be 0).

Input image and reprojection info
---------------------------------
input_filename:          Z:\test\A2000097\A2000097_mos_b5.hdf
output_filename:         Z:\test\A2000097\A2000097_mos_b5_res.tif
input_filetype:          HDF-EOS
output_filetype:         GEOTIFF
input_projection_type:   SIN
input_datum:             WGS84
output_projection_type:  GEO
output_datum:            WGS84
resampling_type:         NN
input projection parameters:  6371007.18 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
output projection parameters: 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00

input image corners (lat/lon):
    UL:  50.00 31.11
    UR:  50.00 62.23
    LL:  40.00 26.11
    LR:  40.00 52.22

input image spatial subset corners (lat/lon):
    UL:  50.00 31.11
    UR:  50.00 62.23
    LL:  40.00 26.11
    LR:  40.00 52.22

     band          select   type lines smpls pixsiz      min    max   fill
  1) 250m 16 days red reflectance   1  INT16  4800  9600 231.6564      0  10000  -1000


SINUSOIDAL PROJECTION PARAMETERS:

   Radius of Sphere:     6371007.181000 meters
   Longitude of Center:     0.000000 degrees
   False Easting:      0.000000 meters
   False Northing:     0.000000 meters

NNResample : processing band 250m 16 days red reflectance

Output image info
-----------------
output image corners (lat/lon):
    UL:  49.999999995507 26.108145782925
    UR:  49.999999995507 62.230145782925
    LL:  39.999999995507 26.108145782925
    LR:  39.999999995507 62.230145782925

output image corners (X-Y projection units):
    UL:  26.108145782925 49.999999995507
    UR:  62.230145782925 49.999999995507
    LL:  26.108145782925 39.999999995507
    LR:  62.230145782925 39.999999995507

     band               type lines smpls pixsiz      min    max   fill
  1) 250m 16 days red reflectance  INT16  5000 18061 0.0020      0  10000  -1000

End Time:  Fri Aug 25 10:16:50 2006

Finished processing!
  • 0

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6

Сообщение Максим Дубинин » 07 сен 2006, 18:05

http://gis-lab.info/qa/modisimport-tool.html
С этим быстрее и проще.
  • 0

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6


Вернуться в Материалы сайта

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 1


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).