New answers tagged convert
0
When using pyproj, note the differences from various releases in how it is used to transform data. Here are a few examples using new/old capabilities based on the question:
Using pyproj >= 2.2.0
import pyproj
print(pyproj.__version__) # 2.4.1
print(pyproj.proj_version_str) # 6.2.1
proj = pyproj.Transformer.from_crs(3857, 4326, always_xy=True)
x1, y1 = (-...
2
The question was missing a little bit of info but with the dataset in hand, I was able to roughly find out what was going on.
First off, each point is actually composed of 6 different points that each contain a single attribute that belongs to the actual point, which is the one in the middle:
You need to isolate the true point with the tool Extract by ...
0
The -ct parameter of ogr2ogr is used in the same way than all the other parameters that takes their value from a string. Just add -ct and then the string to use, enclosed between double quotes. As documented in https://gdal.org/programs/ogr2ogr.html you need GDAL version 3.0 or higher for usint this option.
I tried you pipeline for the states.shp shapefile ...
2
You cannot open directly a DGN file in GRASS GIS. You can import it, but before you have to create a location (see: https://grass.osgeo.org/grass76/manuals/helptext.html) with the EPSG code (projection) of your DGN. You may create a new mapset also (beside the PERNAMENT).
After these steps start your GRASS session. You will see two windows. In the Layer ...
1
Just a sample... You could just divide by 10000000 and that'll work
latitudeE7 = c(205912610, 205913478)
longitudeE7 = c(-1034513483, -1034513949)
coords <- data.frame(latitudeE7, longitudeE7)
coords / 10000000
latitudeE7 longitudeE7
1 20.59126 -103.4513
2 20.59135 -103.4514
2
Proj6 seems to consider that EPSG:4258 and EPSG:4326 mean the same and conversion between them does not change anything.
projinfo -o PROJ -s "EPSG:4258" -t "EPSG:4326"
Candidate operations found: 1
-------------------------------------
Operation n 1:
EPSG:1149, ETRS89 to WGS 84 (1), 1.0 m, Europe - ETRS89
projinfo -k operation EPSG:1149
PROJ string:
+proj=...
0
Some additional information:
To get the code in Ian Turton's answer (GeoTools 22.2, Java 8) to run properly I had to add the following dependencies to the pom.xml file:
<dependency>
<groupId>org.geotools</groupId>
<artifactId>gt-geotiff</artifactId>
<version>${geotools.version}</version>
</dependency&...
0
Let's assume that the Earth is flat and that the meridians are parallel. Let's assume also that BHL-No.3080 point has coordinates Easting = 0 , Northing = 0. We need it to calculate the coordinates of the other points in a planimetric way.
Here we go: We will make a spreadsheet for the calculation of the planimetric coordinates:
The first row has the ...
3
If you have a matching .wld and .prj file then the gt-image module can read and write jpeg and png files.
Once you've read the file in to a GridCoverage you can reproject it and write it out with no difficulty:
AbstractGridFormat format = GridFormatFinder.findFormat(input);
Hints hints = null;
if (format instanceof GeoTiffFormat) {
hints = new Hints(...
-1
I used VeryUtils DWG to Image Converter to convert my 200+ DWG files to JPG files from command line, I wrote a .bat file to do this work, it worked fine for me.
Top 50 recent answers are included
