<div dir="ltr"><div>Hi</div><div><br></div><div>When there are several geoid models for a vertical system, is there a way to force in WKT2 to use a specific one?</div><div>I know in the US I can use the values in the geoid_model table (for instance using the method AuthorityFactory::getGeoidModels), GEOID18 or GEOID12B. Or GGM10 for Mexico, that also uses NAVD88.<br></div><div><span style="font-family:monospace">        GEOIDMODEL["GEOID18"]</span></div><div><br></div><div>However there are other places where those geoid models are not added to the proj.db database, like</div><div><span style="font-family:monospace">fr_ign_RAF18.tif: France, IGN France, RGF93 (EPSG:4965) to NGF-IGN69 height (EPSG:5720)</span></div><div>and <br></div><div><span style="font-family:monospace">fr_ign_RAF09.tif: France, IGN France, RGF93 (EPSG:4965) to NGF-IGN69 height (EPSG:5720)</span></div><div></div><div><br></div><div>Even Rouault mentioned some time ago that the <a href="http://epsg.org">epsg.org</a> purpose of that field was to set not just a name, but the transformation id itself, as shown in <a href="http://docs.opengeospatial.org/is/18-010r7/18-010r7.html">http://docs.opengeospatial.org/is/18-010r7/18-010r7.html</a> , section 10.4<br></div><div><span style="font-family:monospace">        GEOIDMODEL["CGG2013",ID["EPSG",6648]]</span></div><div><br></div><div>So to force the French vertical crs to use one or the other we should write something like <br></div><div><span style="font-family:monospace">        GEOIDMODEL["RAF09",ID["EPSG",8371]]</span></div><div>or <br></div><div><span style="font-family:monospace">        GEOIDMODEL["RAF18",ID["EPSG",8885]]</span></div><div><br></div><div></div><div>Is that right? Is that working in PROJ 9.1.0? I tried in cs2cs, but the output was identical. The log showed that both were using raf18 tif file.<br></div><div>Could it have many entries to choose one automatically? I see in the 
table geoid_model that there are several "operation_code" for the same 
"name", and PROJ chooses one.<br></div><div><br></div><div>I was looking the code trying to understand how the tag GEOIDMODEL is used in a transformation, but didn't get anything.</div><div><br></div><div>Thanks.</div><div><br></div><div>PS I do not know how many vertical systems can use different geoid models, nor what may happen in the future when some countries update their geoid models.<br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div></div></div>