All Questions
Tagged with affine-transformation coordinate-system
60 questions
0
votes
0
answers
29
views
Custom datum transformation pyproj from unknown CRS to WGS84
I have a set of location points that represents real points on earth in Easting and Northing, but these records have lost their information of their original datum, e.g. :
Point
Original X
Original Y
...
0
votes
0
answers
58
views
Rasterio affine transform pixel size meaning
I'm using calculate_default_transform and reproject from rasterio.warp to transform geotiff file from one coordinate reference system (EPSG:25832) to another (EPSG:3857). As a result the pixel size ...
0
votes
0
answers
57
views
Fix for incorrect pixel size post transformation
I have a .tif in WGS 84 and one in EPSG 32632.
EPSG has a resolution of 10m (so the pixel sizes are 10).
When I convert WGS 84 into the EPSG 32632 CRS, even after specifying that the new resolution ...
- 91
1
vote
0
answers
41
views
PostGIS: Custom projection matrix from known points
I have a PNG image of a map (mercator) and known points in a table like this:
coord_map - geometry(Point, 0)
coord_real - geometry(Point, 4326)
010100000000000000000000000000000000000000 (0,0)
...
- 11
0
votes
0
answers
99
views
Convert a homogenous matrix from xyz coordinate system to zyx coordinate system
I want to change my camera pose from xyz coordinate system(where x-forward direction, y left and z up) to zyx coordinate system ( z-forward y- down and x left)
My camera pose is described in an ...
0
votes
1
answer
537
views
Rasterio Affine Transformation - Inverted
I am trying to understand how to properly set up an affine transformation using Rasterio. I have read through guides online and I have used the format provided but for some reason, it does not seem to ...
- 1
0
votes
0
answers
159
views
Correct Format for Transformation Rastio - from_origin
I am attempting to translate an HDF image to a GeoTIFF file and I encountered an issue when working with the transform module from rasterio. The issue comes when defining the transformation using '...
- 1
0
votes
0
answers
52
views
How to create a vector for coordinate transformation?
Very new to geospatial data here. I have a number of satellite images with lower left and upper right corner points in 3857. When I lay these images in QGIS against Google's satellite imagery, I'm ...
- 1
1
vote
1
answer
487
views
Unexpected crsTransform after reducing an ImageCollection: nominal scale shows 111319.49079327357
I have created an Image from an Imagecollection and I want to export it and to define the crsTransform parameter. I have read the documentation and the post on this forum and tried to export it , and ...
- 2,215
1
vote
0
answers
86
views
CRS rotation in Geoserver
I want to add rotation angle for a custom metric CRS, that I created before in the Geoserver. I tryed two methods:
Adding a WKT-stroke of rotated CRS to the Geoserver as usual. The example of such ...
- 53
1
vote
0
answers
91
views
Combining affine transformations and COMPOUNDCRS
I have a need to define a WKT CRS which is both affine-transformed in 2D, and vertically shifted in the height direction. What I would like to say is something like this:
COMPOUNDCRS["Modified ...
- 153
1
vote
0
answers
654
views
Change the origin point of the GeoTIFF files using Rasterio
I created a GeoTIFF file in Python (v3.8) using Rasterio (v1.2.10). The upper left corner is defined as the GeoTIFF file's origin point, which causes some difficulties in subsequent steps. Is there ...
- 11
2
votes
0
answers
232
views
How to remove banding artifact in R after reprojecting sf with st_transform
I have several global map layers I am plotting using ggplot() + geom_sf() + etc. All data is "+proj=longlat +datum=WGS84 +no_defs". My goal is to have Eurasia be a contiguous land mass. This ...
- 21
0
votes
1
answer
1k
views
GDAL coordinate transformation using C/C++ API (OGRCoordinateTransformation)
I'm in a puzzle of figuring out how to create correct coordinate transformations using GDAL C/C++ API.
I have referenced code from tutorial - https://gdal.org/tutorials/osr_api_tut.html#coordinate-...
- 103
2
votes
1
answer
499
views
Visually correcting systematic shift between two geodatasets using QGIS
The equivalent question has been asked wrt ArcMap, but I'm using QGIS and am also an inexperienced user.
I have multiple rasters, provided as GeoTIFFs, that are all projected to the same custom CRS ...
- 21
1
vote
0
answers
162
views
Transforming raster with matrix
I would like to apply a similarity transformation matrix (potentially a homography matrix) on a valid raster (with proper CRS, resolution, bounds, width height, etc).
The 1st image is the original ...
- 258
0
votes
1
answer
384
views
Projected coordinates points rotating horizontally
I have this Sentinel-2 picture that uses the WGS84 UTM zone 23S as projection type. My picture is rescaled to 3413 by 4741. My objective is to transform the X and Y pixels to EPSG 32723 and then from ...
- 53
0
votes
0
answers
143
views
How to convert an affine matrix setup in a metric coordinate system to a pseudo web mercator System
I have an affine matrix (with rotation and translation), which is setup in a local metric coordinate system.
1 0 0 10
M = 0 1 0 15
0 0 1 14
0 0 0 1
webMercPoint = (20000, 1000, 0)
...
- 101
8
votes
1
answer
2k
views
How to find a particular coordinate on a Sentinel image using proj:transform?
Background
I am retrieving images from the Sentinel satellites. I plot these images, using pixel coordinates. I would like to show points specified using latitude and longitude on those images.
Data ...
- 2,692
1
vote
0
answers
104
views
Custom Projection Definition - Excel Formula to ArcGIS Pro
I have been given an Excel formula for converting a custom projection coordinate point to a known defined system (state plane AZ East)
I'm trying to get create a custom projection definition in ArcGIS ...
- 131
2
votes
1
answer
78
views
Wrong transform after interpolating points
Following this post, I'ms trugglig wwith interpolation and projection of set of points. I have used griddata in order to interpolate, and seems like interpolation hapenns, however, the results look ...
- 2,215
2
votes
1
answer
527
views
Custom-derived CRS by affine transformation from EPSG:32633: Setting as layer-CRS works only for current QGis-session. Information lost after restart
I am working with QGIS 3.16, Win10. With help of another SE-user I created several derived CRS by an affine transformation from EPSG:32633 CRS (see below) in order to shift some related layers a ...
- 507
2
votes
1
answer
661
views
How do I find these values for an Affine Transformation (state plane to local grid)?
so I need to perform some kind of transformation. I have a bunch of points like so:
<record id="Localization Points" >
<record id="Point 1" >
<value name="Lat&...
- 143
1
vote
0
answers
125
views
Cannot georeference an image using a custom projection in QGIS
I am trying to georeference a map using a custom projection I created in QGIS. I get the error Transform failed: "Failed to compute GCP transform: Transform is not solvable" Does anyone have ...
- 11
1
vote
2
answers
1k
views
Reading CRS and Affine transformation from a NetCDF with CF conventions
I have a NetCDF with geographical information. I'd like to obtain a GeoTIFF file with rasterio. The problem is how to translate the geographical information from CF standard to something that GeoTIFF ...
- 389
3
votes
1
answer
2k
views
How to use affine transform parameters in a PROJ or WKT2 string
Can anyone show me how I can construct a PROJ or WKT2 string with affine transformation parameters to create a custom CRS in QGIS?
I wish to use the following affine transformation parameters to ...
- 101
3
votes
1
answer
272
views
rasterio reprojection produces no visible changes
I've been trying to perform translation of a raster using the reproject function from rasterio.warp.
I've calculated the translation amount in X and Y from a series of "control points" ...
- 41
2
votes
1
answer
523
views
Visually correcting systematic shift between two geodatasets using ArcMap
I often have two datasets that are obviously not georeferenced correctly. Most of the time a simple translation (shift by x/y) could fix this.
Is there an easy way to do this e.g. by specifying one or ...
- 2,554
0
votes
1
answer
240
views
Rotate, Scale, Translate in OpenLayers 4 and MapServer
Sometimes, I get spatial data (usually DGNs but I am not sure the format matters) that in the coordinate control sheet a 'Project Datum Factor' is noted. Say: 1.000073496. If I add them as is in ...
- 2,360
1
vote
0
answers
437
views
Transform rectangle from decimal degrees to Cartesian coordinates
I'm trying to achieve the following:
Given four coordinates (lat, lng) decimal degrees and a query coordinate within the formed shape, how can I transform these coordinates into Cartesian space where ...
- 111
1
vote
1
answer
2k
views
Apply same coordinate system to raster image and GeoJSON with rasterio
I have a dataset containing raster images and corresponding features in geojson files. I would like to use rasterio.features.rasterize to create a footprint image of the features on the raster image, ...
- 129
3
votes
0
answers
766
views
Coordinate transformation for antenna pointing
I may have bitten off more than I can chew with some coordinate transformation math. As a background, I'm working on a project where we need to mechanically point an antenna to a receiving tower. The ...
- 181
3
votes
1
answer
1k
views
Spatial transformation of latitudes and longitudes that preserve distances between points
I have a set of lat,lon points.
I would like to apply a spatial transformation (rotate, translate) to place these points somewhere else on the earth surface (for anonymization purposes), but preserve ...
0
votes
2
answers
328
views
Difference results from QGIS project on fly and gdalwarp for EPSG:4326 of Sentinel 2 UTM zone 32?
I always thought QGIS and gdalwarp should return same result for projecting a 2D Sentinel 2 image (e.g: T32UQD_20200421T102021_TCI_60m.jp2) from UTM 32 to EPSG:4326.
What I did was:
Open the .jp2 in ...
- 525
1
vote
1
answer
553
views
Creating Affine for 2D Linear Transformation of Fixed Bounding Box
I am currently trying to create an Affine which can transform data from longitude and latitude to a custom CRS that a client is using. Unfortunately the person who originally made this CRS has left ...
- 59
1
vote
1
answer
139
views
Transforming / stretching map or it's part in particular direction
I am trying to model ancient land surfaces we can 'rollback' geological folding and faulting. This is called palinspastic modelling or mapping.
So I have an area that due to folding and faulting was ...
- 471
2
votes
1
answer
868
views
Defining rotated CRS with WKT
I am trying to figure out how to display a rotated grid correctly with GeoServer with using the affine transformation in WKT as described here:
Defining coordinate reference system with rotation in ...
- 41
5
votes
1
answer
3k
views
Defining a coordinate system in WKT or PROJ format that has an Affine transformaiton and bounds
I need to define a local coordinate system via OGC WKT or PROJ form.
I know they have implemented Affine transformations in Proj 5 but I am unsure of how to implement it. The base projection is from ...
2
votes
0
answers
767
views
Reprojection problem with Rasterio
I have already processed so far NDSI map by using Rrasterio, NumPy, os and matplotlib.pyplot for plotting the data.
Now, I would like to eliminate permanent ice from NDSI map. I am facing one major ...
- 57
3
votes
1
answer
685
views
Converting coordinates from a local coordinate reference system
I am trying to define a local CRS in R based on the information in the screenshot image below
My R code is as follows:
I define the local CRS from the screenshot image
I apply it to my dummy ...
- 631
4
votes
1
answer
1k
views
Non-Earth System in QGIS?
What I need is a function for Non Earth Coordinates in QGIS (both projects and layers). They should be mappable to any projected coordinate system (like SWEREF99TM).
Why: In the Mining industry we ...
- 41
1
vote
1
answer
925
views
How to transform a shapefile with a MN03 (LV03) coordinate system into a shapefile with a MN95 (LV95) coordinate system?
Switzerland is currently introducing new coordinates: by the end of 2016, all federal and cantonal reference data (e.g. the Swiss cadastral system) sould be transformed to the new highly-accurate LV95 ...
2
votes
1
answer
3k
views
Transform coordinates from custom projection to lat long PyProj
I am trying to convert a pair of coordinates specified in US survey feet on a custom transverse mercator projection to lat/long using PyProj. I must be doing something wrong because I cannot get the ...
- 692
2
votes
0
answers
109
views
Reprojecting trouble: shapes saved in one EPSG load as another EPSG in QGIS and are 200m off?
I'm in a mixed working environment with some ESRI licenses and several people working with QGIS. For a project, we created data and received data in EPSG 31468 (Gauß-Krüger Zone 4). We now need to ...
- 387
2
votes
0
answers
403
views
Weather data in rotating grid without time dimension or vertical dimension or CRS
I am a beginner using QGIS 2.18.15. I am having some trouble to find the right coordinate reference system for wind speed data of Europe from the German Weather forecast service (DWD). (ftp://ftp-...
- 83
1
vote
0
answers
373
views
QGIS custom CRS from Mapinfo definition
How do I define a custom CRS in QGIS from the Mapinfo grid definition below?
"MIM Grid",1008,33,7,139.48639389,-20.69008981, 1.000063,0,0,7,0.998796124,-0.049054077,12520.6,0.049054077,0.998796124,...
1
vote
0
answers
148
views
Define Custom CRS for working with SVG
I can't seem to figure out a way to define a custom PROJ4 projection that can mirror geometries (e.g. inverse x to -x and y to -y and vice versa). Other than that, I don't need any transformations or ...
- 61
6
votes
2
answers
5k
views
Transforming data from a rotated pole lat-lon grid into regular lat-lon coordinates in QGIS?
I am working with data from a numerical weather model, the output of which is given in grid cells where I have the lat-lon coordinates of each cell's center point. The lat-lon specifications are, ...
- 151
5
votes
3
answers
463
views
Determining coordinates of a SHP file
i have a SHP file of a building converted from a DXF file. This SHP file is unprojected and have no specific coordinates (origin 0,0). I find WGS84 coordinates (lat, long) of this building from Google ...
- 69
2
votes
0
answers
2k
views
Transformation from local projection to WGS84
I'd like to reproject point that is in local coordinate system to the WGS84, using gdal. But exception "No PROJ.4 translation for source SRS, coordinate transformation initialization has failed" ...
