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

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

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

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

[manisha]

С чем может быть связана ошибка 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


Спасибо

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

OUTPUT_PROJECTION_TYPE = GEO
...
OUTPUT_PIXEL_SIZE = 250

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

[manisha]

Спасибо, исправил, осталось от попыток выкинуть в альберса. 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

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

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

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

Содержание моих файлов (они сначала мозаицируются из двух 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!

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

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