| Age | Commit message (Collapse) | Author |
|
Docs: Use +grids and tif file for deformation examples
|
|
|
|
|
|
|
|
|
|
|
|
Related to https://github.com/OSGeo/PROJ-data/pull/22
An entry is added in the ``other_transformation`` table, using a raw PROJ
string. Later, once the deformation model has been registered in EPSG, we'll
have to add code to map the EPSG transformation method and parameters to
+proj=defmodel
In the meantime, this works pretty well:
```
$ src/projinfo -s EPSG:4959 -t EPSG:7907
Candidate operations found: 1
-------------------------------------
Operation No. 1:
PROJ:NZGD2000-20180701, NZGD2000 to ITRF96 deformation model, unknown accuracy, New Zealand
PROJ string:
+proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=axisswap +order=2,1 +step +proj=defmodel +model=nz_linz_nzgd2000-20180701.json +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=rad +xy_out=deg
WKT2:2019 string:
COORDINATEOPERATION["NZGD2000 to ITRF96 deformation model",
VERSION["20180701"],
SOURCECRS[
GEOGCRS["NZGD2000",
DATUM["New Zealand Geodetic Datum 2000",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,3],
AXIS["geodetic latitude (Lat)",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["geodetic longitude (Lon)",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["ellipsoidal height (h)",up,
ORDER[3],
LENGTHUNIT["metre",1]],
ID["EPSG",4959]]],
TARGETCRS[
GEOGCRS["ITRF96",
DATUM["International Terrestrial Reference Frame 1996",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,3],
AXIS["geodetic latitude (Lat)",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["geodetic longitude (Lon)",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["ellipsoidal height (h)",up,
ORDER[3],
LENGTHUNIT["metre",1]],
ID["EPSG",7907]]],
METHOD["PROJ-based operation method: +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=axisswap +order=2,1 +step +proj=defmodel +model=nz_linz_nzgd2000-20180701.json +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=rad +xy_out=deg"],
USAGE[
SCOPE["unknown"],
AREA["New Zealand"],
BBOX[-55.95,160.6,-25.88,-171.2]],
ID["PROJ","NZGD2000-20180701"],
REMARK["New Zealand Deformation Model. Defines the secular model (National Deformation Model) and patches for significant deformation events since 2000"]]
```
```
$ echo "-41 173 0 2016.5" | PROJ_NETWORK=ON src/cct -d 8 +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=axisswap +order=2,1 +step +proj=defmodel +model=nz_linz_nzgd2000-20180701.json +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=rad +xy_out=deg
-40.99999402 172.99999938 0.00130333 2016.5000
```
```
$ echo "-41 173 0 2016.5" | PROJ_NETWORK=ON src/cs2cs -f "%.8f" EPSG:4959 EPSG:7907
-40.99999402 172.99999938 0.00130333 2016.5
```
|
|
|
|
the user writable directory on Windows (#2246)
Fixes #2245
|
|
Implement wink2 inverse by generic inversion of forward method
|
|
- Move the generic method initiated from adams_ws2 to a
pj_generic_inverse_2d() method
- Use it in adams_ws2
- Use it in wink2
Fixes https://github.com/qgis/QGIS/issues/35512
|
|
Fix identification of (one of the) ESRI WKT formulations of EPSG:3035
|
|
Database: add entry for Ranc08_Circe.mnt grid (New Caledonia)
|
|
Fixes https://github.com/qgis/QGIS/issues/36111
|
|
This geoid was added in EPSG 9.8.11 and has been added in PROJ-data
per https://github.com/OSGeo/PROJ-data/pull/21
|
|
Update to EPSG 9.8.11, ESRI 10.8.1 and import scope and remarks for conversion
|
|
|
|
|
|
|
|
This is in preparation for EPSG 9.8.11 import that supersedes a number of
grid transformation for US transformations, but the superseded and replacement
transformations don't operate on the same (source_crs, target_crs), which is
a bit weird. So in that situation, ignores the supersession.
|
|
Limit peirce_q to northern hemisphere, and fix images for adams_hemi, guyou and peirce_q
|
|
channel (#2236)
|
|
Add runtime checking that sqlite3 is >= 3.11
|
|
```
$ LD_LIBRARY_PATH=$HOME/install-sqlite-3.8.2/lib PROJ_DEBUG=1 src/projinfo -s @in.wkt -t @out.wkt --spatial-test intersects >/dev/null
SQLite3 version is 3.8.2, whereas at least 3.11 should be used
```
|
|
- Central lon for zone 2 should be -d10, not d10
- Extra lobe was missing for zone 11
- New figure generated
- New test suite values generated
|
|
Many fixes regarding BoundCRS, CompoundCRS, Geographic3D CRS with non-metre units
|
|
BoundCRS in the general case let to 0 result (contributes to fixes #2232)
|
|
vertical unit change
|
|
case of transformation between Compound of BoundVerticalCRS to GeographicCRS (contributes to fixes #2232)
|
|
a BoundCRS to a GeographicCRS (contributes to fixes #2232)
|
|
and non metre units (contributes to fixes #2232)
|
|
of a datum and the same datum (relates to #2232)
|
|
|
|
Inconsistent chainging of CRS (fixes #2232)
|
|
- The current implementation of the Interrupted Goode Homolosine
projection emphasizes land area. This is a compliment projection
that emphasizes ocean area.
- A value of lon0=-160 produces a reasonable real-world map.
|
|
|
|
Add a +proj=defmodel transformation for multi-component time-based deformation models
|
|
Allow importing WKT1 COMPD_CS with a VERT_DATUM[Ellipsoid,2002]
|
|
|
|
deformation models
Fixes #1001
Co-authored-by: Chris Crook <ccrook@linz.govt.nz>
|
|
coordinates
|
|
|
|
regenerate associated images
|
|
environment variable
|
|
Otherwise it returns junk (negative latitudes are mapped to
the same location as positive latitudes)
I'm a bit confused by PROJ peirce_q implementation. Looking at
the projection of the world, it looks like this matches the diamond
formulation of the right map shown at
https://desktop.arcgis.com/en/arcmap/latest/map/projections/peirce-quincuncial.htm,
but limited to the inner square of this diamond (which corresponds
to the northern hemisphere). We lack the 4 triangles on left, top,
right and bottom for the southern hemisphere. Furthermore, this
formulation of peirce_q does not seem to have the quincuncial
property of the square formulation (left images of the above ESRI
doc), or the one at https://en.wikipedia.org/wiki/Peirce_quincuncial_projection ...
|
|
createOperations(): when converting CompoundCRS<-->Geographic3DCrs, do not use discard change of ellipsoidal height if a Helmert transformation is involved
|
|
|
|
exporting it back as such (on the same object) (fixes #2228)
|
|
|
|
change (geographic3D horizontal)' operation
|