Python getDOY Examples

Example #1

File: project_L3_daily_UTM.py Project: landsat/MODISprocessing

if len(argv) != 2:
    printUsage()
    exit(1)

back_date = argv[1]
if len(back_date) != 8:
    print("****************************")
    print("* date parameter malformed *")
    print("****************************")
    printUsage()
    exit(1)

_year = back_date[:4]
_month = back_date[4:6]
_day = back_date[6:]
_doy = getDOY(_year, _month, _day)
print("Processing date " + back_date + " (DOY = " + str(_doy) + ").")

modisL3_TSMPath = ensureTrailingSlash(
    ensureTrailingSlash(
        ensureTrailingSlash(modisL3_TSMBasePath + _year) + _month) + _day)
modisL3_UTMPath = ensureTrailingSlash(
    ensureTrailingSlash(
        ensureTrailingSlash(modisL3_TSM_UTMPath + _year) + _month) + _day)
modisL3_UTM_QLPath = ensureTrailingSlash(modisL3_TSM_UTM_QLPath + _year)

for _path in [modisL3_UTMPath, modisL3_UTM_QLPath]:
    if not exists(_path):
        print("Making directory: ", _path, " ...")
        makedirs(_path)

Example #2

File: create_MODIS_GEO_L1A_extract.py Project: RobTheOceanographer/MODISprocessing

def printUsage():
    print('Usage: create_MODIS_GEO_L1B.py date')
    print('where date is a string representing the date to process,')
    print('e.g. 20120607 for May 7, 2012.')


argc = len(argv)
if argc < 2:
    printUsage()
    exit(1)

_date  = str(argv[1])
_year  = _date[:4]
_month = _date[4:6]
_day   = _date[6:]
_doy   = getDOY(_year, _month, _day)

print(_year, _month, _day, _doy)

modisL1A_LACPath = ensureTrailingSlash(ensureTrailingSlash(ensureTrailingSlash(modisL1A_LACBasePath + _year) + _month) + _day)
modisL1A_subPath = ensureTrailingSlash(ensureTrailingSlash(ensureTrailingSlash(modisL1A_subBasePath + _year) + _month) + _day)
modisGEOPath     = ensureTrailingSlash(ensureTrailingSlash(ensureTrailingSlash(modisGEOBasePath     + _year) + _month) + _day)
modisGEO_subPath = ensureTrailingSlash(ensureTrailingSlash(ensureTrailingSlash(modisGEO_subBasePath + _year) + _month) + _day)

for _path in [modisL1A_subPath, modisGEOPath, modisGEO_subPath]:
    if not exists(_path):
        print("Making directory: ", _path, " ...")
        makedirs(_path)

try:
    srcList = listdir(modisL1A_LACPath)

Example #3

File: modis_end_to_end_dwnld_proc_backlog.py Project: RobTheOceanographer/MODISprocessing

    printUsage()
    exit(1)

date = argv[1]
if not len(date)==10:
    print("Wrong parameter: ", argv[1])
    printUsage()
    exit(1)

_year  = date[0:4]
_month = date[5:7]
_day   = date[8:10]
_backdate = _year + _month + _day

# get the backday value needed for get_MODIS_LAC.py
DOY_backday = getDOY(_year,_month, _day)
todate = datetime.date.today()
DOY_today = getDOY(todate.year, todate.month, todate.day)
_backday = DOY_today - DOY_backday



print("Processing products of ", _backdate, "...")

dwnldCmd     = 'python /home/uwe/cronjobs/nasa/modis/get_MODIS_LAC.py AQUA' + ' ' + str(_backday) # download L1A_LAC from nasa
l1aCmd       = 'python ' + l1aScript          + ' ' + _backdate     # L1A_LAC => L1A_sub, GEO_sub
l1bCmd       = 'python ' + l1bScript          + ' ' + _backdate     # L1A_sub, GEO_sub => L1B_sub
qlGenCmd     = 'python ' + quicklookScript    + ' ' + _backdate     # L1B_sub => L1B_RGB
l2Cmd        = 'python ' + l2Script           + ' ' + _backdate     # L1B_sub => L2_sub
gptTSMCmd    = 'python ' + gptTSMScript       + ' ' + _backdate     # L2_sub => L2_TSM
l3binningCmd = 'python ' + l3binningScript    + ' ' + _backdate     # L2_TSM => L3_TSM