I have a polygon and a line feature class inside of a file geodatabase. Both the feature classes, and the map reference Washington State Plane North NAD83. After doing a spatial join moving attributes from the polygons to the lines I notice the line shape_length attribute had taken the values of the shape_length from the…

Yesterday we started a utility corridor survey. We have considerable property line locations from our data base that's tied to control based on the local HARN position NAD83(93). The utility...

I've been using and helping others configure and utilize RTK corrections from various services in Idaho, Oregon, and Florida. A couple of the major RTK service providers have switched from NAD83(2011) to ITRF 2014 with a specific Epoch as their reference frame. This is a time-dependent reference frame that accounts for the…

I am working with a large dataset which has datum in NAD27 and NAD83. I needed to convert the northing and easting coordinates with the NAD27 datum to NAD83 datum using Geopandas, which I was able to do successfully. However, many of the measured land surface elevations for the NAD27 data points are way off when we put the data into our LeapFrog model. I do not have a strong background in GIS, however as the resident computer programmer/engineer for my group I have been thrust into this role for the time being on this project. From my understanding, NAD27 and NAD83 are 2-D coordinate reference systems. I have read about some ellipsoidal NAD83 which can be used to calculate vertical elevation, but I don't know if any ellipsoidal datums exists for NAD27. The database also has latitude/longitudes available in WGS 84 coordinate system which I believe is an ellipsoid. Is there a way to estimate the elevation of the latitudes and longitudes using the WGS 84 ellipsoid system? I am in over my head a bit with this problem and am having trouble locating appropriate references and reading to understand this problem.

While trying to run commands in examples of GRASS GIS r.li module, I get the error. This is the output I get from commands in the examples grass_cmd_output(Language:Japanese) (Fri Sep 4 13:12:18 2020) g.gui.rlisetup (Fri Sep 4 13:12:21 2020) コマンド終了 (2 秒) (Fri Sep 4 13:12:24 2020) g.region raster=geology_30m -p projection: 99 (Lambert Conformal Conic) zone: 0 datum: nad83 ellipsoid: a=6378137 es=0.006694380022900787 north: 228500 south: 215000 west: 630000 east: 645000 nsres: 30 ewres: 30 rows: 450 cols: 500 cells: 225000 (Fri Sep 4 13:12:24 2020) コマンド終了 (0 秒) (Fri Sep 4 13:12:27 2020) r.li.patchdensity input=geology_30m conf=my_conf out=patchdens ERROR: Illegal areas disposition (Fri Sep 4 13:12:28 2020) コマンド終了 (0 秒) Following this article, there was the similar(=same?) problem only on Windows OS. I'm using Windows 10 64bit and GRASS 7.8.3 (Language:Japanese). Could you suggest any solutions to fix this error?

I created a polyline using a 30-day trial version of AutoCAD Civil 3D. Before drawing the polyline, I assigned the NAD83 California Zone 5 coordinate system to my drawing and set the units to US Survey Feet. I then exported the file as a 2018 2D AutoCAD DWG. However, ArcGIS Pro 3.5 Standard cannot r...

Hi, this is a question mostly to verify (or counter) my understanding of a very, very small difference in the results from suggested workflows for transforming NAD83 data to NAD83 (2011). I am not trying to split hairs over this question, but just trying to know why. I realize that I am exceeding th...

I am receiving the following error below when running gdalwarp to crop an input orthomosaic in the NAD83(2011)-Alaska Zone 1 CRS. The command still runs and the output cropped TIFF file renders fine when opening in QGIS or ArcGIS. What could be causing this error? I am trying to understand if GDAL does not fully support NAD83-Alaska Zone 1. Also, the input GeoTIFF and shapefile are both in NAD83-Alaska Zone 1. EPSG:6394 info: https://epsg.io/6394 gdalwarp command: gdalwarp -of GTiff -cutline [shapefile] -crop_to_cutline [input-tif] [output-tif] -wo NUM_THREADS=ALL_CPUS -co BIGTIFF=YES -co TILED=YES -co COMPRESS=JPEG -multi gdal version: 3.0.2 Terminal Output: -------------------------------------------- ERROR 6: Cannot find coordinate operations from `EPSG:6394' to `BOUNDCRS[SOURCECRS[GEOGCRS["NAD83(2011)",DATUM["NAD83 (National Spatial Reference System 2011)",ELLIPSOID["GRS 1980",6378137,298.257222101004,LENGTHUNIT["metre",1]],ID["EPSG",1116]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]]],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["latitude",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["longitude",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]]],ABRIDGEDTRANSFORMATION["Transformation from NAD83(2011) 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]]]]' Creating output file that is 185799P x 213951L. Processing rw-11-29a.tif [1/1] : 0ERROR 6: Cannot find coordinate operations from `EPSG:6394' to `BOUNDCRS[SOURCECRS[GEOGCRS["NAD83(2011)",DATUM["NAD83 (National Spatial Reference System 2011)",ELLIPSOID["GRS 1980",6378137,298.257222101004,LENGTHUNIT["metre",1]],ID["EPSG",1116]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]]],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["latitude",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["longitude",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]]],ABRIDGEDTRANSFORMATION["Transformation from NAD83(2011) 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]]]]'

EPSG:6404 Projected coordinate system for United States (USA) - Arizona - counties Coconino; Maricopa; Pima; Pinal; Santa Cruz; Yavapai. State law defines use of International feet (note: not US surve...

I have a point feature class in an Oracle 18c GDB that is projected as NAD83 UTM 17N (SRID = 26917). I would like to extract the X&Y coordinates from the points as WGS84 latitude & longitude decimal degrees (using SQL). The coordinates would be used to generate a Google Street View URL (hence the lat/long). I think this can be done by using SDE.ST_GEOMETRY's ST_TRANSFORM function: ST_Transform takes two-dimensional ST_Geometry data as input and returns values converted into the spatial reference specified by the spatial reference ID (SRID) you provide. However, I'm stuck on the SRID part: sde.st_x(sde.st_transform (shape, ??????)) as longitudex Is there a specific SRID I should be using to convert the UTM points to WGS84 lat/long?

I have been trying to get QGIS on Windows 10 to convert NAD83(2011) LongLatEllipsoidHeight files to NAD83(2011) / UTM zone 12N with NAVD88 orthometric elevations by applying GEOID18. It seems that recent versions of QGIS installed either by using osgeo4w-setup.exe or with an MSI installer no longer get a full set of proj files. Started with Ver 3.40 and did not have much luck changing the vertical; even the predefined CRS "NAD83(2011) + NAVD88 height EPSG:6349" did not do anything to transform the vertical. So following the seemingly logical assumptions that QGIS uses the version of proj that is installed with QGIS and cs2cs also uses proj, I decided to break the process down into steps. Opened the OSGeo4W Shell and entered - echo -109.0000000000 44.000000000 1500.000 | cs2cs +init=epsg:6319 +to +init=epsg:6341+5703 -f "%.3f" and things worked fine with cs2cs: successfully converted from 3D NAD83(2011) LLh to NAD83(2011) / UTM zone 12N easting/northing with NAVD / GEOID18 elevations shown below - 660349.4106 4873817.3333 1510.1279 Then I decided to update to the latest LTR version of QGIS. All the online wisdom indicates that doing that through osgeo4w-setup.exe is the way to go and that is how my Ver 3.40 was installed. For some reason that attempt went sideways and I ended up having to go through Control Panel to remove OSGeo4W and reinstall, again using osgeo4w-setup.exe. Now with Ver 3.44.12 freshly loaded I opened the new OSGeo4W Shell and tried the cs2cs command issued earlier and it no longer worked, it just converted the horizontal information to UTM zone 12N leaving the vertical information the same. Checked the contents of C:\OSGeo4W\share\proj and found there were very few files there - had checked there when 3.40 was installed and there were a bunch of geoid files and other things. It appears that at some time in the past I had downloaded from https://download.osgeo.org/proj/ a file named proj-datumgrid-north-america-1.3.zip and copied the contents to the C:\OSGeo4W\share\proj directory, so I did that again. Now cs2cs works fine again. NOTE: All the geoid files installed this way have a GTX extension. cs2cs must like them, as it resumed working? On a second Windows 10 PC that had QGIS Ver 3.34.11 installed, and the C:\OSGeo4W\share\proj directory has 438 items. There is no indication I downloaded anything extra at the time I installed that one. Vertical transformation does not work here either. NOTE: All the geoid files installed on this PC have TIF extensions. cs2cs also works on this PC, so cs2cs must like TIF files as well? On a third Windows 10 PC I uninstalled the QGIS that had been installed using osgeo4w-setup.exe, then installed Ver 3.44.11 from an MSI installer. With that installation method, the proj files are now at C:\Program Files\QGIS 3.44.11\share\proj. That directory only has 16 items, none of which look like geoid models. In that configuration cs2cs no longer converts the vertical information. On the "main" Windows 10 PC when I attempt to transform to NAD88 via EPSG:6349 using the Predefined CRS the "Select Datum Transformations" screen appears. The when the top line is highlighted the proj line reads - +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +inv +proj=vgridshift +grids=us_noaa_g2018u0.tif +multiplier=1 +step +proj=unitconvert +xy_in=rad +xy_out=deg Which indicates to me that proj is expecting a TIF file. Copied the file named us_noaa_g2018U0.tif from the second Windows 10 PC and put it in C:\OSGeo4W\share\proj on the main PC. Still did not work - no vertical transformation. This post is long because it details what I have tried and possibly eliminated so far. In the process I have waded through oceans of AI hallucination and found very little actual information on how to get this working. Any assistance will be greatly appreciated, from direction on which proj files to use from where on through setting up QGIS to accomplish this should-be-simple task. For "test data" one could simply use the data I sent to and received from cs2cs. Thanks!

I use Star*Net for certain projects, Geolab for others. But something in Star*Net always puzzled me. When setting up a coordinate system, they have "State Planes, NAD83" and "State Planes,...

I use Star*Net for certain projects, Geolab for others. But something in Star*Net always puzzled me. When setting up a coordinate system, they have "State Planes, NAD83" and "State Planes,...

I've been struggling to print my scalebar segments at exactly 1 inch for several hours now. Changing the scalebar width to 25.4mm (1 inch) does not solve the issue; my scalebar segments consistently print a few millimeters short of an inch. I'm wondering if there's a way to force QGIS to register my segment widths in terms of inches instead of mm to troubleshoot this issue. For context, my basemap orthoimage is in NAD83 North Carolina (ft-US). My map scale is "800". My scalebar units are in feet, label multiplier of 1.00. QGIS Version 3.30.2

I'm struggling with how to perform a vertical grid shift. I've been searching through previous postings, and I still don't understand what needs to be done. For example, I would like to reproject a DTM that's in EPSG:26911 (I assume the heights are ellipsoidal; I don't have any additional information) to using EPSG:6340+5703 GEOID18. I'm attempting to do this using GDAL with the following: gdal raster reproject -s "+proj=utm +zone=11 +datum=NAD83 +units=m +no_defs" -d "+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +geoidgrids=g2018u0.gtx" --resolution 1.0,1.0 --dst-nodata -9999 -r bilinear --to ALLOW_BALLPARK=NO --to ONLY_BEST=YES --of COG --co COMPRESS=ZSTD --co LEVEL=16 --co PREDICTOR=3 --co BIGTIFF=YES -i .asc -o .tif When I compare the elevation values before and after, there is about a 25-30 m change in value, which doesn't seem right (I'm working with data from southern Nevada). Can I stay in UTM (or any projected CRS) while doing this, or do I need convert to longitude/latitude, do the vertical grid shift, and then reproject back to a projected CRS? I've also been trying to do something similar with LiDAR data using PDAL, and I'm not sure I've been doing that correctly. I'm sure I'm doing something wrong, but I have no idea what.

I'm completing online course GST 101 for QGIS (my version is 2.18.3) and I'm stuck on the Buffer tool. All of my layers are set to CRS "NAD83(HARN)/New Mexico Central (ftUS) (EPSG:2903)" as instructed in the lab. The goal is to use the query builder to filter the municipal limits in Bernalillo county down to only Albuquerque. Then after filtering create a buffer layer of 5280 ft around it. The problem is I get an empty layer from the buffer tool with no entries in the attribute table. So I tried Buffer on the entire layer (with 0 ft buffer to pinpoint the problem areas) and this is the result. There appears to be a number of features in this layer that aren't recognized because I would have expected the entire layer to be selected. I've also tried the query builder to filter Albuquerque and save to a new layer. Then try to buffer the new layer, still results in an empty attribute table. I've tried the shp file included in the lab as well as downloading the file directly from here: City of Albuquerque GIS Data (Municipal Limits / Various municipal limits in Bernalillo county / SHP File). I'm still fairly new at GIS but I'm following a youtube video and pdf instruction but I'm clearly on a newer version of QGIS. Any suggestions on what the problems may be?

I'm using gdalwarp to clip and warp a large (but not huge) GeoTIFF (6.2 GB total size). The GDAL command I am using is: gdalwarp --config GDAL_CACHEMAX 80% -wm 80% -multi -wo NUM_THREADS=ALL_CPUS -t_srs "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0" -tr 10 10 -r bilinear -cutline [datasource] -cwhere [expression] -crop_to_cutline -cblend 100 [srcfile] [dstfile] I want to make sure I'm fully utilizing all available resources so this goes as fast as possible (hence the --config GDAL_CACHEMAX 80% -wm 80% -multi -wo NUM_THREADS=ALL_CPUS flags). However, when my CPU usage never goes above 100% (i.e. that is only 1 of 8 cores is being used) and memory usage is < 1 GB. Is there any way I can speed up this process?