lonlat2pixline outputs incorrect pix/lin values

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
Post Reply
andrew.meredith
Posts: 38
Joined: Tue Dec 31, 2013 10:29 am America/New_York
Answers: 0

lonlat2pixline outputs incorrect pix/lin values

by andrew.meredith » Thu Dec 07, 2023 12:03 pm America/New_York

Hi,

The OCSSW program lonlat2pixline gives incorrect pixel/line values for some MERIS L1B input files.

I used lonlat2pixline to first get the pixel/line values for a bounding box covering Lake Erie and then used the output to subset the MERIS L1B file to what should be Lake Erie.

lonlat2pixline /data/MERIS/FR_4th/2011/EN1_MDSI_MER_FRS_1P_20110618T161228_20110618T162208_048632_0342_20171231T234956_0100.SEN3/xfdumanifest.xml -83.683300 41.275309 -78.614971 43.006608
# 2857 4329 9138 10163
sline=9138
eline=10163
spixl=2857
epixl=4329

l1bextract_safe_nc -v --spixl=2857 --epixl=4329 --sline=9138 --eline=10163 /data/MERIS/FR_4th/2011/EN1_MDSI_MER_FRS_1P_20110618T161228_20110618T162208_048632_0342_20171231T234956_0100.SEN3 EN1_MDSI_MER_FRS_1P_20110618T161228_20110618T162208_048632_0342_20171231T234956_0100.SEN3.LE3

As you can see in the following image the resulting L1B extract file does not cover Lake Erie.
l1bextract_safe_nc output - NOT LakeErie
l1bextract_safe_nc output - NOT LakeErie
l1bextract_safe_nc_extract_for_LErie.png (477.97 KiB) Not viewed yet
In a small sample, I've noticed the error appears to occur when there's a data gap in the source L1B file .
MERIS L1B source file
MERIS L1B source file
MERIS_l1b_source.png (190.98 KiB) Not viewed yet
I'm running on RHEL linux VM using SeaDAS 8.4 with OCSSW tag V2023.3

Thanks,
Andrew

Tags:

OB.DAAC - SeanBailey
User Services
User Services
Posts: 1465
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Answers: 1
Been thanked: 5 times

Re: lonlat2pixline outputs incorrect pix/lin values

by OB.DAAC - SeanBailey » Mon Dec 11, 2023 9:05 am America/New_York

Andrew,
This is not a lonlat2pixline issue, but rather an issue with the MERIS file. The geolocation is not correct, possibly related to the gap you pointed out. Below are two images created from a subset of the scene you reference. The first is the M13_radiance field, the second shows the same field with a landmask applied based on the file's geolocation information. As you can see, not exactly a good fit :(

This is something to bring up with folks at ESA.

Sean
meris_nolandmask.png
meris_nolandmask.png (213.87 KiB) Not viewed yet
The attachment meris_landmask.png is no longer available
Attachments
meris_landmask.png
meris_landmask.png (28.76 KiB) Not viewed yet

andrew.meredith
Posts: 38
Joined: Tue Dec 31, 2013 10:29 am America/New_York
Answers: 0

Re: lonlat2pixline outputs incorrect pix/lin values

by andrew.meredith » Tue Jan 02, 2024 10:22 am America/New_York

Thanks Sean. I will follow up with ESA.

Andrew

Post Reply