[GRASS-dev] Invalid timestamp file for raster

[Nikos Alexandris]

Hallo Soeren,

I think I well understand the magic behind TGIS by now.

I am doing exactly as you explain: while importing Landsat imagery, the
script (of which I am not pleased as it is now, no tests yet :-/)

- gets the MTL file
  (https://github.com/NikosAlexandris/i.landsat.import/blob/master/i.landsat.import.py#L332)
- gets the timestamp
  (https://github.com/NikosAlexandris/i.landsat.import/blob/master/i.landsat.import.py#L424)
- prints it on-screen as TGIS expects
  (https://github.com/NikosAlexandris/i.landsat.import/blob/master/i.landsat.import.py#L539)
  yet, exports in a file if requested
  (https://github.com/NikosAlexandris/i.landsat.import/blob/master/i.landsat.import.py#L542
  and
  https://github.com/NikosAlexandris/i.landsat.import/blob/master/i.landsat.import.py#L782)

I am doing my job for the moment. Yet the timestamping part needs a complete rewrite.

I note, however, your tip to create a single STRDS for each band.

Thank for the magic,
Nikos

ps- Apologies for not finding the time to keep up with filing tickets in
trac, documenting stuff in our Wiki and be more participative.



* Sören Gebbert <soerengebbert at googlemail.com> [2017-10-24 12:56:43 +0200]:

>Hi Nikos,
>the strategy that i use is to create an input file for t.register that
>contains the raster names and the time stamps for registration:
>
>map_name | start_time | end_time
>...
>...
>
>In case of Landsat scenes, there are not end-time stamps. Create this
>file while importing and fill it with map names and time stamps.
>
>Use this file with t.register and register all maps in a STRDS that
>have been imported. I would suggest to create a single STRDS for each
>band, so you can use map algebra for computation.
>
>Best regards
>Sören
>
>
>
>
>2017-10-24 0:59 GMT+02:00 Nikos Alexandris <nik at nikosalexandris.net>:
>> * Veronica Andreo <veroandreo at gmail.com> [2017-10-23 22:54:03 +0200]:
>>
>>> Hi Nikos,
>>>
>>> Have you tried with t.register?
>>>
>>> You would need to set the start_time="2013-04-13 10:00:09.935626" plus an
>>> increment="10 seconds". I have not tested, but it should work.
>>>
>>> I think it's better than r.timestamp if you want to create time series,
>>> even though it would be great to fix the bug in r.timestamp, anyway.
>>
>>
>> I am r.timestamp-ing Landsat bands while importing them via a script and
>> sourcing meta from MTL files. I noticed that in some cases, the time
>> stamping did not work.
>>
>> Cheers, Nikos
>>
>> _______________________________________________
>> grass-dev mailing list
>> grass-dev at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-dev

-- 
Nikos Alexandris | Remote Sensing & Geomatics
GPG Key Fingerprint 6F9D4506F3CA28380974D31A9053534B693C4FB3 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 228 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20171025/8c313502/attachment.sig>