
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>

</head>

<body dir="ltr">

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

so next steps: I'll contact the owner and pass on the remarks on the documentation and I will ask him to participate in the discussion.</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

regards and thank you for your quick responses</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Karin</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div id="appendonsend"></div>

<hr style="display:inline-block;width:98%" tabindex="-1">

<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Van:</b> Even Rouault <even.rouault@spatialys.com><br>

<b>Verzonden:</b> donderdag 14 februari 2019 16:42<br>

<b>Aan:</b> proj@lists.osgeo.org<br>

<b>CC:</b> Kristian Evers; Praet Karin<br>

<b>Onderwerp:</b> Re: [PROJ] request for integration of optimized transformatiegrid BD72 to ETRS89</font>

<div> </div>

</div>

<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">

<div class="PlainText">On jeudi 14 février 2019 15:00:34 CET Kristian Evers wrote:<br>

> Hello Karin,<br>

> <br>

> As long as your grids are released on a suitable open source compatible<br>

> license, yes it is possible to add them to PROJ. If you are quick (as in<br>

> today) and open a pull request in the proj-datumgrid repository [0] it can<br>

> become part of the coming version 6 of PROJ. I am not sure if<br>

> transformations to and from EPSG:31370 will be using the grids by default<br>

> without also modifying the source code of PROJ. Even Rouault should be able<br>

> to clear that up though.<br>

<br>

This grid is present in the EPSG dataset, and imported in the proj.db of PROJ <br>

master. So this is mostly a matter of having it installed in the PROJ_LIB <br>

directory.<br>

<br>

Adding it in the grid_alternatives database of the table would also be nice so <br>

that it is tagged as being 'known' by PROJ with a hint to which package/URL it <br>

can be downloaded from<br>

<br>

On a master checkout without the grid, projinfo reports:<br>

<br>

<br>

$ PROJ_LIB=data src/projinfo -s EPSG:4313 -t EPSG:4258 -o PROJ<br>

-------------------------------------<br>

Operation n°1:<br>

<br>

EPSG:15928, BD72 to ETRS89 (2), 0.2 m, Belgium - onshore<br>

<br>

PROJ string:<br>

+proj=pipeline +step +proj=push +v_3 +step +proj=axisswap +order=2,1 +step <br>

+proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=cart +ellps=intl +step <br>

+proj=helmert +x=-106.8686 +y=52.2978 +z=-103.7239 +rx=-0.3366 +ry=0.457 <br>

+rz=-1.8422 +s=-1.2747 +convention=coordinate_frame +step +inv +proj=cart <br>

+ellps=GRS80 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step <br>

+proj=axisswap +order=2,1 +step +proj=pop +v_3<br>

<br>

-------------------------------------<br>

Operation n°2:<br>

<br>

unknown id, Null geographic offset from Belge 1972 to ETRS89, unknown <br>

accuracy, World<br>

<br>

PROJ string:<br>

<br>

<br>

-------------------------------------<br>

Operation n°3:<br>

<br>

EPSG:8369, BD72 to ETRS89 (3), 0.01 m, Belgium - onshore<br>

<br>

PROJ string:<br>

+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert <br>

+xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=bd72lb72_etrs89lb08.gsb <br>

+step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1<br>

<br>

<br>

<br>

Now if I create a fake (or true...) grid bd72lb72_etrs89lb08.gsb in data,<br>

<br>

-------------------------------------<br>

Operation n°1:<br>

<br>

EPSG:8369, BD72 to ETRS89 (3), 0.01 m, Belgium - onshore<br>

<br>

PROJ string:<br>

+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert <br>

+xy_in=deg +xy_out=rad +step +proj=hgridshift +grids=bd72lb72_etrs89lb08.gsb <br>

+step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1<br>

<br>

-------------------------------------<br>

Operation n°2:<br>

<br>

EPSG:15928, BD72 to ETRS89 (2), 0.2 m, Belgium - onshore<br>

<br>

PROJ string:<br>

+proj=pipeline +step +proj=push +v_3 +step +proj=axisswap +order=2,1 +step <br>

+proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=cart +ellps=intl +step <br>

+proj=helmert +x=-106.8686 +y=52.2978 +z=-103.7239 +rx=-0.3366 +ry=0.457 <br>

+rz=-1.8422 +s=-1.2747 +convention=coordinate_frame +step +inv +proj=cart <br>

+ellps=GRS80 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step <br>

+proj=axisswap +order=2,1 +step +proj=pop +v_3<br>

<br>

So the grid would be used in priority<br>

<br>

-- <br>

Spatialys - Geospatial professional services<br>

<a href="http://www.spatialys.com">http://www.spatialys.com</a><br>

</div>

</span></font></div>

</body>

</html>