FCI Geolocation

[Glan]

Hi 

I want to obtain pixelwise geolocation of the MTG FCI W_XX-EUMETSAT-Darmstadt,IMG+SAT,MTI1+FCI-1C-RRAD-FDHSI-FD--CHK-BODY--DIS-NC4E_C_EUMT_20250701121004_IDPFI_OPE_20250701120707_20250701120801_N_JLS_O_0073_0030.nc product.

How can i do this.

Glan

[Martin Raspaud]

Hi,

Can you show us what you tried so far, so we can help you find the solution?

Best regards,
Martin
________________________________________
From: pyt...@googlegroups.com <pyt...@googlegroups.com> on behalf of Glan <mani...@gmail.com>
Sent: 31 July 2025 16:05:33
To: pytroll
Subject: [pytroll] FCI Geolocation

--
You received this message because you are subscribed to the Google Groups "pytroll" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pytroll+u...@googlegroups.com<mailto:pytroll+u...@googlegroups.com>.
To view this discussion, visit https://groups.google.com/d/msgid/pytroll/0c2f768f-3c6a-4c0d-805c-09952501eee0n%40googlegroups.com<https://groups.google.com/d/msgid/pytroll/0c2f768f-3c6a-4c0d-805c-09952501eee0n%40googlegroups.com?utm_medium=email&utm_source=footer>.

[David Hoese]

I actually had a separate user email me directly for this exact same problem. In their case they needed the lon/lats in binary files on disk. I'm not familiar with the FCI reader's defaults so it is possible these arrays are flipped, but this is the basic code:

from satpy import Scene

from glob import glob

scn = Scene(reader="fci_l1c_nc", filenames=glob("/data/satellite/fci/uncompressed/RC0067/*BODY*"))
scn.load(["vis_04"])
lons, lats = scn["vis_04"].attrs["area"].get_lonlats()
lons.tofile("fci_vis04_lons_f64.dat")
lats.tofile("fci_vis04_lats_f64.dat")

Dave

[Andrea Meraner]

Hi all, 

indeed, what Dave proposes will do. Just note that if you want the lon-lat arrays to be "upright" instead of flipped (likely the case), you can modify the load command as below to get the flipping done automatically:

scn.load(["vis_04"], upper_right_corner="NE")

Cheers, 
Andrea