Hi Dirk,

there is no such thing as "just exporting". What exactly do you do? What does gdalinfo say about the original file and the exported file?

Cheers, Hannes

Am 05.10.23 um 10:10 schrieb Stulgies Dirk via QGIS-User:

# Unexpected change in pixel values during processing (Metashape/QGIS)

## Software

* QGIS 3.28

* Agisoft Metashape 2.0.1

## Description

I used Agisoft Metashape to create RGB GeoTiff. The source for this GeoTiff are images taken with a UAV (DJI M3M).

When loading the GeoTiff into QGIS and processing the file, the pixel values for the same position differ between the original file and the process result.

My first thought was that I somehow resampled, converted, or transformed the pixel values during the processing.

So I simplified the process step by step to the bare minimum, but the result stayed the same.

Even just exporting the file leads to changed pixel values.

## CRS not recognized

Something I noticed while trying to solve this problem is that QGIS doesn't recognize the GeoTiffs CRS (EPSG:32632) as a predefined CRS but as a custom CRS. This is curious because I use a predefined CRS from Agisoft Metashape to create the GeoTiff.

I tried to reproject the file and also to just change the CRS in the file's QGIS properties, but the pixel values still differ after processing the file.

### CRS defined in Agisoft Metashape:

```

Projected Coordinate System: WGS 84 / UTM zone 32N (EPSG::32632)

Linear Units: metre (EPSG::9001)

Projection Method: Transverse Mercator

    Latitude of natural origin: 0

    Longitude of natural origin: 9

    Scale factor at natural origin: 0.9996

    False easting: 500000

    False northing: 0

Geographic Coordinate System: WGS 84 (EPSG::4326)

Angular Units: degree (EPSG::9102)

Geodetic Datum: World Geodetic System 1984 ensemble (EPSG::6326)

Ellipsoid: WGS 84 (EPSG::7030)

Prime Meridian: Greenwich (EPSG::8901)

```

### CRS regognized in QGIS:

```

BOUNDCRS[

SOURCECRS[

PROJCRS[

"WGS 84 / UTM zone 32N",

BASEGEOGCRS[

"WGS 84",

DATUM[

"World Geodetic System 1984",

ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]

],

PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],

ID["EPSG",4326]

],

CONVERSION[

"UTM zone 32N",

METHOD["Transverse Mercator",ID["EPSG",9807]],

PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],

PARAMETER["Longitude of natural origin",9,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],

PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],

PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],

PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]

],

CS[Cartesian,2],

AXIS["easting",east,ORDER[1],LENGTHUNIT["metre",1]],

AXIS["northing",north,ORDER[2],LENGTHUNIT["metre",1]],

ID["EPSG",32632]

]

                ],

TARGETCRS[

GEOGCRS[

"WGS 84",

DATUM[

"World Geodetic System 1984",

ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]

],

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]],

ID["EPSG",4326]

]

                ],

ABRIDGEDTRANSFORMATION[

"Transformation from WGS 84 to WGS84",

METHOD["Position Vector transformation (geog2D domain)",ID["EPSG",9606]],

PARAMETER["X-axis translation",0,ID["EPSG",8605]],

PARAMETER["Y-axis translation",0,ID["EPSG",8606]],

PARAMETER["Z-axis translation",0,ID["EPSG",8607]],

PARAMETER["X-axis rotation",0,ID["EPSG",8608]],

PARAMETER["Y-axis rotation",0,ID["EPSG",8609]],

PARAMETER["Z-axis rotation",0,ID["EPSG",8610]],

PARAMETER["Scale difference",1,ID["EPSG",8611]]

                ]

]

```

## Question

The questions I have are:

1. What could be the reason for this behavior?

2. How do I prevent the pixel values from changing?

Kind regards,

i.A. Dirk Stulgies

**********************************

NPZ Innovation GmbH

Hohenlieth-Hof, D-24363 Holtsee

Tel.+49-4351-736 277 <tel:%2B49-4351-736%20277>__

Fax +49-4351-736 271 <tel:%2B49-4351-736%20271>

Mail:[email protected]__

Web:www.npz-innovation.de <http://www.npz-innovation.de/>

Firmensitz Hohenlieth, Holtsee

Amtsgericht Kiel, HRB 15342 Kl

Geschäftsführer: Dr. Gunhild Leckband, Thorsten Pommer


------------------------------------------------------------------------
This email was scanned by Bitdefender

_______________________________________________
QGIS-User mailing list
[email protected]
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user

--
Johannes Kröger / GIS-Entwickler/-Berater

---------------------------------------------
Aufwind durch Wissen!
Web-Seminare und Online-Schulungen
bei derwww.foss-academy.com
---------------------------------------------

WhereGroup GmbH
c/o KK03 GmbH
Lange Reihe 29
20099 Hamburg
Germany

Tel: +49 (0)228 / 90 90 38 - 36
Fax: +49 (0)228 / 90 90 38 - 11

[email protected]
www.wheregroup.com
Geschäftsführer:
Olaf Knopp, Peter Stamm
Amtsgericht Bonn, HRB 9885
-------------------------------
_______________________________________________
QGIS-User mailing list
[email protected]
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Reply via email to