search

GFZ / enpt

EnMAP Processing Tool - A Python package for pre-processing of EnMAP Level-1B data
https://git.gfz-potsdam.de/EnMAP/GFZ_Tools_EnMAP_BOX/EnPT
Other
11 stars 2 forks source link

Projection type of the raster output files 'Geographic’ not working #7

Closed msoppa1 closed 9 months ago

msoppa1 commented 10 months ago

Description

When I run the atmospheric correction Polymer using the command line and setting the projection to Geographic, the output is always in UTM. The same happens with Sicor or if I use the GUI.

What I Did

enpt --CPUs 4 --auto_download_ecmwf True --average_elevation 0 --blocksize 100 --deadpix_P_algorithm spectral --deadpix_P_interp_spatial linear --deadpix_P_interp_spectral linear --disable_progress_bars True --drop_bad_bands True --enable_ac True --mode_ac water  --polymer_additional_results True --ortho_resampAlg gauss -od /Volumes/AWIMari/EnMAP_v4_paper/Level2/Polymer/BahiaBlanca/20220928/ -of GTiff --path_l1b_enmap_image /Volumes/AWIMari/EnMAP_v4_paper/Level1/BahiaBlanca/20220928/ENMAP01-____L1B-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z.ZIP --polymer_root /Users/msoppa/Documents/polymer/polymer-v4.16.1 --run_deadpix_P True --scale_factor_boa_ref 10000 --scale_factor_toa_ref 10000 -tgtprj Geographic --threads -1 --vswir_overlap_algorithm vnir_only
danschef commented 10 months ago

@msoppa1 I cannot reproduce this with the test dataset I have here. I guess the dataset you are using was processed with different parameters by the ground segment. So could you just upload your dataset (ENMAP01-____L1B-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z.ZIP) somewhere?

msoppa1 commented 10 months ago

I just sent you the link to download the file per email

danschef commented 10 months ago

I still cannot reproduce this issue with your data. I improved some tests in EnPT (see here) to make sure the parameters set in the GUI or via the command line are correctly parsed by the backend code and also added tests to make sure the output projection is really the desired one but however, all tests are passing and I could not find a problem in the code.

Here is the output of EnPT when I run it with your settings via the GUI (note the line 2023/12/21 12:10:55: Computed common target extent of orthorectified image (xmin, ymin, xmax, ymax in EPSG 4326): (-61.94776475443619, -39.225301828518184, -61.51176497427002, -38.91094067002115)):

QGIS version: 3.32.2-Lima
Qt version: 5.15.8
Python version: 3.11.4
GDAL version: 3.7.2
GEOS version: 3.12.0-CAPI-1.18.0
PROJ version: Rel. 9.2.1, June 1st, 2023
PDAL version: 2.5.6 (git-version: 18a095)
Algorithm started at: 2023-12-21T12:05:30
Algorithm 'EnPT - EnMAP Processing Tool (v0.8.2)' starting…
Input parameters:
{ 'CPUs' : 32, 'auto_download_ecmwf' : False, 'average_elevation' : 0, 'blocksize' : 100, 'deadpix_P_algorithm' : 1, 'deadpix_P_interp_spatial' : 0, 'deadpix_P_interp_spectral' : 0, 'disable_progress_bars' : True, 'drop_bad_bands' : True, 'enable_ac' : True, 'json_config' : '', 'mode_ac' : 1, 'n_lines_to_append' : None, 'ortho_resampAlg' : 2, 'output_dir' : '/home/gfz-fe/EnPT/Output/20231221__2', 'output_format' : 0, 'output_interleave' : 2, 'path_dem' : None, 'path_earthSunDist' : '', 'path_l1b_enmap_image' : '/home/gfz-fe/scheffler/temp/EnPT/20231220_issue_Mariana/ENMAP01-____L1B-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z.ZIP', 'path_l1b_enmap_image_gapfill' : '', 'path_reference_image' : None, 'path_solar_irr' : '', 'polymer_additional_results' : True, 'polymer_root' : '/home/gfz-fe/scheffler/python/polymer-v4.16.1', 'run_deadpix_P' : False, 'scale_factor_boa_ref' : 10000, 'scale_factor_toa_ref' : 10000, 'target_epsg' : None, 'target_projection_type' : 1, 'threads' : -1, 'vswir_overlap_algorithm' : 2, 'working_dir' : '' }

Calling EnPT with the following command:
enpt --CPUs 32 --auto_download_ecmwf False --average_elevation 0 --blocksize 100 --deadpix_P_algorithm spatial --deadpix_P_interp_spatial linear --deadpix_P_interp_spectral linear --disable_progress_bars True --drop_bad_bands True --enable_ac True --mode_ac water --ortho_resampAlg gauss --output_dir /home/gfz-fe/EnPT/Output/20231221__2 --output_format GTiff --output_interleave pixel --path_l1b_enmap_image /home/gfz-fe/scheffler/temp/EnPT/20231220_issue_Mariana/ENMAP01-____L1B-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z.ZIP --polymer_additional_results True --polymer_root /home/gfz-fe/scheffler/python/polymer-v4.16.1 --run_deadpix_P False --scale_factor_boa_ref 10000 --scale_factor_toa_ref 10000 --target_projection_type Geographic --threads -1 --vswir_overlap_algorithm vnir_only

The log messages of the EnMAP processing tool are written to the *.log file in the specified output folder.
/home/gfz-fe/mambaforge/envs/enpt_full_dev/lib/python3.11/site-packages/enpt/cli.py:190: RuntimeWarning: No digital elevation model provided. Note that this may cause uncertainties, e.g., in the atmospheric correction and the orthorectification.
return EnPTConfig(**parsedArgs_to_user_opts(cli_args))
2023/12/21 12:05:33: Reading metadata for VNIR detector...
2023/12/21 12:05:34: Reading metadata for SWIR detector...
2023/12/21 12:05:34: Reading image masks in VNIR sensor geometry.
2023/12/21 12:05:40: Converting DN values to radiance [mW/m^2/sr/nm] for VNIR detector...
2023/12/21 12:05:41: Converting DN values to radiance [mW/m^2/sr/nm] for SWIR detector...
2023/12/21 12:05:43: Computing SNR from VNIR TOA radiance.
2023/12/21 12:05:52: Computing SNR from SWIR TOA radiance.
2023/12/21 12:05:55: Computing VNIR geolayer...
2023/12/21 12:05:57: Computing SWIR geolayer...
2023/12/21 12:06:17: No DEM for the VNIR detector provided. Falling back to an average elevation of 0 meters.
2023/12/21 12:06:17: No DEM for the SWIR detector provided. Falling back to an average elevation of 0 meters.
2023/12/21 12:06:20: Transforming the 'mask_landwater' attribute from VNIR to SWIR sensor geometry.
2023/12/21 12:06:23: Starting atmospheric correction for VNIR and SWIR detector in 'water' mode. Source radiometric unit code is 'TOARad'.
2023/12/21 12:06:23: Running atmospheric correction in 'water' mode, i.e., ACWater/Polymer is applied to water surfaces only. Note that land surfaces will NOT be included in the EnMAP L2A product.
vnir
2023/12/21 12:06:23: Running polymer...
2023/12/21 12:06:23: Currently only tested with ERA5 data, i.e. set True for "auto_download_ecmwf"
Initializing ENMAP product of size (1024, 1000)
2023/12/21 12:06:23: Landmask is EnMAP L1B
Download 2022-09-28 14:00:00 -> /home/gfz-fe/scheffler/python/polymer-v4.16.1/ANCILLARY/ERA5/2022/09/28/era5_20220928_140000.nc
2023-12-21 12:06:24,203 INFO Welcome to the CDS
2023-12-21 12:06:24,203 INFO Sending request to https://cds.climate.copernicus.eu/api/v2/resources/reanalysis-era5-single-levels
2023-12-21 12:06:24,567 INFO Request is queued
2023-12-21 12:06:25,619 INFO Request is running
2023-12-21 12:06:27,172 INFO Request is completed
2023-12-21 12:06:27,172 INFO Downloading https://download-0008-clone.copernicus-climate.eu/cache-compute-0008/cache/data9/adaptor.mars.internal-1703156785.9395874-9703-19-3494c306-8d89-4ab5-871d-35736197b245.nc to /home/gfz-fe/scheffler/python/polymer-v4.16.1/ANCILLARY/ERA5/2022/09/28/era5_20220928_140000.nc.tmp (9.9M)
Download 2022-09-28 15:00:00 -> /home/gfz-fe/scheffler/python/polymer-v4.16.1/ANCILLARY/ERA5/2022/09/28/era5_20220928_150000.nc
2023-12-21 12:06:32,554 INFO Welcome to the CDS
2023-12-21 12:06:32,554 INFO Sending request to https://cds.climate.copernicus.eu/api/v2/resources/reanalysis-era5-single-levels
2023-12-21 12:06:32,893 INFO Request is queued
2023-12-21 12:06:35,499 INFO Request is running
2023-12-21 12:06:37,802 INFO Request is completed
2023-12-21 12:06:37,802 INFO Downloading https://download-0006-clone.copernicus-climate.eu/cache-compute-0006/cache/data3/adaptor.mars.internal-1703156795.6351657-20166-15-7a8f5c57-00dc-4333-9d34-1a899a4ca948.nc to /home/gfz-fe/scheffler/python/polymer-v4.16.1/ANCILLARY/ERA5/2022/09/28/era5_20220928_150000.nc.tmp (9.9M)
Starting processing at 2023-12-21 12:06:41.127358
Processing block: size (100, 1000), offset (0, 0)
Processing block: size (100, 1000), offset (100, 0)
Processing block: size (100, 1000), offset (200, 0)
Processing block: size (100, 1000), offset (300, 0)
Processing block: size (100, 1000), offset (400, 0)
Processing block: size (100, 1000), offset (500, 0)
Processing block: size (100, 1000), offset (600, 0)
Processing block: size (100, 1000), offset (700, 0)
Processing block: size (100, 1000), offset (800, 0)
Processing block: size (100, 1000), offset (900, 0)
Processing block: size (24, 1000), offset (1000, 0)
Done in 0:04:10.474827
2023/12/21 12:10:54: Joining results of atmospheric correction.
/home/gfz-fe/mambaforge/envs/enpt_full_dev/lib/python3.11/site-packages/enpt/processors/atmospheric_correction/atmospheric_correction.py:319: RuntimeWarning: invalid value encountered in cast
in_detector.data = (out_detector * self.cfg.scale_factor_boa_ref).astype(np.int16)
2023/12/21 12:10:55: Starting orthorectification...
2023/12/21 12:10:55: Computed common target extent of orthorectified image (xmin, ymin, xmax, ymax in EPSG 4326): (-61.94776475443619, -39.225301828518184, -61.51176497427002, -38.91094067002115)
2023/12/21 12:10:55: Orthorectifying VNIR data using 'gauss' resampling algorithm...
/home/gfz-fe/mambaforge/envs/enpt_full_dev/lib/python3.11/site-packages/sensormapgeo/transformer_3d.py:258: RuntimeWarning: Switched multiprocessing algorithm from 'bands' to 'tiles' due to a timeout in 'bands' mode. Consider using the LLVM instead of the GNU build of OpenMP to fix this issue (install, e.g., by 'conda install -c conda-forge _openmp_mutex=*=1_llvm'.
warn(msg, RuntimeWarning)
2023/12/21 12:13:37: Orthorectifying SWIR data using 'gauss' resampling algorithm...
/home/gfz-fe/mambaforge/envs/enpt_full_dev/lib/python3.11/site-packages/pyresample/kd_tree.py:259: UserWarning: Possible more than 8 neighbours within 30 m for some data points
get_neighbour_info(source_geo_def,
2023/12/21 12:17:26: Merging VNIR and SWIR data...
2023/12/21 12:17:27: Orthorectifying 'mask_landwater' attribute...
2023/12/21 12:17:29: Orthorectifying 'mask_clouds' attribute...
2023/12/21 12:17:31: Orthorectifying 'mask_cloudshadow' attribute...
2023/12/21 12:17:33: Orthorectifying 'mask_haze' attribute...
2023/12/21 12:17:35: Orthorectifying 'mask_snow' attribute...
2023/12/21 12:17:37: Orthorectifying 'mask_cirrus' attribute...
2023/12/21 12:17:40: Orthorectifying 'polymer_logchl' attribute...
2023/12/21 12:17:43: Orthorectifying 'polymer_logfb' attribute...
2023/12/21 12:17:46: Orthorectifying 'polymer_rgli' attribute...
2023/12/21 12:17:49: Orthorectifying 'polymer_rnir' attribute...
2023/12/21 12:17:53: Orthorectifying 'polymer_bitmask' attribute...
2023/12/21 12:17:55: Setting all pixels to nodata that have values in the VNIR or the SWIR only...
2023/12/21 12:18:04: Generating L2A metadata...
2023/12/21 12:18:20: Write product to: /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z
Writing GeoArray of size (1004, 1392, 209) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-SPECTRAL_IMAGE.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLASSES.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLOUD.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLOUDSHADOW.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_HAZE.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_SNOW.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CIRRUS.TIF.
Writing GeoArray of size (1004, 1392, 3) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_VNIR.TIF.
Writing GeoArray of size (1004, 1392, 3) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_SWIR.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_LOGCHL.TIF.
2023/12/21 12:18:24: Currently, L2A dead pixel masks cannot be saved yet.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_LOGFB.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_RGLI.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_RNIR.TIF.
Writing GeoArray of size (1004, 1392) to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_BITMASK.TIF.
2023/12/21 12:18:24: Currently, the L2A metadata XML file does not contain all relevant keys and contains not updated values!
2023/12/21 12:18:24: Writing metadata to /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-METADATA.XML
2023/12/21 12:18:24: L2A product successfully written!
2023/12/21 12:18:24: Total runtime of the processing chain: 0:12:52.209030

ready.
The output folder '/home/gfz-fe/EnPT/Output/20231221__2' contains:

ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z
ENMAP01-____L1B-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z.log

/home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z
...where the folder 'ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z' contains:

ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_BITMASK.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_LOGCHL.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_LOGFB.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_RGLI.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-ACOUT_POLYMER_RNIR.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-METADATA.XML
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CIRRUS.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLASSES.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLOUD.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_CLOUDSHADOW.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_HAZE.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_QUALITY_SNOW.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_SWIR.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-QL_VNIR.TIF
ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-SPECTRAL_IMAGE.TIF

Execution completed in 775.56 seconds (12 minutes 56 seconds)
Results:
{'outfolder': '/home/gfz-fe/EnPT/Output/20231221__2',
'outraster': '/home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-SPECTRAL_IMAGE.TIF',
'success': True}

Loading resulting layers
Algorithm 'EnPT - EnMAP Processing Tool (v0.8.2)' finished

Here is the output of gdalinfo for the L2A output:

>>> gdalinfo /home/gfz-fe/EnPT/Output/20231221__2/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z/ENMAP01-____L2A-DT0000003967_20220928T143657Z_002_V010400_20231205T013044Z-SPECTRAL_IMAGE.TIF

[...]
Size is 1392, 1004
Coordinate System is:
GEOGCRS["WGS 84",
    ENSEMBLE["World Geodetic System 1984 ensemble",
        MEMBER["World Geodetic System 1984 (Transit)"],
        MEMBER["World Geodetic System 1984 (G730)"],
        MEMBER["World Geodetic System 1984 (G873)"],
        MEMBER["World Geodetic System 1984 (G1150)"],
        MEMBER["World Geodetic System 1984 (G1674)"],
        MEMBER["World Geodetic System 1984 (G1762)"],
        MEMBER["World Geodetic System 1984 (G2139)"],
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]],
        ENSEMBLEACCURACY[2.0]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    USAGE[
        SCOPE["Horizontal component of 3D system."],
        AREA["World."],
        BBOX[-90,-180,90,180]],
    ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (-61.948078423342793,-38.910626936329983)
Pixel Size = (0.000313668906594,-0.000313733691115)
[...]
Corner Coordinates:
Upper Left  ( -61.9480784, -38.9106269) ( 61d56'53.08"W, 38d54'38.26"S)
Lower Left  ( -61.9480784, -39.2256156) ( 61d56'53.08"W, 39d13'32.22"S)
Upper Right ( -61.5114513, -38.9106269) ( 61d30'41.22"W, 38d54'38.26"S)
Lower Right ( -61.5114513, -39.2256156) ( 61d30'41.22"W, 39d13'32.22"S)
Center      ( -61.7297649, -39.0681212) ( 61d43'47.15"W, 39d 4' 5.24"S)
[...]

Also note that I am using version 0.8.2 of the GUI and version 0.19.5 of the backend. However, actually there were no changes since your versions that could have changed something in the output projection.

Could you run EnPT again on your machine and post the log messages here if the output really has the wrong projection?

danschef commented 9 months ago

I close now, just re-open if the issue still persists.