<div dir="ltr">If the 3x3 matrix describes a pure rotation you could probably use proj=helmert. If it includes shear, I think you are short of luck</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den tor. 12. dec. 2019 kl. 13.31 skrev Tammo Jan Dijkema <<a href="mailto:t.j.dijkema@gmail.com">t.j.dijkema@gmail.com</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi all,<br>
<br>
I'm new to proj, so sorry if this is a trivial question (I've read through docs for an hour without finding an answer, so it can't be too trivial).<br>
<br>
I have a custom coordinate system, which is defined through an affine transformation from ETRS89. All coordinates I care about are in xyz Cartesian coordinates. In particular, to get pqr coordinates, I do<br>
<br>
coords_pqr = (coords_etrs - pqrorigin_etrs) * (3x3 matrix)<br>
<br>
Now I'm generating raster images in PQR, which I would like to store in GeoTiff, so if I understand correctly I need to write a proj string into the GeoTiff to describe the projection. More background: I'm writing the PQR-images in python, and want to export them so that I can load them into QGIS, with appropriate coordinates.<br>
<br>
I see that an affine transformation is part of proj6, but I cannot expect users to have that installed. Is there a way to abuse another projection to achieve the same thing?<br>
<br>
If there is another channel to ask this question, please let me know. Thanks in advance!<br>
<br>
Tammo Jan Dijkema<br>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>