New answers tagged

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