|
Hi rk and thanks for your reply. I did try to answer yesterday but I think due to the amount of images I attached to explain your questions, my post must have bounced? Therefore I'll try to answer this time only with less images. Regarding confirmation of the shift, rotate and scale error, we have a series of GCPs (Ground Control Points) around the project which are basically optical targets, which can be clearly seen in the image. These GCPs are coordinated to sub 1cm accuracy in HS2TN15. By comparing the true coordinates of any one of these to the one given by Manifold in the imported image we can compute the shift at this point. By using any 2 of them near to the opposite extremeties of the image, i.e. the longest baseline, we can compute scale and rotation etc. When I checked the projection info given by manifold after import it is saying EPSG:9300. However in saying this, I did change the default settings prior to the last import attempt, so that all imports should be to this transformation? I can't remember what the default setting was, in to try this again but I can tell that it is making no difference as the imported image is in the exact same position. I have confirmed this from some notes that I made on the given Manifold coordinates for some GCPs from the first import. Regarding the same projection info from QGIS (image confirmed to be in the correct position), I am getting the following: WKT BOUNDCRS[ SOURCECRS[ PROJCRS["unknown", BASEGEOGCRS["unknown", DATUM["unknown", ELLIPSOID["unknown",6378137,298.257222932958, LENGTHUNIT["metre",1, ID["EPSG",9001]]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433, ID["EPSG",9122]]]], CONVERSION["unnamed", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",51.9304992, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",-1.0175018, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",1, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",231211.675, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",333385.634, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["easting",east, ORDER[1], LENGTHUNIT["metre",1, ID["EPSG",9001]]], AXIS["northing",north, ORDER[2], LENGTHUNIT["metre",1, ID["EPSG",9001]]]]], 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 unknown 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]]]] Proj4 +proj=tmerc +lat_0=51.9304992 +lon_0=-1.0175018 +k=1 +x_0=231211.675 +y_0=333385.634 +a=6378137 +rf=298.257222932958 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +type=crs Extent Extent not known Hope this helps?
|