4
I think I've found the coordinate system you're looking for.
Your unknown coordinate reference system seems to be referred to in Rwanda as ITRF 2005 / TM Rwanda. It is a projected coordinate system using Transverse Mercator with the following parameters :
False Easting : 500,000
False Northing : 5,000,000
Central meridian : 30
Latitude of origin : 0
Scale ...
4
It looks like those latitudes and longitudes are based on arcminutes somehow. If I divide your numbers by 600,000, I arrive at -20.2806183, -45.6169416, close to the point you mention in Calcario Cazanga mining site in Arcos, Brazil:
If you apply this logic to all your coordinates (divide them by 600,000), and then assign epsg 4326, that might be the ...
2
Make a sample raster:
> r = raster()
> r[] = 10 * (1:ncell(r))
Note the value of the cell is going to be 10 times the cell number...
Which row is latitude -20?
> row = rowFromY(r, -20)
> row
[1] 111
Row 111. Ok. What are the cell numbers for row 111?
> cells = cellFromRow(r, row)
> cells[1:10]
[1] 39601 39602 39603 39604 39605 39606 ...
2
The arcpy.Geometry type supports a projectAs(spatial_reference, {transformation_name}) request, so it should be simple to project a singleton point (PointGeometry) from GeoJSON to the spatial reference of your target layer.
Of course, projecting the point is only the first step in your process, since selection and rendering is done at a layer level, not ...
2
It is some version of State Plane Nevada East with measurements in feet. For example, choosing ESRI:102707 from the EPSG database (http://epsg.io/102707) yields coordinates in that range near 35.16 N lat, -114 W long.
Note that there are several alternative specifications of the State Plane Nevada East projection, so if you don't get the exact one, your ...
2
If you want to write all the values of the ycoord field into a .txt file, you can execute the following code in the Python Console of QGIS (change the parameters with your own):
layer = iface.activeLayer() # your layer
features = layer.getFeatures() # features of your layer
with open('C:/Users/v.bre/Downloads/ycoord.txt','w') as my_file: # Change the path ...
1
You will need some coordinates for your postcodes. The easiest place to get them is by downloading code-point from the Ordnance Survey. You can then join the two layers together to make a map.
1
It seems to be EPSG:3421. (955999, 26397850) translates to Lat 35.1612555° Lon -114.5801733° at epsg.io.
When looking for projected CRS on the basis of given coordinates, I'm using SHPProjFinder tool. This is local Java program, based on similar algorihms as projfinder.com, but with much improved efficieny.
1
If you are interested in masking the data, I would recommend rioxarray.
An example of doing so can be found here.
Here is a targeted example for what you probably want to do:
from shapely.geometry import mapping
from shapely.wkt import loads
import rioxarray
geom = mapping(loads('POLYGON ((5.5937209 52.24012314 0,5.5936411 52.24048512 0,5.59413417 52....
1
Use arcpy.Project_management() to convert from one coordinate system to another (including to/from a projected from/to a geographic CRS).
See the documentation at: Project
You could convert either your shapefile, or your GeoJSON to an in_memory feature class with the same CRS as the other, and then work with that in_memory feature class instead of the ...
Only top voted, non community-wiki answers of a minimum length are eligible
